JASON培训教材-14页解析

1、Jason地学综合研究平Jason软件培训资料Jason软件集合了油气勘探开发不同阶段的储层预测和油气藏描述技术， 它致力于各种资料、各种认识的全面综合，提供符合各种资料、各种认识的储层预测和油气藏描述结果。指导油气藏的勘探和开发，提高钻井成功率，降低风险。Environment Plus 运行环境及分析工具地质框嚮I型 wavelets 子波曲算地農度浦（JiimMwi精细皤层/?6气载就地质總计模拟 EunauonosJ 数据分痕换工M Largo测井a銭处理分折工具弹性反灌 Vielmeua 速度建模主要模块序 号软件模块功能描述1En virome nt- plus运行环境 及分析工具

2、集总了 Jason的运行环境，各种分析工具和 各种辅助性模块。如数据加载显示、（2D、3D ）分析解释、井编辑、合成记录标定、交 绘图、制图等。2EarthModel地质框架模型创建一个由地震坐标描述的、融合了构造（层位、断层）、地质/沉积模式、测井信息 和初始权重分布的，用于储层/油藏定量描述 的地质框架模型。3Wavelets子波估算提供用于合成记录标定和反演的地震子波 估算和理论子波计算的工具，在单井多井、 或无井的情况下，都可以由单道或多道地震 信息估算出最佳子波。4Inv ertrace/ Inv ertrace-plus地震反演在地质框架模型和测井资料控制约束下进 行地震反演并利用

3、反演声阻抗实现储层/油藏的描述和预测。5In verMod测井反演在地质框架模型控制下，用地震资料作约 束，以测井资料为主体的储层/油藏描述和预 测。6StatMod地质统计模拟以地质框架模型、测井和地震资料为基础， 以层为单位，利用储层/油藏参数的空间分布 规律和空间相关性进行随机模拟 /随机反演。 获得一组等概率的储层/油藏参数模型。7Fun cti onMod数据分析变换工具利用系统预定义或用户自定义的函数（表 格）将3维、2维和1维数据变换生成新的 数据。Jason软件是一套综合应用地震、测井和地质等资料解决油气勘探开发不同阶段储层预测和油气藏描述实际问题的综合平台。其中子波估算（Wa

4、velets）和层位标定、地质框架模型（Earthmodel）、地震反演（In vertrace、In vertrace plu* 测井反演（In vermod）、地质统计模拟（Statmod）和数据分析变换（Functionm od） 是主要模块和关键技术。下面根据实际工作步骤来介绍Jas on软件的主要模块和关键技术及应用注意事项。数据加载数据加载顺序为地震7层位7测井7其它（如人文、子波等）；输出可根据 需要有选择性地输出。注意事项：地震数据类型（是2D还是3D）、线道号和XY坐标在SEG-Y道 头中的正确位置、输入数据的字节数（至少为16位）。井数据输入文件的格式与所选的格式模板文件必

5、须一致包括输入 文件本身的声波和密度的单位（us/ft, us/m, g/cm3, kg/m3）、模板 文件中深度的类型（测量深度、TVD等）和单位（m, ft等）。子波估算和层位标定技术这部分工作是通过 Modeling下的 Wavelets 和Analysis下的 Well log editing and seismictie 两个模块完成的。通过子波估算和井曲线编辑的交互迭代，由井旁地震道和井中的阻抗曲线估算出与地震最佳匹配的地震子波。并实现子波估 算、合成记录的制作和层位标定。其技术特点是：同时估算子波的振幅谱和相位谱；子波估算和层位标定同时 完成；方法多样，可处理有井和无井、单井和多

6、井、直井和斜井；质量控制手段 多样。子波估算和层位标定技术的方法如下：1）计算理论子波（如 Ricker） （Wavelets 7dit7Create synthetic wavelet.2) 生成波阻抗(Al)曲线(Well log editing and seismic tie Calculate Al)。3）用理论子波和Al曲线在 Well log editing and seismic tie 中生成合成记录, 用Shift all in MD移动测井曲线（即合成记录）使得合成记录的主要同相轴与地震的对齐。将移动后的曲线存为新井。4）在Wavelets的Input中，选择3）中的新井，

