基于模型测井约束地震反演
地震约束下的密度反演

地震约束下的密度反演摘要:广义线性联合反演对其模型的先验了解依赖性很大,只有在给定初始猜测接近真实地球物理模型时才可以求出有意义的解估计。
地震波旅行时可以认为是目前最精确的地球物理探测量,本文在地震波旅行时的限制和纵波速度的约束下,利用常规最小二乘线性反演进行密度反演。
利用该方法得到的反演结果比较接近真实值,作为先验信息可初始化重震迭代联合反演的模型参数。
并且在速度参数接近准确的情况下,该反演结果比单由速度密度关系式得到的结果更接近真实值。
关键字:密度反演;旅行时限制;速度约束Density Inversion under the Restriction of Seismic DataAbstract: Joint inversion depends on much of the prior understanding of the model. Only when the given initial guess is close to the real physical model we can find out meaningful solution. We think that seismic waves travel time is the most accurate geophysical detecting quantity. Based on the restrictions of seismic waves travel time and the constraints of p-wave velocity, we can finish the density inversion using conventional least square linear inversion. The inversion results we get from the method as the priori information can initialize the parameters of the model used in iterative joint inversion of seismic and gravity data. When the velocity parameters are close to the accurate, the inversion results are close to the real value.Key Words: density inversion; restrictions of travel time; constraints of velocity0 引言地球物理反演是在地球物理学中利用地球表面观测到的物理现象推测地球内部介质物理状态的空间变化及物性结构的一个分支[1]。
利用测井约束反演技术提高王541东地区地震分辨率

利用测井约束反演技术提高王541东地区地震分辨率[摘要]针对地震频带宽度的限制,原有的砂体解释技术已不能满足油田开发的需要,测井约束地震反演是目前具有可能性和可操作性的比较成熟的提高地震分辨率的实用技术,在提高牛庄油田王541砂体分辨率应用取得成功。
[关键词]约束反演地震资料地质储层精细描述稀脉冲中图分类号:p631.81 文献标识码:a 文章编号:1009-914x (2013)11-0002-01在陆相薄储层地质条件下,要提高砂体的解释质量和描述精度,就必须提高地震剖面的分辨率。
用已知地震信息和测井资料作为约束条件,推算出高分辨率的地层波阻抗资料,并结合多学科知识,为王541砂体储层的深度、厚度、物性的精细描述提供可靠的依据。
1.王541东地区的开发地质现状王541东地区位于牛庄洼陷东部、王541开发区东北地区,包括王541,王542 ,王543 ,王53 ,王70等沙三中开发区块及其外围地区,区内含油层系有沙二段分流河道、沙三中浊积岩、沙三下深水浊积扇和沙四段浊积岩,而其中探明储量最多、产能规模大的主要含油层系是以三角洲前缘滑塌浊积岩为主要储层的沙三中岩性油藏,分布于牛庄中部牛25-王70断层两侧,根据目前的开发实践,大部分砂体含油,具有埋藏适中,产量高的优势。
在开发区接合部、外围的零散探井也钻遇沙三中油层,如南部的牛301井沙三中试油获14.2t工业油流,而东南部的牛302、303、牛90等多口井未见油气显示,在东部的通61、王72井区零星发现了沙三中含油砂体,说明区内浊积岩沉积是非常复杂的,从打到沙三中的井的数量来看,在东部地区探井数量少,沙三中属于低勘探区,潜力巨大,目前该区处在多块三位的结合部,影响了油藏描述的精度,目前该区的南坡大连片三维资料维地震资料分辨率低,不能够满足目前砂体描述精度的要求,影响了该区的滚动开发进程。
2.王541东沙三中多井测井约束反演处理2.1 南坡连片地震资料分析本次收集到的地震资料为南坡连片叠前连片处理地震资料,从地震资料情况来看如下图:地震资料的信噪比很高,但是地震资料的分辨率不高,做了该地震资料的频谱分析,显示地震资料的主频在23左右,以作为地震的分辨率,那么h=砂岩速度/(4*主频),该区砂岩的平均速度在3500米,那么分辨率为h=3500/4*23=38米,从该区钻井资料来看,储层一般单层厚度都小于10米,因此地震资料的分辨率是制约该区储层描述的瓶颈之一,必须对地震资料进行再处理,以提高地震资料的分辨率。
地震模型正演

