petrel地震反演
地震反演技术
一、概述
2、正演（Forward Modeling) 正演（ 正演和反演相反， 正演和反演相反，它是对一个假设的地质模 给定某些参数（如速度、层数、厚度） 型，给定某些参数（如速度、层数、厚度）用 理论关系式（数学模型） 理论关系式（数学模型）推导出某种可测量的 如地震波）。 量（如地震波）。 在地震勘探中， 在地震勘探中，正演的一个重要应用就是 制作合成地震记录，进行地震标定。 制作合成地震记录，进行地震标定。另一个重 要应用是进行历史拟合。 要应用是进行历史拟合。
ρi+1Vi +1 − ρiVi Z i+1 − Z i Ri = = ρi+1Vi+1 + ρiVi Z i+1 + Z i
University of Petroleum
1、波阻抗递推公式 对应的波阻抗为： Z = Z (1 + Ri ) i +1 i
1 − Ri
递推公式：
Z n+1
ቤተ መጻሕፍቲ ባይዱ
1 + Rn = Z0 ∏ n =0 1 − R n
University of Petroleum
四、基于模型的反演
1、稀疏脉冲反演方法存在的问题 稀疏脉冲反演方法的输出为矩形波阻抗曲线形式，地 层边界清晰，对厚层碳酸盐岩地区较为合适。然而其致命 的弱点是要求反射系数是稀疏的，而实际上大多数地震道 的反射系数是稠密的。 2、基于模型的反演的基本思路 模型为基础的方法，或简称模型法，首先构造一个地质 模型，并将其与地震资料进行比较，然后利用比较的结果 ，迭代地更新模型，直至其与地震资料资料吻合为止。
University of Petroleum
三、递推反演方法
Petrel软件实例操作流程
Petrel软件实例操作流程Petrel软件实例操作流程目录第一章Petrel简介 (1)一、安装并启动Petrel (1)二、界面介绍 (2)第二章Petrel处理流程介绍 (7)一、数据预备 (7)二、断层建模 (15)三、Pillar Gridding (23)四、Make Horizon (28)六、Layering (35)七、建立几何建模 (36)八、离散化测井曲线 (37)八、数据分析 (38)九、相建模 (43)十、属性建模 (55)十一、体积运算 (64)十二、绘图 (68)十三、井轨迹设计 (70)十四、油藏数值模拟的数据输入和输出 (73)第一章Petrel简介一、安装并启动Petrel把安装盘放入光驱，运行Setup.exe程序，依照提示就能够顺利完成安装，在安装的过程中同时安装DONGLE的驱动程序，安装的过程中不要把DONGLE插入USB插槽，安装完毕，再插入DONGLE，假如LICENSE过期，请和我们技术支持联系。
然后按下面的顺序打开软件。
1. 双击桌面上的Petrel图标启动Petrel。
2. 假如是第一次运行Petrel，在执行Petrel运行前会显现一个Petrel的介绍窗口。
3. 打开Gullfaks_Demo项目。
点击文件>打开项目，从项目名目中选择Gullfaks_2002SE.pet。
二、界面介绍〔一〕、菜单条/ 工具条与大多数PC软件一样，Petrel软件菜单条有标准的〝文件〞、〝编辑〞、〝视图〞、View等下拉菜单，以及一些用于打开、储存project的标准工具，在菜单条下面的工具条里还有更多工具。
在Petrel里，工具条还包含显示工具。
此外在第二个工具条里还有位于Petrel 项目窗口的右端的按钮，它具有附加的Petrel相关的功能。
后面的工具条称为功能条，这些工具是否有效取决于选择进程表中的哪个进程。
操作步骤1.点击上面工具条中的每一项看会显现什么。
Petrel地震地质解释和建模使用技巧
Petrel地震地质解释和建模使用技巧Petrel 合成记录工作流制作合成地震记录，进行层位标定和确定时深关系是地震解释工作中非常重要的环节。
从Petel2009.1.1，开始Petrel里有两个制作合成记录的模块，一个叫Synthetics，一个叫Seismic-Well tie。
这里介绍如何使用Synthetics模块制作合成地震记录。
从Petrel 2007开始Synthetics模块有了很大改进。
最重要的变化是其结果可在Global well logs下有相应的synthetic目录，其相应时深关系可在数据表中显示。
对同一口井可产生多个合成记录，如图1-1，1-2所示。
Synthetics模块制作合成记录工作流主要分为两大步骤：按照已有数据产生合成记录通过welltop 进行时深关系调整（bulkshift或sqeeze/stretch）一、 生成合成记录1. 双击synthetic模块，打开合成记录主界面（如下图），选择create new folder，从界面中well 到well seismic 四个界面对合成记录中所需数据进行选择或创建，如图2所示。
Well：选择要做合成记录的井，可多选，但每口井必须有相应的数据（DT和子波）。
Sonic and time：确定原始输入数据及时深关系。
根据实际数据品质，如果有checkshot，可用来做DT曲线校正；所有井上时深关系以工区井目录，以及每口井的Settings界面里Time界面下设置为准，Synthetics界面里的Overwrite global time log项不启用。
Create synthetic seismogram：创建合成记录选择创建合成记录所需数据：Density、Acoustic Impedence、Reflectiotion coefficients和Wavelet。
如果这些数据都不存在，或者希望修改参数重新创建，则点击黄色星状按钮创建新数据。
Petrel地震地质解释和建模使用技巧2013
Schlumberger Private
Trademark Information Software application names used in this publication are trademarks of Schlumberger. Certain other products and product names are trademarks or registered trademarks of their respective companies or organizations.
Schlumberger Private
ห้องสมุดไป่ตู้
Petrel 地震地质解释和建模使用技巧
I
2.1 Petrel 合成记录工作流
Schlumberger
2.1 Petrel 合成记录工作流
制作合成地震记录，进行层位标定和确定时深关系是地震解释工作中非常重要的环 节。从 Petel2009.1.1，开始 Petrel 里有两个制作合成记录的模块，一个叫 Synthetics，一 个叫 Seismic-Well tie。这里介绍如何使用 Synthetics 模块制作合成地震记录。
Petrel 地震地质解释和建模 使用技巧 2013
斯伦贝谢科技服务(北京)有限公司
Copyright Notice © 2009 Schlumberger. All rights reserved. No part of this manual may be reproduced, stored in a retrieval system, or translated in any form or by any means, electronic or mechanical, including photocopying and recording, without the prior written permission of Schlumberger Information Solutions, 5599 San Felipe, Suite 1700, Houston, TX 77056-2722.
[Petrel2014使用技巧] Geobody体雕刻与神经网络算法划分地震相
![[Petrel2014使用技巧] Geobody体雕刻与神经网络算法划分地震相](https://img.taocdn.com/s3/m/eb5bb9c633d4b14e852468e2.png)
Geobody体雕刻与神经网络算法划分地震相Petrel软件为用户提供了多种划分地震相的方法，用户可以通过Geobody实现多地震属性体的融合，充分利用不同岩相在频率、振幅、相位上的特征，实现地震相（异常体）的划分与识别，为建模提供岩相的初步认识；此外，利用Train Estimation Model模块也可以直接对地震体进行聚类分析。
随着勘探的逐渐深入，井资料不断完善，单井相的认识会成为研究岩相分布的硬数据，如何充分利用井资料，有效结合地震多属性是更为准确的获得相（异常体）空间分布的关键。
利用Geobody体雕刻与神经网络算法划分地震相的方法能够充分利用井的硬数据和地震的软数据，具有更加可靠、可控、精确的特点，不仅可以用来对地震相带的识别，也可以用于气藏、油藏、特殊目标体的空间预测，为有利目标的选取和建模提供更可靠的空间认识（图6）。
1．生成多种地震属性通过Volume Attribute模块提取多种能够反映地震相（异常体）的地震属性体；a. 如针对于河道相，可以选择，Sweetness: 甜点属性Genetic Inversion: 遗传反演Reflection Intensity: 反射强度Graphic Equalizer: 图形均衡Chaos: 去噪属性Structural Smoothing: 构造平滑Local Structural Dip:构造倾角RMS Amplitude: 均方根振幅Filter: 滤波Variance (Edge Method): 方差属性（相干体属性）3D Curvature: 曲率Iso-frequency: 等频体b. 针对于气藏，除上述部分属性外，也可以选择T*Attenuation: 频率衰减相关属性2．Geobody体雕刻Geobody体雕刻技术可以帮助用户实现“所见即所得”，用户可以根据单种或多种地震属性融合的方法有效的识别河道、深水浊积体等典型相（异常体）的特征，立足于对单井的认识，能够很容易的借助体雕刻功能获得某一井区附近特定相的空间展布特征（图1）3.Geobody体提取经过体透视和体提取功能获得的Geobody体可以通过Convert to seismic horizon 的功能获得砂体或者异常体（气藏）的顶底面的高度，进而获得其厚度；此外，还可以将提取出来的Geobody体通过Convert to points转换成对应的散点数据（图2-3）。
[Petrel2014使用技巧]如何解决反演体属性在时深转换过程中数值变化的问题
![[Petrel2014使用技巧]如何解决反演体属性在时深转换过程中数值变化的问题](https://img.taocdn.com/s3/m/7413c314b52acfc789ebc9ec.png)
如何解决反演体属性在时深转换过程中数值变化的问题
存在问题：地震波阻抗反演体在时深转换以后其数值发生变化。
解决办法：已知在时间域中已经做好的地震反演体，其数值仅有0与1，0代表页岩，1代表砂岩。
然后正常进行速度模型建立。
速度模型建立完成后，在Volume attributes中选择Depth conversion-General depth conversion进行地震反演体时深转换。
在Parameters-Interpolation中软件默认选择Smooth 方式。
假如直接进行时深转换，其深度域的地震反演体数值会发生变化，这是由于差值因素以及转换过程中采样间隔变化所致。
将Interpolation方法由Smooth改变为None即可解决该问题。
petrel数据输入格式
Wells中包含的数据1.wellheader数据井名X Y 顶深底深KB（补心海拔）symbol（井别，可自定义，可缺省）文件输入格式：well heads2.井斜井轨迹数据MD AZIM（方位角）INCL（井斜角）--常用输入格式，还有其它输入格式文件输入格式：well path/deviation3.测井曲线文件可直接用测井数据文件，只要将对应的孔渗饱数据所在的列对应好就可以了。
文件输入格式：well logs（las3.0）Well Tops中的数据1.分层数据（必须）X Y Z（TVD垂深）Wellpoint层名井名horizon（正常）default（缺省）（格式不是唯一的，只要和文件中的相关行对应起来就好！）输入文件格式：Petrel Welltops等厚图X Y 厚度属性平面图如果有则新建文件夹输入X Y 属性（孔渗饱，净毛比等）断点数据X Y Z（TVD）地震数据1.SEGY数据体可接受2D,3D地震数据，同时地震反演的数据也可以输入到Petrel中。
2.层位解释线（Seismic line,Surface）。
3.断层解释（Fault stick,polygon）。
X Y 断点深度（TVD）4.速度资料（Surface,Point）输入文件类型为：SEG-Y seismic dataDomain（模版）：Elevation Time其它数据1.Horizons（Time）文件类型：General lines/points，Domain:Elevation Time , Line Type: Seismic 3D Lines2.Faults（Time）文件类型：ZMAP+lines，Domain:Elevation Time ，Line Type: Fault sticks3.Fault polygons（Time）文件类型：ZMAP+lines，Domain:Elevation Time ，Line Type: Fault Polygon4.Surfaces（Time）文件类型：Zmap+grid，Domain:Elevation Time ，Surface:Comfortable5.Velocity文件类型：Zmap+grid，Domain:Velocity。
Petrel中文操作手册(1-3)
Petrel地质建模部分实用操作手册目录第一章Petrel简介 (1)1.1 关于本手册 (5)1.2 Petrel Workflow Tools功能简介 (6)1.3 安装并启动Perel (7)1.4 Perel用户界面简介 (8)1.5 菜单及工具条 (10)第二章数据格式说明及加载 (11)2.1 数据准备 (11)2.2 数据输入 (15)练习2-1 创建一个工区 (16)练习2-2 加载井数据 (16)练习2-3 加载分层数据 (20)练习2-4 加载地震数据 (21)2.3 编辑整理数据 (22)2.3.1 自动生成断层多变形 (22)2.3.2 生成/编辑Polygons (23)2.3.3 生成/编辑Surface (24)2.3.4 面体积计算 (27)2.4 井数据的管理 (27)2.4.1归类引擎（Saved searches） (28)2.4.2井滤波（well filter） (29)第三章聚类分析与井相关 (30)3.1 属性（地震、测井）聚类分析和判断 (30)3.2 井相关（Well Correlation） (37)第四章速度模型与域转换............................................................................. 错误!未定义书签。
4.1 构建速度模型.................................................................................. 错误!未定义书签。
4.2 时间-深度转换 ................................................................................ 错误!未定义书签。
第五章框架模型............................................................................................. 错误!未定义书签。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
Statistical a priori model
AI
Horizontal covariance
Vertical
covariance
0
1
0
time
Mean Mean +/- 2*s.d.
20ms
After: Marion, Robotham and Swaby, TFE
Schlumberger Private
❖copy $sp1/workstation/*.msi to Windows NTFS partition
Petrel things
GF-Petrel connection
GeoFrame part
Basemap -> send -> export to Petrel
Geology Office -> tools -> export well data to Petrel
seismic domain. Possibility of extending beyond the seismic bandwidth.
Inversion Limitations
Estimate the frequencies available within the seismic by wavelet estimation or spectral analysis.
Removes the effects of the wavelet within the seismic bandwidth. Forces well ties to be made and understood. Reservoir properties are separated from the overburden. May provide quantitative predictions on the reservoir properties. Stratigraphic interpretation may be improved. Interpreting in the impedance domain is frequently easier than in the
4 DPsaila 7/5/2007
Geostatistical Inversion
Impedance models are geologically constrained:
• Stratigraphic framework
• 3D a priori model
• 3D heterogeneity model
Benefits of Seismic Inversion Inversion Limitations Major applications on the market Petrel things
What is Seismic Inversion?
What is Seismic Inversion?
Now the normal assumption is that the seismic trace is considered to be a primary only reflectivity model convolved with the seismic wavelet which is summed with some uncorrelated noise and can be represented by the equation:
Petrel things
GF-Petrel connection
Petrel part
❖install GeoFrameDataConnectorForPetrel2009.1.msi
Petrel things
GF-Petrel connection
Petrel part
Petrel things
Normal
Wavelet extraction Time-shift wavelet Sonic calibration curve
Abnormal
Sonic filter/smoothing/blocking Squeeze/stretch Space-shift wavelet Average wavelet
where Vp is p-wave velocity and ρ is density
Seismic reflections are given by the change in the Acoustic Impedance of the layers and are represented by the equation
Initial modeling
Post-stack Inversion steps
Inversion
Post-stack Inversion steps
Post-inversion
Benefits of Seismic Inversion
More geoscientists understand the concept of impedance and geology than the seismic trace. Thus, working in the impedance domain is a great mechanism for integrating with the various disciplines in a multidisciplinary asset team.
Post-stack Inversion steps
Well to seismic calibration
Wavelet extraction
❖For every well ❖Short time window analysis ❖Not care about the shape of wavelet
Background and Theory
Post-stack Inversion steps
Well to seismic calibration Initial modeling Inversion Post-inversion
Post-stack Inversion steps
Well to seismic calibration
There are no facts, only interpretations.
Friedrich Nietzsche
Contents
What is seismic inversion Background and theory Post-stack Inversion steps
Well to seismic calibration Initial modeling Inversion Post-inversion
Petrel things
High Resolution Geostatistical Inversion(SSI)
/display/index.do?id=id3010054
/GTeC/Projects/MMRD/RSConceptStudies/2006-Inversion/
The algorithm is a stochastic simulation approach producing multiple high resolution realizations of acoustic impedance that tie at the well control and are consistent with the seismic data. A prior model is used to constrain the vertical and lateral trends of acoustic impedance and a variogram model is used to constrain the patterns of spatial variability of the inverted acoustic impedance. The method can also be extended to handle AVA inversion of multiple seismic angle stacks.
Mean producing many high resolution
3D AI results
SD
8 DPsaila 7/5/2007
0
AI (km/s . g/cc) 0.7
After: Marion, Robotham and Swaby, TFE
Geostatistical Inversion - Products
Major applications
PrismSeismic – CRYSTAL
Major applications
PrismSeismic – REFRACT
DFN
Petrel things
GF-Petrel connection
GeoFrame part
❖fuse 2.7.4 + ntfs-3g 2009.4.4, NTFS driver for Linux ❖GeoFrame 4.4 + TK-2 patch
Schlumberger Private
Geostatistical Inversion - Products
One realization
• Impedance results are at high
resolution (~2-8 m, ~1-4 ms)
4
AI (km/s . g/cc)
9
• Describes non-uniqueness by
Background and Theory