您的位置:主页 > 成功案例 > 斯伦贝谢软件使用技巧 > GeoFrame操作技巧 >

成功案例

GeoFrame操作技巧>>

39.GeoFrame预测砂岩厚度工作流程

基本思路是,使用地震属性包Seismic Attribute Toolkit(简称SATK)里的算法Cosine Correlation (Freq Indexed)提取分频数据体,进一步在该属性体基础上提取出调谐频率属性。使用调谐频率属性与砂体厚度建立关系,用调谐频率来预测砂体厚度。LPM 和SeisClass 帮助我们分析分频结果并与井数据相关,详细流程见图1。

下面举例说明如何使用该workflow预测砂岩厚度。

一,使用地震属性包SATK里的Spectral Decomposition(CCT) 功能提取分频属性体。

1 打开Seismic Attribute Tookit (SATK);

2 点击Survey 并输入新的Output Class,如图2所示;

3 如图3所示,打开Spectral Decomposition; (CCT)

4 双击Cosine Correlation(Frequency Indexed)设置CCT参数;

5 开时窗在砂岩顶上10ms和砂岩下10ms之间。

6 打开Output Volume,Output to Grid。

Output volume项将以目的层为中心提取分频体,提取的结果是一个以时间为索引的多频体。也就是说,显示在Seis3DV里时垂向上用的是时间,实际上是频率,借助时间来表达频率。每个道上的振幅值还是地震振幅。由于在频率域表达比在常用的时间域要复杂,所以借助时间来表达不同频率所对应的振幅值。

Output to Grid将按照参数输出等频的频率网格。即对所选时窗内的地震属性进行统计,分别输出某一频率对应的振幅变化。

7 设置频率范围从0到75,增量为5。

8 选择Houston_ASAP作为Output Horizon,参数设置如图3。

9 点击run,计算属性。

10 在Geoviz或者Seis3DV里浏览计算得到的分频属性体,如图4所示。

11 在Basemap或者SeisClass Explore里查看不同频率的平面图,如图5所示。

二,提取频率时间。在分频属性体基础上提取调谐频率。

1 对CCT生成的时间索引的频率体计算Time at Max Amplitude属性。该属性得到当振幅达到最大时所对应的时间。

前边已经解释过,分频计算得到的是时间索引的频率体,所以这里得到的时间,实际上对应的是当振幅达到最大时(发生谐振)的频率即调谐频率,不过是以时间的方式表达的,需要进一步换算成频率。

2 设置Survey Input Class为CCT_cosine_freq,在该数据对应的时窗范围,例如0-300ms (0-75HZ)范围内计算调谐时间属性网格(实际是调谐频率)。参数设置如图6。

可以在Basemap里显示该属性网格。

对于现在的结果,我们有两个问题,一是这个调谐时间属性网格有异常值需要去掉。二是调谐时间属性网格必须除以采样率来转换成真正的调谐频率属性网格。

3 去掉异常数据。

1) 从Seismic Catalog里,选择Interpretation Manager > Grid Compution。

2) 打开Equation Operations并设置公式让小于0的为Null。见图7。

3) 保存输出的新网格,参数设置为Null值,并定义好属性的单位。见图8。

4 将调谐时间属性网格转换成调谐频率属性网格。

1)将调谐频率属性网格创建一个公式除以采样率(4ms),将时间转换为频率值,得到真正的调谐频率,见图9。

2)输出结果作为一个新的网格命名为调谐频率,单位为频率,属性为相应频率,见图10。

结果显示如图11:

3) 使用Surface Model Editor把TUNNING_FREQUENCY网格平滑和插值。GeoFrame Basemap > Gridding > Surface Model Editor。

三,计算井点处储层物性参数。

使用Geology Office里的Reservoir Summation (Resum)模块计算井点处砂岩厚度。详细步骤参考Geology office培训手册。

四,井震结合预测砂岩厚度。

使用Log Property Mapping (LPM)分析地震调谐频率属性与井点处砂层厚度之间的关系,进行储层预测。

1 打开LPM,点击Main Selection tab,显示如图12。

2. 打开Data Analysis Tab,选择Net_Thickness物性和Response Frequency属性。如图13所示。

3 点击Edit Intersection,设置默认半径为100feet。

4 计算Regression Coefficient,计算回归系数。

5 点击Regression Coefficient Matrix去显示Net_Thickness和Response Frequency的交绘。

6 显示Regression Line,如图14所示。

虽然相关性十分好,通过频率体(绿线)评估厚度值,有些井可能不十分的匹配,可能是由于有些井的储层很薄引起的。

7 在Data Analysis里,计算Calibration function在Net_Thickness和Response Frequency之间的回归方程。

8 打开Data population界面,LPM根据上一步得到的回归方程,根据井点上的砂岩厚度来预测整个平面上的砂层厚度。

9 误差校正。

点击Calculate Residual Scatter旁边的Compute按钮,计算井点处砂岩厚度值与预测值之间的误差。

点击Calculate Residual Grid旁边的Compute按钮计算误差网格。

点击Apply Residual旁边的Compute按钮计算校正后的砂岩厚度网格。

显示计算结果NET_THICKNESS_MAPPED_WITH_SEISMIC在Basemap,并投上Houstion_ASAP的构造等值线。

这样就可以看到构造有利部位以及储层物性比较好的部位的分布,综合评价油藏特性。

地址:北京市朝阳区酒仙桥路20号颐堤港一座4层

邮编:100015 传真:010-64309502

Copyright © 2002-2024 斯伦贝谢科技服务(北京)有限公司  京ICP备 06043577 号-1