1、叠前地震反演叠前地震反演印兴耀中国石油大学山东东营 反演反演叠前反演叠前反演叠后反演叠后反演EI反演反演AVO反演反演利用叠后资料在井的约束下得到利用叠后资料在井的约束下得到AI利用叠前资料获取利用叠前资料获取AVO属性属性利用叠前资料利用叠前资料及弹性波阻抗获取岩石弹性参数及弹性波阻抗获取岩石弹性参数反演分类反演分类水平叠加水平叠加角道集叠加角道集叠加(8度度)角道集叠加角道集叠加(16度度)角道集叠加角道集叠加(24度度)水平叠加与角度部分叠加水平叠加与角度部分叠加 页岩、含油、气、水砂岩的速度随深度变化关系页岩、含油、气、水砂岩的速度随深度变化关系VpVpVpVp1 1、弹性模量、弹性模
2、量 岩石弹性参数及其相互关系岩石弹性参数及其相互关系 2、纵波速度和横波速度、纵波速度和横波速度3/42KVpsV岩石弹性参数及其相互关系岩石弹性参数及其相互关系 3、泊松比、泊松比222215.0)/(1)/(5.0pspsVVVV岩石弹性参数及其相互关系岩石弹性参数及其相互关系 泊松比的物理意义泊松比的物理意义是流体的体积模量，所以可以叫做流体的不可压缩性度是流体的体积模量，所以可以叫做流体的不可压缩性度量（量（fluid incompressibilityfluid incompressibility）。）。是抵抗剪切应变的度量，它常被称为刚性模量是抵抗剪切应变的度量，它常被称为刚性模量
3、（modulus of rigiditymodulus of rigidity），不可压缩性的度量），不可压缩性的度量（imcompressibilityimcompressibility）或剪切模量（）或剪切模量（shear modulusshear modulus）。(泊松比泊松比)：是反映岩性与含油气性的重要参数，它是用岩石纵向是反映岩性与含油气性的重要参数，它是用岩石纵向拉伸和横向压缩的比值来表示拉伸和横向压缩的比值来表示222215.0)/(1)/(5.0常用弹性参数的意义常用弹性参数的意义n含油气地层一般都具有较低的含油气地层一般都具有较低的和和/n 与流体无关，仅反映岩石骨架的性
4、质与流体无关，仅反映岩石骨架的性质n和和能够灵敏地反映储层属性，并且物理能够灵敏地反映储层属性，并且物理意义明确，在岩性和流体预测方面具有重大意义明确，在岩性和流体预测方面具有重大的意义的意义和和n描述岩石弹性特征的主要地震参数有岩石的弹性模描述岩石弹性特征的主要地震参数有岩石的弹性模量量、密度密度 、纵波速度、纵波速度 和横波速度和横波速度 等。弹性模等。弹性模量反映岩石应力量反映岩石应力应变关系的特征，密度反映岩石应变关系的特征，密度反映岩石的比重，速度则反映地震波在岩石中的传播的特征的比重，速度则反映地震波在岩石中的传播的特征，它是弹性模量的函数，它是弹性模量的函数。)3/42K拉梅常数
5、拉梅常数 剪切模量剪切模量 体积模量，体积模量，k222215.0)/(1)/(5.02.332.843.353.854.364.872.082.593.093.604.114.610.820.991.161.341.511.681.852.032.202.372.54ZvpZvpsvZsvZZvs/Zvp/SH CrossplotDepth Interval:2517.00-3527.00Date:Sun Dec 21 21:02:27 2003 Z Axis:SH0.0010.0020.0030.0040.0050.0060.0070.0080.0090.00100.00Wells:CB3
6、2123.4029.1334.8540.5846.3020.5426.2731.9937.7143.4449.167.2110.5713.9417.3020.6724.0327.4030.7734.1337.5040.86ZlamdarhoZlamdarhoohruimZohruimZZmiurho/Zlamdarho/SH CrossplotDepth Interval:2517.00-3527.00Date:Sun Dec 21 21:02:16 2003 Z Axis:SH0.0010.0020.0030.0040.0050.0060.0070.0080.0090.00100.00Wel
7、ls:CB321VpVsLamdaRhoMiuRho井的砂泥岩纵横波速度交汇图井的砂泥岩纵横波速度交汇图LamdaRho和和MiuRho交汇图交汇图砂砂岩岩泥泥岩岩砂砂岩岩泥泥岩岩0.4 0.8 Vp Vs交会图交会图 交会图交会图/交汇图在在 /交会图中交会图中可以看出可以看出，油层（红色）基本都分布，油层（红色）基本都分布在低泊松比、低在低泊松比、低/的区域。说明该区油层的特点之一是比水的区域。说明该区油层的特点之一是比水层具有更低的泊松比。当泊松比大于层具有更低的泊松比。当泊松比大于0.32时，砂岩为水层。时，砂岩为水层。水层油水同层油层砂岩 泥岩Nml 砂泥岩交会图砂岩 泥岩Nml 砂
8、泥岩交会图砂岩 泥岩Nml 砂泥岩交会图砂岩 泥岩亮点技术的出现，使地震烃类检测能力有了很大提高，但实践表明，该技术也有其局限性，并不是所有的亮点都与气层有关，某些特殊的岩性体也可以形成地震剖面上的“亮点”八十年代初，Ostrander将反射系数随入射角变化应用于亮点型含气砂岩的识别，发现含气砂岩反射振幅随偏移距增加而增大，含水砂岩反射振幅随偏移距增加而减小，这一发现极大地改善了烃类检测的能力AVO分析是一项利用振幅随偏移距变化特征来分析岩性和油气藏的地震勘探技术，其基本思想是振幅随入射角变化(AVA)的规律与岩性参数有关。AVO也有其陷阱，并不是所有气层的AVO特征都呈增加现象，特殊岩性体也
9、会引起AVO增加现象n充分利用地震反射振幅随偏移距不同或入射角不同而充分利用地震反射振幅随偏移距不同或入射角不同而变化的特征，通过对这些数据的纵横波反演技术研究，变化的特征，通过对这些数据的纵横波反演技术研究，能得到能得到纵横波波阻抗，纵横波速度，密度，泊松比纵横波波阻抗，纵横波速度，密度，泊松比等多等多种参数；提供了研究岩性和流体变化规律的更多、更有种参数；提供了研究岩性和流体变化规律的更多、更有效的反演数据体成果效的反演数据体成果n叠前地震反演较叠后反演推进一步叠前地震反演较叠后反演推进一步叠前反演v 从含有丰富地质信息的叠前资料入手，深入研究从含有丰富地质信息的叠前资料入手，深入研究叠前
10、记录的正演模拟、叠前波形反演技术，采用弹性叠前记录的正演模拟、叠前波形反演技术，采用弹性波反演方法由叠前道集中获取地层纵横波速度、密度波反演方法由叠前道集中获取地层纵横波速度、密度等参数；然后利用不同角度的角道集数据得到弹性波等参数；然后利用不同角度的角道集数据得到弹性波阻抗数据体。这些技术相互补充，使叠前地震技术的阻抗数据体。这些技术相互补充，使叠前地震技术的研究进一步得到深入研究进一步得到深入v叠前弹性波反演和叠前弹性波反演和EIEI有机结合，可以建立一套较为有机结合，可以建立一套较为完善的联合反演技术，从中得到的地层速度、密度、完善的联合反演技术，从中得到的地层速度、密度、泊松比、弹性模
11、量等岩性参数数据体，将为研究储层泊松比、弹性模量等岩性参数数据体，将为研究储层的变化规律提供新的手段的变化规律提供新的手段111,222,P垂直入射反射垂直入射反射2221111212,ZZZZZZRPPZoeppritz方程 CBA2112221122111121211222212112221111221122112sin2cos2sin2cos2cos2sin2cos2sinsincossincoscossincossinPSPPPSSPSPSPSSPVVVVVVVVVVVVVVVAPSPPPSPPTTRRB11112cos2sincossinC0R2R1TRPS112211PPPP112
12、21122PPPSPPPPPPPTVVVVVVV当地震波垂直入射到界面上时有o10按斯奈尔定理有：o21210由Zoeppritz方程解得简单模型合成CMP道集，共60道，最大炮检距1500米近偏移距、中偏移距和远偏移距各20道的叠加结果反射2缺失(第二和第三层阻抗相等)反射3微弱某区典型含油砂岩的地震反射道集某区典型含油砂岩的地震反射道集某区典型含水砂岩的地震反射道集某区典型含水砂岩的地震反射道集泥岩泥岩含油砂岩含油砂岩泥岩泥岩含水砂岩含水砂岩含油性研究含油性研究Vp、Vs、叠前弹性波叠前弹性波 波形反演波形反演 弹性参数弹性参数叠前地震叠前地震数据数据角道集数据角道集数据弹性波阻抗弹性波阻
13、抗约约 束束纵横波速度、密度、泊松比、纵横波速度、密度、泊松比、弹性模量等地层弹性参数弹性模量等地层弹性参数叠前弹性波叠前弹性波波阻抗反演波阻抗反演叠前波动方程反演算法基本流程算法基本流程随机产生一族弹性地层模型计算每一地层模型的合成地震道集合成道集与实际观测数据的比较遗传算法修改地层模型收敛？否是计算适应度值保存结果保存结果结束结束叠前弹性波反演叠前弹性波反演叠前弹性波反演 叠前地震反演能够从地震数据中得到叠前地震反演能够从地震数据中得到P P波和波和S S波信息波信息。由于。由于P P波对孔隙流体的变化较敏感，而波对孔隙流体的变化较敏感，而S S波主要与岩石波主要与岩石骨架有关，流体的变化
14、对其影响较小。为了检测流体类骨架有关，流体的变化对其影响较小。为了检测流体类型，需要同时提供岩石的型，需要同时提供岩石的P P波和波和S S波信息波信息在叠前波形反演的文献中，绝大多数是平面波反演，在叠前波形反演的文献中，绝大多数是平面波反演，而实际地震记录是球面波记录而实际地震记录是球面波记录基于声波理论的反演，得不到横波速度基于声波理论的反演，得不到横波速度局部线性化反演方法需要地层参数的导数构成雅可比局部线性化反演方法需要地层参数的导数构成雅可比矩阵，而对地下进行精确模拟需要几百层地层，求导非矩阵，而对地下进行精确模拟需要几百层地层，求导非常困难常困难Why?叠前弹性波反演叠前弹性波反演
15、从左向右分别为初始群体最佳个体的P波速度、S波速度、密度及其叠前合成记录。图中红色曲线为真实井曲线，蓝色曲线为初始估计模型，可见初始模型参数与井数据差别较大。图中中间的黑色细线是地层参数的低频趋势，两侧的黑线是参数的约束范围。叠前弹性波反演叠前弹性波反演初始群体最佳个体的叠前合成记录(b)与井曲线叠前道集(a)的对比，可见二者存在很大差别。(a)(b)从左向右分别为反演的P波速度、S波速度、密度及其合成记录。图中红色曲线为真实井数据，蓝色曲线为反演得到的地层参数，可见反演结果与真实模型有很好的吻合。叠前弹性波反演叠前弹性波反演从残差剖面上看出，两个地震道集有很好的相似性叠前弹性波反演叠前弹性波
16、反演叠前弹性波反演叠前弹性波反演时间(s)P波速度 S波速度密度 Vp Vs 密度（蓝色曲线为实际井数据，红色曲线为叠前波形反演结果，纵轴表示深度，单位：米）v从实际井曲线与反演结果的对比来看，P-波速度和S-波速度匹配的很好，密度曲线匹配的较好v密度曲线能较好地反映地层的变化，含气砂岩的密度较低v与泥岩相比，砂岩表现高P-波速度和高S-波速度。含气后，P-波速度降低，而S-波速度基本不变。充分表明P-波速度的降低是含油气所致而不是由岩性变化引起密度曲线(深度域)合成道集(左)与实际地震记录(右)对LS36-1-1井附近的CDP道集进行了遗传算法反演，下图分别是合成道集与实际地震记录的比较，箭
17、头所示的同相轴同样显示了气砂岩反射引起的振幅随偏移距增大现象反演结果与井曲线的对比，左图是纵波速度曲线，右图是密度曲线；其中红色是反演曲线，兰色是实际井曲线 经过对比认为，反演结果是可靠的，从而为弹性波阻抗反演提供了先验信息。合成的叠前记录合成的叠前记录LS36-1-2LS36-1-2井附近井附近的实际的实际CMPCMP道集道集 根据根据AVO理论，零炮检距（或小炮检距）剖面理论，零炮检距（或小炮检距）剖面可以近似为声阻抗可以近似为声阻抗AI（Acoustic Impedance）的函）的函数，它与岩石的密度和纵波速度有关。为了充分利数，它与岩石的密度和纵波速度有关。为了充分利用大炮检距地震振
18、幅信息，需要利用一种用大炮检距地震振幅信息，需要利用一种与入射角与入射角有关的弹性阻抗有关的弹性阻抗EI（Elastic Impedance），），这种弹这种弹性阻抗形式上性阻抗形式上类似于声阻抗类似于声阻抗AI，不但与岩石的密度不但与岩石的密度和纵波速度，而且与和纵波速度，而且与横波速度、入射角横波速度、入射角有关有关。弹性波阻抗弹性波阻抗(EI)nEI最早出现在最早出现在Connolly于于1999年发表的论文中，随后在年发表的论文中，随后在2000年的年的SEG年会上同时出现了年会上同时出现了4篇论文对篇论文对EI进行了研进行了研究，究，ARCO公司介织了他们申请了专利的弹性波阻抗公司介
19、织了他们申请了专利的弹性波阻抗反演方法，认为在求取的反射系数的稳定性方面要好反演方法，认为在求取的反射系数的稳定性方面要好于于Conno11y方法，而且计算的方法，而且计算的EI和和AI数值在同一个尺数值在同一个尺度下：度下：CGG公司发表论文认为公司发表论文认为EI反演效果要好于反演效果要好于AV0反反演，因为演，因为EI在抗噪能力方面比在抗噪能力方面比AV0方法具有优势：同时方法具有优势：同时BP Amoco公司在会上又提出了扩充弹性波阻抗方法，公司在会上又提出了扩充弹性波阻抗方法，可以用于流体和岩性预测可以用于流体和岩性预测n这些进展说明：弹性波阻抗已经成为波阻抗反演进一这些进展说明：弹
20、性波阻抗已经成为波阻抗反演进一步发展的方向之一，地震反演的发展正在走向步发展的方向之一，地震反演的发展正在走向AI和和EI相结合、相结合、AI和和AVO相结合的道路。相结合的道路。ZoeppritzZoeppritz方程及其近似方程及其近似P-P、P-S 反射和透射反射和透射112PPSVSVP111,222,2ZoeppritzZoeppritz方程及其近似方程及其近似Aki,Richards和Frasier近似cbaR)(第一项包含P-波速度第二项包含密度第三项包含S波速度Shuey近似ZoeppritzZoeppritz方程及其近似方程及其近似)sin(tan2sin)1()(22220
21、ARRRPP31ZoeppritzZoeppritz方程及其近似方程及其近似Smith和Gidlow近似2222tan21sin2221)(Ru用这一近似在地震道集上做加权叠加可得到储层岩石性质的信息2sinGPRPP其中P 垂直入射时的纵波反射系数G 反射系数函数的梯度AVOZoeppritzZoeppritz方程及其近似方程及其近似从中可以得到许多属性剖面依据的依据的Zeoppritz近似方程近似方程SSPSPPPSPVVVVVVVVR2222222sin42secsin4121)(Aki_Richards(1980)近似近似用的就是该近似方程弹性波阻抗反演处理利，本次和、够得到利用该近似
22、通过反演能222222sin2tan21sin4tan121)(PSSSPSPPPVVIIVVIIR222222sec4121sin2sec21sec2141)(PSPSPVVVVRFatti(1994)近似近似 Gray(1999)近似近似 和、够得到利用该近似通过反演能SPII和、够得到利用该近似通过反演能 弹性波阻抗弹性波阻抗Elastic Impedance地震波垂直入射时的褶积模型，对于与角度相关的数据，有：地震波垂直入射时的褶积模型，对于与角度相关的数据，有：S()S()R()R()*W()W()S()S()为角度地震道，为角度地震道，R()R()是与角度有关的反射系数，它是由是与
23、角度有关的反射系数，它是由P P波速度、波速度、S S波速度和密度，根据波速度和密度，根据ZoeppritzZoeppritz方程计算得到的；方程计算得到的；W()W()是角度子波，由角度地震数据得到。就像对反射系数进是角度子波，由角度地震数据得到。就像对反射系数进行积分而得到声阻抗行积分而得到声阻抗(Acoustic Impedance)Acoustic Impedance)一样，角度反射系一样，角度反射系数用来计算弹性波阻抗。数用来计算弹性波阻抗。EIEI并不是一个可以进行物理测量的属性，它是由并不是一个可以进行物理测量的属性，它是由ZoeppritzZoeppritz方方程或其近似式得出
24、的用来解释地震数据的属性，弹性波阻抗反程或其近似式得出的用来解释地震数据的属性，弹性波阻抗反演是目前为解决岩性解释而提出的比较适用的方法演是目前为解决岩性解释而提出的比较适用的方法 。寻求一个函数，使其与声阻抗有相似的性质，以便于寻求一个函数，使其与声阻抗有相似的性质，以便于以任意角入射时反射系数均可以用下面的公式表示：以任意角入射时反射系数均可以用下面的公式表示：)()()()()(11iiiitftftftfR函数函数 称作称作EIEI（弹性波阻抗），当地层弹性参数变化（弹性波阻抗），当地层弹性参数变化不大时，用不大时，用EIEI的对数值来表示反射系数是准确的的对数值来表示反射系数是准确的
25、 )(tf)ln(2121)(EIEIEIR 去掉两边的微分和对数项，得到下式：去掉两边的微分和对数项，得到下式：)sin41(sin8)tan1(222KKSPVVEI该式即为该式即为P波的弹性波阻抗波的弹性波阻抗 弹性波阻抗反演流程弹性波阻抗反演流程n用测井曲线合成横波速度，进而计算井中的用测井曲线合成横波速度，进而计算井中的EI曲曲线，对于同一个线，对于同一个CMP道集来说，道集来说，EI是入射角的函数是入射角的函数n 将偏移距数据转化为角道集数据，将偏移距数据转化为角道集数据，估算每个角度估算每个角度数据所对应的子波数据所对应的子波n 利用与迭后反演相类似的反演算法计算出各个角利用与迭
26、后反演相类似的反演算法计算出各个角度的弹性波阻抗剖面度的弹性波阻抗剖面 n波阻抗反演波阻抗反演n弹性波阻抗反演弹性波阻抗反演子波提取子波提取 在不同角度的道集上进行合成记录标定在不同角度的道集上进行合成记录标定5-15度的度的y551井的井的合成记录层位标合成记录层位标定图，从左至右定图，从左至右分别为：角度子分别为：角度子波、地震道、合波、地震道、合成记录和测井曲成记录和测井曲线。右边的曲线线。右边的曲线分别为：红色的分别为：红色的波阻抗曲线、蓝波阻抗曲线、蓝色的纵波速度曲色的纵波速度曲线线5-15度/y551井5-15度合成记录标定度合成记录标定15-25度的度的y551井井的合成记录层位
27、的合成记录层位标定图，从左至标定图，从左至右分别为：角度右分别为：角度子波、地震道、子波、地震道、合成记录和测井合成记录和测井曲线。右边的曲曲线。右边的曲线分别为：红色线分别为：红色的波阻抗曲线、的波阻抗曲线、蓝色的纵波速度蓝色的纵波速度曲线曲线15-25度/y551井15-25度合成记录标定度合成记录标定25-35度的度的y551井井的合成记录层位的合成记录层位标定图，从左至标定图，从左至右分别为：角度右分别为：角度子波、地震道、子波、地震道、合成记录和测井合成记录和测井曲线。右边的曲曲线。右边的曲线分别为：红色线分别为：红色的波阻抗曲线、的波阻抗曲线、蓝色的纵波速度蓝色的纵波速度曲线曲线2
28、5-35度/y551井25-35度合成记录标定度合成记录标定三、弹性阻抗反演三、弹性阻抗反演n角道集资料角道集资料5-15度角道集剖面度角道集剖面 Ang5-15度15-25度角道集剖面度角道集剖面 Ang15-25度度25-35度角道集剖面度角道集剖面 Ang25-35度度四、弹性阻抗反演四、弹性阻抗反演Ang5-15度 “AVOAVO反演反演”，是从地震数据中抽取，是从地震数据中抽取AVOAVO属性，然后在合适的测井曲属性，然后在合适的测井曲线的约束下将其转化。线的约束下将其转化。“弹性波阻抗弹性波阻抗”反演，首先是在模拟的测反演，首先是在模拟的测井曲线约束下转化一定角度范围的角道集，再抽
29、取所需的弹性参井曲线约束下转化一定角度范围的角道集，再抽取所需的弹性参数。这两种方法只在子波不随偏移距变化时是相同的，当子波由数。这两种方法只在子波不随偏移距变化时是相同的，当子波由于噪音等因素发生变化时，会在不同的于噪音等因素发生变化时，会在不同的AVO属性之间产生属性之间产生“泄漏泄漏”，在通常使用的入射角范围内，这种，在通常使用的入射角范围内，这种“泄漏泄漏”压制了一些抗噪压制了一些抗噪属性。这使得属性。这使得AV0反演不能提供满意的解决办法，然而，弹性波阻反演不能提供满意的解决办法，然而，弹性波阻抗能够避免这个问题，因此弹性波阻抗具有优越性。抗能够避免这个问题，因此弹性波阻抗具有优越性
30、。AVO与与EI反演反演IPRP(P-P)声阻抗IPRP(P-P)弹性波阻抗IPRS(P-S)弹性波阻抗ISRS(S-S)剪切阻抗 角道集叠加过程可以角道集叠加过程可以简单描述为：把固定简单描述为：把固定偏移距的记录转换成偏移距的记录转换成固定入射角固定入射角(或一定或一定角度范围内的叠加角度范围内的叠加)记录，将振幅随偏移记录，将振幅随偏移距变化的关系转换为距变化的关系转换为振幅随入射角变化的振幅随入射角变化的关系关系计算时间计算时间T()jT()j对应的炮检距对应的炮检距X()jX()j在点在点(T()j-1(T()j-1，X()j-1),(X()j-1),(T()jT()j，X()j)X
31、()j)之间建立线性关系之间建立线性关系大角度对应时间大角度对应时间T(min)T(min)小角度对应时间小角度对应时间T(max)T(max)提取各道的提取各道的T(min)T(min)到到T(max)T(max)之间的振幅值做角度之间的振幅值做角度部分叠加部分叠加对应于角度对应于角度(i)(i)的角度道的角度道A()A()处理后的处理后的CMPCMP道集资料道集资料确定确定CMPCMP道集的中心角道集的中心角序列序列(1)(1)、(2)(2)、(i)(i)确定中心角确定中心角(i)(i)的的范围范围(i)min(i)min到到(i)max(i)max计算层速度计算层速度V Vintint均
32、方根速度均方根速度V Vrmsrms角度道数据沿图中曲线拾取。可以看出，一般情况下，同一时间（深度）处远偏移距对应大角度，小偏移距对应小角度；同一偏移距角度随深度减小。零偏移距双程旅行时(ms)05001000150020002500300005001000150020002500炮检距（）炮检距（）零炮检距双程旅行时（）零炮检距双程旅行时（）炮检距、入射角及旅行时关系图炮检距、入射角及旅行时关系图81216202428偏移距(m)偏移距和角度的关系偏移距和角度的关系n角叠加剖面可以为弹性阻抗反演提供资料角叠加剖面可以为弹性阻抗反演提供资料n角度的估算是一个复杂的过程，其准确性受地震资料前角度
33、的估算是一个复杂的过程，其准确性受地震资料前期处理的影响，为了或得更准确的角道集叠加资料，应期处理的影响，为了或得更准确的角道集叠加资料，应加强地震资料的前期处理，同时在估算角度方面应采用加强地震资料的前期处理，同时在估算角度方面应采用更复杂的数学模型更复杂的数学模型n角部分叠加资料比叠后资料包含了更丰富的信息；角部分叠加资料比叠后资料包含了更丰富的信息；AVOAVO和弹性波阻抗反演要求地震资料具有和弹性波阻抗反演要求地震资料具有AVAAVA特征，因此将特征，因此将CMPCMP道集资料转化为角道集部分叠加资料是非常必要的道集资料转化为角道集部分叠加资料是非常必要的 CBSAPVVEI)(弹性波
34、阻抗的近似算法有多种，在此我们用BP近似来计算弹性波阻抗，根据Connolly 在1999年提出的弹性波阻抗近似公式可以得到弹性波阻抗的如下表达式：*)(*)(*)()(CLnVBLnVALnEILnSP进而由上式可以得到如下表达式：对于同一角度来说A，B，C均为常数：2sin8KB2sin41KC2tan1A从弹性波阻抗体中提取岩性参数体从弹性波阻抗体中提取岩性参数体对于多个不同的角度而言有：)()(*)()(*)()(*)(EILnLnCVLnBVLnAspYCX 这是一个典型的型方程组。写成矩阵形式：上面表达式中：左边第一项是一个n3的常数矩阵；第二项是所要求的纵、横波速度及密度的对数形
35、式；右边是一个n1的常数矩阵，它是由反演得到的各个角度的弹性波阻抗体分别取对数得到的。)()()()()()()()()()()()()()()()()()()(21222111nispnnniiiEILnEILnEILnEILnLnVLnVLnCBACBACBACBA从弹性波阻抗体中提取岩性参数体从弹性波阻抗体中提取岩性参数体通过最小平方法得到各样点的)()()(LnVLnVLnsp、spVV可以进一步得到：Ip、Is、泊松比、纵波模量、横波模量及拉梅参数等。对其分别取指数就得到各采样点的岩性参数体：从弹性波阻抗体中提取岩性参数体从弹性波阻抗体中提取岩性参数体参数Lamda*Rho和Mu*R
36、ho的获取由于由于LamdaLamda*RhoRho和和MuMu*RhoRho能够灵敏地反映储层属性，并且物理意义明确，能够灵敏地反映储层属性，并且物理意义明确，在岩性和流体预测方面具有重大的意义，因此我们有必要从反演中拾取这在岩性和流体预测方面具有重大的意义，因此我们有必要从反演中拾取这两个属性。有两种方法可以获取两个属性。有两种方法可以获取LambdaLambda*RhoRho和和MuMu*RhoRho，在此称其为直接，在此称其为直接法和间接法法和间接法 在在ConnollyConnolly（19991999）方程的基础上，采用回归算法求取）方程的基础上，采用回归算法求取AlphaAlph
37、a、BetaBeta和和RhoRho，根据各岩性参数间的关系进一步求得，根据各岩性参数间的关系进一步求得IpIp、IsIs、LambdaLambda*RhoRho和和MuMu*RhoRho 用用IpIp和和IsIs及及LambdaLambda*RhoRho和和MuMu*RhoRho的显式近似方程的显式近似方程FattiFatti近似和近似和GrayGray近似来近似来求取它们，这种方法是直接得到求取它们，这种方法是直接得到IpIp和和IsIs及及LambdaLambda*RhoRho和和MuMu*RhoRho77感谢亲观看此幻灯片，此课件部分内容来源于网络，感谢亲观看此幻灯片，此课件部分内容来源于网络，如有侵权请及时联系我们删除，谢谢配合！如有侵权请及时联系我们删除，谢谢配合！