地震模型正演与反演简介一、地震模型正演(seismic forward modeling)的概念如果我们已知地下的地质模型,它的地震响应如何?地震模型正演就是通过室内模拟得到地质模型对于地震波的响应。
地震模型正演包括物理模拟和数值模拟,数值模拟就是应用相应的地球物理方程和数值计算求解已知的地质模型在假定激发源的作用下的地震相应。
通常,我们针对特定的勘探区块,应用期望或实际的采集参数通过地震正演模拟野外地震采集,得到单炮记录,再通过速度分析、动校正、叠加、偏移等处理得到成像数据。
图1为Marmousi速度模型,图2为正演得到的炮集记录,图3为正演得到的叠加剖面。
图1 Marmousi模型图2正演炮集图3 正演叠加剖面二、数值模型正演方法通常,我们提到的模型正演为数值模拟的模型正演,目前常用的数值模拟地震模型正演方法包括基于射线原理的射线追踪法,以及基于波动方程的有限差分法、有限元法、积分方程法、快速傅里叶变换法和拟谱法等。
射线追踪法主要反映地震波的运动学特征,有限差分、有限元法则适合复杂地质构造的正演模拟,积分方程法涉及复杂的数学推导,快速傅里叶变换法在频率域计算得到正演数据。
三、数值模型正演的步骤数值模拟求解地震模型正演问题的步骤主要包括以下三个方面:1) 地质建模,根据研究对象和问题建立地球物理或地质模型;2) 数学建模,根据应用的物理手段和地球物理模型建立相应的数学模型;3) 模拟计算,选择正演计算方法,编写计算程序进行数值模拟计算。
四、什么是地震反演地震反演技术就是充分利用测井、钻井、地质资料提供的丰富的构造、层位、岩性等信息,从常规的地震剖面推导出地下地层的波阻抗、密度、速度、孔隙度、渗透率、沙泥岩百分比、压力等地球物理信息。
反演就是由地震数据得到地质模型,进行储层、油藏研究。
地震资料反演可分为两部分:1)通过有井(绝对)、无井(相对)波阻抗反演得到波阻抗、速度数据体。
2)利用测井、测试资料结合波阻抗、速度数据进行岩性反演,得到孔隙度、渗透率、砂泥百分比、压力等物理数据。
地震波阻抗反演方法综述

地震波阻抗反演方法综述一、地震反演技术研究现状地震反演方法是一门综合运用数学、物理、计算机科学等学科发展起来的新技术新方法,每当数学方法、物理理论有了新的认识和发展时,就会有新的地震反演技术、方法的提出。
随着计算机技术的不断发展、硬件设施的不断升级,这些方法技术得到了实践验证和提升,反过来地震反演技术运用中出现的新问题、新思路又不断促使数学方法、地球物理学理论的再次发展。
时至今日,地震反演技术仍然是一个不断发展、不断成熟、不断丰富着的领域。
反演是正演的逆过程,在地震勘探中正演是已知地下的地质构造情况、岩性物性分布情况,根据地震波传播规律和适当的数学计算方法模拟地震波在地下传播以及接收地震波传输到地表信息的过程。
地球物理反演就是使用已知的地震波传播规律和计算方法,将地表接收到的地震数据通过逆向运算,预测地下构造情况、岩性物性分布情况的过程。
地震波阻抗正演是对反演的理论基础和实现手段。
1959年美国人Edwin Laurentine Drake在宾夕法尼亚州开凿的第一口钻井揭开了世界石油工业的序幕。
从刚开始的查看地质露头、寻找构造高点寻找石油,到通过地震剖面的亮点技术寻找石油,再到现在运用多种科学技术手段进行油气资源的预测,石油勘探经历了一个飞速的发展历程。
声波阻抗(AI)是介质密度和波在介质中传播速度的乘积,它能够反映地下地质的岩性信息。
声波阻抗反演技术是20世纪70年代加拿大Roy Lindseth博士提出的,通过反演能够将反映地层界面信息的地震数据变为反映岩性变化的波阻抗(或速度)信息。
由于波阻抗与地下岩石的密度、速度等信息紧密联系,又可以直接与已知地质、钻井测井信息对比,因此广泛应用于储层的预测和油藏描述中,深受石油工作者的喜爱。
70年代后期,从地震道提取声波资料的合成声波技术得到了快速发展,以此为基础发展的基于模型的一维有井波阻抗反演技术,提高了反演结果的可靠性。
进入80年代,Cooke等人将数学中的广义线性方法运用于地震资料反演,提出了广义线性地震反演。
地球物理反演总结

