jason地震反演软件培训

《jason地震反演软件培训》由会员分享，可在线阅读，更多相关《jason地震反演软件培训（106页珍藏版）》请在金锄头文库上搜索。

1、地震反演软件培训地震反演软件培训 第一章第一章国内外研究现状国内外研究现状第二章第二章 叠后约束地震反演原理及常用反演方法叠后约束地震反演原理及常用反演方法第三章第三章波阻抗反演影响因素波阻抗反演影响因素第四章第四章复杂岩性储层反演关键技术环节复杂岩性储层反演关键技术环节第五章第五章多种类型复杂岩性储层反演技术对策及应用多种类型复杂岩性储层反演技术对策及应用第六章第六章Jason软件结构软件结构第七章第七章Geoframe系统与系统与JASON系统数据交换系统数据交换第八章第八章反演数据的综合研究反演数据的综合研究第九章第九章叠前地震纵、横波阻抗（速度）约束反演叠前地震纵、横波阻抗（速度）约束

2、反演目录目录国内外研究现状国内外研究现状 隐蔽性油气藏已成为勘探开发重要目标，隐蔽性油气藏已成为勘探开发重要目标，岩性油气藏占有其中很大的比重。对河道砂岩性油气藏占有其中很大的比重。对河道砂体、三角洲浊积砂体、砂砾岩体等复杂岩性体、三角洲浊积砂体、砂砾岩体等复杂岩性储层进行高精度预测和描述，成为重要的研储层进行高精度预测和描述，成为重要的研究问题。究问题。 将具有高横向分辨能力的将具有高横向分辨能力的地震地震资料、高资料、高垂向分辨能力的垂向分辨能力的测井测井资料与资料与地质地质结合起来进结合起来进行地震反演，把界面反射型地震资料反演成行地震反演，把界面反射型地震资料反演成岩层单元型波阻抗等地

3、层物性参数，以岩层岩层单元型波阻抗等地层物性参数，以岩层为目标进行储层解释。为目标进行储层解释。概述概述 融合多种资料和技术所长的约束地震反演处理融合多种资料和技术所长的约束地震反演处理与解释应用技术是研究复杂岩性储层变化规律的十与解释应用技术是研究复杂岩性储层变化规律的十分有效的储层描述新技术。分有效的储层描述新技术。 该项技术通过科学的数学算法和计算机软件，该项技术通过科学的数学算法和计算机软件，综合应用地震、测井、地质等资料，并将它们有机综合应用地震、测井、地质等资料，并将它们有机地结合为一体，把地震资料中包含的丰富岩性、物地结合为一体，把地震资料中包含的丰富岩性、物性、流体信息反演成岩

4、层单元型性、流体信息反演成岩层单元型波阻抗数据波阻抗数据或地层或地层层速度数据以及岩性和孔隙度等地层物性资料。使层速度数据以及岩性和孔隙度等地层物性资料。使其能与钻井、测井资料直接对比分析，以岩层为目其能与钻井、测井资料直接对比分析，以岩层为目标进行地质、储层解释，在研究储层的空间特征和标进行地质、储层解释，在研究储层的空间特征和厚度变化规律方面具有独特的优势。厚度变化规律方面具有独特的优势。 5 5、近几年，从叠后反演发展到叠前反演、近几年，从叠后反演发展到叠前反演 。叠前同步反演能够反演出纵横波阻抗等，直接进行岩性叠前同步反演能够反演出纵横波阻抗等，直接进行岩性和流体识别。和流体识别。Zo

5、eppritzZoeppritz方程组是叠前反演的理论基础，但关系复杂方程组是叠前反演的理论基础，但关系复杂难以直接求解。难以直接求解。ShueyShuey和和AkiAki、RichardsRichards等在二十世纪等在二十世纪八十年代对其进行了简化研究。八十年代对其进行了简化研究。AliAli Tura Tura（19991999）、）、David MDavid M Dolberg Dolberg（20002000）等）等AVOAVO反演、反演、AVAAVA反演方反演方法。法。JASONJASON公司等推出了实用的叠前和叠后地震反演公司等推出了实用的叠前和叠后地震反演软件系统。软件系统。反

6、演方法分类反演方法分类1 1、基于原始地震资料、基于原始地震资料 叠前反演类、叠后反演类叠前反演类、叠后反演类 叠后地震反演使用叠后地震资料，反演叠后地震反演使用叠后地震资料，反演的纵波阻抗能够在一定程度上反映储层的变的纵波阻抗能够在一定程度上反映储层的变化规律，但不能给出横波波阻抗。化规律，但不能给出横波波阻抗。 叠前地震反演使用叠前地震道集资料，叠前地震反演使用叠前地震道集资料，具有良好的保真性和多信息性。能同时反演具有良好的保真性和多信息性。能同时反演纵、横波阻抗等反映地层岩性、物性和流体纵、横波阻抗等反映地层岩性、物性和流体特征的多种参数。特征的多种参数。2 2、基于数学算法、基于数学

7、算法线性反演类、非线性反演类线性反演类、非线性反演类3 3、基于实现方式、基于实现方式基于反褶积的反演方法：包括基于地层反褶积的道积分、基于反褶积的反演方法：包括基于地层反褶积的道积分、递归、广义线性等反演方法；基于最大似然反褶积的稀疏脉冲递归、广义线性等反演方法；基于最大似然反褶积的稀疏脉冲反演方法反演方法; ;基于最大后验概率准则算法反褶积的模拟退火等。基于最大后验概率准则算法反褶积的模拟退火等。基于波动方程的反演：基于波动方程的反演：BornBorn反散射。反散射。基于随机过程的反演：随机反演、随机模拟等。基于随机过程的反演：随机反演、随机模拟等。基于特征分析的反演：特征反演、神经网络反

8、演等。基于特征分析的反演：特征反演、神经网络反演等。其它：混沌反演等其它：混沌反演等反演方法分类反演方法分类4 4、基于地震和测井的相对作用、基于地震和测井的相对作用带限反演：道积分、带限反演：道积分、VlogVlog、SeislogSeislog和和VerilogVerilog* *、递归等递归等测井约束下的宽带反演：广义线性、宽带约束反演、稀测井约束下的宽带反演：广义线性、宽带约束反演、稀疏脉冲反演等。疏脉冲反演等。地震约束下的测井内插外推：随机模拟、随机反演等。地震约束下的测井内插外推：随机模拟、随机反演等。测井测井地震联合反演地震联合反演 ：特征反演、神经网络反演等。：特征反演、神经网