7、给定合适的时窗范围，利用Wavelets 中的 Estimate 下的 Estimate wavelet amplitude and phae spectra 模块, 根据井曲线中的AI曲线和井旁地震道估算子波。4）中估算的子波生成新/压缩等）使得更多的同相轴对齐，将经过编辑的曲线再存为新井，再返回到4）。重复4）和5）直到5) 回至U Well log editing and seismic tie 中，选择 的合成记录。在时间曲线的控制下，编辑测井曲线（拉伸获得最佳的子波和最佳的匹配。这时用于反演和制作合成地震的子波以及井的时 深关系就算确定（如图1）。商mate wavElEt am p

8、l il ude and pNdse spwH匚.Mrl iS 4 THILk± .tfhiR.LuLAB子波估算卜IfifirddnLC * CLi厨MrHida P-Inp«ddfic» FurKtiofis出Ml卽 iirt Oi'jw 城SbH'i 1 呻 却 hbShift nil in MUXH井曲线编辑hLAlgtu'wise朋m砂 朋rt no_l LI . I . Iss:般;V-? UJM”L叫m叫IPIt打jo”UAlhbd*Lh-hIkIhH-EULfaUt atirrJdLui IW l-J*'养Femnc

9、 ucLI topCofjy lot亡JI丄口沪-:.- :I*r:CPrLr binJMPrla- umid.h m rtt kW "IvL11 L-Zft Ftim girt 诸1片空MBSCut 411 Id翠子波CiLz h HhLnbhX Ufl">-F Rcfuvc- sectifln iiJl Jofs)r_!|_ T_ iT亠/ r 1.I T-川UndJiL 曲"IVflbi LIIFIbpL f-T ! 1RM ucll top rww rwll top:bnift veII tcpdti di" Fhh _Tlh Hifl B

10、araXiBT, rsf-a ,-bhiLuu>ru*子波估算tcii/sqiitLzrStretch/呻iH in MP EMLIJ电 gff保存编辑后的井讪吊Ddit'ibtiR 恤ibMk屮, hP4FTP-IfiPMMbwwe * 仇u业L7'd打图1子波估算和井曲线编辑交互迭代过程对于多井，可以用1） 5）的方法估算每口井的子波和确定时深关系，然也可以选一口井孔环后再用平均子波对每口井的标定结果和时深关系进行调整;境和曲线质量较高的井，用 1） 5）的方法估算一个最佳子波，再将这个子波用于其它井并使合成记录与地震记录达到最佳匹配。在这一交互迭代的过程中，要注意以

11、下几点:Wavelets中的时窗选择对子波影响极大，一般所选时窗应包含目的层段,其上下时间应选在地震波形变化较小的位置，时窗长度应为子波长度的3 倍左右。Well log editing and seismic tie 中的曲线拉伸和压缩尽量在较大范围内 进行，拉伸 /压缩不应改变曲线的变化趋势，由时间线进行监视 。估算的子波振幅谱应与地震的尽量接近，相位谱在有效频带内变化尽量 小。地质框架模型（ Earthmodel）地质框架模型是由 Modeling 下的 Earthmodel 模块完成。 它由五个子模块构 成：1）Model builder with TDC2）Model builder

12、 without TDC3）Model generator4）Model Interpolator5）Well Curve generatorModel builder with TDC 和 Model builder without TDC 的作用基本相同，即 利用地震解释层位 /断层，根据 framework table 描述的接触关系形成具有地质意 义的三维密封的框架模型， 同时形成以层为单位的沿层内插测井曲线的权重系 数。其差别是 MB with TDC 在时/深域对模型进行同时描述，即存在时深关系 ； 而 MB without TDC 要么在时间域要么在深度域对模型进行描述，即没有时

13、深关 系。Model generator就是依据Model builder的描述信息，生成具体的测井曲线 模型。 Model Interpolator 就是将已存在的地震或属性模型按 Model builder 的描述 信息进行内插或重新网络化。 Well Curve generator 生成测井曲线。其技术特点是：建立测井分层和地震解释层位的联系并将地震解释层位/断是关键和一行代表2）。层内插成密闭的层位 ；用以层为单位的沿层内插测井曲线的权重系数、 层厚等参 数描述的初始模型；可处理正断层和逆断层；可描述各种地层接触关系（整合、 不整合、河道等） 。在Earthmodel中，编辑描述地质框

