1、数据加载（1）地震数据加载 进入JASON软件主窗口菜单:Datalinks---Seismic/Property data： 进入下一级界面： 该Parameters菜单中每个步骤依次执行
该Parameters菜单中每个步骤依次执行 (a)Creat/edit segy format definition… (b)Select/ edit transport parameter… (c) Desired JGW format: (32bits or Floating point) (d)Existing files….
Transport --- (2)层位数据 Datalinks---Horizons… 出现下一级子菜单 (a)Parameters---Select/edit transport paramenters…. 选择层位，断层等所在目录即可。 (b)Existing files (same to seismic data) （3）井位数据加载 Datalinks--- wells…
Input— (a)Input---Set template file…出窗口（未抓） (b)Input—Select files… Edit---- (a)Edit—Edit template…出窗口
在Generic log –column (select all) 一个接着一个编辑，选中（击黑）点Edit header—出现 LAS edit log header对话框，在JGW type中根据曲线的类型选择相应说明，并修改单位。 深度为Measured depth, 单位为m；声波为P-sonic单位为：数据在200以内为英制，大于200为米制。Edit—Existing files… (same to seismic data) Transport—出窗口 （a）Edit….出现编辑坐标的窗口，输入X，Y坐标 （b）Import wells..加载井数据。 依次直到加载完毕所有井数据。 在底图中可以查看井数据的加载是否正确
第一章数据的加载 jason是目前最常用的反演软件，它操作上的特点是它需要什么数据或参数就给它什么数据或参数。下面是它的主窗口（图1）。 图1 因为作反演之前已经将坨163区块进行了构造解释，所以可以直接从lanmark中将地震、测井、层位数据导入jason中，操作比较简单。步骤如下： 一,选择工区（即一个文件夹） 主窗口——File——Select Project（图2），弹出图3。选择一个工区，ok。 图2 图3 1. 数据的导入 主窗口——Datalinks——Landmark——Landmark Link（2003）（图3），弹出图4。
图3 图4 2. 工区的选择 File——Seisworks project：选地震工区t163，ok。（图5） 图5 File——Openworks project——选SHNEGCAI, 选井列表t163，ok。（图5） 此时，图5 窗口的状态栏将会发生变化，以上选择的工区将会显示。（图6） 图6 3. 地震数据的导入 Select——Import——Seismic/property data（图7）,弹出图8。选cb 3dv（纯波数据，作反演时一定要用纯波数据），ok。
图7 图8 4. 层位数据的导入Select——Import——Horizons，选择反演时需要的层位和断层（图9）。 图9 5. 井数据的导入 Select（图7）——wells，弹出图10。选择需要的井，ok。
图10 E： Transport——Import，以上所选的landmark中的数据将传入jason中。 图11 第二章合成记录的建立 在jason上建合成记录的特点是精度高，但随意性大。建立合成记录的步骤是:井曲线、地震数据、子波的加载，子波的编辑和评价，合成记录的生成和编辑。 1. 井曲线、地震数据、子波的输入 主窗口——Analysis——Well log editing and seismic tie（图1），弹出图2。 图1
Jason反演Wavelets中文手册13．Wavelets Wavelets应用程序的目标是： 利用各种技术评估地震子波： 单独从地震数据 利用来自于一口井或多口井的反射系数 利用地震数据的部分（partial）叠加或multiple partial 叠加 定位井轨迹以优化井和地震数据的关系 横向变化时，进行子波插值 为建模生成合成子波 为得到零相位地震数据计算反褶积 评估Q和产生基于子波的Q 13.1 快速开始 子波估算作为科学，带有艺术性和经验性。应用一个正确的子波是获取满意的储层描述结果的关键。 为了运行Wavelets，必须下面的数据： ◆地震数据 ◆阻抗井或者选择的位置（Impedance well(s) or selected locaitons） ◆选择一个时窗，在时窗内来执行估算 ◆每口井选择一定数量的地震道或者位置用于估算
13.2 菜单导航 13.2.1 File 菜单 13.2.2 Input 菜单 Seismic mode 子波能够通过Zero-offset（ZO）地震数据和A V A非Zero-offset地震数据估算而来。选择你喜欢工作的模式。 Preferred Types 选择喜欢的查询输入测井曲线的数据类型。对于Zero offset模式，它是P-Impedance类型。对于A VA模式，可以是P-Sonic，S-Sonic和密度类型。Seismic Seismic 菜单有两部分。在第一部分中，指定了地震数据文件。你必须提供一个包含2D或3D数据的地震数据（.mod）。地震数据必须代表真正振幅偏移反射信息。在Zero-offset模式中，合成地震记录是通过子波与零偏移反射系数曲线褶积而来的，其中反射系数曲线来源于输入的P-Impedance曲线。在AV A模式中，角度依赖的反射系数来源于P-Sonic，S-Sonic和密度井曲线，所使用的角度是是在输入的A VA地震数据文件中指定的。 在第二部分中，需要设置地震数据的A V A参数。这个对话框要在其他程序中使用，所以子波参数是被灰化的。
电阻率法和电磁法联合反演进展 地球物理反演过程是一个将观测数据转换为地质-地球物理模型的过程，其首要解决的问题之一就是尽可能的减少地球物理场的多解性，从而得到一个可靠的地质模型。然而，地球物理方法所采集到的数据集通常受到地表干扰或方法自身缺点的影响，使得反演结果难以全面的认识地下地质概况。直流电阻率法( DirectCurrent Ｒesistivity，DC) 在浅层勘探中具有较高的灵敏度，但是由于装置展布的限制和深部传导电流受到极大的削弱，因此很难获取到深部的有效信息。电磁法，受到地形起伏和近地表横向电性分布不均匀，易产生静态效应，极易混淆真假异常，增加地质解释的难度。由此可见，单一地球物理方法的反演容易造成地下地质解释的模糊性。在综合地球物理中同时运用多种地球物理数据进行同步、顺序、剥离、伸展等方式计算同一地质体的物性特征和几何展布称为联合反演，联合反演是综合地球物理工作中不可或缺的一种重要的定性和定量解释工具。结合多种地球物理数据进行联合反演能够有效的减少模型解空间。这主要是由于: (1) 不同地球物理方法获得的解空间不尽相同。利用不同地球物理方法的优势，一种方法中的零空间可以通过联合反演在另一种方法中得到补充。 (2) 不同地球物理方法测量的物性参数不尽相同。地下岩矿石包含了多种物理属性，通过在同一区域对不同物性参数的测量，从不同的侧面提高对该区域岩性及其范围的识别。 (3) 不同地球物理方法所受的干扰因素不尽相同。某种方法部分受到强干扰的数据可以用另一种方法的数据经行校正，有时比单一方法采集更多数据更为有
效。 自1975年V ozoff和Jupp首次提出进行MT和DC数据联合反演(V ozoffandJupp，1975)以来，中外学者对联合反演从理论到应用都做了许多研究:Sasaki(1989)、Sharma(SharmaandKaikkonen，1999；SharmaandVerma，2011)等对大地电磁和直流数据进行了联合反演研究。Zeyen和Pous(1993)进行了带有先验模型的磁法和重力联合反演研究。Hering等(1995)详细论述了DC与浅层地震的联合反演算法。Bosch等(2006)运用蒙特卡洛方法进行了三维重磁的联合反演研究。Moorkamp等(2011)提出了一种适用于地震、MT和重力的三维联合反演框架。国内学者对联合反演也做了许多贡献，于鹏等(2007)利用改进的全局寻优的快速模拟退火算法，实现了重力和地震资料的约束同步联合反演。万玲等(2013)提出地面磁共振MＲS与瞬变电磁(TEM)联合反演方法，提高了解释结果的准确度。彭淼等(2013)研究了基于交叉梯度耦合约束的大地电磁与地震走时资料的三维联合反演算法。陈晓等(2010)采用非线性模拟退火方法实现了加入有效模型约束的大地电磁与地震的同步联合反演，使反演的解更实际更稳定。陈华根等(2012)在实际资料处理中应用了MT和重力的模拟退火联合反演并取得了较好的效果。刘彦等(2012)在对国内外大地电磁与地震数据联合反演的研究现状分析的基础上，总结了电震联合反演算法的类型。 回顾联合反演的发展历史，联合反演是综合地球物理定量解释的重要工具，是未来地球物理学的一个发展方向。从总体上来说国内的联合反演主要集中在重磁电磁与地震方法的结合，而更具有合理性的基于相同物性基础的联合反演较为少见。
RockTrace 同时AVA 约束稀疏脉冲反演 自从2000 年秋季推出以来，RockTrace 就对行业应用和集成PSTM 地震数据的方法产生了重大影响。它是唯一的定量集成测井曲线弹性岩石物性和AVA 地震数据的技术，可以生成标定的定量岩石物性三维数据体。 RockTrace 以InverTrace Plus技术为基础，并将该技术扩展到了AVO 域。在InverTrace Plus中，应用的约束以波阻抗(Z p) 为依据。在RockTrace 模块中，目标是在波阻抗之外求解出横波阻抗(Z s) 和密度，因此，对这三个参数分别地设置约束条件。在并行处理过程中，能够生成的弹性参数类型组合为： ?纵波阻抗、横波阻抗和密度 ?纵波阻抗、纵/横波速度比和密度 ?纵波时差、横波时差和密度 ?纵波速度、横波速度和密度 和InverTrace Plus一样，应用全局模式时，一个空间控制项被加到目标函数上，同时，一个相当大的地震道数据体被整体转换。RockTrace反演应用多个不同角度叠加道集的地震数据体，并且生成三个弹性参数数据体。该算法是InverTrace Plus模块整体多道反演算法的扩展。它是以一个地震数据体作为输入信息，并且只生成一个阻抗数据体（上述数据体之一）。通常，RockTrace 算法可以产生三个弹性参数数据体，并且保留了许多其它模块的原理和约束算法。 这是业界一项独特的技术，具有以下技术优势： ?反演得到的弹性参数是岩石的真实属性，与储层属性有关。 ?当采用Knott-Zoeppritz 方程式选项时，无任何物理近似假设。 ?允许振幅和相位随偏移距变化。通过对每隔输入部分叠加数据体，计算唯一的子波来实现反演。 ?在反演过程中，弹性参数可以直接进行各自约束。 ?岩石物理关系可用于约束弹性参数对。 ?由于所有输入数据必须和单一的输出模型相一致，降低了噪音的影响。 ?最终的弹性参数模型可重构输入地震数据，这也是反演优化算法的一部分。 ?由于处理过程是整合了所有数据并与单一模型保持一致，所以质量控制分析更简单，而且结果更加机关内却。 RockTrace 可同时求解多个褶积方程，包括纵波速度、横波速度和密度。此外，S 代表地震角度叠加道集，KZ 表示Knott-Zoeppritz 反射系数方程，W 表示角度叠加道集的相应子波。由于RockTrace 使用完全的KZ 方程，因此可支持P-P 波形数据、P-S 转换波形数据，或者二者的结合。RockTrace 内的组件 RockTrace 由三个模块构成： ?P-P 和P-S 弹性阻抗曲线生成器 ?同时AVA 约束稀疏脉冲反演 ?垂直数据对齐 另外，在Wavelets模块中，具有RockTrace 软件许可，就可以运行多角度叠加数据子波估算工具。
Jason软件培训资料 Jason软件集合了油气勘探开发不同阶段的储层预测和油气藏描述技术，它致力于各种资料、各种认识的全面综合，提供符合各种资料、各种认识的储层预测和油气藏描述结果。指导油气藏的勘探和开发，提高钻井成功率，降低风险。 主要模块
Jason软件是一套综合应用地震、测井和地质等资料解决油气勘探开发不同阶段储层预测和油气藏描述实际问题的综合平台。其中子波估算（Wavelets）和层位标定、地质框架模型（Earthmodel）、地震反演（Invertrace、Invertraceplus）、测井反演（Invermod）、地质统计模拟（Statmod）和数据分析变换（Functionmod）是主要模块和关键技术。下面根据实际工作步骤来介绍Jason软件的主要模块和关键技术及应用注意事项。 一、数据加载 数据加载顺序为地震→层位→测井→其它（如人文、子波等）；输出可根据需要有选择性地输出。 注意事项： ●地震数据类型（是2D还是3D）、线道号和XY坐标在SEG-Y道 头中的正确位置、输入数据的字节数（至少为16位）。 ●井数据输入文件的格式与所选的格式模板文件必须一致包括输入 文件本身的声波和密度的单位（us/ft，us/m，g/cm3，kg/m3）、模板 文件中深度的类型（测量深度、TVD等）和单位(m，ft等)。 二、子波估算和层位标定技术 这部分工作是通过Modeling下的Wavelets…和Analysis下的Well log editing and seismic tie…两个模块完成的。通过子波估算和井曲线编辑的交互迭代，由井旁地震道和井中的阻抗曲线估算出与地震最佳匹配的地震子波。并实现子波估算、合成记录的制作和层位标定。 其技术特点是：同时估算子波的振幅谱和相位谱；子波估算和层位标定同时完成；方法多样，可处理有井和无井、单井和多井、直井和斜井；质量控制手段多样。 子波估算和层位标定技术的方法如下： 1）计算理论子波（如Ricker）(Wavelets…→Edit→Create synthetic wavelet... )。
第四篇约束稀疏脉冲波阻抗反演 1、概述： 约束稀疏脉冲波阻抗反演方法（简称CSSI ）是目前应用比较广泛的一种方法，它是以趋势约束稀疏脉冲反褶积算法为基础。对地下的反射吸收序列作如下假设：它的分布是稀疏分布的，而不是像其它的反演方法那样假设其是连续分布的。在反演过程中，它是通过寻找一个使目标函数最小的反射系数脉冲数目来得到波阻抗数据的，此时合成的地震记录可以和实际地震记录匹配的最佳。 2、基本假设和方法原理 约束稀疏脉冲反演（CSSI ）是一种递推地震波阻抗反演方法，基于范数反褶积和最大似然反褶积。它假设 地层的反射系数序列是由一系列稀疏且不连续的大反射系数和高斯分布的小反射系数叠加而成。从地质角都解释，大反射系数代表地下不连续界面、岩性分界面。首先由约束稀疏脉冲的基本假设，得到目标函数： )1ln()(2)ln(2)(1 )(1 1 22 1 22 λλ----+= ∑∑==M L M K n N K r R J L K L K （1） 2R 为反射系数均方值，2N 为噪音均方值；第K 个采样点的反射系数为r(K),n(K)为噪 音；M 为反射层数；L 为总采样数；λ为既定反射系数似然值；可经迭代计算出反射系数。其次用最大似然反褶积反演得到反射系数。采用递推算法结合初始波阻抗趋势，得到波阻抗模型。 ) () (1) 1()(i R i R i Z i Z +-= （2） Z(i)为第i 层的波阻抗值；R(i)为第i 层的反射系数。 最后约束稀疏脉冲反演依据目标函数逐一对地震道计算的初始波阻抗进行调整， 并调整反射系数，得到目标函数： Z L a d s L r L F q p ?+-+=-11)()(λ （3） r 是反射系数，Z ?是阻抗趋势差；d 是地震道；s 是合成地震道；λ是残差权 重因子；a 是趋势权重因子；p,q 是L 模因子。式（3）右边分别是反射系数的绝对值和，合
前言 JASON地学综合研究平台(JASON GEOSCIENCE WORKBENCH)为用户提供的跨越地震、地质、测井资料的综合分析研究工具，它可满足油气勘探开发不同阶段对储层的油气藏定量研究的需求。JASON把不同学科的有效信息的融合作为客观存它的基石，最大限度地利用不同信息的优势，为用户提供符合不同学科信息的客观可靠的油气藏参数模型。 JASON软件的重要特点就是随着越来越多的非地震信息（测井、测试、地质）的引入，由地震数据推演的油气藏参数模型的分辩率和细节会得到不断的改善。用户可根据需要，由JASON的模块构建自己的研究流程。 主要模块及功能如下： Enviaonment-Plus 运行环境及分析工具 数据输入与输出（地震segy格式、测井las格式和层位ASCII格式） 各种数据，各种方式的显示（井、层位、地震等的2D/3D显示） 合成记录标定 2D/3D解释（地震体与属性体解释） 交会图/直方图分析 三维立体显示与三维（地质/储集）体自动解释 沿层、层间属性提取、沿层属性切片 层位数据计算（平滑、加/减、拟合等） 地层异常检测 处理工具包（重采样、滤波、互相关等） 等值线 Wavelets 子波估算 用多种技术估算地震子波（无井估算地震子波、单井或多井估算地震子波等） 空变子波 理论子波 估算子波的振幅谱与相位谱 计算平均子波 VelMod 速度建模 建立三维速度模型（用均方根速度，逐层的层速度编辑与平滑） 时深转换，深时转换 提供阻抗的低频模型 EarthModel 地质框架模型 构建以层为基础的地层框架模型 生成以地层框架模型为基础的测井曲线内插模型 提供用于地震反演的低频模型 InverTrace 地震反演的储层与油气藏描述 提供可靠的地震反演声阻抗数据体 在地震反演声阻抗数据体上解释，可提高解释的精度和可靠性 预测产层有效厚度和平均孔隙度 InverMod 基于测井的精细储层油气藏描述
中国矿业大学（北京） 研究生课程考试试卷 考试科目油气藏评价与开发技术 考试时间2017年2月26号 学号TSP1600201017 姓名李泽辰 所属学院地球科学与测绘工程学院 类别（硕士、博士、进修生）硕士
评分及评语（由阅卷人填写）： 阅卷人签名：任课教师签名：时间：时间：
测井约束地震反演进行预测的原理和过程 摘要:测井资料具有很高的纵向分辨率, 但在横向上探测范围却很小;地震资料在横向上可连续追踪地层信息，但其纵向分辨率却受到记录频带的限制而远远低于测井信息。采用普通地震资料进行反演，通常受地震频带宽度的限制，其精度和分辨率都不能满足油田开发对薄砂层的要求。测井约束地震反演从地质模型出发将测井与三维地震资料有机结合，以地震解释的层位或层序为基础，从井点出发内插外推，形成初始波阻抗模型，再通过不断修改、更新模型，使模型正演合成地震资料与实际地震数据吻合度最高，则该最终模型的数据便是反演结果。关键词测井约束地震反演储层预测 从20世纪80年代初开始,广泛利用波阻抗反演技术进行储层预测收到了很好的效果; 90年代,在波阻抗反演基础上发展起来的测井约束地震反演技术应用日趋广泛,并有很多成功的实例。测井约束反演作为一种储集层横向预测技术,以测井资料丰富的高频信息和完整的低频成分补充地震资料有限带宽的不足,综合地质认识等信息作为约束条件,可以得到高精度的波阻抗资料,弥补了普通的地震剖面不能分辨薄层砂体的不足。 1 井约束条件下储层反演的基本原理 井约束条件下的储层反演是基于褶积模型的地震道反演, 从测井曲线出发, 用井的波阻抗模型作为反演的初始模型,进而求取反射系数序列, 从井旁地震道提取地震子波,用所求得的地震子波与反射系数序列褶积生成合成记录, 比较合成记录与地震道的误差。当误差达到最小时,其波阻抗模型就作为井旁道反演后的波阻抗曲线, 按地震解释层位建立模型, 把测井声阻抗曲线沿此模型横向外推, 充分利用测井资料丰富的高频信息和完整的低频成份补充地震带宽的不足,
反演技术 前言 一. 反演的概念、目的 二. 反演的发展历史及趋势 三. 反演的基本方法 四. 地震反演难题的解决方案 五. 反演的实质 六. 反演的基本流程 七. AVO反演处理简介 地震、测井、钻井是石油工作者认识地下地质构造、地层、岩性、物性、含油气性的最重要的信息来源。虽然测井、钻井仅能提供井孔附近的有关信息，尤其是有关岩性、物性、含油气性的信息，但是这些信息往往具有很高的分辨率，可信度、准确性，能确切地指出含油气层的位置，定量化分析与储层、油藏有关的参数。然而一个油气田勘探、开发方案的设计、实施、调整仅靠测井、钻井资料是远远不够的，
必须与地震资料相结合进行综合分析才能取得良好效果。 地震资料的分辨率虽然远远不及测井、钻井，但是随着地震勘探技术的发展，从光电记录、模拟记录到数字记录，从二维到三维，地震资料的信噪比、分辨率、成像的准确性都获得了极大的提高，由于地震资料包含大量地下地质信息，覆盖面积广，具有三维特性，所以这项技术的使用越来越受到石油工作者的重视，如何利用地震资料研究地下地质构造、地层？如何进行储层预测、油藏描述？如何进行油藏、含油气层的预测？ 这些问题促使地球物理学家、地质学家开发应用了一系列地震资料特殊处理技术，如地震资料反演技术、地震属性分析技术、AVO 分析技术，这些技术充分利用测井、钻井、地震的长处，使人们对地下储层、油藏的研究从点到面、从二维到三维、从三维可视化研究到油藏动态监测、从定性研究到定量化研究，大大提高了钻探成功率，有效地指导了油田开发，为提高油田最终采收率起到了积极的作用，因此地震技术被列为二十一世纪石油工业发展的首要技术，相信地震资料特殊处理技术（地震资料反演技术、地震属性分析技术、AVO分析技术）也必将在我国油田勘探、开发中起到越来越重要的作用。 一. 反演的概念、目的 地震资料反演技术就是充分利用测井、钻井、地质资料提供的丰富的构造、层位、岩性等信息，从常规的地震剖面推导出地下地层的波阻抗、密度、速度、孔隙度、渗透率、沙泥岩百分比、压力等信息。那么如何理解这个概念？还是让我们看看什么是正演吧！ 1.正演的概念 如果我们已知地下的地质模型，它的地震响应如何？通过模拟野外地震采集，得到单炮记录，再通过速度分析、动校正、叠加、偏移得到合成剖面这一过程就是正演。
几种常用的反演方法综述 一、稀疏脉冲反演(C onstrained Sparse Spike Inversion) 1、原理： ①首先假设地下地层的波阻抗模型所对应的反射系数序列模型是稀疏的，即由起主导作用的强反射系数序列和具有高斯背景的弱反射系数序列叠加而成。 ②将地震记录与子波进行稀疏脉冲反褶积得到地层反射系数，一般是使用最大似然反褶积求得一具有稀疏特性的反射系数序列Ri。 根据①的假设可以导出最小目标函数： R（K）为第一个采样点的反射系数，M 为反射层数, N为噪音变量的平方根，L 为采样 总数，? 根据目标函数，对每一道，从上到下推测反射系数的位置点，判断反射系数的幅值大小。如此反复迭代修改每个反射系数的位置和幅度，使最后的修改误差最小符合似然比值的判别标准即可，这样就完成了一道的反褶积，得到该道的反射系数的分布。 ③通过最大似然反演导出波阻抗Zi 反演公式为Zi=Zi-1*[(Ri+1)/Ri]. 具体的计算方法是稀疏脉冲序列每次建立的反射系数为一个脉冲，然后在地震资料中提取子波与初始反射系数进行褶积，得到一个初始合成地震记录，并用此合成地震记录与实际地震纪录作对比得到他们之间的残差，利用这个残差的大小来修改反射序列中脉冲的个数再次进行褶积运算，得到新的合成地震记录，再与实际地震资料对比，就这样循环迭代，直到残差达到最小，最后得到一个与实际地震资料达到最佳逼近的合成地震记录，获得宽频带的反射系数。
图1 稀疏脉冲反演每次建立反射序列为一个脉冲，增加脉冲进行循环迭代约束稀疏脉冲反演采用的是一个快速约束趋势的反演算法，约束条件主要是波阻抗趋势和地质控制，而波阻抗趋势又是由解释层位和断层来控制的，从而可以把地质模式融入进去得到一个宽带的结果，恢复地质信息中缺少的低频和高频成分。 约束稀疏脉冲反演的最小误差函数是： 第二项为原始地震道与合成地震道的均方差的总和； 第三项为趋势协调的补偿 i 是地震道样点号；di是原始地震道；Si是合成地震记录；ri 为地震道采样 点的反射系数；ti是波阻抗趋势；Zi是地震道采样点的波阻抗值，介于井约束的 最大和最小波阻抗之间；ɑ是趋势最小匹配加权因子，一般情况下ɑ=1；p、q是L 模 因子，一般情况下p =1，q=2是调节或平衡因子，与信噪比大小有关。 根据目标函数，我们可以看出：在约束稀疏脉冲反演中，反射系数的稀疏、原始地震道 与合成记录的残差最小这两项是相互矛盾的。?值的大小反映了合成地震道与实际地震道匹 配程度的好坏。若?值太大，着重强调残差最小，即过于追求合成地震记录与原始地震记 录的吻合程度，导致会把一些噪音当作有效信息出现在反演剖面中，同时，由于反射系数的
StatMod MC入门手册 Chapter 1.工作流程 Chapter 2.基本的输入输出数据 输入数据 输出数据 岩性实现岩性概率体 属性实现 地质统计学参数 岩石物理分析 地层网格模型地震数据 测井曲线 ……………………...5% ...….………………..15% ..……………………..5% ……………………...50% ……………………...10% ….………………….15% 百分数表示每个步骤所用时 间占整个项目时间的百分比Stage 4:反演 Stage 2:地质统计学参数分析 Stage 3:模拟 Stage 5:协模拟 Stage 1: 项目准备 Stage 6: 不确定性分析与风险评估
Chapter 3.详细操作步骤 操作步骤以StatMod MC培训数据为例 第一步.首先完成一个高质量的叠后CSSI反演 这一步的目的是为地质统计学提供一个好的研究基础, 这个“好”主要体现在： (1)好的井震标定, 目标区的相关值达到0.85以上; (2)好的叠后反演结果, 用来质控地质统计学模拟和反演结果, 是地质统计学反演结果横向预测准确 度的参照物; (3)利用叠后反演结果进行砂体雕刻, 对目标区的岩性展布、比例有一个总体上正确的把握, 这些认 识都是地质统计学的初始输入。 （说明：在提供的培训数据中已经为用户做了以上准备，用户可以从主界面中打开该培训数据所在工区, 然后用Map View看工区底图，用Section View查看地震数据、叠后CSSI反演数据、地质框架模型与 层位数据以及井数据与子波 , 并用Well Editor检查井震标定情况） 第二步. 数据准备 ●●井曲线重采样 这一步将测井数据重采样至地质微层采样间隔，具体操作为： (1)JGW主界面→ Analysis→ Processing toolkit； (2)Input→ Data selection→ Data type：选Well, 点击Input file(s)右边List选择任意井(可以选多井), 然后在弹出的界面Select logs中选择任意井曲线(可以多选),点击OK退出; (3)Parameters→ Resample log, 在弹出界面Processing toolkit中填写重采样间隔(注意s 与ms单位), 点击OK退出; (4)Output→ Define process, 从Select from中选择Resample log, 点击??输入到右边的Process里面; (5)Output→ Generate, 在弹出的界面中填写输出路径和输出文件名,然后点击Generate,开始计算重 采样的曲线。 ●●计算岩性曲线：这里说的岩性曲线和测井解释的岩性不同, 而是指能通过弹性属性(主要指纵波 阻抗) 区分出来的岩性, 所以一般可分两到三种岩性。练习里通过交会图的方式教授如何划分岩性。 ●●地震数据：线、道间距尽量一致，如果不一致需要进行跳道处理，避免往反演结果中引入各向 异性。 ●●地质框架模型：微层采样间隔与需要分辨的薄层厚度要匹配。 第三步. 地质统计学参数分析 这里说的地质统计学参数主要指三个参数：概率密度函数(probability density function, 简称pdf，描述某一属性在空间的概率分布情况)、变差函数（描述某一属性随距离的变化，是距离的函数）、云变换（描述两个属性之间的相关关系）。下面先介绍前两个参数: ●●概率密度函数 (1)JGW主界面→ Modeling→ StatMod MC, 弹出StatMod MC主界面，主要填写内容如下：
Jason 约束稀疏脉冲反演在苏里格气田 下石盒子组砂岩储层预测中的应用 X 刘盼喜,时义刚 (中国石化中原油田分公司物探研究院,河南濮阳 457001) 摘 要:苏里格气田为鄂尔多斯盆地岩性油气藏,充分利用高分辨率地震资料,在经过测井曲线标准化后,在储层精细标定基础上进行约束稀疏脉冲反演(C SSI)对下盒子组砂岩进行储层预测,优选滚动探井,钻遇32.8m/L 储层并获高产气流,取得了滚动勘探的成功,为同类岩性油气藏的勘探开发提供了参考和借鉴。 关键词:苏里格气田;岩性气藏;滚动勘探;约束稀疏脉冲反演 中图分类号:P 631.4 文献标识码:A 文章编号:1006—7981(2012)01—0135—02 苏里格气田位于长庆靖边气田西北侧的苏里格 庙地区。鄂尔多斯盆地陕北斜坡的西北侧.目前已 发现上古生界二叠系下石盒子组盒7、盒8、山西组山 1、山2及下古生界奥陶系马家沟组马五5、马五6等 多套含气层段。苏里格气田是一个低孔、低渗、低压、 低丰度、大面积分布的岩性气藏(图1) 。图1 研究区位置图 1 资料解释 1.1 合成记录及层位标定 层位标定是储层横向预测的关键。它将地震资 料与钻井资料相互沟通,使二者之间建立一个准确 的对应关系。为后续的地震地质层位追踪解释、储层 预测等提供可靠目标。利用Welledit 模块作合成记 录,针对目的层段从波形、频率、相位等特征与地震 资料对比,对合成记录进行微小的压缩、拉伸,使其 与地震相关性最佳(图2)。 1.2 层位解释 利用Discovery 解释系统,以储层标定为基础, 对目的层,进行对比追踪精细解释,获得所需的地震 层位数据,使其作为目标反演中保持地震数据能够 在同一低频背景下运算的约束条件,从而提高初始模型精度。分别编制了盒8底与山西组底的深度构造图(图3) 。 2 约束稀疏脉冲波阻抗反演(C SSI)2.1 测井响应特征盒8砂岩储层主要电性特征(表1)为“三低一大”即低自然伽玛、低密度、低声波时差;大幅度自然电位异常。储层速度特征表现为泥岩和致密砂岩有较明显的速度差异,砂岩具高阻抗特征,泥岩阻抗值较低(表)。 135 2012年第1期 内蒙古石油化工X 收稿日期5 作者简介刘盼喜(6),汉族,工程师,年毕业于石油大学勘探系,现工作于中原油田物探研究院,从事构造解 释、储层预测工作。1:2011-11-1:199-1992
第五期用户通讯 随机模拟与随机反演 地震 属性随机模拟、岩性指示随机模拟与随机反演： 1. 需要输入的数据： 1、wells (logs and tops) 2、seismic.mod (?=4ms …) 3、horizons 4、wavelets (?=4ms …) 5、aitm.mod (?=4ms …) 6、solid model (with TDC) 2. 主要工作流程： 1、建立岩性类型与单位 2、利用 density 与 gamma 曲线建立岩性曲线 (lithology) 3、建立岩性概率模型 (shale_prob.mod, sand_prob.mod) 4、直方图与变异图的分析 并通过交会图确定两类数据的相关系数 5、随机模拟 6、计算多次孔隙度随机模拟结果的平均值 7、计算多次模拟的标准偏差 8、计算小层中的砂岩累计厚度 9、随机反演 （一）建立岩性类型与单位： Utilities → Project management → Project parameters Edit → Type and units … → → Lothologic data type → Data type :
Ok (存入文件 usertypelist.txt) （二）利用 density 与 gamma 曲线建立岩性曲线 (深度域、lithology)：例如：工区有六口井：1.wll 2.wll 3.wll 4.wll 5.wll 6.wll 经过分析：当 GR<85, density<2300kg/m^3 时为 sand 当 GR>85, density.2300kg/m^3 时为 shale 1、建立一个函数（利用 density 与 gamma 曲线建立岩性曲线）： FunctionMod → Edit → Functions → New 键入 Output type: Output unit: Parameters : Parameters kind 选 Parameters name Parameters type 选 Parameters unit Ok Parameters kind 选 Parameters name Parameters type 选 Parameters unit Ok Use domain range definitions →
稀疏脉冲法 包括最大似然反褶积、L1范数反褶积、最小熵反褶积、最大熵反褶积、同态反褶积等，稀疏脉冲反演是基于脉冲反褶积基础上的递推反演方法,其基本假设是地层的强反射系数是稀疏分布的。从地震道中根据稀疏的原则提取反射系数,与子波褶积后生成合成地震记录;利用合成地震记录与原始地震道残差的大小修改参与褶积的反射系数个数,再作合成地震记录;如此迭代,最终得到一个能最佳逼近原始地震道的反射系数序列。该方法适用于井数较少的地区,其主要优点是能够获得宽频带的反射系数,较好地解决地震反演的多解性问题,从而使反演结果更趋于真实。 约束稀疏脉冲反演采用一个快速的趋势约束脉冲反演算法,用解释层位和井约束控制波阻抗的趋势和幅值范围,脉冲算法产生了宽带结果,恢复了缺失的低频和高频成分;同时,再加入根据井的波阻抗的趋势约束。约束稀疏脉冲反演最小误差函数是J=∑(ri)p+λq∑(di-si)q++α2∑(ti-Zi)2(1) 式中:ri为样点的反射系数;zi为样点的波阻抗;di是原始地震道;si 是合成地震道;Zi介于井约束的最大和最小波阻抗之间;ti是用户提供的波阻抗趋势;α为趋势最小匹配加权因子;p,q为L模因子;i是地震道样点序号;λ为数据不匹配加权因子。 如果从最大似然反褶积中求反射系数r(t),则在上述过程中为了得到可靠的反射系数估计值,可以单独输入波阻抗信息作为约束条件,从而求得最合理的波阻抗模型 Z(t)=Z(t-1)（1+r(t)）/（1-r(t)）(2) 稀疏脉冲法假设反射系数是稀疏的、离散的，利用测井资料可以得到井旁道的准确反射系数，通过上述反褶积方法，在测井资料、地质模型的约束下，逐道递推子波、反射系数，从而反演出波阻抗、速度等数据。 常规递推法与稀疏脉冲反演法主要是利用反褶积方法来恢复反射系数序列，由经过标定的反射系数序列递推出相对波阻抗，然后加上从声波测井和地质模型中得到的低频分量，最终得到反演波阻抗。这两类方法的主要缺陷是选择可靠低频信息较为困难，由反射系数递推波阻抗过程中误差积累快，当反射系数存在较大误差时，递推出来的波阻抗剖面会面貌全非。 此外，经过反褶积处理的结果，并不代表真正的反射系数序列，
周一下午：地震数据加载、管理及质量控制 Types of files created I n JGW : Grid：Run.rapxy.rap这两个文件必须时刻保存。 Seismic :seismic.mod seismic.minseismic.mind三个文件一起拷贝有效 Horizons: Horizons.horHorizons.hinHorizons.hind时间域、深度域和属性都是以.hor格式加载Wells：well.wil这个文件包含多种信息，加载的层位等都在这个文件中 Other optional data:相、子波等 Data classes and types 一个里面classes包含很多types，单位要一致。 Daturms : 其他软件是指海平面为基准的，海平面以上的海拔补心高，JASON是以地震采集的基准面为基准，这个值不能修改，只能修改测井的数据得到准确的值。加载测井曲线时要注意。 自动抓的图片保存在JASON_SNAPSHCTS文件夹里 标准工区建立：一级目录:年份_油田名_工区名 打开软件的界面： File-creat project
点击OK 填写上表：一般只定义单位即可：
点击OK：工区建好，在盘里可以见到工区 加载地震数据：Datalinks–sesmic-SEG-Y-disk SEG-Y import
首先定义一些参数、格式，然后再加载地震数据体. 定义格式Parameters 里面第一行 首先是选择标准SEG-Y,点击右下角，就出来了。 CDP的选择：根据实际数据，如果选择specify的话format要选择16或者32字节，不能选8。 点击Qickverify of settings files name，选中地震数据所在的文件，点击地震数据文件，OK
StatMod MC 入门手册 Chapter 1. 工作流程 Chapter 2. 基本的输入输出数据 输入数据 输出数据 岩性实现 岩性概率体 属性实现 地质统计学参数 岩石物理分析 地层网格模型 地震数据 测井曲线 …………………… ...5% ...….………………..15% ..……………………..5% ……………………...50% ……………………...10% ….………………….15% 百分数表示每个步骤所用时 间占整个项目时间的百分比
Chapter 3.详细操作步骤 操作步骤以StatMod MC培训数据为例 第一步.首先完成一个高质量的叠后CSSI反演 这一步的目的是为地质统计学提供一个好的研究基础, 这个“好”主要体现在： (1)好的井震标定, 目标区的相关值达到0.85以上; (2)好的叠后反演结果, 用来质控地质统计学模拟和反演结果, 是地质统计学反演结果横向预测准确 度的参照物; (3)利用叠后反演结果进行砂体雕刻, 对目标区的岩性展布、比例有一个总体上正确的把握, 这些认 识都是地质统计学的初始输入。 （说明：在提供的培训数据中已经为用户做了以上准备，用户可以从主界面中打开该培训数据所在工区, 然后用Map View看工区底图，用Section View查看地震数据、叠后CSSI反演数据、地质框架模型与 层位数据以及井数据与子波 , 并用Well Editor检查井震标定情况） 第二步. 数据准备 ●●井曲线重采样 这一步将测井数据重采样至地质微层采样间隔，具体操作为： (1)JGW主界面→ Analysis→ Processing toolkit； (2)Input→ Data selection→ Data type：选Well, 点击Input file(s)右边List选择任意井(可以选多井), 然后在弹出的界面Select logs中选择任意井曲线(可以多选),点击OK退出; (3)Parameters→ Resample log, 在弹出界面Processing toolkit中填写重采样间隔(注意s 与ms单位), 点击OK退出; (4)Output→ Define process, 从Select from中选择Resample log, 点击??输入到右边的Process里面; (5)Output→ Generate, 在弹出的界面中填写输出路径和输出文件名,然后点击Generate,开始计算重 采样的曲线。 ●●计算岩性曲线：这里说的岩性曲线和测井解释的岩性不同, 而是指能通过弹性属性(主要指纵波 阻抗) 区分出来的岩性, 所以一般可分两到三种岩性。练习里通过交会图的方式教授如何划分岩性。 ●●地震数据：线、道间距尽量一致，如果不一致需要进行跳道处理，避免往反演结果中引入各向 异性。 ●●地质框架模型：微层采样间隔与需要分辨的薄层厚度要匹配。 第三步. 地质统计学参数分析 这里说的地质统计学参数主要指三个参数：概率密度函数(probability density function, 简称pdf，描述某一属性在空间的概率分布情况)、变差函数（描述某一属性随距离的变化，是距离的函数）、云变换（描述两个属性之间的相关关系）。下面先介绍前两个参数: ●●概率密度函数 (1)JGW主界面→ Modeling→ StatMod MC, 弹出StatMod MC主界面，主要填写内容如下：
一、Jason综合地震反演和储层预测 地质框架模型（Earthmodel） 子波提取和分析（wavelets） 约束稀疏脉冲反演（invertrace） 多参数岩性特征反演（invermod） 地质统计随机模型与随机反演（statmod） (一)地质框架模型（Earthmodel） 为储层和油气藏定量描述创建一个由地震坐标描述的地质框架模型。这个模型融合了构造（层位、断层）、地质沉积模式、测井资料和初始权重分布等信息。地质框架模型是整个Jason地震反演和储层、油气藏定量描述的基础。
1.目的 1)构造以层为基础的、用于Jason其它模块的参数描述模型，即地质框架模型。 2)生成以地质框架模型为基础的测井曲线内插模型。 3)提供用于地震反演的低频模型。 4)生成平滑、封闭的层位顶、底、厚度平面图。
2.模块及功能 （1）Model builder with/without TDC(模型建造器) 用构造、地质、沉积、测井等资料形成一个参数化的时间、深度域的三维封闭模型。这个参数化的模型包括了层位、断层、地层接触关系（整合、不整合、河道、礁等），测井曲线和基于层位的权重系数。 （2）Model generator（模型生成器） 由模型建造器形成的参数化三维封闭模型创建不同测井曲线的三维属性模型（如波阻抗、声波速度、孔隙度、SP等） （3）Model interpolator（模型内插器） 在参数化的三维封闭模型控制下生成不同网格密度的三维属性体。 （4）Well curve generator（测井曲线生成器） 从三维属性体中抽取任意位置上的测井曲线。 3.特点 1)地质框架模型含有生成属性模型所有参数。 2)地质框架模型包括地震、地质、沉积、测井等资料和信息。 3)提供多种内插算法，供用户选择。 (二)子波提取和分析（wavelets） 提供用于合成记录与反演的地震子波估算或理论子波计算工具。在单井、多井或无井的情况下，都可以由单道或多道地震信息估算出最佳地震子波。 1.目的 1)用各种各样的技术估算地震子波。包括地震资料（无井）估算地震子波、用单井或多井和井旁地震道估算地震子波、用单井或多井和不同偏移距部 分叠加资料估算相应的子波。 2)确保合成记录与井旁道的最佳匹配。 3)如果需要，可内插形成空变子波体。 4)估算使地震资料零相位化的反褶积算子。 2.算法与工具 1)以地震资料的自相关为基础，估算子波的振幅谱。 2)估算输入子波的常振幅谱。 3)同时估算子波的振幅谱和相位谱。 4)估算用于Rocktrace的不同角道集的子波。 5)子波特征编辑（如极性、坐标位置、角度范围等） 6)标定子波振幅。