弹性波阻抗理论和实现方法
基于柯西先验分布的叠前弹性阻抗反演方法及应用
基于柯西先验分布的叠前弹性阻抗反演方法及应用郝前勇;宗兆云【摘要】以贝叶斯反演为代表的概率化反演方法既能考虑观测数据的不确定性,又可以考虑待求解参数的先验信息,在实际地震反演中备受青睐.经研究表明,在柯西先验信息下获取的反演结果更具有稀疏性,且具有高分辨特征.叠前弹性阻抗反演是一种基于多角度部份叠加剖面的叠前地震反演方法,信息量丰富,计算效率高.这里在贝叶斯框架下,实现了基于柯西先验的叠前弹性阻抗反演方法,并提取了对储层流体敏感的弹性参数.实际资料应用表明,基于柯西先验的弹性阻抗反演方法合理可靠,具有较高的分辨能力,且提取的弹性参数能够较好地吻合实际钻遇结果.%Inversion based on probability scheme represented by Bayesian inversion has been widely utilized in seismic data interpretation. On one hand, it contains the uncertainty in observed data, on the other hand, it also contains the prior information for the model parameters to be estimated. Cauchy prior probability density information leads to more sparse inversion result with high resolution. Elastic impedance inversion as a kind of pre-stack seismic inversion methods with partial angle-stack has also been utilized in practice widespread for its high efficiency and abundant information. In this paper, a kind of elastic impedance Bayesian inversion with Cauchy prior information is proposed, and with the inverted elastic impedances in different angle, Lame parameters are extracted. Real data test shows that the proposed method is reliable and has high resolution. The extracted lame parameters match with drilling result well.【期刊名称】《物探化探计算技术》【年(卷),期】2012(034)006【总页数】6页(P717-722)【关键词】柯西先验;弹性阻抗;参数提取;贝叶斯框架【作者】郝前勇;宗兆云【作者单位】中国石油大学地球科学与技术学院,山东青岛266555;中国石油大学地球科学与技术学院,山东青岛266555【正文语种】中文【中图分类】P631.3+40 前言概率化反演是求解不适定反演问题的有效方法之一。
ept软件模块培训教材_弹性波阻抗反演
EPT软件/功能模块培训系列教材GMAX v1.0 – Elastic Impedance Inversion弹性波阻抗反演1.模块功能2.原理和方法3.参数和使用说明4.应用注意事项EPT公司, 1.模块功能弹性波阻抗反演（EI）是叠前地震反演重要方法之一。
基于流体置换模型技术，应用纵波声波时差、密度、泥质含量、孔隙度、含水饱和度和骨架、流体的各种弹性参量，反演井中横波速度。
根据井中纵波速度、横波速度和密度计算井中弹性波阻抗，在复杂构造框架和多种储层沉积模式的约束下，采用地震分形插值技术建立可保留复杂构造和地层沉积学特征的弹性波阻抗模型，使反演结果符合研究区的构造、沉积和异常体特征。
采用广义线性反演技术反演各个角度的地震子波，得到与入射角有关的地震子波。
在每一个角道集上，采用宽带约束反演方法反演弹性波阻抗，得到与入射角有关的弹性波阻抗。
最后对不同角度的弹性波阻抗反演纵横波阻抗，进而获得泊松比等弹性参数, 对储层的几何、物性和含流体特性进行精细描述。
叠前地震弹性参数反演的关键技术包括:◆基于流体替换模型的井中横波速度反演技术◆与偏移距有关的子波反演◆复杂地质构造情况下弹性阻抗建模◆纵横波阻抗、泊松比、拉梅系数和剪切模量反演2.原理和方法地震反射振幅不仅与分界面两侧介质的地震弹性参数有关，而且随入射角变化而变化。
叠前弹性波阻抗反演技术利用不同炮检距的地震数据及横波、纵波、密度等测井资料，联合反演出与岩性、含油气性相关的多种弹性参数，综合判别储层物性及含油气性。
正是由于叠前弹性波阻抗反演利用了大量地震及测井信息，所以进行多参数分析的结果较叠后声阻抗反演在可信度方面有很大提高，可对含油气性进行半定量—定量描述。
传统的A VO 和岩石物理分析是提取和分析纵横波速度的异常变化来确定孔隙流体和岩性的变化。
纵横波速度和密度对反射系数的重要性，可以从平面波的Zoeppritz方程中看出。
但是，在波动方程中，Md2U/dX2= d2U/dX2，（U是位移），其表达式并不与地震波速度直接相关，而与岩石密度和弹性模量相关。
波阻抗反演基本原理及Jason反演软件介绍
二、反演技术的发展历程
3.90年代早中期，测井约束反演诞生和发展时期。 突破传统地震频带的限制，具有比直接反演更高的分辨率。 商品化的软件如strata、BCI、Jason等。 4.90年代中后期 ，非线性反演理论为基础的各种反演算法 出现和发展时期。如地震特征反演、随机反演、模拟退火反 演、概率法神经网络反演、遗传算法、小波反演等。
子波与反射系数的褶积得到地震记录
S(t)=W(t)*R(t)
一、概述
地质模型
低速层 高速层
3.地震褶积模型
反射系数 分步褶积 地震响应
1 2
低速层 3 高速层 低速层 高速层 6 更高速层 4 5
一、概述
4.结论
1.地震数据的生成基于褶积原理，因此地震剖面绝不 是地质剖面，简单的把地震剖面当地质剖面解释会产 生错误，特别是薄互层沉积时，错误会更多。 2.由于地震反演基于地震数据，而地震数据的分辨率 低，识别不了薄互层沉积时，地震反演的结果也就识 别不了薄互层。 3.为了得到好的反演结果，仅仅用地震数据是完成不 了的，因此许多新技术应运而生，通过不同的数学方 法，把地震数据与测井数据结合，试图得到高分辨率 的反演结果，识别薄互层，来指导或直接找油找气。
• Principal component analysis （主因子分析） • Model estimation （模型估算） • Model generator （模型生成器）
四、Jason6.2平台反演模块的构成
Invermod
Inv er c t ra e
Ro
ck tra ce
Statmod
一、概述
1.前言
地震反演的目的: 根据地震资料，反推出地下介质的 波阻抗、速度和密度等岩石地球物理参 数的分布，估算储层参数，并进行储层 预测，以便为油气田的勘探和开发提供 可靠的基础资料。
弹性波的基本理论
第一章 弹性波的基本理论第一节 弹性理论概述 一、弹性介质的概念地震勘探中将地层叫做介质。
由于地震勘探是研究人工激发的地震波在岩层中的传播规律来探测地下地质体的的存在和确定岩土物理力学参数的地球物理方法，它的地球物理前提是岩矿石间的弹性差异，因此需要研究介质的弹性性质。
人工激震后，岩石附近的质点发生破碎，介质产生的是塑性形变；远离震源的介质质点会发生振动，发生体积和形状的变化，但由于受到的作用力极小，且作用时间极短，随着外力的消失而消失，岩层的这种随外力消失而恢复原形的形变称为弹性形变。
产生弹性形变的介质叫弹性介质。
在弹性介质内传播的地震波称地震弹性波。
研究地震弹性波可用弹性波理论，如虎克定律等。
（一）各向同性介质和各向异性介质 对某一特定岩层，如果沿不同方向测定的物理性质均相同，称各向同性介质，否则是各向异性介质。
（二）均匀介质、层状介质若介质的弹性性质不仅与测定方向无关，而且与坐标位置无关，就称为均匀介质；如速度v=c （常数）。
非均匀介质中，介质的性质表现出成层性，称这种介质为层状介质；其中每一层是均匀介质；不同介质层的分界处称界面（平面或曲面）；两个界面之间的间隔称为该层的厚度。
将速度v 是空间连续变化函数的介质定义为连续介质。
连续介质是层状介质的一种极限情况。
即当层状介质的层数无限增加，每层厚度无限减小，层状介质就过渡为连续介质，如 v=v 0 (1+βz) 叫线性连续介质。
（四）单相介质和双相介质只考虑单一相态的介质称单相介质，即把组成地层的岩石都视为单一固体相由于岩石往往由两部分组成，一部分是构成岩体的骨架，称基质，另一部分是由各种流体或气体充填的孔隙，由于地震波经过岩石基质和流体孔隙传播的速度不同，因此从波传播来说，这种岩石是由两种相态组成，称这种岩石为双相介质。
二、弹性模量 （一）应力与应变1. 应力：弹性体受力后产生的恢复原来形状的内力称内应力，简称为应力。
应力和外力相抗衡，阻止弹性体的形变。
叠前弹性阻抗反演技术(李国发)
叠前弹性阻抗反演技术李国发刘洋王濮（中国石油大学资源与信息学院信息科学与地球物理系，北京102249）摘要常规的叠后波阻抗反演方法不能全面、准确地提供储层的的岩性和物性信息。
本文介绍了一种新的叠前反演方法——叠前弹性阻抗反演，与常规叠后反演和叠前AVO反演相比，叠前弹性阻抗反演方法能够更好地提供储层岩性和流体性质的信息。
国内外应用实例表明：弹性阻抗反演技术是继叠后反演和叠前A VO反演之后，地震反演技术在地球物理勘探领域的又一重要进展。
关键词声波阻抗，弹性阻抗，叠前反演，弹性参数随着油田勘探和开发工作的不断深入，目前的处理方法、解释结果和储层分析手段已经不能很好地满足精细储层预测地要求。
尽管传统的迭后反演方法在精细储层描述中发挥了重要的作用，但是迭后信息对地震波动力学特性的模糊以及反演结果的多解性限制了迭后反演方法在储层预测中作用。
迭前A VO反演技术虽然在一定程度上利用了地震资料的迭前信息，但是其理论局限以及对地震资料品质的要求使得该方法在实际应用中存在诸多“陷阱”，成功的实例鲜有报道。
特别是目前对迭后反演技术和A VO技术的人为修饰和“改进”，脱离理论、不切实际的渲染和套用，加剧了这两项技术在应用中的混论。
迭后反演技术很好地利用测井信息，但丢失了地震信号的迭前信息；A VO技术利用了地震资料的迭前信息，但是不能很好地利用测井信息。
迭后反演具有较好的稳定性和分辨率，但不能得到泊松比、流体特性等信息；A VO反演技术能够得到有关岩性和流体的信息，但是稳定性和分辨率较低。
随着与角度有关的弹性阻抗概念的提出和建立，迭前弹性阻抗反演技术很好的融合了迭后反演和A VO反演的优点，弹性阻抗中既有纵波信息也有横波信息，进而可以得到纵波速度、横波速度、密度、泊松比、弹性模量等诸多岩石物性参数，这些参数不同于迭后反演方法由波阻抗派生得到的信息，而是通过反演得到的独立信息。
因此迭前弹性阻抗反演方法在很大程度上提高了利用反演技术进行储层预测和储层描述的能力，是地震反演技术新的发展方向。
应用地球物理学原理第二章04弹性波的特征
03
弹性波在地壳中的传播
地壳的分层结构
地壳是地球最外层的硬壳，由 岩石和土壤组成，具有明显的 分层结构。
地球的地壳分为多个板块，板 块之间的相互作用可以产生地 震波。
地壳的分层结构对弹性波的传 播具有重要影响，不同层中的 波速和传播方向可能不同。
弹性波在不同介质中的传播
弹性波在固体、液体和气体中传播时具有不同的特征。
地下结构的不确定性可能导致弹性波传播模型的 误差，从而影响解释结果的准确性。
需要对地下结构进行详细调查和建模，以获得更 准确的弹性波传播特征。
数据处理与解释的复杂性
01
02
03
弹性波数据的处理涉及 多种算法和技术，如滤 波、反演、成像等，处
理过程较为复杂。
弹性波数据的解释需要 丰富的专业知识和经验 ，对解释人员的素质要
应用地球物理学原理第二章 04弹性波的特征
目录
• 弹性波的基本概念 • 弹性波的物理特性 • 弹性波在地壳中的传播 • 弹性波的应用 • 弹性波的局限性
01
弹性波的基本概念
弹性波的定义
弹性波
在弹性介质中传播的波动现象，由于介质的弹性性质，当 受到外力作用时，介质发生形变并产生恢复力，这种恢复 力会以波动的形式在介质中传播。
资源开发规划
通过分析地下岩层的弹性波特征，评 估资源的可开采性和开发风险，为资 源开发提供科学依据。
环境保护监测
利用弹性波技术监测环境变化，如土 壤污染、地下水污染等，为环境保护 提供技术支持。
05
弹性波的局限性
对地下结构的依赖性
弹性波的传播特性与地下结构密切相关，不同的 地下介质对弹性波的传播有显著影响。
弹性波的传播方式
弹性波可以通过反射、折射、散射等方式传播， 其传播路径和速度受到介质的不均匀性和边界条 件的影响。
地震波在分层介质中的反射、透射与...
第1章 绪论1.1 弹性波场论概述在普通物理的力学部分，我们曾经着重讨论过物体在外力作用下的机械运动规律。
在讨论时，由于物体变形影响很小，我们将其忽略，而将物体视为刚体或简化为质点，这是完全正确的。
然而，实际上任何物体在外力作用下不仅会产生机械运动，而且会产生变形。
由于变形物体内部将相互作用，产生内力、应力和应变。
当应力或应变达到一定极限时，物体就会破坏，这一点在研究材料和工程力学中尤其要考虑，地球介质也不例外，地壳运动或地震都会产生地质体的应力或应变。
在弹性力学中，主要讨论对物体作用时的变形效应，物体不再假定为刚体，而是弹性体、塑性体，应当视为可变形体，我们研究的视角也从外部整体过渡到内部局部。
长期的生产实际和科学实验均已表明，几乎所有的物体都具有弹性和塑性。
所谓的弹性是指物体的变形随外力的撤除而完全消失的这种属性。
所谓的塑性是指物体的变形在外力的撤除后仍部分残留的这种属性。
物体的弹性和塑性受诸多因素影响而发生改变，并在一定的条件下相互转化。
因此，确切地，应当说成物体处于弹性状态或塑性状态，而非简单地说物体是弹性体或塑性体。
在弹性力学中，只讨论物体处于弹性状态下的有关力学问题，这时物体可称为弹性体。
由上所述，弹性力学又称弹性理论，研究的对象是弹性体，其任务是研究弹性体在外界因素（包括外力，温度等）作用下的应力、应变和位移规律。
简单地说，弹性力学就是研究弹性体的应力、应变和位移规律的一门学科。
弹性力学是固体力学中很重要的一个分支。
而固体力学是从宏观观点研究固体在外力作用下的力学响应的科学，它主要研究固体由于受外力作用所引起的内力（应力）、变形（应变）以及与变形有直接关系的位移的分布规律及其随时间变化的规律。
可见，应力、应变和位移是空间和时间的函数。
与固体力学对应的还有流体力学等。
固体力学还包括材料力学，断裂力学等等。
弹性力学本身又分为弹性静力学（Elasticity Statics ）和弹性动力学（Elasticity Dynamics ）。
弹性波
斯通利波
在两种不同介质的半空间体的交界面上传播的波称为斯通利波，因斯通利首先发现并研究这种波而得名。它是一种波速与两个介质的性质有关的变态瑞利波。斯通利波的存在与介质的弹性拉梅常数和介质密度有关。在两个介质的拉梅常数λ1、G1和λ2、G2满足λ1/G1=λ2/G2=1的情况下，存在条件如图所示，如果两个介质的密度ρ1和ρ2之比ρ1/ρ2和G1/G2在图示坐标系中对应的点落在曲线A和曲线B之间，斯通利波就存在。在地震学中，理论上已证明斯通利波是存在的，但尚未观测到。
式中为拉普拉斯算符；α和β分别为纵波波速和横波波速;嗞=嗞(x,y,z,t)为标量势；ψx=ψx(x,y,z,t)、ψy=ψy(x,y,z,t)、ψz=ψz(x,y,z,t)为矢量势φ(x,y,z,t)的三个分量。ψx、ψy、ψz统称为波函数，它们和嗞同坐标系中的三个位移分量u、v、w的关系为：
上述波动方程是根据下面的假设导出的：①弹性介质中各质点间的相对位移为无穷小量；②介质是完全线弹性的，即应力和应变之间呈均匀线性关系，服从胡克定律；③介质是各向同性的；④不计外力（如重力、体积力、摩擦力等）。
在精确理论发展的同时，近似解理论也得到发展。有限差分方法先被用于解决短杆中弹性波的传播问题，后被推广到一些复杂结构中波的传播问题。有限元法逐步用于研究弹性波问题，开始用于分析细杆中弹性波的传播，后用于分析各种结构（柱、板、壳体）中的波的传播以及层状介质、正交异性介质中的波的传播等。非线性弹性波的传播问题的研究也取得初步成果。
弹性波理论
地震波交错网格高阶差分数值模拟研究摘要: 地震波数值模拟技术是勘探地球物理学中的重要组成部分，研究通过弹性波一阶速度——应力方程，采用交错网格高阶有限差分法实现了地震波在各向同性介质中的高精度的数值模拟，并采用完全匹配层( PML) 吸收边界来消除边界反射，可取得较好的效果。
通过模型的正演计算和复杂模型的处理结果表明，交错网格高阶有限差分法数值模拟是一种快速有效的地震波数值模拟方法。
关键词: 地震勘探; 交错网格; 有限差分; 数值模拟引言地震数值模拟是模拟地震波在介质中传播的一种数值模拟技术，随着地震波理论在天然地震和地震勘探中的应用，地震模拟技术便应运而生，并随着地震波理论和计算机技术的发展，地震数值模拟技术自20世纪60年代以来也得到了飞速发展，形成了目前具有有限差分法、有限元法、虚谱法和积分方程法等各种数值模拟方法的现代地震数值模拟技术。
有限差分法是偏微分方程的主要数值解法之一。
在各种地震数值模拟方法中，最早出现的数值模拟方法是有限差分法。
Ａlterman和Karal(1968)首先将有限差分法应用于层状介质弹性波传播的数值模拟中。
此后，Boore(1972)又将有限差分法用于非均匀介质地震波传播的模拟。
Alford等（1974）研究了声波方程有限差分法模拟的精确性。
Kelly等（1976）研究了用有限差分法制作人工合成地震记录的方法。
Virieux(1986)提出了应用速度——应力一阶方程交错网格有限差分法模拟P——SV波在非均匀介质中的传播。
交错网格方法提高了地震模拟的精度和稳定性，并消除了部分假想。
有限元法也是偏微分方程的数值解法之一。
Lysmer和Drake（1972）最早将有限元法应用于地震数值模拟。
Marfurt(1984)研究对比了模拟弹性波传播的有限差分法和有限元法的精度。
Seron等（1990,1996）给出了弹性波传播有限元模拟方法。
Padovani等（1994）研究了地震波模拟的低阶和高阶有限元法。
叠前弹性波阻抗反演在储层预测中的应用
４３ 期
朱 波 ，等 ：叠 前 弹 性 波 阻 抗 反 演 在 储 层 预 测 中 的 应 用
４６ ９
图 １ 纵 波 阻 抗 直 方 图 Ｆｉｇ．１ Ｐ－ｗａｖｅｉｍｐｅｄａｎｃｅｈｉｓｔｏｇｒａｍ
其中：犞Ｐ 为纵波阻抗；犞Ｓ 为零 入 射 角 横 波 阻 抗 反 射 系数；ρ为密度；ΙΡ 为纵波阻抗；ΙＳ 为横波阻抗；α为 入射角。利用 Ｆａｔｔｉ方 程 对 多 个 角 道 集 数 据 体 实 施 联 合 井 约 束 稀 疏 脉 冲 反 演 ，可 同 时 得 到 纵 波 阻 抗 、横 波阻抗、密度三 个 数 据 体。 联 合 井 约 束 稀 疏 脉 冲 反 演是对单个叠加数据体井约束稀疏脉冲反演的一种 推广，将地震波数据 体 转 换 成 垂 直 入 射 的 声 波 阻 抗 数据体 。 ［３－４］ 这种转换是通过 核 心 目 标 函 数 的 最 小 化实现的，它包括稀 疏 的垂直反射系 数的 Ｌ１ 模、合 成记录与原始地震道 之 差 的 Ｌ１ 模 和 其 他 约 束 条 件 函数，由反演出的 Ｐ 波阻 抗 和 Ｓ 波 阻 抗 可 以 计 算 出 岩石弹性参数。
（ ） 犚（α）≈
（１＋ｔａｎ２α）１ ２
ΔΙΡ ΙΡ
－８
犞Ｓ 犞Ｐ
２ ２ΔΙΙＳＳｓｉｎ２α＋
［（ ） ］ ４
犞Ｓ 犞Ｐ
２
ｓｉｎ２α－ｔａｎ２α
Δρ ２ρ
（１）
收 稿 日 期 ：２０１９０２２５ 基 金 项 目 ： 中 国 石 油 天 然 气 股 份 有 限 公 司 重 大 科 技 专 项 （２０１６Ｅ－０１） 第一作者：朱波（１９８３－），男，高 级 工 程 师，研 究 方 向 为 地 震 采 集、处 理、解 释 及 综 合 研 究，Ｅｍａｉｌ：ｚｈｕｂｏｑｈ＠ｐｅｔｒｏｃｈｉｎａ．ｃｏｍ．
弹性波基础理论
地震波基础知识
人工地震—由人为活动引起的振动。
工业爆破、地下核爆炸造成的振动
炸药震源—深部勘探
机械震源—浅层勘探
地震波基础知识
弹性介质：若某物体在外力作用下产生形变；当 外力去掉之后，物体能迅速恢复到受力前的形态和 大小，物体的这种性质称为弹性。具有这种性质的 物质，称为弹性介质 塑性：如果外力超过物体的弹性极限，或外力 作用时间太长，当外力消失时，物体不能恢复原 状，物体的这种性质被称为塑性。
5个弹性参数之间的关系：
Vs2 (3V p2 4Vs2 ) E 2 2 V V p s V p2 2Vs2 2(V p2 Vs2 ) Vs2 K (V p2 4 Vs2 ) 3 2 (V 2V 2 ) p s
地震波基础知识
在某一时刻t0开始在介质中激起 波源的振动。过了一段时间，到了 时刻△t(△ t> t0)，波源的振动 可能停止，而t1=t0+△t时的波正 在振动，这个时候，介质中分几个 区域，分界面S上，介质中的各点 刚刚开始振动，这一曲面S叫波在 时刻t的波前；在分界面S’上，介 质中的各质点刚刚停止了振动，这 一曲面S’叫波在△t时刻的波后 （波尾）。
地震波基础知识
P波
P波早于S波到达地面 S波破坏性更大
S波
地震波基础知识
面波 纵波和横波都在介质内部传播，统称为体波。 根据弹性力学理论，还有两种仅仅沿弹性介质表 面传播，离开表面而深入介质内部就会衰减。 一种是沿介质与大气接触的自由表面传播的面波， 称为瑞雷面波 。 另外一种则是沿 两弹性介质之间 的传播的面波 ， 称为勒夫面波。
地震波基础知识
横波 又叫做切变波或 S波。它是由旋转力作用， 弹性介质产生形状形变，这种形变引起的振动称 为横波。该波的传播方向与质点的振动方向相垂 直。质点振动在水平平面中的横波分量称为 SH 波，在垂直平面中的横波分量称为SV波。
光弹性原理及实验方法
§3.1 光测弹性法简介
光测弹性法的定义：
光弹（光测弹性力学）：利用光学方法进行应力分析的全场测试技术。
光测弹性法的基本原理：
采用具有所谓“暂时双折射效应”的高分子材料，如：环氧树脂、聚碳酸脂等 ，制成与实物形状和尺寸几何相似的模型（模型制作）；
并且给模型施加与实物相似的载荷，在特定的偏振光场中观察，模型中出现与 应力有关的干涉条纹（条纹获取）；
电测法回顾
实际测量考虑：
z 温度补偿、防湿 z 栅长的选择 z 粘贴方向的影响 z 横向效应系数的修正 z 长导线的影响和修正
力学信号（应变）
（广义）胡克定律等
应力
z 单向应力 z 平面应力状态
（主应力方向已知） z 平面应力状态 （主应力方向已知）
电阻应变仪 （惠斯登电桥）
电桥的加减特性：
ΔU BD
光程
10
§3.2 光学基本知识
(2) 光程差与相位差
若初始位相相同的两束相干光，其中一束通过折射率为n1，厚度为h的介质。另
一束通过折射率为n2，厚度为h的介质，出介质时两束光的光程差与相位差为：
Δ
=
v
(t1
−
t2
)
=
v
⎛ ⎜ ⎝
h v1
−
h v2
⎞ ⎟ ⎠
=
h
( n1
−
n2
)
位相差： δ = 2π Δ λ
h
E1 n1
E2 n2
Q1:如n1>n2,问（1）离开介质时两束光的光矢量振动方程；（2）那束光位相滞后；（3） 位相滞后多少?
E1
=
a sin [ω(t
−
t1 )
弹性波阻抗反演
数据体中提取声波阻抗； 从大角度数据体中提取 ＥＩ。 从 而 改 进 了 孔 隙 度 和 岩 性 的 预 测 精 度 。
根据反演数据， 发现砂岩储层与某些页岩具有 相同的 ＡＩ 范围， 但 ＥＩ 坐标上有些砂岩的值却始终 比对应的页岩低。这就意味着人们虽然不能够用 ＡＩ 直接从砂岩中区分页岩， 但 如果同时使 用 ＡＩ 和 ＥＩ， 就 可 以 划 出 一 个 趋 势 线 ， 它 有 助 于 从 页 岩 中 区 别 砂 岩 。这 种 综 合 研 究 储 层 的 新 方 法 有 助 于 降 低 钻 井开发的风险和成本。 ２．２ 纵波与横波同步阻抗反演
正的测井曲线构建低频数据分量； ④选择反演的约 束参数， 该反演能分别最佳的匹配油井的阻抗和密 度。在反演算法中， 分别对每一个角度进行子波估 算， 这样能够减少因调谐效应产生的振幅变化。计 算结果产生一个 ＡＩ－ＳＩ 交会图。
同步反演优于单一反演或单一地震振幅提取 的映射。在该油田东部靠 近 ６＃ 井的地方 ， 远、近偏 移距叠加的高振幅区域被认为在这个构造隆起处 中 具 有 较 高 的 气—油—比（ ＧＯＲ） ， 与 此 区 域 相 对 低的 ＡＶＯ 梯度与低Ｎ／Ｇ估算相一致。
ＢＰ Ａｍｏｃｏ 公司的 Ｃｏｎｎｏｌｌｙ 于 １９９９ 年正式发 表了弹性波阻抗（ Ｅｌａｓｔｉｃ Ｉｍｐｅｄａｎｃｅ， 简称 ＥＩ） 反演 方 法 的 经 典 论 文［１］， 标 志 着 反 演 领 域 中 最 重 要 的 物 理意义的提高； ２０００ 年 ＳＥＧ 年会上出现了 ４ 篇 ＥＩ 反演研究的论文， ２００１ 又发表了 １０ 篇； 在 ２００１— ２００２ 年的《Ｔｈｅ Ｌｅａｄｉｎｇ Ｅｄｇｅ》和《Ｇｅｏｐｈｙｓｉｃｓ》中 也 有 多 篇 关 于 弹 性 波 阻 抗 的 论 文 。这 些 论 文 都 认 为
波阻抗和阻抗-概述说明以及解释
波阻抗和阻抗-概述说明以及解释1.引言1.1 概述波阻抗和阻抗是电路中一个非常重要的概念，它们在电子电路设计和分析中扮演着关键的角色。
波阻抗是指在电磁波传输过程中，描述电磁波在不同介质之间传播时发生反射和折射的特性。
而阻抗则是描述电路中对电流的阻碍程度的物理量，是电路分析中一个基本的参数。
本文将首先介绍波阻抗的概念，包括其定义和性质，随后将探讨波阻抗的计算方法，以及在电路中的具体应用。
最后，我们将总结波阻抗的重要性，并展望其在未来的发展方向。
通过对波阻抗和阻抗的深入研究，可以帮助读者更好地理解电路中的电磁波传输过程，以及优化电路设计和性能。
1.2 文章结构:本文将首先介绍波阻抗的概念，包括其定义、特点和作用，帮助读者全面了解波阻抗的基本知识。
接着将详细讨论波阻抗的计算方法，包括静态波阻抗和动态波阻抗的计算公式及应用。
最后，我们将探讨波阻抗在电路中的应用，并结合实例分析波阻抗在不同电路中的作用和影响。
通过这些内容的呈现，读者将深入了解波阻抗在电路领域的重要性以及未来可能的发展方向。
1.3 目的：本文旨在深入探讨波阻抗的概念、计算方法以及在电路中的应用。
通过对波阻抗的研究，我们可以更好地理解电路中的信号传输和阻抗匹配问题，提高电路性能和效率。
同时，我们也将探讨波阻抗在未来的发展趋势，为相关领域的研究提供一定的参考和启示。
通过本文的阐述，读者将能够全面了解波阻抗的重要性，并对其应用和发展方向有更深入的认识。
2.正文2.1 波阻抗的概念波阻抗是电磁波在传播过程中所遇到的阻力，也可以理解为介质对电磁波传播的阻碍程度。
波阻抗与介质的电导率、磁导率以及波的频率有关。
波阻抗通常用符号Z表示，单位为欧姆(Ω)。
在电磁波传播中，波阻抗是一个非常重要的参数，它影响着电磁波在不同介质中的传播速度和能量损失。
波阻抗的大小可以描述介质对电磁波的吸收、反射和透射情况，是电磁波传播过程中的关键因素之一。
波阻抗的概念在电磁学和通信领域中被广泛应用。
弹性力学 弹性波
u 0 (x,t) w 0
从而 e 0
而 2u 0
2
2
x 2
2w 0
代入不计体力的运动微分方程，可见其第一、第三式成为恒 等式，第二式简化为：
2
t 2
c22
2
x 2
E
c2 2(1 )
c2 为横波在弹性体中的传播速度。由于横波的体积应变
24
e 0 ，故横波为等容波。
横波的波动方程的通解为：
(x,t) f1(x c2t) f2(x c2t)
显然，整个通解表示朝相反两个方向传播的两个波，它的 位移沿着y方向，而传播方向是沿着x方向，传播速度等于
常量 c2 。
25
§11-4 球面波
如果弹性体具有圆球形的孔洞或具有圆球形的外表面， 则在圆球形孔洞或圆球形外表面上受到球对称的动力作用时 ，由孔洞向外传播或由外表面向内传播的弹性波，称为球面 波。
球面波是球对称的。利用球对称的基本微分方程：
E(1 ) (1 )(1 2)
(
d 2ur dr2
2 r2
dur dr
2 r2
ur )
kr
0
此时，ur ur (r, t) ，而不计体力时，用径向惯性力
2tu2r代替 kr ，
26
即得：
(1
E(1 )(1
) 2
)
(
2ur r 2
2 r
ur r
17
§11-3 纵波与横波
一、纵波 [定义] 弹性体的质点运动方向平行弹性波的传播方向（图示）
纵波的传播形式
18
将x轴取为波的传播方向，则弹性体内任取一点的位 移分量都有：
u u(x, t) 0 w 0
弹性波阻抗公式精度及参数敏感性分析
弹性波阻抗公式精度及参数敏感性分析冯清源【摘要】从表达形式、理论基础、假设条件出发,系统地介绍了7种纵波弹性阻抗的概念,比较了各弹性阻抗对4类AVO砂岩模型的近似精度,并分析各表达式对背景参数变化的敏感程度.研究结果表明,Verwest弹性波阻抗精度最高,尤其当入射角较大时,且与背景参数无关,稳定性好.【期刊名称】《中州煤炭》【年(卷),期】2019(041)005【总页数】7页(P95-101)【关键词】反射系数;弹性波阻抗;精度分析;敏感性分析【作者】冯清源【作者单位】神华地质勘查有限责任公司,北京 102209【正文语种】中文【中图分类】TD7120 引言传统波阻抗反演基于地震波垂直入射的假设条件，利用叠加道集资料实现弹性参数的求解。
由于地下储层特征的复杂多样，叠后波阻抗反演提供的有限信息，往往造成解释结果的严重不确定性。
作为传统波阻抗反演的延伸和扩展，弹性阻抗反演充分融合常规声阻抗反演和叠前AVO反演的技术特点，利用不同入射角的地震叠加数据提取弹性参数信息，更有利于储层流体的识别和储层物性参数预测[1]。
Connolly[2](1999)首次提出了弹性阻抗EI(Elastic Impedance)的概念，并给出了相应的计算公式。
该方法存在2大缺陷：①反演结果对背景参数K值(K=(VS/VP)2)的依赖性较高，K值的设定是否合适将直接影响实际反演结果的准确程度；②弹性波阻抗随着入射角的变化而变化，且在量纲上与声波阻抗存在差异，不便于两者的对比。
为了克服Connolly表达式量纲不一致的问题，Whitcombe[3](2002)对弹性阻抗进行了归一化处理(Elastic Impedance Normalization)。
Whitcombe等[4](2002)利用学者前期在流体因子和岩性预测方面的研究成果，发现体积模量、剪切模量、纵横波速度比、拉梅常数等弹性参数的变化与Shuey[5]定义的AVO参数A、B和C存在明显的关联性，为了建立各弹性参数变化特征与弹性波阻抗差异之间的对应关系，扩大弹性波阻抗的应用范围，提出了扩展弹性波阻抗EEI(Extended Elastic Impedance)。
几种地震波阻抗反演方法的比较分析与综合应用
Ｅ｛（! ｘ；ｚ０） ｚｉ，ｉ＝１，…，Ｎ｝＝Ｐｒｏｂ｛（ｚ ｘ）≤ｚ０
ｚｉ，ｉ＝１，…，Ｎ｝
（ ４）
通过相应的指示条件期望的估算可以得到相
应的条件概率分布值。根据条件数据，利用指示克
Ａｂｓｔｒ ａｃｔ：Ｓｅｉｓｍｉｃ ｉｎｖｅｒｓｉｏｎ ｉｓ ｔｈｅ ｅｓｓｅｎｔｉａｌ ｆｉｒｓｔ－ｓｔｅｐ ｔｏ ｅｎｓｕｒｅ ｔｈｅ ｐｒｏｐｅｒ ｉｎｔｅｇｒａｔｉｏｎ ｏｆ ｓｅｉｓｍｉｃ ｄａｔａ， ｗｅｌｌ－ｌｏｇｓ，ｖｅｌｏｃｉｔｙ ａｎｄ ｇｅｏｌｏｇｉｃａｌ ｉｎｆｏｒｍａｔｉｏｎ． Ｔｈｅ ａｕｔｈｏｒｓ ｍａｄｅ ｃｏｍｐａｒｉｓｏｎ ａｍｏｎｇ ｓｏｍｅ ｓｅｉｓｍｉｃ ｉｍｐｅｄａｎｃｅ ｉｎｖｅｒｓｉｏｎ ｍｅｔｈｏｄｓ ｕｓｕａｌｌｙ ａｐｐｌｉｅｄ ｂｙ ｏｉｌ ｃｏｍｍｕｎｉｔｙ ｉｎ ｐｒｉｎｃｉｐｌｅｓ，ｋｅｙ ｗａｙｓ，ｗｏｒｋ ｆｌｏｗｓ，ａｄｖａｎｔａｇｅｓ， ｄｉｓａｄｖａｎｔａｇｅｓ ａｎｄ ａｐｐｌｙｉｎｇ ｃｏｎｄｉｔｉｏｎｓ． Ｔｈｅ ｓｅｉｓｍｉｃ ｉｎｖｅｒｓｉｏｎ ｐｒａｃｔｉｃｅ ｉｎｄｉｃａｔｅｓ ｔｈａｔ ｇｏｏｄ ｒｅｓｅｒｖｏｉｒ ｐｒｅｄｉｃｔｉｏｎ ｒｅｓｕｌｔｓ ｃｏｕｌｄ ｂｅ ｇａｉｎｅｄ ｆｏｌｌｏｗｉｎｇ ｔｈｅ ｒｕｌｅｓ ａｎｄ ｗｏｒｋ ｆｌｏｗｓ，ｗｉｔｈ ｒｅｄｕｃｔｉｏｎ ｏｆ ｔｈｅ Ｅ ＆ Ｐ ｒｉｓｋ．
! ! ＯＢＪＦ＝ （（ ｒｉ）ｐ＋"（ｑ ｄｉ－ｓｉ）ｑ）＋#２ （ ｔｉ－Ｚｉ）（２ ３）
式中 ＯＢＪＦ 为目标函数，其物理意义匹配误差， 第一项为反射系数绝对值的和，第二项为地震记录 与合成道之差，第三项为与波阻抗趋势之差平方 和。目标函数最小时，反射系数在测井数据约束下， 合成记录与实际地震数据匹配误差最小。ｒｉ 为反射 系数，λ为地震匹配系数，ｄｉ 为地震数据，ｓｉ 合成地 震记录，α趋势匹配系数，ｔｉ 为用户定义的趋势，主 要是根据测井数据编辑的一个粗略的波阻抗体，Ｚｉ 为用户定义的控制范围内的波阻抗值，ｐ 为反射系 数因子，一般取 ２，ｑ 地震偏差因子，一般取 ２。在迭 代运算过程中，先使用较小的脉冲个数，产生一个 初始模型，然后修改模型，使目标函数达到最小。之 后，不断增加脉冲个数，重复进行迭代，直到反演结 果没有大的改进时就停止迭代，输出反演结果。
交流阻抗的原理与应用
交流阻抗的原理及应用-测聚苯胺修饰电极的电化学性能一、实验目的(1)掌握交流阻抗法（EIS)的实验原理及方法。
(2）了解Nyquist图和Bode图的意义。
（3）学会用Zsimpwin软件对实验数据进行拟合.二、实验原理交流阻抗法（alternating current impedance，AC impedance）阻抗测量原本是电学中研究线性电路网络频率响应特性的一种方法,引用到研究电极过程，成为电化学研究中的一种实验方法。
控制通过电化学系统的电流或电势在小振幅的条件下随时间按正弦规律变化，同时测量相应的系统电势或电流随时间的变化，此时电极系统的频响函数就是电化学阻抗。
通过阻抗可以分析电化学系统的反应机理、计算系统的相关参数。
交流阻抗法是一种以小振幅的正弦波电位(或电流）为扰动信号，益加在外加直流电压上，并作用于电解池，通过测童系统在较宽频率范围的阻抗谱，获得研究体系相关动力学信息及电极界面结构信息的电化学测量方法。
对于一个电解池系统，当在电极两端施加一定电压时,阴阳极会构成一个回路，在这个回路中，电子和离子的传递受到一定的阻力的作用，包括：溶液的阻力，电极的阻力。
而这些阻力正好可以用电阻R进行表征。
再者,在电极和溶液界面上，两相中的剩余电荷会引起静电相互作用，以及电极表面与溶液中的各种粒子（溶剂分子、溶剂化了的离子和分子等）的相互作用.复数阻抗的测量是以复数形式给出电极在一系列频率下的阻抗，不仅能给出阻抗的绝对值，还可给出相位角，可为研究电极提供较丰富的信息。
对于一个纯粹电化学控制的电极体系，可等效成如图2一1所示的电路。
图2一1测试电池的等效电路图2一1中，R e 为溶液电阻，C P 为电极/溶液的双电层电容，R P 为电极电阻。
此等效电路的总阻抗为:2p2p 2222p 2p 2e 1jw -1R C R C R C RP R Z PP ωω+++=其中，实部是2p2p 2pe 1R C R R Z ω++=，虚部是2p2p 2p 2p ,,R C 2ω1R j ωωZ -+=对于每一个w 值,都有相应的Z ’与Z ’’,在复数阻抗平面内表示为一个点连接各w 的阻抗点，得到一条曲线，成为复数阻抗曲线，如图2一2所示。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
假设 K 为常数, 可以将上式写为:
l n E I l n V P 1 t a n 2 l n V S 8 K s i n 2 l n 1 4 K s i n 2 ln V P (1tan2)V S 8K sin2 14K sin2
2020/5/29
2020/5/29
-
弹性阻抗理论
Aki-Richards 等式可写为(三项式) :
R A B s2 i n C s2 itn a 2 n
whe:rA e12VVPP, B1 2 V V PP4V VP S2 V V SS2V VP S2
and: C1VP. 2 VP
Connolly (1999)提出, 类似于声阻抗, 我们定义弹性阻抗(EI) 为:
在谈论弹性阻抗反演之前, 简单回顾反演方法在地震勘探 运用的历史.
2020/5/29
-
地震岩性估算
CDP/CRP 道集
迭加 反演
估算 Z= VP
传统的地震岩性估算是在地震数据叠加偏移后数据(叠后数据)道进行 反演. 这只能估算声阻抗, 不足于推测流体成分.
2020/5/29
ห้องสมุดไป่ตู้
-
道集
AVO 分析
属性 1
2020/5/29
-
Applying Elastic Impedance to the Colony Data
In this exercise, we will start with the 2D Colony dataset which we analyzed earlier and apply the Elastic Impedance analysis to it.
-
然后进行积分和取幂, 可得到以下弹性阻抗 EI等式:
E I() V P (1 ta n 2 )V S ( 8 K s in 2 ) (1 -4 K s in 2)
where K VVSP22
如果 = 0o, EI 弹性阻抗就还原为声阻抗(AI), 等式为:
E0 IA I V P
2020/5/29
-
以上等式使用的是Aki-Richards所有三项等式. 对于入射角大 于300, 该等式将不是线性的拟合. 对于大角度 (大偏移距), 如果只用等式的前两 项,等式可写为:
E I() V P (1 s in 2 )V S ( 8 K s in 2 ) (1 4 K s in 2)
where
弹性阻抗理论
利用Hampson-Russell软件 实现弹性波阻抗反演
介绍
先回顾弹性阻抗反演的原理,然后利用Hampson-Russell 一组软件分析弹性阻抗反演一个实际例子.
弹性阻抗Elastic Impedance (EI) 是由BP 公司Patrick Connolly 提出.
实现弹性阻抗反演可用AVO和STRATA软件自然组合来完 成.
属性 2
估算 VP, VS, 和
叠加
反演
估算 Z= VP
AVO 方法 同时使用多属性估算VP, VS, 和 , 可估算推测流体和
岩性.
2020/5/29
-
AVO分析可能用的属性
但是哪两种(或更多属性)可用来最好地预测流体和岩性这些参数?
目前AVO可提取和使用的一系列不同属性:
偏移距或角度段限(范围限制)叠加 弹性阻抗反演 截距和梯度 提取的纵横波反射系数 RP,RS 由 RP 和 RS 反演到纵横波阻抗 ZP, ZS 纵横波阻抗 ZP, ZS 转换Lambda-Mu-Rho 分析
2020/5/29
-
从范围限制叠加到弹性阻抗
范围限制叠加使用的是常偏移距或常相位道集,方法很成熟并可 避免同相轴未拉平问题. 但它究竟有何含义? Patrick Connolly(BP), 提出新颖的方法来解释范围叠加的含义 , 称之为弹性阻抗. 弹性阻抗是基于Aki-Richards 等式, 下几张图片解释弹性阻抗 的概念和由来.
-
角度叠加进行交绘图
这是近远偏移 距(角度)的交 绘图,高震幅 区块用不同颜 色突出显示.
2020/5/29
-
Top GAS Base GAS Coal
这是交绘图中的突出区块对应在剖面上的位置. 可看出气砂岩已被很 好地描绘.
2020/5/29
-
角度范围叠加
(a)
(b)
(a) 近角度叠加 (0-15o)时间切片, (b)远角度叠加 (15-30o)时间切片(三维河道砂). 对时间切片震幅包络求得是10MS窗口内的平均值并运用了Z-SCORE转换. 异常 河道砂得以清楚显示.
R1E I 1lnEI
2EI 2
2020/5/29
-
假设
K
V V
S P
2
注意到
si2n ta2n ta2n si2n ,
我们可以把 Aki-Richards 等式重写为:
lE n 1 2 I V V P P 1 ta 2 n V V S S 8 K s2 i n 1 4 K s2 i n
K
VVPS
2
当 = 0o, 得到:
E0 IA I V P
2020/5/29
-
弹性阻抗 – 油饱和度的影响
将声阻抗(AI)转换入射角为30° 的弹性阻抗通常曲线类似但阻抗绝对值变低. 角度增加其视声波阻抗降低. 油砂岩比页岩视声波阻抗降低幅度更大.如下图
2020/5/29
Connolly 1999
-
弹性阻抗 – 实例
下图 显示声阻抗和弹性阻抗(30度角)的曲线比较:
2020/5/29
-
重叠显示更好地区分弹性阻抗和声阻抗在含油砂岩处的差别:
2020/5/29
-
EI 反演步骤
道集
AVO 分析
近角度叠加1 远角度叠加2
反演弹性阻抗 反演弹性阻抗
EI(1)
EI(2)
2020/5/29
-
Case Study from the Colony Sand Play in Alberta, Canada
本讨论先使用前面两个属性的方法. 后面再讨论其它的属性方法.
2020/5/29
-
范围限制叠加
道集
AVO 分析
近叠加
远叠加
以上是范围限制叠加的简单流程. 利用常偏移距或常角度的范围叠加已经非常 成熟. 问题是如何解释其结果?
2020/5/29
-
(a)
(b)
2020/5/29
气砂岩范围限制叠加
(a) 近角度 (0o-15o) 叠加 (b) 远角度 (15o-30o)叠加 数据来自加拿大Colony 区 块. 在630MS有一个“亮点” 震幅在远角度 叠加上比近角 度叠加要强的多. 从井资料 分析这是气砂岩导致 “亮点”.