下载文档原格式
PETREL操作流程1.前期数据准备地震数据体，断层线FAULT LINS OR 断层棍FAULT STICKS，FAULTPOL YGONS，数字化的等值线。
工区内各井的坐标，顶深，海拔，底深（完钻井深），东西偏移，方位角，倾角，砂岩分层数据，砂层等厚图，测井曲线（公制单位），单井相，各层沉积相图，砂岩顶面构造图，单井岩性划分，测井解释成果表，含油面积图。
（在编辑数据的过程中，命名文件时最好数据文件名都和井名一致）2.数据加载①加载井口数据（WELL HEADERS）WELL_NAME X Y KB TOP BOTTOM SYMBOL井名X坐标Y坐标海拔顶深底深（完钻）井的类型②加载井斜数据（WELL PATH）第一种数据格式MD TVD DX DY AZIM INCL斜深垂深东西偏移南北偏移方位角倾角第二种数据格式MD INCL AZIM第三种数据格式TVD DX DY（单井用WELL LOGS，多井加井斜可用PRODUCTION LOGS）③加载分层数据（WELL TOPS）（包括断点数据）MD WELLPOINT 层名WELL NAME-1500 HORIZON Nm31 NP1-1600 FAULT Nm32 NP1以WELL TOPS加载之后删除系统的缺省项，新建4项，对应输入数据的列，名称进行编辑，Sub-sea Z values must be negative!(低于海平面的Z值都为负)，该选项在编辑时不要选中④加载测井曲线（WELL LOGS）LAS格式文件MD RESIS AC SP GR曲线采用0.125m的点数据（1m8个点数据），注意有的曲线单位要由英制转换为公制，如：AC 英制单位μs/in要换成工制单位μs/m，再用转换程序转换为LAS格式文件进行输入，以提高数据的加载速度。
如果有孔渗饱数据，按相同格式依次排列即可。
在/INPUT DATA中设置数据的排列顺序，曲线内容较多，系统缺省项只有MD，所以要用SPECIFY TO BE LOADED定义新的曲线，对应加载数据的列数，名称和属性进行编辑。
反演问题的数值解法研究第一章引言反演问题是指通过观测数据得到模型参数或物理参数的过程。
在许多领域中，反演问题都是非常重要的，如地球物理学、医学成像、无损检测等。
由此带来的数值计算问题也是非常重要的，因为反演问题涉及到从离散的观测数据中推断出连续的参数，需要依赖数值方法来求解。
本文主要呈现了一些常用的反演问题的数值解法的研究，包括线性反演问题和非线性反演问题。
我们将对各种反演问题的数值解法进行介绍，包括正则化方法、Bayesian方法、梯度下降等。
第二章线性反演问题线性反演问题是指观测数据与模型参数之间的函数关系是线性的反演问题。
我们通常将这种问题表示为$A\mathbf{x}=\mathbf{b}$，其中$A$是线性算子，$\mathbf{x}$是模型参数，$\mathbf{b}$是观测数据。
线性反演问题的数值解法可以使用奇异值分解（SVD）或者正则化方法。
其中，SVD可以将线性反演问题转换为一个完全指定和完全可逆的问题，可以得到唯一的解。
但是，由于数值算法的限制和观测数据误差的影响，SVD不一定是最好的解决方案。
为了解决这个问题，我们可以使用正则化方法。
正则化方法是一种通过增加稳定性约束条件来处理不适定反演问题的技术。
这些约束条件可以有效地减少反演问题的不确定性。
常用的正则化方法包括Tikhonov正则化和阻尼最小二乘法。
Tikhonov正则化是通过加入二次惩罚项来限制解的大小，从而使得解更加平滑。
阻尼最小二乘法是通过同时加入观测数据误差和模型误差的项来解决线性反演问题。
这两种方法都可以通过基于SVD的方法求解。
需要注意的是，对于线性反演问题，只有当观测数据是无误的时候才能得到正确的解。
这是因为线性反演问题的解非常敏感，即使存在微小的误差，也会导致解的失真。
第三章非线性反演问题与线性反演问题不同，非线性反演问题的观测数据与模型参数之间具有非线性关系。
常见的非线性反演问题包括逆时偏移（RTM）、全波形反演（FWI）和电磁成像等。
第一章反演理论第一节基本概念一．反演和正演1．反演反演是一个很广的概念，根据地震波场、地球自由振荡、交变电磁场、重力场以及热学等地球物理观测数据去推测地球内部的结构形态及物质成分，来定量计算各种有关的物理参数，这些都可以归结为反演问题。
在地震勘探中，反演的一个重要应用就是由地震记录得到波阻抗。
有反演，还有正演。
要正确理解反演问题，还要知道正演的概念。
2．正演正演和反演相反，它是对一个假设的地质模型，给定某些参数（如速度、层数、厚度）用理论关系式（数学模型）推导出某种可测量的量（如地震波）。
在地震勘探中，正演的一个重要应用就是制作合成地震记录。
3．例子考虑地球内部的温度分布，假定地球内部的温度随深度线性增加，其关系式可表示成：T(z)=a+bz正演：给定a和b，求不同深度z的对应温度T(z)反演：已经在不同点z测得T(z)，求a和b。
二．反演问题描述和公式表达的几个重要问题1．应用哪种参数化方式——离散的还是连续的？2．地球物理数据的性质是什么？观测中的误差是什么？3．问题能不能作为数学问题提出，如果能够，它是不是适定的？4．对问题有无物理约束？5．能获得什么类型的解，达到什么精度？要求得到近似解、解的范围、还是精确解？6．问题是线性的还是非线性的？7．问题是欠定的、超定的、还是适定的？8．什么是问题的最好解法？9．解的置信界限是什么？能否用其它方法来评价？第二节反演的数学基础一．解超定线性反问题1．简单线性回归可利用最小平方法确定参数a 、b 使误差的平方和最小。
⎪⎪⎩⎪⎪⎨⎧∑-∑∑∑-∑=-=∑∑-=22)()(x x n y x xy n b x b y n x b y a （1-2-1） 拟合公式为：bx a y+=ˆ （1-2-2） 该方法的公式原来只适用于解超定问题，但同样适用于欠定问题，当我们有多个参数时，称为多元回归，在地球物理领域广泛采用这种方法。
此过程用矩阵形式表示，则称为广义最小平方法矩阵方演。
石油勘探中地震反演技术的使用方法与数据处理流程解析地震反演技术是石油勘探中非常重要的一项技术，它通过分析地震波在地下介质中传播过程中的特性，来推断地下油气储层的位置、性质和形态等信息。
本文将对地震反演技术的使用方法与数据处理流程进行详细解析。
地震反演技术的使用方法分为几个基本步骤，包括数据采集、数据处理、模型构建与反演、结果解释等。
下面将对每个步骤进行详细介绍。
首先是数据采集。
石油勘探中常用的地震勘探方法有地面反射法、地震层析法以及地震井间反射法等。
在数据采集阶段，需要布设地震检波器，并记录地震数据。
通常会采集多个不同位置的地震记录，以获得更全面的地下信息。
数据采集完成后，接下来是数据处理。
数据处理的目标是将采集到的原始地震数据转化为可以进行反演的有效信息。
数据处理中的一项重要任务是去除噪声，包括地震数据中的随机噪声和仪器本底噪声等。
此外，还需要对数据进行质量控制和预处理，以确保后续的反演和解释工作的准确性和可靠性。
在数据处理完成后，就可以进行模型构建与反演了。
模型构建的目标是根据采集到的地震数据，构建地下介质模型，并将模型中的参数用于地震波传播方程的计算。
地震反演则是根据模型构建的结果，将地震数据与模拟合成数据进行比较，以推断地下构造的性质、位置、形态等信息。
地震反演中最常用的方法是基于反演算法的全波形反演技术，该技术通过求解正问题和反问题的数值算法，来精确地恢复地下介质的特征。
最后是结果解释。
经过地震反演后，得到的结果通常是一系列地下介质的参数分布，包括速度、密度、衰减等。
结果解释则需要借助地震学理论、地质学和地球物理学等方法，将这些参数与地下岩石的性质进行对比和解释。
通过分析地下构造的特征，可以推断油气储层的位置、性质、形态等信息，并为后续的油气勘探提供重要的依据。
除了使用方法外，地震反演技术的数据处理流程也需要进行详细解析。
通常的数据处理流程包括数据预处理、地震成像、速度模型更新等步骤。
![[Petrel2014使用技巧]地震解释窗口中Ghost的应用](https://img.taocdn.com/s1/m/74fad470783e0912a2162ae2.png)
地震解释窗口中Ghost的应用
应用一：叠合显示多种属性
1、点击Ghost图标，画矩形圈定范围，确认在Windows栏Ghost 被激活，如图1；
2、在Input勾选另一属性体（注意不是Inline或Xline），如图1显示的为Envelope属性。
此时如果移动窗口或改变窗口大小，可以点击右键更新对应数据。
3、Windows栏双击Ghost，在settings里可以设置Ghost窗口的透明度。
右键菜单可以删除Ghost。
应用二：辅助地震解释
Ghost还可以用来进行断层上下盘同相轴的对比，与邻近剖面的对比，帮助确定解释的层位。
此时，Ctrl+Shift可以旋转Ghost窗口，Shift+RMB恢复垂直状态，如图2。
应用三：
在Ghost窗口里显示井曲线，通过压缩和拉伸Ghost窗口，可以与地震剖面之间进行初步的对比，如图3所示。
断层建模中如何较好的处理复杂接触关系的断层在做构造模型的过程中，遇到断层发育的油藏，断层交接关系又很复杂，即使勉强处理的断层接触关系，在网格化之后，网格的质量很难保证，往往容易形成畸形的网格，给数模的网格质量也很不理想。
在处理这种问题的时候，petrel2012给我们提供了一个好的解决办法，避免了上述的诸多问题。
工作流程如下1、在Structural framework下，双击Geometry definition进程（图1），创建一个新的几何体，从地震数据体得到工区范围，选择对应的深度域或时间域，点击ok（图2）。
2、双击Fault framework modeling进程，弹出窗口如下（图3），从input中输入所有的解释断层，设置相应的参数，点击ok，则自动生成含有所有断层相互连接关系的断层格架（图4）。
在3D窗口中查看断层之间的交接关系，对接触关系复杂又处理不好的断层，记下它们的断层名字。
3、展开Corner point gridding文件夹，双击Fault model from structural framework进程(图5)，在这个进程界面下，对第二步中记下的断层Closing FaultEast、Closing Fault South，在include标签下对应的断层名字处，不打勾，让它们不参与计算，设置相应的参数，点击ok（图6）。
即对处理不好的复杂断层先不让它参与形成断层模型，以保证整体断层模型的网格质量，然后按照常规的构造建模流程，进行pillar gridding<make horizons<make="" zones<layering，各个阶段都要进行相应的质量检查与控制。
一系列的小层构造模型建完之后，开始添加复杂的断层进入到成型的构造框架中。
采用阶梯状网格算法，保证复杂断层附近的网格不发生扭曲变形。
双击Stair-step faulting进程（图7），弹出界面窗口中即为没有参与断层模型建立的复杂断层（图8），在Stair-step标签下勾选，点击ok。
如何解决反演体属性在时深转换过程中数值变化的问题
存在问题：地震波阻抗反演体在时深转换以后其数值发生变化。
解决办法：已知在时间域中已经做好的地震反演体，其数值仅有0与1，0代表页岩，1代表砂岩。
然后正常进行速度模型建立。
速度模型建立完成后，在Volume attributes中选择Depth conversion-General depth conversion进行地震反演体时深转换。
在Parameters-Interpolation中软件默认选择Smooth 方式。
假如直接进行时深转换，其深度域的地震反演体数值会发生变化，这是由于差值因素以及转换过程中采样间隔变化所致。
将Interpolation方法由Smooth改变为None即可解决该问题。
如何解决反演体属性在时深转换过程中数值变化的问题
存在问题：地震波阻抗反演体在时深转换以后其数值发生变化。
解决办法：已知在时间域中已经做好的地震反演体，其数值仅有0与1，0代表页岩，1代表砂岩。
然后正常进行速度模型建立。
速度模型建立完成后，在Volume attributes中选择Depth conversion-General depth conversion进行地震反演体时深转换。
在Parameters-Interpolation中软件默认选择Smooth 方式。
假如直接进行时深转换，其深度域的地震反演体数值会发生变化，这是由于差值因素以及转换过程中采样间隔变化所致。
将Interpolation方法由Smooth改变为None即可解决该问题。