9、络反演等。反演方法分类反演方法分类第一章第一章国内外研究现状国内外研究现状第二章第二章 叠后约束地震反演原理及常用反演方法叠后约束地震反演原理及常用反演方法第三章第三章波阻抗反演影响因素波阻抗反演影响因素第四章第四章复杂岩性储层反演关键技术环节复杂岩性储层反演关键技术环节第五章第五章多种类型复杂岩性储层反演技术对策及应用多种类型复杂岩性储层反演技术对策及应用第六章第六章Jason软件结构软件结构第七章第七章Geoframe系统与系统与JASON系统数据交换系统数据交换第八章第八章反演数据的综合研究反演数据的综合研究第九章第九章叠前地震纵、横波阻抗（速度）约束反演叠前地震纵、横波阻抗（速度）约束

10、反演目录目录第一节第一节 稀疏脉冲地震波阻抗反演稀疏脉冲地震波阻抗反演第二节第二节随机模拟地震波阻抗反演随机模拟地震波阻抗反演第三节第三节协模拟岩性及孔隙度反演协模拟岩性及孔隙度反演第四节第四节 常用反演流程常用反演流程第二章第二章 叠后约束地震反演原理及常用反演方法叠后约束地震反演原理及常用反演方法JASONJASONEnvironmentEarthModelStatModInverModelLargoWaveletsInverTraceRockTraceFunctionModVelModFunctionMod（一）稀疏脉冲地震波阻抗约束反演方法（一）稀疏脉冲地震波阻抗约束反演方法 1、方法

11、原理、方法原理 F=Lp(r)+Lq(s-d)+-1L1(Z)+in(L1(Z-Zin)+xL1(Z-Zx) 其中： Lp(r)+Lq(s-d)是是基本地震道反演目标函数基本地震道反演目标函数基本地震道反演目标函数基本地震道反演目标函数L1(Z)纵向软趋势约束纵向软趋势约束纵向软趋势约束纵向软趋势约束L1(Z-Zin)+L1(Z-Zx)横向软空间约束横向软空间约束横向软空间约束横向软空间约束稀疏脉冲反演方法是目前使用较多的一种稀疏脉冲反演方法是目前使用较多的一种重要重要的的地震波阻抗反演方法。在地震波阻抗反演方法。在最大似然反褶积最大似然反褶积和和L1范数范数反褶积反褶积研究的基础上，形成了稀

12、疏脉冲地震波阻抗反研究的基础上，形成了稀疏脉冲地震波阻抗反演方法。演方法。其其基本假设基本假设是，地下地层的反射系数序列是由一系是，地下地层的反射系数序列是由一系列服从高斯分布的大反射系数和小反射系数背景迭合列服从高斯分布的大反射系数和小反射系数背景迭合而成的。从地质意义上讲，大反射系数代表的是地下而成的。从地质意义上讲，大反射系数代表的是地下不连续界面和岩性分界面。不连续界面和岩性分界面。在反演过程中，通过计算地震脉冲的均方根值和在反演过程中，通过计算地震脉冲的均方根值和噪音的均方根值，估算所给定的采样有反射的似然值，噪音的均方根值，估算所给定的采样有反射的似然值，迭代优化求解稀疏脉冲反射系

13、数模型所代表的地下波迭代优化求解稀疏脉冲反射系数模型所代表的地下波阻抗模型。阻抗模型。2、技术关键、技术关键（1）硬约束和软约束条件硬约束和软约束条件建立建立。 （2）值等值等关键参数关键参数试验。试验。 （3）频带频带补偿补偿、拓宽、拓宽。3、实际资料使用特点、实际资料使用特点使用地震资料品质较好。可用测井速度使用地震资料品质较好。可用测井速度或地震处理速度进行约束。或地震处理速度进行约束。4、反演结果特点、反演结果特点体现地震资料振幅、频率、相位特征，体现地震资料振幅、频率、相位特征，适用范围较广适用范围较广。但反演结果分辨率不太高。但反演结果分辨率不太高。Y84Y118-1Y118-1Y

14、8410mreservoirdelineation稀疏脉冲反演波阻抗剖面与对应的地震剖面对比稀疏脉冲反演波阻抗剖面与对应的地震剖面对比（一）（一） 稀疏脉冲地震波阻抗约束反演方法稀疏脉冲地震波阻抗约束反演方法 The Menu of Constrained sparse spike inversionThe Menu of Constrained sparse spike inversion稀疏脉冲反演稀疏脉冲反演趋势控制趋势控制TheTrendeditorMenuTheTrendeditorMenu 稀疏脉冲反演参数选择及质量检查稀疏脉冲反演参数选择及质量检查地震记录与合成记录相关性讯噪比误

15、差入稀疏脉冲反演稀疏脉冲反演质量控制参数质量控制参数可可可可可可以以以以以以通通通通通通过过过过过过质质质质质质量量量量量量控控控控控控制制制制制制加加加加加加强强强强强强反反反反反反演演演演演演结结结结结结果果果果果果的的的的的的可可可可可可靠靠靠靠靠靠性性性性性性和和和和和和准准准准准准确确确确确确性性性性性性数据道合并滤波器显示数据道合并滤波器显示稀疏脉冲反演稀疏脉冲反演频带补偿频带补偿(道合并道合并)补充地震缺少的低频成份从算法上，稀疏脉冲反演的计算流程可分为三步：从算法上，稀疏脉冲反演的计算流程可分为三步：、反射系数反演、反射系数反演采采用用最最大大似似然然反反褶褶积积进进行行反反射

16、射系系数数的的反反演演，最最大大似似然然反反褶褶积积对对地地层层的的假假设设认认为为：地地层层的的反反射射系系数数是是由由较较大大的的反反射射界界面面的的反反射射和和具具有有高高斯斯背景的小反射叠加组合而成，导出一个最小目标函数：背景的小反射叠加组合而成，导出一个最小目标函数：式式中中，R2和和N2分分别别为为反反射射系系数数和和噪噪音音的的均均方方值值，r(K) 和和n(K)表表示示第第K个个采采样样点点的的反反射射系系数数和和噪噪音音，M表表示示反反射射层层数数，L表表示示采采样样总总数数，表表示给定反射系数的似然值。示给定反射系数的似然值。通过多次迭代，求取反射系数。通过多次迭代，求取反

17、射系数。、根根据据反反射射系系数数的的反反演演结结果果结结合合阻阻抗抗趋趋势势计计算算一一个个初初始的波阻抗始的波阻抗根根据据最最大大似似然然反反褶褶积积计计算算得得到到的的反反射射系系数数，结结合合初初始始阻阻抗抗模模型型，采用递推算法，反演得到初始的波阻抗模型：采用递推算法，反演得到初始的波阻抗模型：式中，式中，Z（i）为第为第i层的波阻抗值，层的波阻抗值，R(i)为第为第i层的反射系数。层的反射系数。、结合井的约束条件进行波阻抗反演、结合井的约束条件进行波阻抗反演约约束束稀稀疏疏脉脉冲冲反反演演对对每每一一道道依依据据目目标标函函数数对对计计算算出出的的初初始始波波阻阻抗抗进进行行调整，

18、包括对反射系数的调整。目标优化函数为：调整，包括对反射系数的调整。目标优化函数为：式式中中，r为为反反射射系系数数序序列列，z为为与与阻阻抗抗趋趋势势的的差差序序列列，d为为地地震震道道序序列列，s为为合合成成地地震震道道序序列列，为为残残差差权权重重因因子子，为为趋趋势势权权重重因因子子，p、q为为L模模因因子子。具具体体的的，右右式式第第一一项项反反映映了了反反射射系系数数的的绝绝对对值值和和，第第二二项项反反映映了了合成声波记录与原始地震数据的差值，第三项为趋势约束项。合成声波记录与原始地震数据的差值，第三项为趋势约束项。约束稀疏脉冲反演是基于道的反演，它的实质就是在阻抗趋势的约束约束稀

19、疏脉冲反演是基于道的反演，它的实质就是在阻抗趋势的约束下，用最少数目的反射系数脉冲达到合成记录与地震道的最佳匹配。下，用最少数目的反射系数脉冲达到合成记录与地震道的最佳匹配。稀疏脉冲反演稀疏脉冲反演处理流程处理流程数据加载数据加载初始模型初始模型分析和预处理、极性分析分析和预处理、极性分析井旁子波井旁子波初标定初标定相关性分析相关性分析综合标定综合标定选择反演井、范围、子波等选择反演井、范围、子波等QCQC控制选择参数控制选择参数设置趋势和边界等约束条件设置趋势和边界等约束条件残差分析残差分析生成相对阻抗体生成相对阻抗体设置滤波器设置滤波器高频补偿高频补偿低频补偿低频补偿输出反演成果输出反演成

20、果色标色标子波优选子波优选校正校正层位解释层位解释断层解释及趋势面分析断层解释及趋势面分析地层结构分析地层结构分析框架模型框架模型Earth modelConstraint sparseTrace merge约束稀疏脉冲波阻抗反演中关键的处理环节约束稀疏脉冲波阻抗反演中关键的处理环节 (1)(1)、做好地层、构造解释、建好初始模型、做好地层、构造解释、建好初始模型a)合成记录标定；合成记录标定；b)完完成成精精细细的的地地层层和和构构造造解解释释，要要求求除除了了包包括括顶顶底底反反演演控控制制层层外外，对对地地层层结结构构和和接接触触关关系系存存在在差差异异的的层层系系必须进行加密层控解释，