14、架结构表 （即framework table） 难点。Framework table用来描述地层的接触关系以及微层的内插模式， 一个地层，并由老到新依次定义地层之间的接触关系和内插模式（如图Framework table的建立是通过图中 Edit下的 Append row on top or insert row 菜单 实现。两者均弹出Edit row的界面（如图2）。Framework table的编辑界面（如图2）。其中In dex层位或断层的标识号，每个层位或断层都有唯一的标识号，且从下往上依次为0, 1, 2,。Primary fit In terface (PEI)定义层的顶界面并用

15、来指导层内的内插 模式Stratigra phy定义由Primary fit In terface确定地层内微层的内插模 式以及层顶和它下面的地层关系。定义模式包括：proportional to topand base, prop orti onal to top, channel 和 portional to base, truncation, reef,图2地质框架结构表及所描述的地质现象Base In dex用来确定削截位置和断层位置，也就是说削截和错段 只发生在Base Index所指定的层位上面。Trend fit interface(TFI)选择项，如果选择TF

16、I，由PFI定义的地层 内插模式将由平行于TFI的模式所取代。Seco ndary fit In terface(SFI)选择项，在顶面部分被剥蚀的情况， 可以指定该顶面下部的未被剥蚀的层面，作为SFI，以完善层面的几何形态并指导内插。编辑地质框架结构表时要注意以下几点：从底层开始向上逐层编辑。先建断层下盘的地层，后建断层上盘的地层。被断层切割的层不能作为 datum。四、地震反演(Invertrace 禾口 Invertrace plus)Jason中地震反演 以约束稀疏脉 冲反演 Constrained sparse spik In versio n CSSI）为主，另外还提供了递归道反演

17、（Recursive trace In version RTI） 和道合并（Trace mergingAITM ）即加低频等模块。实际中大量使用的是约束 稀疏脉冲反演和道合并。其优势在于：趋势和约束可为常数（1D）、沿层变化（2D） 或空变（3D）;提供反演算法中各种参数选择的质量控制图版；低频分量遵从地 质模型。适用于勘探开发的各个阶段。约束稀疏脉冲反演（CSSI）的关键是：编辑趋势（Edit trend）约束稀疏脉冲反演中的趋势是指研究区块的波阻抗变化背景（如图3）,它 是 CSSI 的初始模型。In vertrace 可选 con sta nt 或者 Laterally varyi ng

18、 , Co nsta nt 是 指反演所用的趋势在每个层上为常数；Laterally varying是指反演所用的趋势在每个层上随XY变化。In vertrace plus要选择Model file作为趋势。它表明趋势在 空间每一点都可能发生变化。对于Constant,可以用手编辑，也可以由井计算，用From logs即可；Laterally varying，要用波阻抗的层位文件；Model file,要用波阻抗的*.mod文 件。编辑约束（Edit hard trend constraints ）约束限定CSSI反演中波阻的变化范围。因此编辑约束时必须使得约束的最 大最小值为相应层位的波阻抗

19、的最大最小值（如图3）。Constraints is:其选项和含义与 Edit trend中的trend is相同。图3编辑趋势和约束口nrtiTjinI- cdMorijiHJt CLirl«Vp-lrrl二1 Jontrairt hdilnr LKjflLbtn:H. Mil"iakKtL_ C«fiitr_ Ie:JLpw ufiiLr-ILIrUflh-iL Lil. All呼皿Tvw K-J LrChkiIkbikplambda 值Lambda是反射系数（阻抗）和地震数据匹配的平衡因子。低lambda值将导 致较少的反射脉冲，阻抗模型就缺少变化细节。太小

20、的lambda值导致非常差的地震数据匹配；太大的lambda值可能导致虚假的反射脉冲或变化细节。所以要 选择合适的lambda值使得反射系数的匹配和地震的匹配都比较好。Trace Merging就是将来自于Earthmodel （测井资料）或Velmod（地震速度） 的阻抗模型中的低频成份与来自 CSSI反演的阻抗模型合并形成一个绝对阻抗模 型，便于岩性、孔隙度等方面的研究。在滤波器设计中要确定绝对阻抗模型（合 并模型）的哪一部分频率成份来自低频模型，哪一部分频率来自CSSI反演模型。五、 测井反演（Invermod）Invermod的目的是利用测井和地震资料建立基于地质框架模型的高分辨率储层