地球物理反演总结一、名词解释(30)二、简答(30)三、综述(40)第一章正问题:给定一个问题,寻找答案反问题:给定一个答案,寻找问题适定性问题:解一定存在;解的唯一性;问题发生一些小的变动仅导致问题的解发生小的变动非适定性问题:解不一定存在;解可能不唯一;问题中小的变动导致问题解较大变动正演问题(正问题):已知模型m,求解数据d的过程反演问题(反问题):已知数据d,求解模型参数m的过程地震反演(SeismicInversion):把常规的界面型反射剖面转换成岩层型的测井剖面,将地震资料变成可与测井资料直接对比的形式,实现这种转换的处理过程叫地震反演。
地震反演在石油勘探开发中的应用:1、微构造识别;2、岩性预测;3、储层参数评价4、流体识别(烃类检测)地球物理反演:根据各种位场、地震波、地球自由振荡、交变电磁场、以及热学的地球物理观测数据去推测地球内部的结构形态及物质成分,定量计算各种有关的物理参数地震勘探中应用最广的反演问题是地震波阻抗反演,地震波阻抗反演是储层地球物理研究的最基本的处理技术之一,通过地震波阻抗反演把常规的界面型地震反射剖面转换成岩层型的伪测井剖面,因而使地震资料转变成可以与钻井资料直接对比的剖面形式,可以说波阻抗反演是地震资料处理的最终处理结果。
地震反演的目的:根据地震资料,反推出地下介质的波阻抗、速度和密度等岩石地球物理参数的分布,估算储层参数,并进行储层预测,以便为油气田的勘探和开发提供可靠的基础资料。
第二章地震正演定义:地震正演是根据设计的地质模型,选择速度、密度、波松比等地层参数,用某种方法求得地震响应,通过与实际剖面对比,合理解释复杂的地质现象。
地质模型:物理模型、数值模型算法原理:褶积模型、绕射迭加模型、射线追踪模型、波动方程模型地震物理模拟:在实验室内将野外的地质构造和地质体按照一定的模拟相似比制作成物理模型,并用超声波或激光超声波等方法对野外地震勘探方法进行模拟的一种地震模拟方法。
地震反演

