论文部分内容阅读
[摘 要]砂巖作为一种有利的储层在工区中如何展布这一问题十分重要。本文以吉林油田大南洼陷为例，从井出发，统计工区中每口井相应深度的砂岩厚度，依据从反演波阻抗数据体中提取的层间均方根振幅属性，利用序贯高斯模拟方法对井间的砂岩厚度进行预测，达到了良好的效果。通过对不同层段砂岩厚度的计算，对工区内扇体的展布规律和发育过程有了一定的了解。
[关键词]砂岩厚度；波阻抗数据体；层间属性；序贯高斯模拟
中图分类号：P588.21+2.3 文献标识码：A 文章编号：1009-914X（2014）48-0011-01
引言
大南洼陷位于伊通盆地中部，西部受西北缘断层遮挡，北部受二号断层遮挡，大段泥岩发育，埋深条件有利，是良好的油气供给区。受边界断层影响，盆地边缘较陡，主要发育冲积扇、水下扇等沉积类型，扇体中的砂体成为有利的储层，因此研究工区中砂体的展布特征就显得尤为重要。
随机模拟的方法有很多，各有其优缺点和适用范围。其中序贯高斯模拟算法（Sequential Gaussian Simulation，SGSIM）是一种连续性变量的随机模拟方法，其应用十分广泛。SGSIM被认为是模拟连续性变量的首选方法，SGSIM的模拟结果比确定性算法更符合地质规律[1]。
1 砂岩厚度平面预测
本次研究的思路是首先统计工区中不同层段内每口井的砂岩厚度，作为井资料。然后对反演数据体提取层间的RMS（均方根振幅）属性，作为地震资料，再使用序贯高斯模拟法对井间的砂岩厚度进行预测。
1.1 统计井点处的砂岩厚度
大南洼陷主要使用的地震层位有Tf、Tfp、Th、Tg，分别对应奢岭组一段，双阳组三段、双阳组二段和双阳组一段（图1），依据各井的卡奔图，使用excel表格对工区中的51口井分别进行了统计（图2），除个别井资料不全外，其余井都有第一手的数据。
Tf到Tfp层段内各井点处的砂岩厚度如下（图3）
2.2 从波阻抗数据体中提取层间属性
使用波阻抗数据体是考虑到波阻抗相比原始地震数据有更丰富的地质含义，通过使用纵波阻抗也可以在一定程度上区分岩性。一般来说砂岩的阻抗高，泥岩的阻抗低。从地震属性的基本定义出发，Brown将地震属性分为4类：时间属性、振幅属性、频率属性和吸收衰减属性。源于振幅的属性可能提供地层和储层信息[2]。而均方根振幅属性对振幅变化非常敏感，可以用来识别振幅异常或描述层序。
2.3 序贯高斯模拟法
序贯高斯模拟有着计算速度快，方法简单的特点，适合于模拟一些中间值很连续、极端值很分散的物性参数。下图是序贯高斯模拟的基本流程（图5）[1]。
通过序贯高斯模拟我们得到Tf-Tfp层段砂体厚度平面图（图6），从图中可以看出序贯高斯模拟以井点处为硬约束，地震属性为软约束，非线性的对井间砂岩的厚度进行预测，根据色标值的深浅，我们可以很明显的看出紫色的部分砂岩较厚，有从边界向中心展布的趋势，符合扇体发育的基本特征。
2 砂岩厚度平面分析
依据工区的地层表，我们分别对Tf-Tfp、Tfp-Th、Th-Tg层段的砂岩厚度进行了预测，Th-Tg层段最早沉积，从下图（图7）可以看出在双阳组一段时，砂体的厚度范围在0到105米之间，冲积扇或水下扇携带着泥沙进入盆地沉积，该时期扇体普遍发育，主要方向为东北向、东南向，西南向，如箭头所示。
到了双阳组二段时（图8），砂岩的厚度范围增大，在0到228米之间。显示出该时期水流携带的泥沙量增多，搬运作用增强，把泥沙搬运到了盆地的中央地区。相比双一段时扇体发育的方向，双二段扇体主要发育在西北向、西南向，如箭头所示。
在双阳组三段时（图9），砂岩的厚度范围减少，在0到92米之间。各个方向均有扇体发育，东北向、东南向的水流重新活跃，但是水流能量较弱，扇体的展布未能到达盆地的中心。
3 结论
大南洼陷各个时期的沉积过程比较杂乱，扇体叠置严重，以往我们对此问题认识不清。本次研究的思路就是从平面出发，直观的反映各个时期砂体的分布规律，为我们刻画扇体和了解各个时期扇体的发育情况创造了条件。通过序贯高斯模拟这一方法，我们可以较为准确的预测井间的砂体厚度，为接下来的分析工作奠定了坚实的基础。
参考文献
[1] 张晓坤，随机模拟在三维地质建模中的应用：[硕士学位论文].北京：中国地质大学，2009.
[2] 杨占龙，彭立才，等.地震属性分析与岩性油气藏勘探[J].石油物探，2007，46（2）：131-136.
[关键词]砂岩厚度；波阻抗数据体；层间属性；序贯高斯模拟
中图分类号：P588.21+2.3 文献标识码：A 文章编号：1009-914X（2014）48-0011-01
引言
大南洼陷位于伊通盆地中部，西部受西北缘断层遮挡，北部受二号断层遮挡，大段泥岩发育，埋深条件有利，是良好的油气供给区。受边界断层影响，盆地边缘较陡，主要发育冲积扇、水下扇等沉积类型，扇体中的砂体成为有利的储层，因此研究工区中砂体的展布特征就显得尤为重要。
随机模拟的方法有很多，各有其优缺点和适用范围。其中序贯高斯模拟算法（Sequential Gaussian Simulation，SGSIM）是一种连续性变量的随机模拟方法，其应用十分广泛。SGSIM被认为是模拟连续性变量的首选方法，SGSIM的模拟结果比确定性算法更符合地质规律[1]。
1 砂岩厚度平面预测
本次研究的思路是首先统计工区中不同层段内每口井的砂岩厚度，作为井资料。然后对反演数据体提取层间的RMS（均方根振幅）属性，作为地震资料，再使用序贯高斯模拟法对井间的砂岩厚度进行预测。
1.1 统计井点处的砂岩厚度
大南洼陷主要使用的地震层位有Tf、Tfp、Th、Tg，分别对应奢岭组一段，双阳组三段、双阳组二段和双阳组一段（图1），依据各井的卡奔图，使用excel表格对工区中的51口井分别进行了统计（图2），除个别井资料不全外，其余井都有第一手的数据。
Tf到Tfp层段内各井点处的砂岩厚度如下（图3）
2.2 从波阻抗数据体中提取层间属性
使用波阻抗数据体是考虑到波阻抗相比原始地震数据有更丰富的地质含义，通过使用纵波阻抗也可以在一定程度上区分岩性。一般来说砂岩的阻抗高，泥岩的阻抗低。从地震属性的基本定义出发，Brown将地震属性分为4类：时间属性、振幅属性、频率属性和吸收衰减属性。源于振幅的属性可能提供地层和储层信息[2]。而均方根振幅属性对振幅变化非常敏感，可以用来识别振幅异常或描述层序。
2.3 序贯高斯模拟法
序贯高斯模拟有着计算速度快，方法简单的特点，适合于模拟一些中间值很连续、极端值很分散的物性参数。下图是序贯高斯模拟的基本流程（图5）[1]。
通过序贯高斯模拟我们得到Tf-Tfp层段砂体厚度平面图（图6），从图中可以看出序贯高斯模拟以井点处为硬约束，地震属性为软约束，非线性的对井间砂岩的厚度进行预测，根据色标值的深浅，我们可以很明显的看出紫色的部分砂岩较厚，有从边界向中心展布的趋势，符合扇体发育的基本特征。
2 砂岩厚度平面分析
依据工区的地层表，我们分别对Tf-Tfp、Tfp-Th、Th-Tg层段的砂岩厚度进行了预测，Th-Tg层段最早沉积，从下图（图7）可以看出在双阳组一段时，砂体的厚度范围在0到105米之间，冲积扇或水下扇携带着泥沙进入盆地沉积，该时期扇体普遍发育，主要方向为东北向、东南向，西南向，如箭头所示。
到了双阳组二段时（图8），砂岩的厚度范围增大，在0到228米之间。显示出该时期水流携带的泥沙量增多，搬运作用增强，把泥沙搬运到了盆地的中央地区。相比双一段时扇体发育的方向，双二段扇体主要发育在西北向、西南向，如箭头所示。
在双阳组三段时（图9），砂岩的厚度范围减少，在0到92米之间。各个方向均有扇体发育，东北向、东南向的水流重新活跃，但是水流能量较弱，扇体的展布未能到达盆地的中心。
3 结论
大南洼陷各个时期的沉积过程比较杂乱，扇体叠置严重，以往我们对此问题认识不清。本次研究的思路就是从平面出发，直观的反映各个时期砂体的分布规律，为我们刻画扇体和了解各个时期扇体的发育情况创造了条件。通过序贯高斯模拟这一方法，我们可以较为准确的预测井间的砂体厚度，为接下来的分析工作奠定了坚实的基础。
参考文献
[1] 张晓坤，随机模拟在三维地质建模中的应用：[硕士学位论文].北京：中国地质大学，2009.
[2] 杨占龙，彭立才，等.地震属性分析与岩性油气藏勘探[J].石油物探，2007，46（2）：131-136.