21、/油气藏模型。基本算法包括 Geologic In version GI和Seismic characterIn version SCI。GI和SCI也可以联合使用 （Both ）。基本思路是通过井孔位 置的信息（波阻抗或地震道）的线性组合形成合成地震道，比较合成地震道与实 测地震道，修改权重系数，使二者达到最佳匹配，再利用这个权重系数 内插测井 曲线形成相应的反映储层/油气藏特征的数据模型。其特点是：遵从地质框架模型、分辨率高； 地震资料对井间的变化起一定的控制作用。 适用于地质情况较为 简单、测井资料质量好、井较多的滚动勘探开发阶段。GI 和 SCI 的差别就在形成合成地震道的方法不同。前

22、者利用阻抗曲线(井 点)和子波生成合成地震道，后者直接利用井点位置的地震道作为合成地震道。Invermod 由四部分组成：1)2)3)4)Principal component analysisModel estimati onModel gen erator Well curve gen erator其中 Model generator和 well curve generator与 Earthmodel 中基本相同。六、 随机模拟 /随机反演( Statmod)Statmod 是 Jason 中 的 随 机 模 拟 软 件 包 ， 有 Statmod-Analysis 和 Statmod-M

23、odeling 两部分组成。随机模拟技术充分利用地震资料横向分辨率高 (采样密集)、测井资料纵向分辨率高的特点， 解决地质上参数分布范围重叠或 多值对应关系的问题，得到高分辨率的参数模型 。Statmod-Analysis 用来分析数据，形成基于数据的直方图和变换(histograms and transformS 以及变异函数(Variograms)。Statmod-Modeling 进行随机模拟和随机反演。 随机模拟的算法包括： 克 里金(kinging),协克里金(Co-Kriging),序贯高斯模拟(SGS Sequential Gaussian Simulation) ， 序 贯 高

24、 斯 协 模 拟 (SGCCSSequential Gaussion Collocated Co-Simulation)，序贯 门槛值协模拟(STISSequential Threshold Indicator Simulation)，序贯指示模拟(SIS Sequential Indicator Simalation)，带趋势序贯 指示模拟 (SIS With trendSequential Indicator Simulation with trend)，随机反演中包括最快速下降法(greedy)和模拟退火法(Simulated Annealing)。Seque ntial Gaussia

25、 n Simulati on21.9*aCon diti oning data17.126.2o11.221.9*Ran dom locati onbdc.5.Estimate the local con diti onal p robability distributio n (lc pd) at that locati onDraw at ran dom a si ngle value from the IcpdIn clude the n ewly simulated value in the set of con diti oning dataRe

26、p eat ste ps 1 through 4 un til all grid no des have a simulated valueChoose at ran dom a grid node at which we have not yet simulated a value希细解择测井曲统岩性曲线储层物弊模型直方動 变异函数地甑演1 统计特征分析(Statmod-Analysis)Statmod的所有模拟技术(除指示类 SIS外)都假定输入数据为正态分布。 然而，实际数据并非都为正态分布。因此Statmod-Analysis首先要分析输入数据 的分布特征(histogram)，找到一

27、种变换将非正态分布的数据映射为正态分布 (Transform)。然后再分析数据的空间相关性，用合适的函数拟合得到变异函数 (Variogram)。Statmod-Analysis就是要根据实际数据的分布特征(直方图)和 离散变异图寻找合适的变换(transform)和变异函数(variogram)。数据分布直方图和变换( Data for histograms and tramsforms) 根据随机模拟要模拟的储层参数和约束参数选择数据并产生相应的数据分 布直方图。 再选择合适的变换类型拟合数据分布直方图。 提供三种变换类型 (即 Gaussian, Log-Gaussion和Table)。

28、随机模拟时，用拟合数据的分布函数把实际 数据变换成正态分布。数据离散变异图和变异函数( Data for variogram sampling and modeling) 数据选择与Data for histograms and tramsforms相同。注意所选数据必须有变 换，即已作了 Data for histograms and tramsforms数据空间相关性是由变异函数(Variogram)描述的。变异函数相当于相关函数的倒数，即变异函数的值越大相关性越差 。变异函数能同时描述区域变量的随机性和结构性， 能从数学上对区域变量进行严格分析， 是空间变异规律和空间结构分析的有 效工具

29、。变异函数是通过合适的连续函数拟合数据的离散变异图而得到。 Jason中可用的连续函数包括高斯函数(Gaussian)、球函数(Spherical),指数函数(Exponential)和块金值(Nugget)。变异函数中的程高(Sill)和变程(range) 反映了数据的离散程度和空间相关性。分析确定数据分布(histogram)和变异函数(variogram)时，要注意以下 几点。分层段进行。选择较小的数据采样间隔( 1ms)。选择合适的 分析间距 (histogram: Nr of interval)和相关间距 (variogram： Nr of lags， interval)。选择合适的

30、程高(Sill)和变程(Range)。程高选在离散变异图散 点变平的位置(如图4)。变程选为程高所对应的滞后距离(Lagdistance) (如图 4)。SillNugget Evponeniiall Sp herical GaussianLag distance地震、地质等)和由3D储层/油气藏模型。可<Range图4变异函数2 随机模拟(StatMod-Modeling)StatMod-Modeling就是要利用确定性信息(测井、 StatMod-Analysis导出的不确定性地质统计信息来创建 用的技术包括：1)2)3)4)或者它们的组合。要强调的是，Krigi ngSeque n