21、必须进行加密层控解释，c)要求断层解释完全闭合；要求断层解释完全闭合；d)根根据据地地震震反反射射特特征征建建立立准准确确的的地地层层接接触触关关系系，建建立立框架模型框架模型c)根据测井信息、断层面和地层框架模型建立准确根据测井信息、断层面和地层框架模型建立准确的初始波阻抗模型的初始波阻抗模型. .T4中1中2中3中4中5中6中7T6T6中8中中5中中6中4中中4中中3深入研究三角洲的九期次深入研究三角洲的九期次沉积，解释九套地质层位。沉积，解释九套地质层位。精细地质层位解释精细地质层位解释T1T1T2T2 综合标定综合标定反演的迭代次数参数对反演效果的影响反演的迭代次数参数对反演效果的影响

22、最大反射系数反演迭代次数最大反射系数反演迭代次数Max.#iterationreflectivityinversion波阻抗反演迭代次数波阻抗反演迭代次数Max.#iterationsP-impedanceinversion稀疏迭代稀疏迭代次数次数Sparsity目标函数优化迭代次数目标函数优化迭代次数Objectivefunctionminization稀疏迭代稀疏迭代次数次数Sparsity目标函数优化迭代次数目标函数优化迭代次数Objectivefunctionminization做好质量控制做好质量控制合理使用井的约束作用合理使用井的约束作用滨南地区沙三段多井约束稀疏脉冲反演波阻抗剖面

23、滨南地区沙三段多井约束稀疏脉冲反演波阻抗剖面高速火成岩高速火成岩WashFanY118-53Y118-53Y132-5Y132-5断层清楚断层清楚河道砂体横截面河道砂体横截面Y132-5Y132-5 Y941Y941 Y942 Y942河道砂体纵切面河道砂体纵切面渤南地区沙三段稀疏脉冲反演波阻抗剖面渤南地区沙三段稀疏脉冲反演波阻抗剖面过牛110井的反演剖面Xin154Xin154-1Xin154-2Niu109-45反反演演效效果果分分析析预预测测井井效效果果分分析析辛辛154-2辛辛154第一节第一节稀疏脉冲地震波阻抗反演稀疏脉冲地震波阻抗反演第二节第二节 随机模拟地震波阻抗反演随机模拟地震

24、波阻抗反演第三节第三节协模拟岩性及孔隙度反演协模拟岩性及孔隙度反演第四节第四节 常用反演流程常用反演流程第二章第二章 叠后约束地震反演原理及常用反演方法叠后约束地震反演原理及常用反演方法 （二）、随机模拟地震波阻抗反演方法（二）、随机模拟地震波阻抗反演方法 1、方法原理、方法原理 将地质统计模拟与地震反演紧密结合在一将地质统计模拟与地震反演紧密结合在一起。起。2、技术关键、技术关键 变差函数分析拟合变差函数分析拟合 约束模型精细建立约束模型精细建立 等概率反演结果筛选等概率反演结果筛选3、实际资料使用特点、实际资料使用特点 适用于钻井资料丰富的三维工区：适用于钻井资料丰富的三维工区： 变差函数

25、与储层空间变化特征结合：变差函数与储层空间变化特征结合： 与其它反演成果相互结合与其它反演成果相互结合4、反演结果特点、反演结果特点 反演结果分辨率高。反演结果分辨率高。在钻井较少地区，在钻井较少地区，会因统计规律不强而产生假象。会因统计规律不强而产生假象。StatModMenudisplayCalculationandCalculationandAnaylsisAnaylsisofTransformofTransformCalculationandCalculationandAnaylsisAnaylsisofofVariogramVariogramStatModModelingThe ob

