弹性波阻抗理论和实现方法
基于柯西先验分布的叠前弹性阻抗反演方法及应用
基于柯西先验分布的叠前弹性阻抗反演方法及应用郝前勇;宗兆云【摘要】以贝叶斯反演为代表的概率化反演方法既能考虑观测数据的不确定性,又可以考虑待求解参数的先验信息,在实际地震反演中备受青睐.经研究表明,在柯西先验信息下获取的反演结果更具有稀疏性,且具有高分辨特征.叠前弹性阻抗反演是一种基于多角度部份叠加剖面的叠前地震反演方法,信息量丰富,计算效率高.这里在贝叶斯框架下,实现了基于柯西先验的叠前弹性阻抗反演方法,并提取了对储层流体敏感的弹性参数.实际资料应用表明,基于柯西先验的弹性阻抗反演方法合理可靠,具有较高的分辨能力,且提取的弹性参数能够较好地吻合实际钻遇结果.%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 前言概率化反演是求解不适定反演问题的有效方法之一。
高应变
摩阻力为R（x），冲击脉冲峰经过传感器位置的时间
为t1，则脉冲到达x的时间为t1+x/c，则土阻力R（x）
产生的上行压缩波到达传感器位置的时间为t1+ 2x/c，
由前面的分析可知， R（x）到达桩顶后使得桩顶测
点处的力增加R(x)/2,而速度则减少R (x)/（2Z)，这样，
桩侧土阻力R (x) 使得桩顶的力与速度和刚度的乘积
V1dቤተ መጻሕፍቲ ባይዱV1u V 2d V 2u
由于Pu=-Z·Vu，Pd=Z·Vd ，则有
（5-24）
P1u Z 2 Z1 P1d 2Z1 P2u
Z2 Z1
Z2 Z1
（5-25）
P2d 2Z 2 P1d Z1 Z 2 P2u
Z2 Z1
Z2 Z1
（5-26）
假定介质是均质（弹模E和波速C不变）的，则 波阻抗Z（EA/C）的变化是由截面积A的变化引 起的，上式变为：
身轴向应变和桩身加速度的时程曲线，通过适当的计
算得到桩顶的力及速度的时程曲线，通过对桩土进行
单元划分，设定每一单元的桩土参数，以一条实测曲 线作为输入的边界条件，用波动方程求解桩顶的另外
• 一条曲线，通过不断地调整桩土单元的各种参数（阻 力，阻尼，刚度等），使计算的曲线与实测的曲线充
分的吻合，而最后得到的各单元的桩土参数可以认为
的差值F-ZV增加了一个与R(x)相同大小的值。由于土
阻力一经激发后便维持不变，所以t1+2x/c时刻桩顶测
• 的快捷、廉价相比，高应变法检测桩身完整性虽然是附带性的， 但由于其激励能量和检测有效深度大的优点，特别在判定桩身 水平整合型缝隙、预制
• 桩接头等缺陷时，能够在查明这些“缺陷”是否影响竖向抗压 承载力的基础上，合理判定缺陷程度。因此可以作为低应变法 检测这类缺陷桩的一种补充验证手段。
地震波阻抗反演方法综述
地震波阻抗反演方法综述一、地震反演技术研究现状地震反演方法是一门综合运用数学、物理、计算机科学等学科发展起来的新技术新方法，每当数学方法、物理理论有了新的认识和发展时，就会有新的地震反演技术、方法的提出。
随着计算机技术的不断发展、硬件设施的不断升级，这些方法技术得到了实践验证和提升，反过来地震反演技术运用中出现的新问题、新思路又不断促使数学方法、地球物理学理论的再次发展。
时至今日，地震反演技术仍然是一个不断发展、不断成熟、不断丰富着的领域。
反演是正演的逆过程，在地震勘探中正演是已知地下的地质构造情况、岩性物性分布情况，根据地震波传播规律和适当的数学计算方法模拟地震波在地下传播以及接收地震波传输到地表信息的过程。
地球物理反演就是使用已知的地震波传播规律和计算方法，将地表接收到的地震数据通过逆向运算，预测地下构造情况、岩性物性分布情况的过程。
地震波阻抗正演是对反演的理论基础和实现手段。
1959年美国人Edwin Laurentine Drake在宾夕法尼亚州开凿的第一口钻井揭开了世界石油工业的序幕。
从刚开始的查看地质露头、寻找构造高点寻找石油，到通过地震剖面的亮点技术寻找石油，再到现在运用多种科学技术手段进行油气资源的预测，石油勘探经历了一个飞速的发展历程。
声波阻抗（AI）是介质密度和波在介质中传播速度的乘积，它能够反映地下地质的岩性信息。
声波阻抗反演技术是20世纪70年代加拿大Roy Lindseth博士提出的，通过反演能够将反映地层界面信息的地震数据变为反映岩性变化的波阻抗（或速度）信息。
由于波阻抗与地下岩石的密度、速度等信息紧密联系，又可以直接与已知地质、钻井测井信息对比，因此广泛应用于储层的预测和油藏描述中，深受石油工作者的喜爱。
70年代后期，从地震道提取声波资料的合成声波技术得到了快速发展，以此为基础发展的基于模型的一维有井波阻抗反演技术，提高了反演结果的可靠性。
进入80年代，Cooke等人将数学中的广义线性方法运用于地震资料反演，提出了广义线性地震反演。
反射波法测桩原理及方法
反射波法测桩基本原理
假定桩为一维弹性杆件，其材料由介质密度ρ确定。如果 用C 表示在杆件中传播的应力波（纵波）波速，则C 与杆件 的弹性模量E 和介质密度ρ的关系用以下公式表示：
C=
E
ρ
用A 表示杆件任一截面的面积，则定义波阻抗为：
Z = ρ CA =
EA C
在应力波作用下，波阻抗综合反映了杆件材料的物理性 质及其几何尺寸的变化。
V I + V R = VT F I + F R = FT
FI Z1
−
FR Z1
=
FT Z2
Z 1 (V I − V R ) = Z 2V T
Байду номын сангаас
联立求解可得：
质点作用力反射系数 质点振动速度反射系数
R R
F
= = = =
FR Z = FI Z V V
R I
2 2
− Z + Z
2 2
1 1 1 1
V
= − = =
RV
VR Z 2 − Z1 = = − VI Z 2 + Z1
桩底支承介质较桩身混凝土软以及桩身缩径、断裂、 离析、夹泥等均会出现波阻抗减小的现象。此时，反 射系数Rv＞0，在速度检测器接收的曲线上，反射波相 位与入射波相位相同。由于RF＜0，则反射波为上行拉 伸波，与入射波（下行压缩波）传播方向一致。 如果桩完全断裂，Z2≈0，发生全反射，质点应力为 0，质点振动速度加倍。应力波在断裂位置以上多次反 射，形成震荡曲线。检测不到断裂位置以下信息。
反射波时域曲线实例
钢筋混凝土灌注桩： 旋挖成孔，灌注C35 商品混凝土。 上图为重锤敲击， 波形呈现低频大振 幅衰减振动，表明 桩头浅部有严重缺 陷； 下图为轻锤敲击， 缺陷位置明显。
地球物理前沿读书报告
浙江大学地球科学学院地球物理前沿读书报告弹性波理论与岩石物理学[作者姓名]目录摘要 (1)第1章基本的对应关系[1] (2)1.2 骨架颗粒特征、分选及充填型式与速度 (2)1.3 岩石密度与速度 (2)1.4 孔隙度与速度 (2)1.5 孔隙流体和波阻抗 (3)第2章以QHD326油田为例[2] (4)2.1 纵横波速度（泊松比）与砂岩孔隙含流体类型之间的关系 (4)2.2 拉梅常数λ和密度ρ乘积与砂岩孔隙含流体类型之间的关系 (4)2.2 小结 (4)第3章以辽东湾凹陷为例[3] (5)3.1 岩性与纵横波速度的关系 (5)3.2 纵横波速比和波阻抗与岩石物性的关系 (5)结论与认识 (6)弹性波理论和岩石物理之间的关系摘要岩石物理学是一个多学科交叉、融合的学科，包含了物理学、地球物理学、测井、地质学等学科；它是从理论和实验上研究岩石的物理特性，搞清流体类型、流体饱和度、岩比、孔隙率、裂隙等各种因素对岩石电磁波、弹性波的影响，以及它们在地球物理和岩石物理测井数据中的响应。
通过综合研究地球物理、岩石物理测井等资料，找出这些响应与储集层之间的关系。
从其定义中我们可以知道，既然岩石的物理特性能够影响弹性波，那么反过来，弹性波的不同特征就可以作为其对地下不同物理性质的岩石的响应。
因此，在弹性波理论和岩石物理学之间存在某种联系，具体的讲就是弹性波对地下各个界面或对不同岩体的响应特征与地下分界面或者岩石物理性质之间的联系，这种联系可能是一对多、多对一或者多对多。
实现的基本思路是利用各种地震弹性参数的空间变化预测和描述地下储集岩石某种物理参数的变化，如岩石骨架、孔隙、流体甚至温度、压力等，进而实现岩性预测和油气检测的目标。
根据上述分析，显对地下地质构造和岩性探测可以用弹性波（人工地震），这也就是在地球物理中弹性波理论的重要性所在。
第1章基本的对应关系[1]1.2 骨架颗粒特征、分选及充填型式与速度，当组成岩石的颗粒分选变差时，颗粒间排列充填得更加紧密，一方面会减小粒间孔隙空间，另一方面使得骨架变得更加坚硬，导致岩石速度增加。
0-叠后反演-Jason-讲
四、Jason6.2平台反演模块的构成
Invermod
Invertrace
Velmod Environment Plus s
Statmod
（地震反演）
Functionmod
• Recursive trace inversion (RTI 递归反演) • Constrained Sparse Spike Inversion (CSSI 约束稀疏脉冲反 演) • Trace merging (AITM 道合并) • Netpay and porosity estimation (NPPE 有效厚度和孔隙度估 算) 依据绝对阻抗与储集岩和孔隙度的关系计算有效厚度和孔隙度。
三、反演方法分类
1、基于原始地震资料 叠前反演类、叠后反演类 2、基于数学算法 线性反演类、非线性反演类 3、基于实现方式 基于反褶积的反演方法：包括基于地层反褶积的道积分、 递归、广义线性等反演方法； 基于最大似然反褶积的稀疏脉冲反演方法; 基于最大后验概率准则算法反褶积的模拟退火等。 基于波动方程的反演：Born反散射。 基于随机过程的反演：随机反演、随机模拟等。 基于特征分析的反演：特征反演、神经网络反演等。 其它：混沌反演等
• Principal component analysis （主因子分析） • Model estimation （模型估算） • Model generator （模型生成器）
四、Jason6.2平台反演模块的构成
Invermod
Rocktrace
Velmod Environment Plus Earthmodel
Velmod
Environment Plus Earthmodel
Statmod
Functionmod
三种弹性波阻抗公式比较
三种弹性波阻抗公式比较
张奎;倪逸
【期刊名称】《石油地球物理勘探》
【年(卷),期】2006(041)0z1
【摘要】弹性波阻抗反演是一种结合了常规声阻抗反演和叠前AVO反演的地震反演技术.本文通过对不同的模型进行试算,分析比较了Connolly,whitcombe等和Verwest等提出的三种弹性波阻抗公式.研究结果表明,由verwest EI公式计算得到的反射系数比其他两种EI公式求得的反射系数精度高.与其他两种弹性波阻抗相比,Verwest EI可以作为更好的岩性指示器指导岩性和流体预测.
【总页数】4页(P7-10)
【作者】张奎;倪逸
【作者单位】东方地球物理公司物探技术研究中心;东方地球物理公司物探技术研究中心
【正文语种】中文
【中图分类】P61
【相关文献】
1.三种弹性波阻抗公式比较 [J], 张奎;倪逸
2.半无限弹性体动态位移的计算公式及其与布希涅斯克公式的比较 [J], 徐玉华;杨维汉
3.弹性波阻抗公式精度及参数敏感性分析 [J], 冯清源
4.弹性波阻抗公式精度及参数敏感性分析 [J], 冯清源
5.Whitcombe扩展弹性波阻抗公式的改进 [J], 苑春方;彭苏萍;吕焱;王晓伟;蔡利文
因版权原因，仅展示原文概要，查看原文内容请购买。
基于Fatti近似的弹性阻抗方程及反演
Ｆ ｔ 近 似 更 准 确 的反 射 系数 ， 后 对 这 个公 式进 行 了标 准 化 以 实 现 不 同 角 度 的 弹性 阻抗 间 量 纲 的 统 一 ， 后 用 标 准 ａｔ ｉ 然 最 化 后 的 公 式 进 行 了反 演 ， 反 演得 到 的 不 同 角度 的 弹 性 阻抗 数 据 体 中 可 直 接 提 取 得 到 纵 横 波 阻 抗 数 据 体 ． 用 实 例 从 应
王保丽 ， 印兴耀 ， 张繁昌 ， 李爱山
（ ． 国 石 油 大 学 （ 东 ） 球 资 源 与 信 息 学 院 ， 营 ２７ ６ ； ２ 中海 石 油 研 究 中心 国际 研 究 部 ， 京 １０ ２ ） １中 华 地 东 ５０ １ ． 北 ００ ７
摘
要 用 Ｃ ｎ ｏｌ 弹 性 阻 抗 （ Ｉ ｅ ｓｉ ｉ ｅ ａ ｃ ） 式 进 行 反 演 只 能 直接 得 到 纵 、 波 速 度 和 密度 的信 息 ， 后 ｏ ｎｌ ｙ的 Ｅ ，ｌ ｔ ａ ｃ ｍｐ ｄ ｎ ｅ 公 横 然
才 可 间接 计 算 得 到 纵 横 波 阻抗 等 其 它 的 参 数 数 据 体 ， 样 便 增 加 了一 步 误 差使 数 据 的 准 确 性 降低 ． 文 首 先 针 对 该 这 本
方 法 的 这 些 缺 点 和 不足 ， 出 了以 Ｚ ｅ ｐ ｉｚ方 程 的 Ｆａｔ 近 似 为 基 础 一 种 新 的 弹性 阻抗 公 式 ， 该 公 式 可 得 到 比 用 提 ｏ ｐ ｒｔ ｔｉ 由
ａ ｒ ｘｍ ａ ｉｎ ａ ｄ ｉｖ ｒ ｉｎ ｐｐ ｏ ｉ ｔｏ ｎ ｎ ｅ ｓｏ
Ｗ ＡＮＧ ｏ ｌ ＹＩ Ｘｉ — ａ ， ＺＨＡＮＧ ｎ ｃ ａ ， ＬｉＡ ｉｓ ｎ Ｂａ — ｉ ， Ｎ ｎｇ ｙ ｏ Ｆａ — ｈ ｎｇ — ｈａ