31、tial Simulati onSeque ntial In dicator Simulatio nGeostatistical Inv ersi on在运行 StatMod-Modeling 之前，必须要用 StatMod-Analysis获得要模拟层的属性、岩性的变换和变异函数。关于每种方法所需的输入参数如 下表。模拟方法确定性信息不确定性信息Gaussia n Simulati onsolid model、layers、primary data lithology masks(可选)、trace gateTran sform for p rimary dataVariogram for

32、p rimary dataGaussia n co-Simulati onSolid model、layers、primary data、Secondary data、lithology mask(可选)、trace gateTran sform for p rimary dataTran sform for sec on dary dataVariogram for p rimary dataVariogram for sec on dary dataGaussia n collocated co-simulati on同 Gaussian co-simulationTran sform f

33、or p rimary data Tran sform for sec on dary data Vario ngram for p rimary dataCorrelati on coeffice nt of p rimary and sec on dary dataTreshold In dicator simulatio nsolid model、layers、primary data、(应为门槛区间分类数据)lithology masks(可选)、trace gateThreshold variogram/variogramIn dicator simulatio nsolid mod

34、el、layers、lithology data、trace gateIn dicator variogramLithology p robability(c on sta nt value)In dicator simulatio nwith trend同上In dicator variogramLithology p robability(horiz on or model)Run Inversionsolid model、layers、primary data、 lithology masks(可选)、seismic data、 Wavelet、conversion table、trac

35、egate、 Time gateTran sform for p rimary dataVariogram for p rimary data注：solid model 地质框架模型，来自 model builder。layers 由 model builder 中 framework table 所定义。只对所选岩性中的参数（如砂岩中的孔隙度）进行模拟。 指与primary data相关的参数，但分布密集，多为地震反primary data 要模拟的参数（如孔隙度），但分布稀疏，多为井数据。lithology mask sec on dary data指岩性曲线0演波阻抗数据。Litholo

36、gy data Lithology Simulation (Indicator Simulation with trend)Define the trend models which define the occurrence p robability of sand and shale based on AI.Determine the Variogram function to fit the real dataUse Indicator Simulation to get equal p robability lithology models.Use FunctionMod to gen

37、erate the confident model based on the equal p robability lithology models.Use different cutoff values to get lithology models with different confidence.? Define the trend models which define the occurrence p robability of sand and shale based on AI1> Determine the max. / min. Al values to define

38、 clean sand or shale (An alysis - Cross pl ots and Histogram).2> Determine the probability of sand/shale occurrenee in overlap zone (Create or edit table)3> Use FunctionMod to tran sfer Al model to p robability model.? Determine the Variogram function to fit the real data1> Check the well log to make sure there are the lithology logs in all available well. If not, please create or input (FunctionMod or DataLinks)2> Use StatMod - Analysis Module to generate Variogram function.a. Determine the layer you are interested. Select reas