26、jectives of The objectives of StatModStatMod Modeling are to build 3D reservoir models or Modeling are to build 3D reservoir models or reservoir maps, based on reservoir maps, based on geostatisticalgeostatistical information derived in information derived in StatModStatMod Analysis Analysis. 随机模拟、随

27、机反演的基本原理随机模拟、随机反演的基本原理: :以以地地质质框框架架模模型型、测测井井和和地地震震资资料料为为基基础础，以以层层为为单单位位，利利用用储储层层/ /油油气气藏藏参参数数的的空空间间分分布布规规律律和和空空间间相相关关性性进进行行随随机机模模拟拟/ /随随机机反反演演，获获得得一一组组等等概概率率的储层的储层/ /油气藏参数模型。油气藏参数模型。目的：目的： 生生成成既既满满足足测测井井资资料料和和地地质质统统计计特特性性，又又满满足足地地震震资资料料的的储储层层/ /油油气气藏藏参数模型。参数模型。 更准确地估算各种参数的不确定性，提供参数模拟的可靠性评价。更准确地估算各种参

28、数的不确定性，提供参数模拟的可靠性评价。模块及功能模块及功能: : StatmodStatmod-Analysis-Analysis 利用已知资料分析确定储层利用已知资料分析确定储层/ /油气藏参数的空间分布规律（直方图）和空油气藏参数的空间分布规律（直方图）和空间相关性间相关性( (变异函数变异函数) )。 StatmodStatmod- -ModellingModelling 用不同的地质统计模拟技术实现储层用不同的地质统计模拟技术实现储层/ /油气藏参数的随机模拟油气藏参数的随机模拟 随机反演。随机反演。在在三三维维空空间间，可可以以通通过过储储层层变变量量的的一一系系列列数数值值，模模

29、拟拟得得到到其其它它未未知知空空间间点点所所具具有有的的可可能能的的储储层层参参数数值值。克克里里金金方方法法( (KrigingKriging) )是是一一种种很很好好的的地地质质统统计计模模拟拟方方法法。应应用用变变差差函函数数模模型型所所提提供供的的空空间间结结构构信信息息，通通过过求求解解克克里里金金方方程程组组计计算算局局部部估估计计的的加加权权因因子子即即克克里里金金系系数数进进行行加加权权线线性性估估计计，充充分分考虑了空间数据的结构性和随机性。考虑了空间数据的结构性和随机性。大大量量的的随随机机模模拟拟过过程程实实现现所所得得到到的的空空间间储储层层参参数数体体在在统统计计特特

30、性性上上具具有有相相同同的的概概率率可可能能性性，并并且且与已有的实测数据结果具有同样的吻合程度。与已有的实测数据结果具有同样的吻合程度。假假定定仅仅用用井井中中测测量量值值通通过过随随机机模模拟拟能能勾勾画画出出声声阻阻抗抗曲曲线线。在在理理论论上上我我们们可可通通过过对对比比模模型型的的地地震震响响应应和和声声阻阻抗抗反反演演结结果果来来验验证证模模拟拟结结果果与与3D地震数据的一致性。地震数据的一致性。随随机机模模拟拟地地震震反反演演方方法法正正是是基基于于这这种种思思想想，以以测测井井、地地震震、地地质质资资料料为为基基础础，将将地地质质统统计计模模拟拟与与地地震震反反演演紧紧密密结结

31、合合在在一一起起，反反演演得得到到高高分分辨辨能能力力的的波波阻阻抗抗结结果果。随随机机地地震震反反演演从从随随机机建建模模产产生生的的一一系系列列储储层层模模型型中中，优优选选出出与与地地震震数数据据最最佳佳匹匹配配的的储储层层模模型型，是是通通过过波波阻阻抗抗将将储储层层特特性性和和地地震震记记录录相相联联系系来来直直接接估估计计储储层层参参数数特特征的一个完整的反演过程。征的一个完整的反演过程。反演前资料预处理子波估算地震、地质标定精细构造解释精细模型建立随机地震反演处理储层综合描述随机地震反演方法处理流程随机地震反演方法处理流程F50-1F50-2F8-1F8-2F50-1F50-2樊

32、樊东东地地区区随随机机地地震震反反演演波波阻阻抗抗连连井井剖剖面面煤层煤层砂层砂层泥岩泥岩随机地震波阻抗反演剖面随机地震波阻抗反演剖面（过大（过大11、大、大13井）井）反演阻抗与测井阻抗有良好的对应关系奥陶灰岩太2段海相灰岩精细波阻抗反演精细波阻抗反演反演优点：反演优点：反演分辨率高，纵向与测井波阻抗有很好的吻合，横向变化自然，与三维地震的变化相匹配。可以清楚的分清低速煤层和高速灰岩。横向尖灭点2 2、随机地震反演、随机地震反演随机模拟的基本原理和计算方法基础概念基础概念: :（1 1）、地质统计学）、地质统计学（2 2）、数据分析与概率密度函数）、数据分析与概率密度函数（3 3）、空间结构

33、分析与变差函数）、空间结构分析与变差函数（4 4）、变差函数与克里金插值）、变差函数与克里金插值地质统计学：地质统计学：以区域化变量理论为以区域化变量理论为基础，以变差函数为基本工具，研究基础，以变差函数为基本工具，研究那些在空间上即具有随机性又具有结那些在空间上即具有随机性又具有结构性的自然现象的科学。构性的自然现象的科学。随机模拟的实现方法建立在地质统随机模拟的实现方法建立在地质统计学之上的，所以又称地质统计模拟。计学之上的，所以又称地质统计模拟。随机模拟的基本原理和计算方法基础概念基础概念: :（1 1）、地质统计学）、地质统计学（2 2）、数据分析与概率密度函数）、数据分析与概率密度函

34、数（3 3）、空间结构分析与变差函数）、空间结构分析与变差函数（4 4）、变差函数与克里金插值）、变差函数与克里金插值数据分析采用直方图和散点图等形式分析变量的分布特征不同的分布具有不同概率密度ccdf，知道了方差和数学期望就可以明确一个随机变量的概率密度特征。正态分布:标准正态分布:对数正态分布:不同分布的概率密度公式随机模拟的基本原理和计算方法基础概念基础概念: :（1 1）、地质统计学）、地质统计学（2 2）、数据分析与概率密度函数）、数据分析与概率密度函数（3 3）、空间结构分析与变差函数）、空间结构分析与变差函数（4 4）、变差函数与克里金插值）、变差函数与克里金插值空间结构分析与变

