波阻抗剖面技术是地震岩性勘探中的重要方法之一。它是将地震资料与测井资用线性正、反演技术求取地下地层的反射系数、波阻抗、拟速度等重要于波阻抗剖面可以提供众多地下信息，故该方法已被地球物理学界所注自建的PDP－11－45机的波阻抗剖面系列，在合成记录的研究中，提供的H一是一个成功的方法，达到了引进的CYBER一720所做的合成记录水平；在常型范畴内实现了在时间域中求取反射系数、相对波阻和绝对波阻的方法H－T0转换方法和最佳反褶积算子求取法的成功，保证了地层反褶积取其误差约为1％。这些试验结果，对实际处理是有意义的。为找油气、和岩石性质，有较大的实用价值。 2100433B

地震波阻抗反演是储层预测的有效手段之一，波阻抗与含油气储层有很好的对应性。

波阻抗反演实际上是从地震剖面上消除子波影响,留下反射系数,再由反射系数计算出能反映地层物性变化的物理参数波阻抗。

波阻抗反演的方法很多，其中最常用的是广义线性反演，该方法有较高的反演精度，但是在反演过程中，高频噪声对解病态方程组影响很大,即使加了正则化因子,也无法消除高频噪声的影响。

根据反射系数序列，并给定波阻抗初始值，利用递推波阻抗的反演方法就可以算出相应的波阻抗序列。该方法的核心是反演反子波，但是反演出反子波局部化特征并不理想，因为两翼都有振荡，并不是标准的δ(t)函数，而是带有剩余子波的δ(t)函数，其频率过低，计算出的反射系数频率也较低，相当于带有剩余子波的反射系数。递推波阻抗的反演方法对误差及初值 Z0 均十分敏感，稳定性也差。可以用递推波阻抗反演法反演出的波阻抗作为迭代反演波阻抗的初值。

因此，针对广义线性反演，重点讨论了迭代解方程的波阻抗算法、解方程的正则化方法、波阻抗反演中的分辨率、改善解的二次圆滑滤波、用于控制反演质量的统计量和均方差及波阻抗的约束反演等关键技术问题;针对递推波阻抗的反演,重点讨论了递推波阻抗的反演方法、子波和反子波提取、误差对提取反子波的影响等关键技术问题。同时做了大量的模型试算。2100433B

1 波阻抗计算

对有一定预应力的岩石进行循环冲击试验的装置为自行研制的动静组合加载装置。对于有轴压时的冲击过程中，杆件与试样仍然满足装置赖以存在的一维应力波传输理论。应力波从一种介质传播到另一种波阻抗不同的介质时，入射波在2 种介质的界面处会产生反射波和透射波。由于其值比岩石的波阻抗大许多，认为在整个冲击过程中其大小保持恒定；

在整个冲击过程中并不是每一瞬时的反射波或透射波都夹杂多次反射和透射波的影响，当时，测试所得的反射波中没有 2 次和更高次反射的影响。同理，在透射波采样时间内的区段，试验所得的透射波中也没有夹杂多次透射波的影响。因此，可以利用SHPB 试验系统采集的入射波和反射波或透射波在各自采样区段内相对应的值计算岩石的波阻抗。而其他时间区段的波阻抗不可由试验数据直接计算。

2 计算波阻抗时计算式的选择

在实际试验时，反射波和入射波是用入射杆上的同一个应变片测量的，特别是有轴压的冲击试验时，反射波的起跳点受入射波的影响较大，其对应的时间点和大小不易确定，极易引起误判。而透射波不在这种情况，起跳点相对容易判断，因此计算岩石的波阻抗较为理想。

由于岩石与弹性杆间的波阻抗大小不同，每一入射波经过透反射后必有一反射波和透射波与之相对应；冲击过程中的不同时刻岩石所受的冲击动载大小不同，岩石试件内部裂隙的闭合程度或扩展程度亦不相同，即岩石试件具有不同波阻抗。换言之，在一次冲击过程中，入射波和反射波或透射波与岩石的波阻抗一一对应，不同时刻波阻抗的大小很可能不一样，入射应力大小、透射应力大小和时间点3 个变量必须确定其中2 个不变。根据一维应力波理论可知，在异型冲头形状、长度和冲击速度相同的前提下，在区间固定一时间点，则每次冲击得到的理论上应为定值，对照时间确定出对应的透射波的值。具体选其中哪一个时间点视情况而定，它们之间的差别微乎其微。要强调的是对同一岩石试件的循环冲击过程的损伤进行分析时，应选择同一时间点。

另外，利用声波定义循环冲击过程岩石的累积损伤，也是在冲击载荷作用后测量岩石的波速值。本文的方法与此相似，每次测量（除第1 次）的值都是前1 次冲击作用后岩石在压密后的波阻抗值，这个值表征岩石在前1 次冲击后内部的累积损伤，因此，在测量程序上也与声波定义损伤的一致 。