Jason约束稀疏脉冲反演流程

地震反演是储层预测和油藏描述中非常关键的技术。

地震反演技术将地震数据、地质构造模型、井点测井资料有机地结合为一体；把地震资料的常规界面反射型剖面，反演成能与测井资料直接对比的岩层型波阻抗剖面,以及岩性剖面和孔隙度等地层物性剖面，反演结果分辨储层的能力也有了明显提高。

Jason稀疏脉冲波阻抗反演流程(附流程图)1 子波提取子波长度的选取要适宜，一般为80ms左右，在浅层，地震频带较宽，子波可短些；深层，地震频带窄，子波可略长些；对提取的子波，要从子波的波形、振幅谱、相位谱等方面进行子波质量的判断，要求子波的波形稳定、频域内单峰平滑、有效频带内相位稳定。

2 合成地震记录标定反演过程中，可以利用反演波阻抗数据，对地震地质标定进行精细微调，提高标定的质量。

子波提取与合成地震记录标定有着密切的关系，两者是循环反复进行的。

先用初始理论子波进行初步标定，标定的时深关系基本确定后，再用井旁地震道和测井反射系数重新估算子波。

完成最终的地震地质标定之后，再重新提取更精确合理的子波。

3地质层位解释解释的层位要尽量靠近反演目的层，要十分精细，反映出目的储层的空间变化特点来。

同时，使用时还要对层位中的缺失、奇异点进行修补、剔除，保证反演所用层位的正确性。

4 地震、地质、测井约束模型的建立建立尽可能接近实际地层沉积情况的波阻抗约束模型，是减少反演最终结果多解性的十分重要的一环。

建立波阻抗模型的过程实际上就是把地震界面信息与测井波阻抗信息正确结合起来的过程。

5 反演参数选择选择合理的反演参数是做好反演处理的前提，稀疏脉冲反演中的主要参数有稀疏约束项、拟合误差约束项以及模型约束项等等。

针对不同的地质情况，有目的性地选取合适的反演参数，同时利用相应的质量监控手段，加强对反演信息的反馈处理，可以有效地改进反演结果的质量，起到事半功倍的效果。

6反演结果检查处理过程中的反演数据体反复检查、各种资料对比分析，发现异常、寻找原因、反复反演，知道获得满意的结果为止。

叠后约束稀疏脉冲反演中文培训教程一、叠后约束稀疏脉冲反演的基本原理叠后约束稀疏脉冲反演是一种地震反演方法，它利用地震波传播和反射的物理过程，通过对地震数据进行逆向建模，来获取地下介质的信息。

该方法的核心思想是通过寻找一种最优的模型，使得通过该模型建立的地震波场与观测到的地震数据尽可能吻合。

在叠后约束稀疏脉冲反演中，通常将地下介质的性质用参数化的方式表示，将地震波传播和反射的物理过程用数学模型描述，通过优化方法来求解最优的模型参数。

这些模型参数可以包括地下介质的速度、密度、衰减系数等参数，通过对这些参数的反演，可以获取到地下结构的详细信息。

叠后约束稀疏脉冲反演的关键之处在于它采用了稀疏性约束和地震波形约束。

稀疏性约束是指地下介质的参数在某种表示下是稀疏的，即它们可以被较少的非零系数表示。

地震波形约束是指通过对地震波传播和反射的物理过程进行数学建模，来约束模型参数的取值范围，使得反演得到的模型更符合地震数据的实际观测。

二、叠后约束稀疏脉冲反演的算法流程叠后约束稀疏脉冲反演的算法流程通常包括以下几个关键步骤：1. 数据预处理：对地震数据进行预处理，包括数据去噪、时域频域滤波等操作，以减小噪声对反演结果的影响。

2. 构建反演模型：将地下介质的性质参数化表示，并建立地震波传播和反射的数学模型。

3. 求解反演问题：利用最优化方法，求解使得通过模型建立的地震波场与实际观测数据吻合的最优模型参数。

4. 结果评估：对反演结果进行评估，包括地下结构的成像、模型参数的稳定性等。

5. 反演结果应用：将反演得到的地下结构信息应用于油气勘探、地质勘探等领域。

三、叠后约束稀疏脉冲反演的实际应用叠后约束稀疏脉冲反演在油气勘探、地质灾害预测、地下水资源开发等领域具有重要的应用价值。

在油气勘探中，通过叠后约束稀疏脉冲反演可以获取到油气藏的地下结构信息，帮助勘探人员准确地定位油气藏的位置、形态及大小。

在地质灾害预测中，叠后约束稀疏脉冲反演可以提供地下结构信息，帮助人们预测地震、泥石流等自然灾害的发生概率以及可能的影响范围。

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地学平台是一个将各学科有用信息综合在一起，即为用户提供丰富的、从油气勘探、开发到生产阶段将地震、测井和地质信息有机结合，进行反演、模拟和预测不同阶段储层岩石物性和流体变化的综合研究工具。

您根据可用信息和经济条件，您可以明智的选择应用哪个“工作流程”更好。

1、Enveriment 运行环境（工区管理地震解释和接口）与其它解释系统紧密相连、数据格式灵活、方便。

2、EarthMod 地质框架模型，包括模型建造器，模型生成器和测井曲线生成器。

把地质、地震、测井曲线（斜井轨迹）、沉积模式（整合，不整合，礁体及河道等）建立以地震道为坐标的模型上，成为储层预测的基础，并赋予地质意义。

3、VelMod 速度模型建立工具VelMod 除具有建立速度模型的基本功能外，且有其独特特征：•地质引导：结合地质信息进行插值运算，能把稀疏的速度数据拓展为三维速度场。

•基于地质概念的速度调整•真三维插值4 、Wavelets 二维和三维地震子波分析●无井、单井、多井提子波和空变子波分析●从斜井中进行真三维分析●估算信噪比谱●为使结果稳定利用模型驱动5、 Modtrace 完整的地震反演系统( 包括对大倾角地层)●所用约束稀疏脉冲反演是基于道的递归反演，产生地震带宽内的波阻抗数据。

●其算法的唯一特征应用解释控制，为获得可靠的低频信息提供工具。

在EarthMod 的基础上建立空变约束条件，其结果突出各向异性，提高了反演分辨率。

●多种方式提取子波，包括无井、单井多道、多井多道、直井、斜井提子波和在EarthMod 的基础上提空变子波。

●由于大量地震信息的反演是基于EarthMod 的基础，所以适用与多条复杂断层的地质情况，使其反演结果更接近真实的地质模型。

●给出产层有效厚度图和孔隙度分布图６、InverMod 精细储层描述（多井）•采用一欧洲石油公司的专利，专门致力于薄层预测和精细描述。

对油田滚动勘探开发十分有用。

第四篇 约束稀疏脉冲波阻抗反演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）右边分别是反射系数的绝对值和，合成声波记录与原始地震数据的差值，趋势约束项。

约束稀疏脉冲反演基于地震数据，用波阻抗趋势进行约束，以期用最少的反射系数脉冲，对合成记录与地震道进行最佳匹配。

约束稀疏脉冲反演方法在识别煤层顶板砂岩中的应用说明了稀疏脉冲反演方法的原理、分析了反演过程中需要准备的基础资料和需要注意的关键环节，采用该方法对研究区煤层顶板砂岩进行预测，获得了砂岩的分布规律，经钻井验证可知效果较好。

标签：煤层；砂岩；反演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是目前最常用的反演软件，它操作上的特点是它需要什么数据或参数就给它什么数据或参数。

下面是它的主窗口（图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。