35、差函数 空间结构分析就是分析观测数据的性质及其所空间结构分析就是分析观测数据的性质及其所表征的变量的空间分布特征表征的变量的空间分布特征 （空间的连续性、各向（空间的连续性、各向异性等）。异性等）。 手段：手段：协方差函数、相关函数、变差函数协方差函数、相关函数、变差函数等；等； 变差函数（变异函数变差函数（变异函数VariogramVariogram function function 、结结构函数构函数structure function structure function ）是最能够反映区域化是最能够反映区域化变量的空间结构的函数。一般用变差曲线来表示。变量的空间结构的函数。一般用变差

36、曲线来表示。变差函数的定义：设设Z(x)Z(x)是一个随机函数是一个随机函数 ， Z(x+h)Z(x+h)为沿为沿a a方向上方向上h h距离距离的随机函数值。的随机函数值。变差函数的定义：变差函数是在任一方向变差函数的定义：变差函数是在任一方向a a 相距相距h h的的两个区域化变量两个区域化变量Z(x)Z(x)及及Z(x+h)Z(x+h)增量的方差，它是增量的方差，它是h h和和a a 的的函数：函数：离散的情况下：离散的情况下：变差函数图变程：变程：指区域化变量在空间上具有相关性的范围。在变程范围之内，数据指区域化变量在空间上具有相关性的范围。在变程范围之内，数据具有相关性；而在变程之外

37、，数据之间互不相关，即在变程以外的观测值不具有相关性；而在变程之外，数据之间互不相关，即在变程以外的观测值不对估计结果产生影响。变程的大小反映了变量空间相关性的大小，变程相对对估计结果产生影响。变程的大小反映了变量空间相关性的大小，变程相对较大意味着该方向的观测数据在较大范围内相关，反之，则相关性较小较大意味着该方向的观测数据在较大范围内相关，反之，则相关性较小 。 具有不同变程的克里金插值图像（变程越大，变量的空间相关性越大）不同的地质体由于连续性不同，变程差异是很大的块金值：块金值：变差函数如果在原点间断，这在地质统计学中被称为变差函数如果在原点间断，这在地质统计学中被称为“块金效块金效应

38、应”，表现为在很短的距离内有较大的空间变异性，它可以由测量误差引，表现为在很短的距离内有较大的空间变异性，它可以由测量误差引起，也可以来自矿化现象的微观变异性。在取得有效数据的尺度上，这种起，也可以来自矿化现象的微观变异性。在取得有效数据的尺度上，这种微观变异性是不可得到的。在数学上块金值相当于变量纯随机性的部分。微观变异性是不可得到的。在数学上块金值相当于变量纯随机性的部分。如果无论多么小，两个随机变量都不相关，这种情况称为纯块金效应。如果无论多么小，两个随机变量都不相关，这种情况称为纯块金效应。基台值：基台值：代表变量在空间上的总变异性大小代表变量在空间上的总变异性大小, ,其为块金值和拱

39、高之和。所其为块金值和拱高之和。所谓拱高，为在取得有效数据的尺度上，可观测得到的变异性幅度大小。当谓拱高，为在取得有效数据的尺度上，可观测得到的变异性幅度大小。当块金值等于块金值等于0 0时，基台值即为拱高。时，基台值即为拱高。 变差函数的拟合变差函数的拟合h变差理论变差函数模型理论变差函数模型球状模型：球状模型：由一个真实变程由一个真实变程“和正的方差贡献或基台值和正的方差贡献或基台值c c来确定。来确定。 式中式中 cc基台值；基台值；aa变程；变程；h h 滞后距。滞后距。指数模型：指数模型：由一个真实变程由一个真实变程o(o(有效变程有效变程a a3)3)和正的方差贡献和正的方差贡献c

40、 c来确定。来确定。 变变差差函函数数渐渐近近地地逼逼近近基基台台值值。在在实实际际变变程程“处处，变变差差函函数数为为0.950.95c c。模模型型在在原原点点处为直线。处为直线。高斯模型：高斯模型：由一个真实变程。和正的方差贡献由一个真实变程。和正的方差贡献c c来确定。来确定。变差函数惭近地逼近基台值。在实际变程变差函数惭近地逼近基台值。在实际变程a a处，变差函数为处，变差函数为0.950.95c c。模模型在原点处为抛物线。为一种连续性较好但稳定性较差的模型。型在原点处为抛物线。为一种连续性较好但稳定性较差的模型。幂幂模模型型：由由一一个个幂幂值值和和正正的的斜斜率率c c来确定。

41、来确定。幂幂模模型型为为一一种种无无基基台台值值的的变变差差函函数数模模型型。这这是是一一种种特特殊殊的的模模型型。当当参参数数w w改改变变时时，它它可以表示原点附近的各种形状。可以表示原点附近的各种形状。当当w w1 1时时，变变差差函函数数为为一一直直线线，即即为为线线性性模模型型，这这一一模模型型即即为为著著名名的的布布朗朗运运动动( (其其随随机机函函数数的的理理论论模模型型为为随随机机行行走走过过程程) )的的变变差差函函数数模模型型；当当 时时，变变差差函函数数为为抛抛物物线线形形状，为分数布朗运动状，为分数布朗运动( (fBmfBm) )的变差函数模型。的变差函数模型。三维变差

42、分析三维变差分析变差函数曲线拟合变差函数曲线拟合随机模拟的基本原理和计算方法基础概念基础概念: :（1 1）、地质统计学）、地质统计学（2 2）、数据分析与概率密度函数）、数据分析与概率密度函数（3 3）、空间结构分析与变差函数）、空间结构分析与变差函数（4 4）、变差函数与克里金插值）、变差函数与克里金插值克里金估计克里金估计 克里金估计是一种进行局部估计的方法。它所提供的克里金估计是一种进行局部估计的方法。它所提供的是区域化变量在一个局部区域的平均值的最佳估计量，即是区域化变量在一个局部区域的平均值的最佳估计量，即最优（估计方差最小）、无偏最优（估计方差最小）、无偏( (估计误差的数学期望

43、为估计误差的数学期望为0)0)的估计。的估计。 克里金估计所利用的信息，通常为一组实测数据及其克里金估计所利用的信息，通常为一组实测数据及其相应的空间结构信息。应用变差函数模型所提供的空间结相应的空间结构信息。应用变差函数模型所提供的空间结构信息，通过求解克里金方程组计算局部估计的加权因子构信息，通过求解克里金方程组计算局部估计的加权因子即克里金系数，然后进行加权线性估计。这样，采用克里即克里金系数，然后进行加权线性估计。这样，采用克里金系数进行的局部估计就充分考虑了空间数据的结构性和金系数进行的局部估计就充分考虑了空间数据的结构性和随机性，从而使克里金方法优越于其它的一些传统的统计随机性，从

