加强交流增进友谊!地震反演原理、方法地震反演的实际应用JASON反演的工作流程地震反演软件培训2006年11月物探研究院第一章国内外研究现状第二章叠后约束地震反演原理及常用反演方法第三章波阻抗反演影响因素第四章复杂岩性储层反演关键技术环节第五章多种类型复杂岩性储层反演技术对策及应用第六章Jason软件结构第七章Geoframe系统与JASON系统数据交换第八章反演数据的综合研究第九章叠前地震纵、横波阻抗(速度)约束反演目录国内外研究现状隐蔽性油气藏已成为勘探开发重要目标,岩性油气藏占有其中很大的比重。对河道砂体、三角洲浊积砂体、砂砾岩体等复杂岩性储层进行高精度预测和描述,成为重要的研究问题。将具有高横向分辨能力的地震资料、高垂向分辨能力的测井资料与地质结合起来进行地震反演,把界面反射型地震资料反演成岩层单元型波阻抗等地层物性参数,以岩层为目标进行储层解释。概述融合多种资料和技术所长的约束地震反演处理与解释应用技术是研究复杂岩性储层变化规律的十分有效的储层描述新技术。该项技术通过科学的数学算法和计算机软件,综合应用地震、测井、地质等资料,并将它们有机地结合为一体,把地震资料中包含的丰富岩性、物性、流体信息反演成岩层单元型波阻抗数据或地层层速度数据以及岩性和孔隙度等地层物性资料。使其能与钻井、测井资料直接对比分析,以岩层为目标进行地质、储层解释,在研究储层的空间特征和厚度变化规律方面具有独特的优势。反演技术的发展历程1、1979年(Lindseth)把波阻抗改写为反射系数的积分或求和形式,也可以表达成相位谱900旋转的一种褶积滤波器即递推反演。其商品化软件有Vlog、Dlog、Seislog等。2、Cooke和Schneider在1983年首次提出基于模型的反演。商品化软件有Glog、SLIM等。1990年,Debeye和VanRiel提出稀疏脉冲反演。上述方法在反演过程中都没有利用测井资料的高频信息,因此分辨率较低。3、90年代早中期,测井约束反演诞生和发展时期。突破传统地震频带的限制,具有比直接反演更高的分辨率。商品化的软件如strata、BCI、Jason等。4、90年代中后期,非线性反演理论为基础的各种反演算法出现和发展时期。如地震特征反演、随机反演、模拟退火反演、概率法神经网络反演、遗传算法、小波反演等。5、近几年,从叠后反演发展到叠前反演。叠前同步反演能够反演出纵横波阻抗等,直接进行岩性和流体识别。Zoeppritz方程组是叠前反演的理论基础,但关系复杂难以直接求解。Shuey和Aki、Richards等在二十世纪八十年代对其进行了简化研究。AliTura(1999)、DavidMDolberg(2000)等AVO反演、AVA反演方法。JASON公司等推出了实用的叠前和叠后地震反演软件系统。反演方法分类1、基于原始地震资料叠前反演类、叠后反演类叠后地震反演使用叠后地震资料,反演的纵波阻抗能够在一定程度上反映储层的变化规律,但不能给出横波波阻抗。叠前地震反演使用叠前地震道集资料,具有良好的保真性和多信息性。能同时反演纵、横波阻抗等反映地层岩性、物性和流体特征的多种参数。2、基于数学算法线性反演类、非线性反演类3、基于实现方式基于反褶积的反演方法:包括基于地层反褶积的道积分、递归、广义线性等反演方法;基于最大似然反褶积的稀疏脉冲反演方法;基于最大后验概率准则算法反褶积的模拟退火等。基于波动方程的反演:Born反散射。基于随机过程的反演:随机反演、随机模拟等。基于特征分析的反演:特征反演、神经网络反演等。其它:混沌反演等反演方法分类4、基于地震和测井的相对作用带限反演:道积分、Vlog、Seislog和Verilog*、递归等测井约束下的宽带反演:广义线性、宽带约束反演、稀疏脉冲反演等。地震约束下的测井内插外推:随机模拟、随机反演等。测井—地震联合反演:特征反演、神经网络反演等。反演方法分类第一章国内外研究现状第二章叠后约束地震反演原理及常用反演方法第三章波阻抗反演影响因素第四章复杂岩性储层反演关键技术环节第五章多种类型复杂岩性储层反演技术对策及应用第六章Jason软件结构第七章Geoframe系统与JASON系统数据交换第八章反演数据的综合研究第九章叠前地震纵、横波阻抗(速度)约束反演目录第一节稀疏脉冲地震波阻抗反演第二节随机模拟地震波阻抗反演第三节协模拟岩性及孔隙度反演第四节常用反演流程第二章叠后约束地震反演原理及常用反演方法JASONEnvironmentEarthModelStatModInverModelLargoWaveletsInverTraceRockTraceFunctionModVelModFunctionMod(一)稀疏脉冲地震波阻抗约束反演方法1、方法原理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范数反褶积研究的基础上,形成了稀疏脉冲地震波阻抗反演方法。其基本假设是,地下地层的反射系数序列是由一系列服从高斯分布的大反射系数和小反射系数背景迭合而成的。从地质意义上讲,大反射系数代表的是地下不连续界面和岩性分界面。在反演过程中,通过计算地震脉冲的均方根值和噪音的均方根值,估算所给定的采样有反射的似然值,迭代优化求解稀疏脉冲反射系数模型所代表的地下波阻抗模型。2、技术关键(1)硬约束和软约束条件建立。(2)λ值等关键参数试验。(3)频带补偿、拓宽。3、实际资料使用特点使用地震资料品质较好。可用测井速度或地震处理速度进行约束。4、反演结果特点体现地震资料振幅、频率、相位特征,适用范围较广。但反演结果分辨率不太高。Y84Y118-1Y118-1Y8410mreservoirdelineation稀疏脉冲反演波阻抗剖面与对应的地震剖面对比(一)稀疏脉冲地震波阻抗约束反演方法TheMenuofConstrainedsparsespikeinversion稀疏脉冲反演—趋势控制TheTrendeditorMenu稀疏脉冲反演参数选择及QC质量检查地震记录与合成记录相关性讯噪比误差入稀疏脉冲反演—质量控制参数可以通过质量控制加强反演结果的可靠性和准确性数据道合并滤波器显示稀疏脉冲反演—频带补偿(道合并)补充地震缺少的低频成份从算法上,稀疏脉冲反演的计算流程可分为三步:①、反射系数反演采用最大似然反褶积进行反射系数的反演,最大似然反褶积对地层的假设认为:地层的反射系数是由较大的反射界面的反射和具有高斯背景的小反射叠加组合而成,导出一个最小目标函数:式中,R2和N2分别为反射系数和噪音的均方值,r(K)和n(K)表示第K个采样点的反射系数和噪音,M表示反射层数,L表示采样总数,λ表示给定反射系数的似然值。通过多次迭代,求取反射系数。)1ln()(2)ln(2)(1)(1122122MLMKKJLKLKnNrR②、根据反射系数的反演结果结合阻抗趋势计算一个初始的波阻抗根据最大似然反褶积计算得到的反射系数,结合初始阻抗模型,采用递推算法,反演得到初始的波阻抗模型:式中,Z(i)为第i层的波阻抗值,R(i)为第i层的反射系数。)(1)(1)1()(iRiRiZiZ③、结合井的约束条件进行波阻抗反演约束稀疏脉冲反演对每一道依据目标函数对计算出的初始波阻抗进行调整,包括对反射系数的调整。目标优化函数为:式中,r为反射系数序列,Δz为与阻抗趋势的差序列,d为地震道序列,s为合成地震道序列,λ为残差权重因子,α为趋势权重因子,p、q为L模因子。具体的,右式第一项反映了反射系数的绝对值和,第二项反映了合成声波记录与原始地震数据的差值,第三项为趋势约束项。ZLdsLrLFqp11)()(约束稀疏脉冲反演是基于道的反演,它的实质就是在阻抗趋势的约束下,用最少数目的反射系数脉冲达到合成记录与地震道的最佳匹配。稀疏脉冲反演处理流程数据加载初始模型分析和预处理、极性分析井旁子波初标定相关性分析综合标定选择反演井、范围、子波等QC控制选择参数设置趋势和边界等约束条件残差分析生成相对阻抗体设置滤波器高频补偿低频补偿输出反演成果色标子波优选校正层位解释断层解释及趋势面分析地层结构分析框架模型EarthmodelConstraintsparse…Tracemerge约束稀疏脉冲波阻抗反演中关键的处理环节(1)、做好地层、构造解释、建好初始模型a)合成记录标定;b)完成精细的地层和构造解释,要求除了包括顶底反演控制层外,对地层结构和接触关系存在差异的层系必须进行加密层控解释,c)要求断层解释完全闭合;d)根据地震反射特征建立准确的地层接触关系,建立框架模型c)根据测井信息、断层面和地层框架模型建立准确的初始波阻抗模型.T4中1中2中3中4中5中6中7T6T6中8中5中6中4中4中3深入研究三角洲的九期次沉积,解释九套地质层位。精细地质层位解释T1’T1T2’T2综合标定反演的迭代次数参数对反演效果的影响最大反射系数反演迭代次数Max.#iterationreflectivityinversion波阻抗反演迭代次数Max.#iterationsP-impedanceinversion稀疏迭代次数Sparsity目标函数优化迭代次数Objectivefunctionminization稀疏迭代次数Sparsity目标函数优化迭代次数Objectivefunctionminization做好质量控制合理使用井的约束作用反演过程中井的约束控制约束强度:1、无约束或宽松的约束2、强的井约束)()()(11ZtrendLdsLrLFqp约束方法:无约束(disableanyconstraint)软约束(softconstraint):α-1,硬约束(hardconstraint):趋势+边界条件趋势约束方式:1、固定值Constant2、沿层变化lateral3、模型约束model滨南地区沙三段多井约束稀疏脉冲反演波阻抗剖面高速火成岩WashFanY118-53Y132-5断层清楚河道砂体横截面Y132-5Y941Y942河道砂体纵切面渤南地区沙三段稀疏脉冲反演波阻抗剖面过牛110井的反演剖面Xin154Xin154-1Xin154-2Niu109-45反演效果分析—预测井效果分析辛154-2辛154第一节稀疏脉冲地震波阻抗反演第二节随机模拟地震波阻抗反演第三节协模拟岩性及孔隙度反演第四节常用反演流程第二章叠后约束地震反演原理及常用反演方法(二)、随机模拟地震波阻抗反演方法1、方法原理将地质统计模拟与地震反演紧密结合在一起。2、技术关键●变差函数分析拟合●约束模型精细建立●等概率反演结果筛选3、实际资料使用特点●适用于钻井资料丰富的三维工区:●变差函数与储层空间变化特征结合:●与其它反演成果相互结合4、反演结果特点反演结果分辨率高。在钻井较少地区,会因统计规律不强而产生假象。StatModMenudisplayCalculationandAnaylsisofTransformCalculationandAnaylsisofVariogramStatModModelingTheobjectivesofStatModModelingaretobuild3Dreservoirmodelsorreservoirmaps,basedongeostatisticalinformationderivedinStatModAnalysis.随机模拟、随机反演的基本原理:以地质框架模型、测井和地震资料为基础,以层为单位,利用储层/油气藏参数的空间分布规律和空间相关性进行随机模拟/随机反演,获得一组等概率的储层/油气藏参数模型。目的:·生成既满足测井资料和地质统计特性,又满足地震资料的储层/油气藏参数模型。•更准