地震反演姓名:李雪松班级:油气田s101 学号:201071059一、地震反演的基本定义通俗的讲就是由地震为基础加上其他条件为约束推测出地层岩性构造的过程叫地震反演。
把常规的界面型反射剖面转换成岩层型的测井剖面,将地震资料变成可与测井资料直接对比的形式,实现这种转换的处理过程叫地震反演。
地震反演:地震反演是利用地表观测地震资料,以已知地质规律和钻井、测井资料为约束,对地下岩层空间结构和物理性质进行成像(求解)的过程,广义的地震反演包含了地震处理解释的整个内容。
地震多解性和粗略性地震反演多解性是指同一地震资料可对应用不同的岩层结构,粗略性是指推断的参数少,分辨率低,前者可能导致地下模型的错误,后者影响模型的精度。
理论基础:鲁宾逊褶积模型基础。
其实现过程是:(1)子波提取;(2)逐道修改波阻抗模型道,与子波褶积合成地震记录,使之与相对应地震记录相似度最大化(相关系数最大或绝对误差最小),逐道外推,直到完成全剖面的波阻抗反演。
叠偏地震记录X(t)可表示为:X(t)=K(t)*W(t) (1)式中W(t)为地震子波,K(t)为反射系数。
从井出发优化波阻抗模型,并正演合成地震记录,使之与相对应地震道相关度最大化,形成反演成果道,选择反演成果道中相关系数达到标准的阻抗模型,以此为基础点建立下一道的初始波阻抗模型,并进行上述优化,直至完成全剖面的反演。
技术特点:1.采用逐道相关外推建立(优选)初始阻抗模型。
2.反演成果纵向分辨薄层的能力较强。
前提条件:要有地震偏移资料,构造沉积解释层位,标准化后的声波和密度测井曲线,如有其它相关资料更好。
优点:逐道外推波阻反演对井的依赖较小,单井时通常也有较高的精度,整体建模反演,适应于岩性剧烈变化的地带,井多时反演精度较高。
缺点:逐道外推波阻抗反演在地震资料较差,岩性剧烈变化地带适应性较差,要调整参数进行试验。
整体建模波阻抗反演井少时反演精度不够高。
求取的孔隙度、渗透率和饱和度参数,可信度相对较低。
油藏工程名词解释
1.1油气藏工程名词解释1.1.1 开发地质名词1 火成岩:由地壳、地幔中形成的岩浆在侵入或喷出的情况下冷凝而成的岩石。
2 变质岩:岩浆岩或沉积岩在温度、压力的影响下改变了组织结构而形成的岩石。
3 沉积岩:地表或接近地表的岩石遭受风化(机械或化学分解)、再经搬运沉积后经成岩作用(压实、胶结、再结晶)而形成的岩石。
沉积岩在陆地表面占岩石总分布面积的75%。
沉积岩与石油的生成、储集有密切关系。
它是石油地质工作的主要对象。
4 碎屑沉积岩:在机械力(风力、水力)的破坏作用下,原来岩石破坏后的碎屑经过搬运和沉积而成的岩石。
例如砂岩、黄土等。
火山碎屑岩则是火山喷发的碎屑直接沉积形成的岩石。
5 化学沉积岩:各种物质由于化学作用(溶解、沉淀化学反应)沉积形成的岩石。
如岩盐、石膏等。
6 岩石结构:包括岩石的颗粒、杂基及胶结物之间的关系。
7 岩石构造:指组成岩石的颗粒彼此相互排列的关系。
8 岩层:由成分基本一致,较大区域内分布基本稳定的岩石组成的岩体。
9 层里:受许多平行面限制的岩石组成的沉积岩层状构造。
10 水平层里:层面相互平行且水平的层里。
水平层里表示沉积环境相当稳定。
如深湖沉积。
11 波状层里:层面象波浪一样起伏。
海岸或湖岸地带由于水的波浪击拍形成的层面。
12 交错层里:一系列交替层的层面相交成各种角度的层里。
由于沉积环境的水流或水动力方向改变形成的层里。
13 沉积旋回:岩石的粒度在垂直向上重复出现的一种组合。
14 正旋回:岩石自下而上由粗变细的岩石结构。
例如自下而上为砾岩、砂岩、粉砂岩、泥岩的组合。
15 反旋回:岩石自下而上由细变粗的岩石结构。
例如自下而上为泥岩、粉砂岩、砂岩、砾岩的组合。
16 复合旋回:中部粗顶底部细的沉积组合。
如顶底为泥岩中部为砂岩。
17 沉积韵律:岩层的成份、结构或颜色等有规律重复出现的现象。
18 沉积相:是指在特定的沉积环境形成的特定的岩石组合。
例如河流相、湖相等。
沉积单元级别划分是相对的。
X851井区测井约束反演方法及应用
3 应 用及 效果 分析
首先 , x 5 对 8 1井 区 的两 口井 进 行单 井 解释 。利 用 测井 在纵 向上 的高分 辨率 , 确 划分 储 层 , 准 对须 四
以孔 隙型 为 主 , 四段 孔 隙低 于 1 , 须 0 含水 饱和 度 值
2 1 年第 1 期 00 5
内蒙 古 石 油 化 工
15 4
X. 区 测 井 约 束 反 演方 法 及 应 用 8 1井 5
赵诚 亮 , 苏媛 媛
( 都理工大学 , 川 成都 成 四 6 05 ) 1 0 9
摘
要: 该文论 述 了基 于模 型 的测井 约 束反 演方 法 的基本 原 理及技 术流程 , 对测 井 约束 地 震反 演 中
提取 , 子波相位的判断 , 子波能谱的稳定性 )③层位 ; 对 比分析 ( 合成地震记录与井旁地震道的相关性 ) ; ④ 模 型 的 建立 ( 利用 丰富 的 地 质 信 息 来 约 束 反演 ) ;
⑤ 地震 资 料极性 的确定 及 参与 约 束反 演 的井 的数 量
等。
新场气田位于 四川盆地川西坳陷中段的一个北
外 推 , 生 初始 波 阻抗模 型 ; 后反 复 对初 始 波阻 抗 产 然 模 型进 行 修 改 , 得模 型 的合 成 记 录 不 断 地 逼 近 实 使 际地 震 记 录 , 时 的 波 阻抗 模 型 便 是 最 终 的 反 演模 此
型。
测 井约 束 反演 技 术是 近十 年发 展 起来 的一 种高
A N i t M 坤/ f c gc / m
% 和 《l 。 度} O 结果 一m
1 4 1 1 — 0-5. -. - 0 - 35 1 1 0 9 5 5 26 77 12 7 25 3 3 — — 00 3 0 2 1 8 o 0 —. 层 — o . 4气
处理与反演
16512000 T89-469 16510000 T87-467 16508000 T89-465 ha2
ha3
16506000
16504000 T93-460n
16502000
T93-457 T89-460
T93-693
16500000
16498000
T93-689 T90-692 T89-689 T88-691 4775000 4780000 4785000 4790000
5 4 3
16504000 T93-460n
16502000
T93-457 T89-460
T93-693
2 1 0
16500000
16498000
T93-689 T90-692 T89-689 T88-691 4775000 4780000 4785000 4790000
16496000
郝家沟组下段II油组砂岩岩芯孔隙度平面图
郝家沟组上段I油组砂岩速度的平面分布图
16516000
16514000
16512000 T89-469 16510000 T87-467 16508000 T89-465 ha2
ha3
16506000
16504000 T93-460n
16502000
T93-457 T89-460
T93-693
16500000
地震子波
时深关系
哈2井地震地质层位标定
哈3井地震地质层位标定
主要内容
•工区概况
•地震高分辨率处理 •地震地质层位标定
•测井约束反演
•储层横向预测
测井约束反演
第一步、地震地质层位标定及子波反演 第二步、利用地震构造解释资料,基于沉积模式建立 地质模型,构造波阻抗的低频分量。在建模过程中, 主要用构造解释资料和测井资料,相当于井上的测井 资料,按构造模式进行插值和外推,得到每条测线的 波阻抗、速度等模型,波阻抗模型可作为反演时的初 始波阻抗和约束条件。 第三步、在时间域进行测井约束波阻抗反演,将波阻 抗转换成速度。通过测井数据波阻抗和速度的交绘图, 得到波阻抗和速度之间的关系,可将波阻抗剖面转为 速度剖面。 第四步、将时间域的速度、波阻抗转换到深度域。
测井约束地震波形反演的同伦摄动法
测井约束地震波形反演的同伦摄动法傅红笋;曹莉;韩波【期刊名称】《地球物理学报》【年(卷),期】2012(055)009【摘要】测井数据和地震数据是地震勘探中两种最重要的资料.测井约束地震波形反演是在非线性波形反演的基础上,利用已知测井资料详细的垂直分辨能力和地震资料均匀密集的水平采样特点,通过迭代反演来求取一个具有较高分辨率的速度参数.本文建立了测井约束反演模型,研究了测井约束下地震波形反演的同伦摄动求解方法.同伦摄动法作为一种新的、求解数学物理中各种非线性问题的有效方法,具有计算速度快、计算精度高的优点.这对于提高反演的精度和效率是十分有益的.为了表征该方法的有效性和稳定性,分别对水平层状介质模型和逆冲断层带模型进行了数值模拟,并与Landweber迭代法相对比,结果表明该算法具有更好的收敛性,能够取得更为满意的反演效果.%The well logging data and seismic data are two most important data in the seismic prospecting. Well log constrained seismic waveform inversion is based on the nonlinear waveform inversion techniques. It makes full use of detailed vertical resolution of well log data and horizontal dense sampling of seismic data. Its core is to obtain an approximate velocity parameter of higher resolution by iterative inversion. This paper established the well log constraint inversion model, and studied homotopy perturbation method for solving this problem. Homotopy perturbation method is a novel and effective method, and can solve various nonlinear equations of various branches in mathematics andphysics. This method has the advantages of fast computation and high precision. This improves the accuracy and stability in inversion. In order to show the effectiveness and stability, we carried out numerical simulations of the horizontally stratified medium and the thrust fault belt model, and compared with Landweber iteration. Numerical results show that the new method improves the convergence, and can produce more satisfactory effect.【总页数】7页(P3173-3179)【作者】傅红笋;曹莉;韩波【作者单位】大连海事大学数学系,大连 116026;哈尔滨工业大学数学系,哈尔滨150001;哈尔滨工业大学数学系,哈尔滨 150001【正文语种】中文【中图分类】P631【相关文献】1.地震波形反演与测井联合的三维建模方法 [J], 罗亚能;黄捍东;王玉梅;慎国强;张达;杨亚迪2.基于模型约束非线性反演与测井随机约束反演技术的联合应用分析 [J], 高哲红;陈秀娟;徐爱丽3.二维波动方程测井约束反演的自适应同伦共轭梯度法 [J], 耿淑彤;冯国峰;王迎;李壮4.二维波动方程速度的正则化-同伦-测井约束反演 [J], 傅红笋;韩波5.地震波形约束的蒙特卡洛——马尔科夫链随机反演方法 [J], 周爽爽;印兴耀;裴松;杨亚明因版权原因,仅展示原文概要,查看原文内容请购买。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
基于模型测井约束地震反演刘凤;张瑞香;张广旭【期刊名称】《内江科技》【年(卷),期】2017(038)011【总页数】3页(P35-37)【作者】刘凤;张瑞香;张广旭【作者单位】中国石油大学胜利学院;中国石油大学胜利学院;中国科学院海洋研究所【正文语种】中文地震数据的生成是建立在褶积原理上,但地震剖面绝不能直接用来当作地质剖面,直接将地震剖面当地质剖面来解释会产生许多错误,基于模型的测井约束反演就是建立在褶积模型基础上,应用波动理论进行地质建模运算,结合测井资料所具有的视高频以及甚低频信息,大幅提高地震资料频带宽度,通过不同的内插外推方法进行运算,可以最大限度的提取地震剖面里蕴含的地质细节,直接用于反映地下地层分布和岩相类型。
精确的基于模型测井约束地震反演可以有效的结合地震横向高分辨率、测井详细的纵向分辨率以及丰富完整的地质模型信息,在储层预测与地震资料综合解释上具有非常良好的应用前景。
随着油气勘探开发的深入,储层预测的难度进一步加大,尤其是在砂岩保护层和碳酸盐岩复杂储层预测方面,随着现阶段数字处理技术的发展,建立在地震与地质结合基础上的地震反演就应运而生了。
由Lindseth先生早在20多年首次提出地震反演技术以来,该项技术获得了长足的进展,新方法和新软件不断涌现。
地震反演技术是综合多种勘探地球物理资料(地震资料、测井资料以及岩性等资料)的一门综合处理分析技术,通过特定的处理手段完成多种资料的有机融合,解决复杂构造条件下的地层地质解释问题。
依靠地震资料的约束外推,结合测井资料详细的垂向信息,利用地震与测井数据的互补性质,以测井约束地震进行联合反演已成为目前“广泛”应用的地球物理技术。
此方法可以更好的处理解释薄层,为钻探之前更多地了解地下地层分布和岩相类型,综合地震地球物理数据以及地质岩心等数据,从而促进了地震与地质两大学科间的联合,推动了相关边缘学科的发展,拓展了储层预测与地震资料综合解释手段,丰富了地震地层学的研究方法。
地震反演从使用的地震资料来分,地震反演方法主要分两种:一种是叠前反演,主要有旅行时反演(Tomography)和振幅反演(波动方程、AVO反演、弹性波阻抗反演);一种是叠后反演,叠后反演主要有振幅反演(波阻抗反演)和波场反演(构造成像)。
从应用的数学方法来分,可分为最优化拟合反演、模拟退火反演、蒙特卡罗反演、基于神经网络的反演、遗传算法反演以及统计随机反演等。
再者按照反演的结果来分,大致可分为:构造反演、波阻抗反演和地震属性的岩石物性多参数反演。
从反演方法上来分,地震反演又可以概括为基于波动理论方程的反演和基于地震褶积模型的间接反演。
目前已经发展较为成熟的是基于地震褶积模型的反演。
从地震波信息研究分析地震反演,叠前反演主要有旅行时反演和振幅反演,叠前应用价值高,但是基于成本等方面,在现场生产中并没有大规模应用;在叠后地震反演的振幅反演和波场反演中,目前应用作为广泛的波阻抗反演就基于叠后振幅反演,由于叠后地震发展程度高,理论方法较为成熟,运算成本低等优点,在储层地震学领域得到了广泛的应用。
本文就叠后反演中的基于模型反演为实例,分析地震反演过程中的重要参量及方法的局限性,以期为储层预测做出有效指示。
地震反演是建立在假定地震波垂直入射假设基础上,综合应用勘探地球物理数据进行的反问题处理,通过已知的测量数据进行实际地质模型反问题的处理,地震资料反演是一项十分复杂的理论与实际问题。
地震波的固有特性决定了所记录的地震波仅有近似的相对振幅信息,地震反演的对象是不完整的离散时间序列, 是一个多解的非线性问题。
叠后波阻抗反演技术发展多年已经成为一套完备可靠的储层预测的手段,并在薄层地质解释以及隐蔽油气藏勘探中发挥了重要作用。
反演是正演的反问题,由已知资料通过反演的算法与地质模型的选择,并将两者通过一定的系数比来匹配,匹配性有很大的主观性,且数学模型的最优化不等于实际地质情况的最优刻画,因此通过对数据进行的反演数学方法的提高只是反演运算的一部分,还要不过实际地质资料的输入进行地震资料多解性的限制,反演的结果就存在随机性、多解性。
地震反演问题的多解性的存在是必然的,通过反演得到的地震响应地质模型也不是唯一的,因此我们要解决地震反演的多解性问题需要利用其他的可以观测数据资料,比如通过测井获得的综合地层速度密度以及物性资料(孔隙度、渗透率、饱和度、孔隙压力等数据)并结合地质学认识(岩心资料、构造认识等等),通过多参数的输入可以相对有效的解决反演的多解性问题,使得反演的解决尽量的接近实际地质问题。
基于模型的测井约束地质反演就是这样一套方法,主要通过基于声波波动方程的地震资料(理论零炮检距叠后地震资料),并结合测井资料的高频信息,拓宽地震信息的频带宽度,使之部分达到井中约束资料的频率范围;综合利用了多种资料,通过合理的参数设置,使得多种资料有效融合,最终更精细的描述地下地层的沉积状况,帮助我们直接找油找气。
基于模型反演皆是建立在地震记录褶积模型基础上,地震子波与反射系数的褶积得到地震记录,地震反射波S等于反射系数R与地震子波W的褶积,而某界面的垂直反射系数R就等于该界面上下介质的波阻抗差与波阻抗和之比,地震波阻抗就是速度与密度的乘积。
其中:是s(t)为地震记录,r(t)为反射系数,w(t)为地震子波。
第i个界面的反射系数定义如下:基于模型的测井约束反演就是建立在褶积模型基础上的,地震数据的生成基于褶积原理,由于地下不同地层存在密度以及速度差,体现在这里就是波阻抗差异,当地震波传播到有波阻抗差的地层界面,就是发生反射形成地震反射波,形成地震剖面上的地震层位,地震剖面绝不能当作地质剖面来对待,如果把地震剖面直接当作地质剖面解释会产生许多错误,因此结合地震反演理论进行运算可以有效的解决地质解释中的假象以及隐藏的地质问题。
(1)地震数据要求。
作为反演的主体数据地震数据,要求必须要是高保真度的资料,不是通常用于解释的增益地震数据需是经过精细处理的零偏移距的叠后纯波数据,地震数据是经过偏移处理的没有多次波的数据,每道地震只与下方地层的反射序列相关,没有AVO变化,通常假设波阻抗反演输入的地震数据其振幅信息反映了地下波阻抗变化情况,地震剖面没有多次波和绕射波的噪声分量,或者说噪声是随机的与地震无关的白噪音。
在地震资料处理时就应该考虑反演的高保真数据的要求,在反褶积、噪声切除,尤其是多次波的处理方面,最终得到更加真实反映地质反射特征的真振幅剖面,尽量做到高分辨率、高信噪比,使得地震更加有效的反映地下波阻抗变化,输入地震资料的品质直接影响了基于模型地震反演的分辨率精度,对真实反映地层的精度影响最大,地震资料的有效频带范围基本确定了地震反演所能识别储层厚度的能力。
(2)测井数据要求(用于反演的测井数据处理和要求)。
用于基于模型测井约束反演的测井资料最主要的是声波测井资料,在进行反演时对输入的测井数据要求必须要进行归全区归一化或者统一标准处理,需要经过环境校正、深度校正等预处理,要去除测井资料在不同采集时期以及不同测井工具之间的系统误差,以及采集过程中的泥浆性质以及井眼状态等多方面因素的影响。
通过前期的预处理,有效的排除系统误差、人为因素以及其他非地质因素的影响,得到建立初始模型的基本测井资料,测井资料是反演高频信息输入最主要的影响因素,因此高分辨的高保真的测井数据是进行反演的最关键因素之一,有效的保证基于模型反演薄层解释准确性的关键,进一步保证了反演储层地质参数的可靠性。
基于模型测井约束地震反演的关键步骤可概括为:①地震解释层位及解释目的层的标定,建立初始地质模型;②提取子波,试验多种方法提取地震子波;③必要的特殊地震数据处理(剩余地震相位校正处理以及反褶积处理);④建立单井波阻抗模型精化及取低频,结合地震解释成果建立整体反演初始模型;⑤多井数据输入以及参数设置,迭代运算;⑥加低频信息等输出反演成果。
其中最重要的关键步骤是初始地质模型的建立以及井震输入影响因子系数。
地质模型的建立是利用层序地层学解释完成的地质格架和利用不同插值方法(用克里金、平方反比等插值方法)等不断调整系数产生的井震联合响应模型,完成合成记录的建立,并与实际地震记录对比,反复调整以减小阻抗模型与实际地震反射的差距,最终完成数据反演运算,产生相应的高分辨率反演成果(波阻抗剖面等)。
首先基于模型测井约束地震反演最主体的数据是地震数据,地震的分辨率直接决定了反演的精度,如果地震的分辨率低,反演结果也不会无中生有大幅提高分辨率,尤其是复杂地层的薄互层沉积识别时,地震的分辨率是决定反演结果的最大影响因素。
另外通过结合新技术新数学方法,在地震数据分辨率一定的情况下,将地震与测井数据通过更合理先进的数据方法有机结合,可以从一定程度上大幅提高反演结果的分辨率,在复杂层状沉积地层中发挥重要作用,可用于指导油气,新方法新技术以及精确的测井资料输入也是影响反演的重要因素。
地震反射是通过s(t)与w(t)在时间域的褶积,就相当于他们在频率域中,反射系数相位谱r和子波相位谱w的乘积。
地震资料中主要包含相对丰富的低频信息,缺少视高频和更低频段的信息,通过基于模型测井约束地震反演,可以通过反演算法将测井的高频信息恢复到反演解释剖面上,恢复这些地震剖面丢失的信息。
无论多么成熟的反演方法均有其局限性。
因此利用高分辨率的地震资料、精确处理的测井数据,并输入多种因素控制的精细的地质模型是取得准确高分辨率反演结果的关键因素。
(1)获取高质量的地震数据。
基于模型地震反演的分辨率首先取决于地震资料本身的分辨率,提高反演结果的分辨率,首先考虑的是通过新技术新手段提高地震资料分辨率,可以采用一些新的采集方法及处理技术,如采用高分辨数据采集、分频或拓频处理、针对复杂地层进行的地震主频的处理技术等技术,以及采用一些新的采集手段提高资料的信噪比,更大的获取有效信号,采集更宽频带信号,以及更好的数据处理分析方法。
地震反演处理过程中无法对噪声部分进行衰减处理,因此用于反演的数据必须保证有较高的信噪比,这样才能保证数据反演结果的高质量,因此采集高质量的有效信号,采用先进的处理方法是获得高质量反演结果的基础。
(2)收集高质量的测井资料。
基于模型测井约束地震反演,高分辨率经过准确校正及归一化处理的测井资料也是反演的关键,测井资料可以提供地震资料本身所不具备的甚低频和视高频段的信息,拟补地震资料在低频以及高频方面的先天缺陷。
对进行反演所采用的测井资料要求参数完成齐全,数据质量要求较高,且最好首先保证有声波时差测井资料,如声波时差数据较差,可采用电阻率数据以及其他测井数据进行反演。
对于采用的反演井资料要进行归一化处理,保证所有数据可靠前提下,才可以应用到基于模型反演的运算中。