44、而使克里金方法优越于其它的一些传统的统计方法如距离平方反比加权和三次样条等插值方法。方法如距离平方反比加权和三次样条等插值方法。 设设 ， 为为区区域域上上的的一一系系列列观观测测点点， ， 为为相相应应观观测测点点处处的的随随机机变变量量。区区域域化化变变量量Z(Z(x x) )在在 处处的的随随机变量机变量 可采用一个线性组合来估计：可采用一个线性组合来估计：式中，式中， 为权系数。为权系数。 从从上上式式可可知知，求求取取的的关关键键是是利利用用统统计计模模型型确确定定的的值值。无偏性和估计方差最小被作为选取无偏性和估计方差最小被作为选取 的标准，即的标准，即 从这两个关系式可推导出求取

45、的克里金方程组。从这两个关系式可推导出求取的克里金方程组。普通克里金估计普通克里金估计 为变差函数，通过求解上述方程，得到为变差函数，通过求解上述方程，得到的解，计的解，计算出差值。算出差值。 同时可以计算出方差值：同时可以计算出方差值：利用拉格朗日乘数法利用拉格朗日乘数法 克里金估计方法克里金估计方法随机模拟技术建立在克立金差值方法的基础之上的，不同的随机模拟技术往往采用了不同的克立金技术。需要了解各种克立金技术的特点。1、简单克里金、简单克里金(SimpleKriging) 为已知常数为已知常数2、普通克立金、普通克立金(OrdinaryKriging) 普通克里金是简单克里金普通克里金是

46、简单克里金( (SK)SK)的最常用的变化形式。与简单克里金不的最常用的变化形式。与简单克里金不同的是，普通同的是，普通 克里金的为未知常数。通过限制所有的权值克里金的为未知常数。通过限制所有的权值之和为之和为l l，它将均值从估计值中过滤掉了，因而无需考虑平稳均值的它将均值从估计值中过滤掉了，因而无需考虑平稳均值的先验知识。普通克里金随机函数先验知识。普通克里金随机函数 服从二阶平稳假设。服从二阶平稳假设。3、协克里金、协克里金(Co-Kriging) 协同克里金的估计方法利用几个变量之间的空间相关性，对其中的协同克里金的估计方法利用几个变量之间的空间相关性，对其中的一个或几个变量进行空间估

47、计，可以提高估计的精度。采样点的数目不一个或几个变量进行空间估计，可以提高估计的精度。采样点的数目不足的情况在油藏描述中是经常遇到的，井比较少，资料不全不准都是造足的情况在油藏描述中是经常遇到的，井比较少，资料不全不准都是造成这种情况的原因。在被估计变量的观察数据较少的情况下，可利用协成这种情况的原因。在被估计变量的观察数据较少的情况下，可利用协同克里金的方法用其相关变量的信息进行弥补，以保证其估计精度。同克里金的方法用其相关变量的信息进行弥补，以保证其估计精度。 如如孔隙度渗透率。孔隙度渗透率。协同克里金估计的初始变量和二级变量的线性组合形式如下：协同克里金估计的初始变量和二级变量的线性组合

48、形式如下： 式中式中 为随机变量为随机变量Z Z在位置在位置0 0处的估计值；处的估计值； 分别为初分别为初始变量的始变量的n n个样本数据；个样本数据； 分别为二级变量的分别为二级变量的m m个样本数据；个样本数据； 和和 为需要确定的协同克里金加权系数。为需要确定的协同克里金加权系数。协克里金估计系统的建立和其它克里金系统的建立方法是大同小异协克里金估计系统的建立和其它克里金系统的建立方法是大同小异的。利用克里金估计的无偏性和最小二乘法可推导出传统的普通协的。利用克里金估计的无偏性和最小二乘法可推导出传统的普通协同克里金估计的方程组如下：同克里金估计的方程组如下：式中，式中， 和和 为拉格

49、朗日因子；为拉格朗日因子； 为协方差。为协方差。 4、同位协克里金、同位协克里金(CollocatedCo-Kriging) 同位协同克里金是协同克里金的一种简化形式，即如果二级变量密同位协同克里金是协同克里金的一种简化形式，即如果二级变量密集取样时，只保留与估计点同位的二级变量。如二级变量为地震属性。集取样时，只保留与估计点同位的二级变量。如二级变量为地震属性。同位协同克里金的估计值为：同位协同克里金的估计值为： 对应的协同克里金方程组只要求知道对应的协同克里金方程组只要求知道Z Z协方差函数协方差函数和和ZYZY互互协协方方差差函函， 后后者者可可以以通通过过以以下下的的模模型型来来近近似

50、：似：其其中中 和和 是是Z Z和和Y Y的的方方差差函函数数，是同位的是同位的Z Zx x数据的线性相关系数。数据的线性相关系数。5、指示克里金、指示克里金(-Kriging) 将变量值表示成将变量值表示成0 0或或1 1的形式（如砂岩的形式（如砂岩1 1、泥岩、泥岩0 0），进行类型变量的），进行类型变量的内插；是可以剔除高的特异值的方法，如电阻率。内插；是可以剔除高的特异值的方法，如电阻率。随机建模的基本原理随机建模的基本原理 随随机机模模拟拟是是指指以以已已知知的的信信息息为为基基础础，以以随随机机函函数数为为理理论论，应应用用随随机机模拟方法，产生可选的、等概率的储层模型的方法。模拟

51、方法，产生可选的、等概率的储层模型的方法。 随随机机函函数数由由一一个个区区域域化化变变量量的的分分布布函函数数和和协协方方差差函函数数( (或或变变差差函函数数) )来表征。来表征。 随随机机模模拟拟的的基基本本思思想想是是从从一一个个随随机机函函数数Z(u)Z(u)中中抽抽取取多多个个可可能能的的实实现现，即即人人工工合合成成反反映映Z(u)Z(u)空空间间分分布布的的可可供供选选择择的的、等等概概率率的的高高分分辨辨率率实实现现，记为记为代代表表变变量量Z(u)Z(u)在在非非均均质质场场A A中中空空间间分分布布的的L L个个可可能能的的实实现现。若若用用观观测测的的实实验验数数据据对

