Jason约束稀疏脉冲反演流程
Jason反演的流程和实例+基础教程
Depth sample interval : 2 m
1ms
310
0.25
(2) Edit data for seismic modeling : Wavelet interpolation …→ 选内插算法
Time gate …→
Use horizons
Set horizons file … *Tinterface.hor
(3) Output solid model :
notdc
Output directory: Test163gz
(4) Interface & layers : Merged
Separate
产生文件： ../ Test163gz /Tinterface.hor, Tthickness.hor, Ttopbase.hor 4、File → Save & Exit 5、剖面显示，检查时间域的 井曲线和分层数据
11
图3
图4 编辑趋势 （曲线中心）
Select well … (逐口井检查趋势，反复做) Save as …→ *.atm ⇐ Ok 同样编辑曲线的带宽，edit——edit hard constraint(图 5)
编辑带宽 （曲线范围） Select well … (逐口井检查趋势，反复做) Save as …→ *.atm ⇐ Ok (2) Edit constraints … (用井的阻抗曲线，约束地震反演的数据范围)
多井），此处选 t712_ed.wll, t714_ed.wll, t717-1_ed.wll, t710_ed.wll (5) Logs …→ 选多种井曲线 (只选三种: p_sonic,
density,
p_impedance) (6) Trace gate …→ 选道范围???
Jason稀疏脉冲波阻抗反演流程
地震反演是储层预测和油藏描述中非常关键的技术。
地震反演技术将地震数据、地质构造模型、井点测井资料有机地结合为一体；把地震资料的常规界面反射型剖面，反演成能与测井资料直接对比的岩层型波阻抗剖面,以及岩性剖面和孔隙度等地层物性剖面，反演结果分辨储层的能力也有了明显提高。
Jason稀疏脉冲波阻抗反演流程(附流程图)1 子波提取子波长度的选取要适宜，一般为80ms左右，在浅层，地震频带较宽，子波可短些；深层，地震频带窄，子波可略长些；对提取的子波，要从子波的波形、振幅谱、相位谱等方面进行子波质量的判断，要求子波的波形稳定、频域内单峰平滑、有效频带内相位稳定。
2 合成地震记录标定反演过程中，可以利用反演波阻抗数据，对地震地质标定进行精细微调，提高标定的质量。
子波提取与合成地震记录标定有着密切的关系，两者是循环反复进行的。
先用初始理论子波进行初步标定，标定的时深关系基本确定后，再用井旁地震道和测井反射系数重新估算子波。
完成最终的地震地质标定之后，再重新提取更精确合理的子波。
3地质层位解释解释的层位要尽量靠近反演目的层，要十分精细，反映出目的储层的空间变化特点来。
同时，使用时还要对层位中的缺失、奇异点进行修补、剔除，保证反演所用层位的正确性。
4 地震、地质、测井约束模型的建立建立尽可能接近实际地层沉积情况的波阻抗约束模型，是减少反演最终结果多解性的十分重要的一环。
建立波阻抗模型的过程实际上就是把地震界面信息与测井波阻抗信息正确结合起来的过程。
5 反演参数选择选择合理的反演参数是做好反演处理的前提，稀疏脉冲反演中的主要参数有稀疏约束项、拟合误差约束项以及模型约束项等等。
针对不同的地质情况，有目的性地选取合适的反演参数，同时利用相应的质量监控手段，加强对反演信息的反馈处理，可以有效地改进反演结果的质量，起到事半功倍的效果。
6反演结果检查处理过程中的反演数据体反复检查、各种资料对比分析，发现异常、寻找原因、反复反演，知道获得满意的结果为止。
叠后约束稀疏脉冲反演中文培训教程
叠后约束稀疏脉冲反演中文培训教程一、叠后约束稀疏脉冲反演的基本原理叠后约束稀疏脉冲反演是一种地震反演方法，它利用地震波传播和反射的物理过程，通过对地震数据进行逆向建模，来获取地下介质的信息。
该方法的核心思想是通过寻找一种最优的模型，使得通过该模型建立的地震波场与观测到的地震数据尽可能吻合。
在叠后约束稀疏脉冲反演中，通常将地下介质的性质用参数化的方式表示，将地震波传播和反射的物理过程用数学模型描述，通过优化方法来求解最优的模型参数。
这些模型参数可以包括地下介质的速度、密度、衰减系数等参数，通过对这些参数的反演，可以获取到地下结构的详细信息。
叠后约束稀疏脉冲反演的关键之处在于它采用了稀疏性约束和地震波形约束。
稀疏性约束是指地下介质的参数在某种表示下是稀疏的，即它们可以被较少的非零系数表示。
地震波形约束是指通过对地震波传播和反射的物理过程进行数学建模，来约束模型参数的取值范围，使得反演得到的模型更符合地震数据的实际观测。
二、叠后约束稀疏脉冲反演的算法流程叠后约束稀疏脉冲反演的算法流程通常包括以下几个关键步骤：1. 数据预处理：对地震数据进行预处理，包括数据去噪、时域频域滤波等操作，以减小噪声对反演结果的影响。
2. 构建反演模型：将地下介质的性质参数化表示，并建立地震波传播和反射的数学模型。
3. 求解反演问题：利用最优化方法，求解使得通过模型建立的地震波场与实际观测数据吻合的最优模型参数。
4. 结果评估：对反演结果进行评估，包括地下结构的成像、模型参数的稳定性等。
5. 反演结果应用：将反演得到的地下结构信息应用于油气勘探、地质勘探等领域。
三、叠后约束稀疏脉冲反演的实际应用叠后约束稀疏脉冲反演在油气勘探、地质灾害预测、地下水资源开发等领域具有重要的应用价值。
在油气勘探中，通过叠后约束稀疏脉冲反演可以获取到油气藏的地下结构信息，帮助勘探人员准确地定位油气藏的位置、形态及大小。
在地质灾害预测中，叠后约束稀疏脉冲反演可以提供地下结构信息，帮助人们预测地震、泥石流等自然灾害的发生概率以及可能的影响范围。
迭后约束稀疏脉冲反演中文培训教程
Data_loading 工区
该数据只是用于数据加载练习，在后面的反演和解释练习中并没有用到该数据。
Wedge2d 工区
该数据是合成的二维数据工区，该工区让您充分理解利用波阻抗数据来开展工作的好 处。
Big_gom 工区
该工区是更大范围的 GOM 工区，您将通过它学习如何解释 InverTracePlus 生成的反演 结果，在该练习中我们会用到工区内的 6 口井，在工区中有一个 InverTracePlus 反演数据体 aitm_ai.mod 文件。
第五章：建立低频模型............................................................................42 5.1 构造格架建立 ............................................................................42 5.2 低频模型的产生.........................................................................44
启动 Map View 模块
在 JGW 主菜单点击 File > Select project 并选择 GOM 工区，点击 Analysis > MAP VIEW 启动该模块。
数据的选择显示
通常有五种方法可以选择您需要加载的数据： －在 MAP VIEW 窗口顶部通过 Input 菜单。 －通过菜单下面的分类快捷按钮 , , , 进行数据加载。 －通过 Plot List，点击显示窗口右边的快捷按钮 可以进入，进行数据加载。
第二章：工区的建立与数据的加载..........................................................24 2.1 工区的建立 ................................................................................24 2.2 数据的加载 ................................................................................25
Jason反演Wavelets中文手册
Jason反演Wavelets中文手册13．WaveletsWavelets应用程序的目标是：利用各种技术评估地震子波：单独从地震数据利用来自于一口井或多口井的反射系数利用地震数据的部分（partial）叠加或multiple partial 叠加定位井轨迹以优化井和地震数据的关系横向变化时，进行子波插值为建模生成合成子波为得到零相位地震数据计算反褶积评估Q和产生基于子波的Q13.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和密度类型。
SeismicSeismic 菜单有两部分。
在第一部分中，指定了地震数据文件。
你必须提供一个包含2D或3D数据的地震数据（.mod）。
地震数据必须代表真正振幅偏移反射信息。
在Zero-offset模式中，合成地震记录是通过子波与零偏移反射系数曲线褶积而来的，其中反射系数曲线来源于输入的P-Impedance曲线。
在AV A模式中，角度依赖的反射系数来源于P-Sonic，S-Sonic和密度井曲线，所使用的角度是是在输入的A VA地震数据文件中指定的。
在第二部分中，需要设置地震数据的A V A参数。
Jason 软件详细介绍与基本操作手册
Jason模块功能简单描述Jason地学平台是一个将各学科有用信息综合在一起，即为用户提供丰富的、从油气勘探、开发到生产阶段将地震、测井和地质信息有机结合，进行反演、模拟和预测不同阶段储层岩石物性和流体变化的综合研究工具。
您根据可用信息和经济条件，您可以明智的选择应用哪个“工作流程”更好。
1、Enveriment 运行环境（工区管理地震解释和接口）与其它解释系统紧密相连、数据格式灵活、方便。
2、EarthMod 地质框架模型，包括模型建造器，模型生成器和测井曲线生成器。
把地质、地震、测井曲线（斜井轨迹）、沉积模式（整合，不整合，礁体及河道等）建立以地震道为坐标的模型上，成为储层预测的基础，并赋予地质意义。
3、VelMod 速度模型建立工具VelMod 除具有建立速度模型的基本功能外，且有其独特特征：•地质引导：结合地质信息进行插值运算，能把稀疏的速度数据拓展为三维速度场。
•基于地质概念的速度调整•真三维插值4 、Wavelets 二维和三维地震子波分析●无井、单井、多井提子波和空变子波分析●从斜井中进行真三维分析●估算信噪比谱●为使结果稳定利用模型驱动5、 Modtrace 完整的地震反演系统( 包括对大倾角地层)●所用约束稀疏脉冲反演是基于道的递归反演，产生地震带宽内的波阻抗数据。
●其算法的唯一特征应用解释控制，为获得可靠的低频信息提供工具。
在EarthMod 的基础上建立空变约束条件，其结果突出各向异性，提高了反演分辨率。
●多种方式提取子波，包括无井、单井多道、多井多道、直井、斜井提子波和在EarthMod 的基础上提空变子波。
●由于大量地震信息的反演是基于EarthMod 的基础，所以适用与多条复杂断层的地质情况，使其反演结果更接近真实的地质模型。
●给出产层有效厚度图和孔隙度分布图６、InverMod 精细储层描述（多井）•采用一欧洲石油公司的专利，专门致力于薄层预测和精细描述。
对油田滚动勘探开发十分有用。
稀疏脉冲反演方法及Jason软件实现
第四篇 约束稀疏脉冲波阻抗反演1、概述：约束稀疏脉冲波阻抗反演方法（简称CSSI ）是目前应用比较广泛的一种方法，它是以趋势约束稀疏脉冲反褶积算法为基础。
对地下的反射吸收序列作如下假设：它的分布是稀疏分布的，而不是像其它的反演方法那样假设其是连续分布的。
在反演过程中，它是通过寻找一个使目标函数最小的反射系数脉冲数目来得到波阻抗数据的，此时合成的地震记录可以和实际地震记录匹配的最佳。
2、基本假设和方法原理约束稀疏脉冲反演（CSSI ）是一种递推地震波阻抗反演方法，基于范数反褶积和最大似然反褶积。
它假设地层的反射系数序列是由一系列稀疏且不连续的大反射系数和高斯分布的小反射系数叠加而成。
从地质角都解释，大反射系数代表地下不连续界面、岩性分界面。
首先由约束稀疏脉冲的基本假设，得到目标函数：)1ln()(2)ln(2)(1)(1122122λλ----+=∑∑==M L M K n NK r RJ LKLK（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反演的流程和实例+基础教程
Generate residual data Ok 4、File → Save & Exit 显示反演的进程，完成后可以在 Analysis——Section View and Interpretion,查看结果。
约束稀疏脉冲反演
1、Input :
(1) Seismic data and wavelets …→
Seismic → list …→ mig.mod Select Wavelet …→ list …⇐ 选一个子波 ：Create_ Wavelet25.mtr
(Hale Waihona Puke ) Trace gate …选整个体
(3) Time gate …→
多井），此处选 t712_ed.wll, t714_ed.wll, t717-1_ed.wll, t710_ed.wll (5) Logs …→ 选多种井曲线 (只选三种: p_sonic,
density,
p_impedance) (6) Trace gate …→ 选道范围???
2、Edit：
(1) Edit framework … -------关键 注意：
Top horizon
Bottom horizon
3、Output → Generate … （选 *.mod）
9
Run now 任选输出内容: Timpedance_from_sonic.mod
Timpedance_from_sonic_and_density.mod Timpedance_from_sonic_gardner.mod Tporosity.mod Tsonic.mod Tdensity.mod *Ok 产生文件： ../test163mod—genetate/*.mod 4、File → Save & Exit
约束稀疏脉冲反演方法在识别煤层顶板砂岩中的应用
约束稀疏脉冲反演方法在识别煤层顶板砂岩中的应用说明了稀疏脉冲反演方法的原理、分析了反演过程中需要准备的基础资料和需要注意的关键环节，采用该方法对研究区煤层顶板砂岩进行预测，获得了砂岩的分布规律，经钻井验证可知效果较好。
标签：煤层；砂岩；反演Abstract：this paper explains the principle of sparse pulse inversion method，analyzes the basic data to be prepared and the key steps to be used in the inversion process，and uses this method to predict the coal seam roof sandstone in the research area，and obtains the distribution of sandstone. The results are consistent with drilling data.Keywords：coal；Sandstone；inversion序言煤田地震勘探主要目标为煤层的展布、断层的发育程度、顶底板含水性等[1-3]。
顶板的含水性与砂质含量有直接的关系[4-5]，因此，寻找顶板砂岩的分布规律是煤田地震勘探的一个重点，本文针对煤层顶板砂岩的地质特点提出采用约束稀疏脈冲反演的方法来预测。
1 研究区概况研究区位于陕北黄土高原北端，毛乌素沙漠东南缘，被第四系风积沙和风沙滩地所覆盖，以风蚀风积沙漠丘陵地貌为主，表、浅层地震地质条件较差。
含煤地层为延安组，共9层煤层，可采煤层为2号、3-1号、4-2号、5号煤层，局部可采煤层为7号、8号、9号煤层。
其中2号煤层赋存于延安组第四段顶部，为区内最厚的主采煤层，煤层厚度为2.01～4.26m，局部夹一层粉砂岩夹矸。
煤层顶板岩性主要为细粒砂岩、粉砂岩，次为泥岩或中粒砂岩。
Jason反演入门培训教程
第一章数据的加载jason是目前最常用的反演软件，它操作上的特点是它需要什么数据或参数就给它什么数据或参数。
下面是它的主窗口（图1）。
图1因为作反演之前已经将坨163区块进行了构造解释，所以可以直接从lanmark中将地震、测井、层位数据导入jason中，操作比较简单。
步骤如下：1、选择工区（即一个文件夹）主窗口——File——Select Project（图2），弹出图3。
选择一个工区，ok。
图2图32、数据的导入主窗口——Datalinks——Landmark——Landmark Link（2003）（图3），弹出图4。
图3图4A 工区的选择File——Seisworks project：选地震工区t163，ok。
（图5）图5File——Openworks project——选SHNEGCAI, 选井列表t163，ok。
（图5）此时，图5 窗口的状态栏将会发生变化，以上选择的工区将会显示。
（图6）图6B 地震数据的导入Select——Import——Seismic/property data（图7）,弹出图8。
选cb 3dv（纯波数据，作反演时一定要用纯波数据），ok。
图7图8C 层位数据的导入Select——Import——Horizons，选择反演时需要的层位和断层（图9）。
图9D：井数据的导入Select（图7）——wells，弹出图10。
选择需要的井，ok。
图10E：数据的传输Transport——Import，以上所选的landmark中的数据将传入jason中。
图11第二章合成记录的建立在jason上建合成记录的特点是精度高，但随意性大。
建立合成记录的步骤是:井曲线、地震数据、子波的加载，子波的编辑和评价，合成记录的生成和编辑。
1、井曲线、地震数据、子波的输入主窗口——Analysis——Well log editing and seismic tie（图1），弹出图2。
图1图2Input——well，选择要输入的井，例如：T714_e。
jason地质统计学反演手册
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,开始计算重采样的曲线。
反演方法综述
④从概率密度函数中随机抽取一个值，计算反射系数并与子波进行褶积得到合成地震记录；
若这个值能使合成地震记录与实际地震记录的匹配程度增加则接受此值，若不增加则以一定的概率接收此值，接受的概率分布由波兹曼函数决定。若拒绝则返回上一步；
降低模拟退火的温度；
2、关键环节：
（1）分析储层地球物理特征
测井资料，尤其是声波和密度测井资料，是初始模型建立的基础和地质解释的基本依据，但是一般情况下声波测井都会受到井口环境例如井壁垮塌、泥浆浸泡等的影响而产生误差，同一井口的不同层段，不同井口的同一层段的误差都不尽相同。因此，用于制作初始波阻抗模型的测井资料必须经过环境校正。声波资料是唯一与地震发生联系的资料，储层与围岩的声波特征不同是进行测井约束反演的先决条件。但是由于储层的固有结构和钻井过程中的工程因素，造成目的层段和围岩声波测井上无明显差异。这就要求在仔细分析测井资料的基础上，对声波测井进行合理的校正，这就是储层地球物理特征重构。
2、技术关键和基本流程
地质统计学反演过程中的关键问题是如何使下一个模拟值(从pdf中抽取值)比前一个模拟值更快地达到规定的与地震数据的匹配程度,以避免陷入局部极小或由于大量的迭代次数而耗费机时。目前多采用模拟退火算法解决这个问题。
（1）模拟退火算法的原理：模拟退火(simulated annealing)算法是一类被称为蒙特卡罗法的随机张弛法,它允许目标函数在增加的方向上作随机的变化,因此能跳出局部极小值,找到全局或逼近全局的最优解。在退火过程中系统的能量服从波尔兹曼概率分布,系统依概率
通过最大似然反演导出波阻抗Zi反演公式为Zi=Zi-1*[(Ri+1)/Ri].
具体的计算方法是稀疏脉冲序列每次建立的反射系数为一个脉冲，然后在地震资料中提取子波与初始反射系数进行褶积，得到一个初始合成地震记录，并用此合成地震记录与实际地震纪录作对比得到他们之间的残差，利用这个残差的大小来修改反射序列中脉冲的个数再次进行褶积运算，得到新的合成地震记录，再与实际地震资料对比，就这样循环迭代，直到残差达到最小，最后得到一个与实际地震资料达到最佳逼近的合成地震记录，获得宽频带的反射系数。
反演及Jason软件介绍
预处理 输入地震、测井、地质等资料
建立地质模型（EarthMod）
子波估算 进行主成分分析PCA 进行模型估算Model Estination 生成反演数据体 Model Generator (Al、KTH、Rt、Sp…..) 用Function Mod生成相应的平面图
InverMod反演流程图
七、地质统计与随机模拟系统
（Statmod）
基于地质统计和随机模拟的储层预测反演 工具。它包括三维序贯高斯条件模拟 （SGCS），三维序贯指示模拟（STS）， 三维克里金及三维协克里金（Krigging CoKrigging），它即可以得到连续变量（孔 隙度、渗透率等）的分布图，又可以获得 非连续变量（岩性等）的分区图。它所得 的结果及满足各种参数的随机性又满足各 种参数的空间分布的相关性。
八、通用的数据参数转换工具
（FunctionMod）
用户根据定义的函数将JGW数据（Jason 内部格式）转换成新的JGN数据。函数 可以是单变量或多变量的。
Jason精细油藏描述软件综合利用多学科技术进行 精细储层研究和描述工具，其中不仅包含有常规 的反演技术，更重要的是包含了独创性的随机反 演技术和综合应用技术。为复杂油气田的勘探开 发提供了独特的有效手段。
四、建立速度模型系统 （Velmod）
基于地质模型的速度模型建立工具。 在无井或有井的情况下，利用迭加速 度为你提供了一个与地质模型一致的 可靠的速度模型。该模型可用于一切 需要速度模型的地震资料处理和解释 工作。
五、地震反演系统（InverTrace）
基于地震道的声阻抗反演和储层参数预 测工具。InverTrace包括递归反演和约束 稀疏脉冲反演。反演参数的预选和约束 的可视化编辑功能及空变约束，既保证 了反演的质量，又提高了效率和效益。 灵活方便（可视化）的低频，带限频率 的合并方法，确保绝对阻抗的合理可靠， 进而获得与声阻抗密切相关的产层有效 厚度和平均孔隙度。
波阻抗反演基本原理及Jason反演软件介绍
Ro ck tra ce
Velmod
Wa ve let s
Environment Plus Earthmodel
Statmod
r go La
Functionmod Jason 的重要特点就是随着越来越多的非地震信息（测井，测试，地 质）的引入，由地震数据推演的油气藏参数模型的分辨率和细节会得到不 断的改善。用户可根据需要由Jason 的模块构建自己的研究流程。
四、Jason6.2平台反演模块的构成
Invermod
Inv r ert a ce
Ro ck tr
Velmod （速度建模）
ac e
Velmod
Wa ve let s
Environment Plus Earthmodel
Statmod
o arg L
在地质框架模型的控制下，用地震迭加 速度建立三维速度模型。适用于无井，少井 或速度横向变化大的研究项目。 • 建立三维速度模型（无井）或提高速度模 型的精度（少井或速度横向变化大） • 压力预测。 时深转换，深时转换。 • 提供阻抗的低频模型。
Rocktrace (弹性反演）
Velmod
Wa ve let s
Environment Plus Earthmodel
rgo La
基于模型的 测井属性反演
地质统计随机模 拟与随机反演
道积分 LandMark
G-LOG VELOG PIVT SEISLOG DELOG RM（GeoQuest） Paradigm(叠前) Strata Jason
ISIS Jason
SLIM BCI(宽带约束反演) 广义逆波阻抗反演 PARM Strata Jason(InverMod)
Jason约束稀疏脉冲反演在苏里格气田下石盒子组砂岩储层预测中的应用
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)为“三低一大”即低自然伽玛、低密度、低声波时差;大幅度自然电位异常。
L1-L2范数联合约束稀疏脉冲反演的应用
L1-L2范数联合约束稀疏脉冲反演的应用
脉冲反演技术是地球物理勘探中常用的一种非常有效的成像手段。
在这个过程中，我们通过对电磁波或者声波进行观测，来逆推出地下储层的物理特征。
然而，由于现实中存在的噪声，以及观测数据的不完善性，使得脉冲反演过程中往往需要引入约束条件，来规避掉问题的局部极小值和过拟合等问题。
其中，L1和L2范数联合约束稀疏脉冲反演技术是近年来备受
关注的一种方法。
具体来说，该方法旨在在保留主要特征的同时，对噪声进行有效的消除，同时加速反演计算。
具体的实现方式是在反演过程中，引入L1和L2范数两个维度来进行约束，这样就可以同步优化反演结果的精度和稀疏性。
在应用方面，该方法的适用范围非常广泛。
例如，在油气勘探中，我们可以将其应用于将地面观测数据转化为地下储层模型；在医学成像中，可以将其用于对人体组织进行精细成像等等。
在这些应用中，L1-L2范数联合约束稀疏脉冲反演技术都可以
起到非常好的优化效果，提高成像分辨率和精度，缩短反演计算时间，降低成本。
综上所述，L1-L2范数联合约束稀疏脉冲反演技术是一种非常
有效的地球物理成像手段。
在今后的应用中，我们可以把它应用于更广泛的领域，不断完善其理论和算法，为地球物理勘探事业的发展进步贡献力量。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
12
JGW_8v2_workflow_SYB
【ASCII Horizon Import-(Untitled Session)】 ParametersExisting files，选择 Append and overwrite overlap Existing XY coordinates，选择 Do not read new XY coordinates。
在外部打开层位文件，查看内容。 第二列为 Inline，第四列为 Xline，第五列为 X 坐标，第六列为 Y 坐标，第七列为 Z 值。 【Edit ascii horizon generic colum format definition】 选中 Column 2， 点击 Edit， 弹出窗口 【Edit column definition】 ， Column definition 设置为 Line number。
11
JGW_8v2_workflow_SYB
【Import ASCII horizon files】 点击 Add，选择要加载的层位和断层文件，如果有深度域的层位，需要点击 Edit 将其设置 为 TVD 格式； 在 JGW horizon file 填写内部文件名称，如 demo。OK。
【ASCII Horizon Import-(Untitled Session)】 设置完成。
弹出窗口【Configure project “Jason_demo”】 选择米制单位：Metric units，勾选或不勾选 Add extended data type list，OK。
2
JGW_8v2_workflow_SYB
新建工区完成，主窗口变成【Jason-Jason_demo】
3
JGW_8v2_workflow_SYB
1
JGW_8v2_workflow_SYB
1. 新建工区
主窗口【Jason-***（工区名） 】 FileCreate project
弹出新工区命名窗口【Create JGW project】 在 Directory 填入工区名称（完整路径名） ：/data/jason_project/Jason_demo，OK。 以后加载的数据及运算产生的数据均保存在此文件夹中。 /data/jason_project 是存放 Jason 工区的目录，Jason_demo 是新建工区的名称。
【Disk SEG-Y Import-(Untitiled transaction)】 ParametersSelect/edit transport parameters 弹出窗口【Select/edit transport parameters】 选择 SEG-Y format 为 Standard disk SEG-Y； 选择地震数据类型：Time/Depth。如果选择 Time，单位默认为 ms； 选择地震资料维度：2D/3D； 点击 Add，选择要加载的地震数据； 在 JGW file 填入内部文件名，例如 demo；OK。
10
JGW_8v2_workflow_SYB
同样方法，分别将第四列、第五列、第六列、第七列设置为 CDP number、X coordinate、 Ycoordinate、Horizon value。OK。 注：第五列、第六列 XY 坐标可以选择为 Skidefinition】 层文件格式设置完成。OK。
2. 数据加载
2.1 地震数据加载
主窗口【Jason-Jason_demo】 DatalinksSeismicSEG-YDisk SEG-Y Import
弹出窗口【Disk SEG-Y Import-(Untitiled transaction)】 ParametersCreate/edit SEG-Y format definition
JGW_8v2_workflow_SYB
△根据 Fugro-Jason 公司 2011 年****Jason 软件培训讲座整理而成，适用于 JGW8v2。 △demo 数据由 Fugro-Jason 公司提供。
目录
1. 新建工区 .............................................................................................................................. 2 2. 数据加载 .............................................................................................................................. 4 2.1 地震数据加载 ................................................................................................................ 4 2.2 层位、断层数据加载 .................................................................................................... 8 2.3 井曲线加载 .................................................................................................................. 14 2.4 井分层加载 .................................................................................................................. 17 3. 子波估算和井震标定 ........................................................................................................ 18 3.1 计算时深关系曲线 ...................................................................................................... 18 3.2 初步估算子波、制作合成记录 .................................................................................. 21 3.3 确定子波极性 .............................................................................................................. 29 3.4 调整合成记录 .............................................................................................................. 30 3.5 子波估算 ....................................................................................................................... 31 3.5.1 选取时窗 ............................................................................................................... 31 3.5.2 子波估算 ................................................................................................................ 32 3.5.3 提取多井子波 ....................................................................................................... 35 4. 建立低频模型 .................................................................................................................... 38 5. 约束稀疏脉冲反演 ............................................................................................................ 45
8
JGW_8v2_workflow_SYB
弹出窗口【Import ASCII horizon files】
【Import ASCII horizon files】 点击 Set，弹出窗口【Select horizon definition】
9
JGW_8v2_workflow_SYB
【Select horizon definition】 一般情况下，不选择软件提供的格式 Horizon file format。 点击 Test file 后的 list，选择要加载的某个层位； 点击 Create column format，弹出窗口【Edit ascii horizon generic column format definition】 窗口左下角显示层位文件共分为 7 列。
【Edit SEG-Y format definition】 选择地震资料维度，2D/3D； 选择 CDP、Line number、X coordinate、Y coordinate 道头位置，一般为 21、17、73、77。 选择 XY unit 为 m。 选择要加载的地震数据文件。 OK。 为了确保道头位置正确，可点击 Show binary headers 查看道头信息。