52、对模模拟拟过过程程进进行行条条件件限限制制，使使得得采采样样点点的的模模拟拟值值与与实实测测值值相相同同( (即忠实于硬数据即忠实于硬数据) )，就称为条件模拟；否则为非条件模拟。，就称为条件模拟；否则为非条件模拟。插值与模拟插值与模拟 克克里里金金插插值值为为局局部部估估计计方方法法，对对估估计计值值的的整整体体空空间间相相关关性性考考虑虑不不够够，它保证了数据的估计局部最优，却不能保证数据的总体最优。它保证了数据的估计局部最优，却不能保证数据的总体最优。克克里里金金插插值值法法为为光光滑滑内内插插方方法法，为为减减少少估估计计方方差差而而对对真真实实观观测测数数据据的的离离散散性性进进行行

53、了了平平滑滑处处理理，但但有有一一些些有有意意义义的的异异常常也也被被光光滑滑作作用用而而光光滑滑掉了。掉了。克里金插值法得到的是确定性的解克里金插值法得到的是确定性的解模模拟拟是是全全局局最最优优估估计计方方法法，在在保保证证全全区区最最优优的的情情况况下下，做做到到局局部部最最优，也就是能够准从全局的分布特征和空间结构特征。优，也就是能够准从全局的分布特征和空间结构特征。模拟得到的是不确定性的解模拟得到的是不确定性的解随机模拟的计算流程序贯模拟（序贯模拟（Si）的的计算步骤计算步骤：（1 1）数据统计分析，统计分布特征；）数据统计分析，统计分布特征；（2 2）空间结构分析，确定变差函数；）

54、空间结构分析，确定变差函数；（3 3）划分模拟网格；）划分模拟网格；（4 4）随机地选择一个待模拟的网格节点；）随机地选择一个待模拟的网格节点；（5 5）根据变差函数和分布特征估计该节点的累积条件分布函数）根据变差函数和分布特征估计该节点的累积条件分布函数( (ccdfccdf) )；（6 6）随机地从随机地从ccdfccdf中提取一个分位数作为该节点的模拟值；中提取一个分位数作为该节点的模拟值；（7 7）将该新模拟值加到条件数据组中；）将该新模拟值加到条件数据组中；（8 8）重复）重复4747步，直到所有节点都被模拟到为止，从而得到一个模拟实现。步，直到所有节点都被模拟到为止，从而得到一个模

55、拟实现。（9 9）重复上述步骤，得到多个模拟实现以供选择。）重复上述步骤，得到多个模拟实现以供选择。Jason包含的随机模拟技术(1)(1)序序贯贯法法：一一步步一一步步进进行行模模拟拟；根根据据条条件件处处理理资资料料和前步和前步( (次次) )模拟的元素模拟每个模型元素；模拟的元素模拟每个模型元素； 序贯高斯模拟（序贯高斯模拟（SGSSGS） 序贯高斯协模拟（序贯高斯协模拟（SGCSSGCS） 序贯高斯配置协模拟（序贯高斯配置协模拟（CGCCSCGCCS） 序贯阀值指示模拟（序贯阀值指示模拟（CTISCTIS） 序贯指示模拟（序贯指示模拟（SISSIS） 带趋势的序贯指示模拟（带趋势的带趋

56、势的序贯指示模拟（带趋势的SISSIS）。）。 (2) (2) 迭迭代代法法：该该方方法法连连续续( (不不断断) )修修正正初初始始成成像像使使得得目目标标函函数数达达到到最最小小，可可以以考考虑虑几几种种类类型型的的约约束束条条件件( (如如MetropolisMetropolis法或模拟退火法法或模拟退火法) )。序贯高斯模拟（序贯高斯模拟（SGSSGS） 高斯随机域是最经典的随机函数。该模型的最大特征是随机变高斯随机域是最经典的随机函数。该模型的最大特征是随机变量符合高斯分布量符合高斯分布( (正态分布正态分布) )。该方法主要用于连续变量。该方法主要用于连续变量( (如孔隙度如孔隙度

57、) )的随机模拟。对于高斯模拟来说，的随机模拟。对于高斯模拟来说，ccdfccdf可用简单克里金来求取。由可用简单克里金来求取。由于于ccdfccdf的正态性，因此，整个模拟过程被极大地简化，序贯地确定的正态性，因此，整个模拟过程被极大地简化，序贯地确定一系列的一系列的ccdfccdf就被简化为解一系列克里金方程组。就被简化为解一系列克里金方程组。序序 贯贯 高高 斯斯 协协 模模 拟拟 ( (Sequential Sequential Gaussian Gaussian Co-Co-Simulation)Simulation)方法方法 序序贯贯高高斯斯协协模模拟拟（SGCSSGCS）方方法法

58、是是应应用用协协克克里里金金求求取取ccdfccdf的的序序贯贯模模拟拟方方法法，其其步步骤骤与与序序贯贯高高斯斯协协模模拟拟相相似似，不不同同的的是是采采用用协协克克里里金金和和正正态态得得分分的的变变差差函函数数模模型型来来确确定定该该节节点点处处随随机机函函数数Y(u)Y(u)的的ccdfccdf函数的参数函数的参数( (均值和方差均值和方差) )，并求取，并求取ccdfccdf。 可以进行波阻抗以外的变量的模拟。可以进行波阻抗以外的变量的模拟。序贯高斯配置协模拟（序贯高斯配置协模拟（CGCCSCGCCS） 采采用用配配置置协协克克里里金金（Collocated Collocated C

59、o-Co-KrigingKriging）和和正正态态得得分分的的变变差差函函数数模模型型来来确确定定该该节节点点处处随随机机函函数数Y(u)Y(u)的的ccdfccdf函函数数的的参参数数( (均均值值和和方差方差) )，并求取，并求取ccdfccdf。 用于不同采样密度的变量的模拟。用于不同采样密度的变量的模拟。序贯指示模拟（序贯指示模拟（SISSIS） 序序贯贯指指示示模模拟拟( (SIS)SIS)既既可可用用于于类类型型变变量量，又又可可用用于于离离散散化化的的连连续续变变量量类类别别的的随随机机模模拟拟。该该方方法法无无需需假假设设原原始始样样本本服服从从正正态态分分布布，而而是是通通

60、过过给给出出一一系系列列的的门门槛槛值值，估估计计某某一一类类型型变变量量或或离离散散化化连连续续变变量量低低于于某某一一门门槛槛值值的的概概率率，以以此此确确定定随随机机变变量量的的分分布布。该该方方法法实实际际上上是是应应用用指指示示克克里里金金求求取取ccdfccdf的的序序贯贯模模拟拟方方法法，其其主主要要特特点点是变量的指示变换、指示克里金和序贯模拟算法。是变量的指示变换、指示克里金和序贯模拟算法。 带趋势的序贯指示模拟（带趋势的带趋势的序贯指示模拟（带趋势的SISTRSISTR） 序序贯贯指指示示模模拟拟为为经经典典的的SISSIS方方法法，其其中中各各相相的的比比例例是是面面向向

61、全全部部网网格格的的，也也就就是是说说，在在模模拟拟不不同同的的网网格格节节点点时时各各相相比比例例是是不不变变的的，显显然然，这这种种算算法法具具有有局局限限性性。当当局局部部变变化化时时，可可应应用用具具有有趋趋势势的的序序贯贯指指示示模模拟拟方方法法( (SISTR)SISTR)，以以对对经经典典SISSIS方方法法进进行行改改进进。SISTRSISTR方方法法通通过过从从地地震震数数据据和和其其它它定定性性和和定定量量数数据据中中提提取取信信息息，从从而而为为每每一一个个模模拟节点提供一个局部相比例，各相比例满足：拟节点提供一个局部相比例，各相比例满足： 相应地，考虑了局部相比例的指示

62、克里金估计公式变成如下形式：相应地，考虑了局部相比例的指示克里金估计公式变成如下形式： 序贯指示同位协同模拟（序贯指示同位协同模拟（SICoSimSICoSim） 该该方方法法是是经经典典SISSIS方方法法的的又又一一个个扩扩展展方方法法，它它通通过过用用同同协协同同指指示示克克里里金金方方法法代代替替指指示示克克里里金金方方法法从从而而成成功功地地将将二二级级变变量量( (如如地地震震属属性性) )引引入入SISSIS算法中。同位协同指示克里金估计方法公式可表述如下：算法中。同位协同指示克里金估计方法公式可表述如下： 埕北埕北35地区随机地震反演波阻抗剖面地区随机地震反演波阻抗剖面薄储层分

63、辨薄储层分辨清楚清楚厚储层识别厚储层识别准确准确测井吻合测井吻合程度高程度高断层清晰断层清晰第一节第一节稀疏脉冲地震波阻抗反演稀疏脉冲地震波阻抗反演第二节第二节随机模拟地震波阻抗反演随机模拟地震波阻抗反演第三节第三节 协模拟岩性及孔隙度反演协模拟岩性及孔隙度反演第四节第四节 常用反演流程常用反演流程第二章第二章 叠后约束地震反演原理及常用反演方法叠后约束地震反演原理及常用反演方法（三（三) ) 协模拟岩性及孔隙度反演方法协模拟岩性及孔隙度反演方法 1、方法原理、方法原理 在波阻抗反演基础上进行的地质统计协在波阻抗反演基础上进行的地质统计协模拟岩性或孔隙度反演。模拟岩性或孔隙度反演。 2、技术关

64、键、技术关键 基础数据关系的仔细统计分析基础数据关系的仔细统计分析 直方和变差函数拟合直方和变差函数拟合 等概率参数体的优选等概率参数体的优选 3、实际资料使用特点、实际资料使用特点 需要有大量实际岩性或孔隙度资料作为需要有大量实际岩性或孔隙度资料作为分析统计基础。分析统计基础。4、反演结果特点、反演结果特点 剔除了波阻抗表达储层与非储层的重剔除了波阻抗表达储层与非储层的重叠值所带来的多解性，能更为直观地识别叠值所带来的多解性，能更为直观地识别目的储层。目的储层。反演前资料预处理子波估算地震、地质标定精细构造解释精细模型建立随机地震反演处理储层综合描述随机岩性模拟反演方法处理流程随机岩性模拟反

65、演方法处理流程随机岩性模拟反演处理王王王王7878井随机反演波阻抗结果井随机反演波阻抗结果井随机反演波阻抗结果井随机反演波阻抗结果王王王王7878井岩性反演结果井岩性反演结果井岩性反演结果井岩性反演结果( (三）协模拟岩性及孔隙度反演方法三）协模拟岩性及孔隙度反演方法埕北埕北3535地区高斯配置协模拟泥质含量反演剖面地区高斯配置协模拟泥质含量反演剖面埕北埕北3535地区高斯配置协模拟泥质含量反演剖面地区高斯配置协模拟泥质含量反演剖面砂岩概率体砂泥岩砂泥岩概率密概率密度分布度分布函数函数砂泥岩概率体求取泥岩概率体高斯模拟求得的砂泥岩概率体剖面阻抗阻抗概概率率大11大13大开11砂泥岩剖面砂泥岩剖

66、面岩性反演与井上岩性吻合程度高（岩性剖面）岩性反演与井上岩性吻合程度高（岩性剖面）分辨率高岩性反演结果分析上部连续性差下部连续性好波阻抗剖面与岩性剖面对比图波阻抗剖面与岩性剖面对比图岩性剖面岩性剖面波阻抗剖面波阻抗剖面砂岩泥岩大11大13大开11岩性反演结果分析大11大13大开11孔隙度反演剖面（过大孔隙度反演剖面（过大11井南北向剖面）井南北向剖面）气层气层气层盒盒3段曲流河段曲流河孔隙度三维立体显示与综合解释盒盒1段网状河段网状河立体立体第一节第一节稀疏脉冲地震波阻抗反演稀疏脉冲地震波阻抗反演第二节第二节随机模拟地震波阻抗反演随机模拟地震波阻抗反演第三节第三节协模拟岩性及孔隙度反演协模拟岩

67、性及孔隙度反演第四节第四节 常用反演流程常用反演流程第二章第二章 叠后约束地震反演原理及常用反演方法叠后约束地震反演原理及常用反演方法反反演演基基础础处处理理分分析析反反演演处处理理运运算算及及应应用用基础资料准备基础资料准备基础资料预处理基础资料预处理稀疏脉冲反演稀疏脉冲反演随机地震反演随机地震反演岩性反演岩性反演反演基础准备反演基础准备储层综合描述储层综合描述地震资料地震资料测井资料测井资料地质资料地质资料地震资料地震资料目标处理及分析目标处理及分析测井资料标准化测井资料标准化统计分析及转换统计分析及转换地震地质标定子波求取精细构造解释约束模型建立反演处理反演处理频带补偿频带补偿全频带结果输出全频带结果输出测井阻抗测井阻抗直方分析直方分析及函数拟合及函数拟合测井阻抗测井阻抗变差分析变差分析及函数拟合及函数拟合随机反演处理随机反演处理测井岩性测井岩性直方分析直方分析及函数拟合及函数拟合测井岩性测井岩性变差分析变差分析及函数拟合及函数拟合阻抗体阻抗体直方分析直方分析及函数拟合及函数拟合岩性反演处理岩性反演处理反演结果应用反演结果应用储层综合描述储层综合描述