论文部分内容阅读
宝浪油田地理位置处于新疆维吾尔自治区巴音郭楞蒙古自治州焉耆回族自治县境内。构造上位于焉耆盆地博湖坳陷北部凹陷宝浪苏木构造带。主要产油层为三工河组。沉积相以上辫状河三角洲前缘水下分流河道和下三角洲前缘水下分流河道为主。岩性包括砂砾岩、砾岩、粗砂岩、细砂和粉砂岩等。岩石成分成熟度和结构成熟度较低。
储层属辫状河三角洲上发育起来的粗粒低渗储层类型。单层砂体厚度从1m至6m不等，常见多层叠加。河道频繁迁移加剧了储层的非均质性。不同砂体内平均孔隙度介于12％-14％之间，平均渗透率介于2.972×10-3μm2-27.3×10-3μm2之间。储层属性和流体分布复杂，砂体厚度大(Ⅲ1油组接近20m)，同时，独特的高陡背斜构造及裂缝的发育，使该地区储层成为典型的非常规储层类型：低密度、低孔隙度、低渗透率、粗粒、高粘土矿物、裂缝和复杂构造。
开发13年以来，一些与地质有关的问题仍然影响该地区储层的有效开发。例如：砂体的三维分布；流动单元的分布；裂缝的特征及其对流体渗流的影响；含水率上升的控制因素，等。
本文主要利用地质统计学、随机模拟和数值模拟等方法，首次提出并开展了针对低渗含裂缝储层进行了综合表征和模拟的方法和实现；首次采用小尺度和大尺度数值模拟相结合的方法系统评价了不同裂缝及不同流动单元组合情况下的地质模型的开发特征及参数。
具体内容主要包括以下五个方面：
第一，流动单元综合研究：在地层划分、岩心分析、岩心测试、储层特征参数和岩石物理相的基础上，采用聚类方法划分了流动单元类型。通过多元回归方法建立了流动单元的预测模型。对比砂体的空间分布，分析了流动单元在平面和剖面上的分布规律。
第二，裂缝单元综合表征：通过岩心观察测定了岩心裂缝的密度、倾角、走向等参数；并通过测井和构造曲率方向分别建立了低角度裂缝和高角度裂缝的预测模型，预测并绘制了不同裂缝的裂缝密度平面分布图。
第三，对建立随机离散裂缝分布模型的方法和参数进行了研究，并对不同裂缝密度的随机离散分布模型进行了比较。
第四，通过对预测的低角度和高角度裂缝的平面分布图进行数字化，并对裂缝的计数密度进行了面积密度的转化，最后建立了研究区的低角度和高角度裂缝随机离散分布模型。并采用多界面(多级界面主要包括油组界面、小层界面和单砂体界面)及流动单元的约束方法建立了低渗基质的地质模型。
第五，对岩心规模及小尺度规模的不同流动单元及不同裂缝产状的模型进行了数值模拟和比较。并在储层地质模型的粗化基础上，对Ⅲ1油组进行了数值模拟，在历史拟合的基础上，分析了剩余油分布和剩余可采储量的平面分布。
通过研究，取得的主要成果和认识如下：
(1)通过FZI=4μm、FZI=2.72μm、FZI=1.85μm和FZI=1.36μm.可以把本地区流动单元类型可以分为5类，流动单元A、流动单元B、流动单元C、流动单元D和流动单元E。在流动单元基础上建立的储层孔渗测井解释模型预测的储层物性与岩心分析的物性的相关系数大于0.7，能满足精细储层建模的要求。
(2)研究区裂缝类型按成因主要分为两类：构造成因和构造-沉积复合成因。复合成因的层状裂缝最发育，断层有关的裂缝次之。将近60％裂缝属于复合成因，40％属于构造成因。裂缝按倾角主要分为两种类型，低角度裂缝和高角度裂缝。高角度裂缝走向以东北、南西和东南方向为主，平均为156°。其倾角分布范围为15°-85°之间，大多数在40°-75°之间分布，平均为54.2°。背斜的不同部位倾角不同，西北部大于东南部。低角度裂缝的倾角在0°-20°之间不等，常呈“薄饼状”，一些裂缝的宽度小于5cm。
(3)三孔隙度方法能够很好的预测该地区的低角度裂缝的密度，临界值为0.3。利用SPSS软件建立了低角度裂缝密度预测的多元回归方程，利用该方法预测的裂缝密度和岩性测量的裂缝之间的相关系数为0.716。在模型建立的基础上，编写了VB代码自动计算了除部分没有补偿中子、声波和密度测井的所有井，并分析了低角度裂缝的密度平面分布。高角度裂缝的预测模型通过分析裂缝发育指数和岩心裂缝密度的关系来建立。
(4)流动单元组合下的岩心尺度数值模拟结果是不同的流动单元储层单独开采下的采收率最高，其次是流动单元A和流动单元B组合，其他组合均低于这两种组合。流动单元组合下的小尺度模拟表明，在相对注入速度相等的情况下，井距对综合采出程度的影响较小；而注入速度对综合采出程度影响较大。
(5)四种不同倾角的裂缝模拟结果显示，在含水量达到98％时，流动单元A、流动单元B和流动单元C中采收率的关系是135°＜45°＜0°＜90°＜无裂缝模型。
(6)和没有裂缝的模型比较，具有裂缝的岩心当含水饱和度达到98％时，流动单元A和流动单元B的含水率减少，流动单元C的含水率增加。由于串通裂缝的层间流动影响，流动单元A的采收率增加，流动单元B和流动单元C的减少，总的采收率降低。
(7)离散裂缝等效模拟结果显示，裂缝的等效处理网格大小对基质渗透率较低的模型模拟结果影响较大。当垂向网格大小等于模型厚度时，模拟的含水率比离散裂缝模型的结果低20％。合理的垂向网格大小和基质渗透率成正比关系。
(8)多界面和流动单元共同控制储层地质建模方法能很好的控制砂体形态和砂体属性参数分布。流动单元A、B、C、D、E和泥岩的最大主变程分别为635m、429m、520m、558m、556m和783m。孔隙度在流动单元B、C、D和E中的最大主变程分别为483m、676m、516m和628m。孔隙度在流动单元B、C、D和E中的最小次变程分别为248m、200m、264m和232m。流动单元A、B、C、D和E的渗透率控制范围分别为52.5～162.6×10-3μm2、18.3～97×10-3μm2、6.6～51.6×10-3μm2、2.5～17.8×10-3μm2、1.4～16.8×10-3μm2。最后，建立了各种储层属性分布模型。
(9)运用随机模拟方法建立了不同裂缝参数条件下的裂缝分布模型，探讨了裂缝密度参数与裂缝模型结果参数的关系。当裂缝面积密度为0.05m2/m3、0.1m2/m3、0.5m2/m3、1m2/m3、2m2/m3和10m2/m3时，裂缝模型的平均孔隙度分别为0.0000、0.0002、0.0019、0.0028、0.0087和0.0461；平均渗透率分别为1×10-3μm2、10×10-3μm2、79×10-3μm2、119×10-3μm2、369×10-3μm2、1260×10-3μm2和1941×10-3μm2；裂缝形状系数平均值分别为3984m-2、371m-2、1476m-2、3×106m-2、8×105m-2、和1×105m-2。粗化的孔隙度和渗透率值与裂缝面密度成直线关系。
(10)利用数值模拟方法对(9)的模型在基质结合不同流动单元情况下的渗流特征进行了对比。发现不同模型的采收率和累计采油量与裂缝面积密度没有直线关系。采收率在基质为流动单元A、流动单元B和流动单元C时的拐点分别为2m2/m3、2m2/m3和1m2/m3；累计产油量在裂缝面密度大于0.5 m2/m3时递增，小于0.5 m2/m时递减。注采井连线与裂缝走向平行时的采收率和累计产油量均小于相应的在注采井连线与裂缝走向垂直下的模拟结果。而含水饱和度达到98％的时间关系为：当基质为流动单元A时，垂直模型晚于平行模型，当基质为流动单元B和C时，垂直模型早于平行模型。
(11)运用随机模拟方法建立了Ⅲ1油组的高角度裂缝和低角度裂缝离散分布模型。并利用Oda方法对该离散裂缝网络模型进行了粗化，得出了裂缝孔隙度、渗透率和裂缝形状指数。模型的平均孔隙度为0.4％，I、J和K方向的平均渗透率分别为11×10-3μm2、22×10-3μm2和13×10-3μm2。I、J和K方向的平均裂缝间距分别为2.78m、4.19m和0.58m。整个模型的裂缝形状指数为18836/m2.在该模型的基础上，利用数值模拟方法对剩余油分布进行了研究。剩余油分布主要受注采井网不完善、次级断层、砂体分布和流动单元储层非均质性等因素影响控制和影响。
储层属辫状河三角洲上发育起来的粗粒低渗储层类型。单层砂体厚度从1m至6m不等，常见多层叠加。河道频繁迁移加剧了储层的非均质性。不同砂体内平均孔隙度介于12％-14％之间，平均渗透率介于2.972×10-3μm2-27.3×10-3μm2之间。储层属性和流体分布复杂，砂体厚度大(Ⅲ1油组接近20m)，同时，独特的高陡背斜构造及裂缝的发育，使该地区储层成为典型的非常规储层类型：低密度、低孔隙度、低渗透率、粗粒、高粘土矿物、裂缝和复杂构造。
开发13年以来，一些与地质有关的问题仍然影响该地区储层的有效开发。例如：砂体的三维分布；流动单元的分布；裂缝的特征及其对流体渗流的影响；含水率上升的控制因素，等。
本文主要利用地质统计学、随机模拟和数值模拟等方法，首次提出并开展了针对低渗含裂缝储层进行了综合表征和模拟的方法和实现；首次采用小尺度和大尺度数值模拟相结合的方法系统评价了不同裂缝及不同流动单元组合情况下的地质模型的开发特征及参数。
具体内容主要包括以下五个方面：
第一，流动单元综合研究：在地层划分、岩心分析、岩心测试、储层特征参数和岩石物理相的基础上，采用聚类方法划分了流动单元类型。通过多元回归方法建立了流动单元的预测模型。对比砂体的空间分布，分析了流动单元在平面和剖面上的分布规律。
第二，裂缝单元综合表征：通过岩心观察测定了岩心裂缝的密度、倾角、走向等参数；并通过测井和构造曲率方向分别建立了低角度裂缝和高角度裂缝的预测模型，预测并绘制了不同裂缝的裂缝密度平面分布图。
第三，对建立随机离散裂缝分布模型的方法和参数进行了研究，并对不同裂缝密度的随机离散分布模型进行了比较。
第四，通过对预测的低角度和高角度裂缝的平面分布图进行数字化，并对裂缝的计数密度进行了面积密度的转化，最后建立了研究区的低角度和高角度裂缝随机离散分布模型。并采用多界面(多级界面主要包括油组界面、小层界面和单砂体界面)及流动单元的约束方法建立了低渗基质的地质模型。
第五，对岩心规模及小尺度规模的不同流动单元及不同裂缝产状的模型进行了数值模拟和比较。并在储层地质模型的粗化基础上，对Ⅲ1油组进行了数值模拟，在历史拟合的基础上，分析了剩余油分布和剩余可采储量的平面分布。
通过研究，取得的主要成果和认识如下：
(1)通过FZI=4μm、FZI=2.72μm、FZI=1.85μm和FZI=1.36μm.可以把本地区流动单元类型可以分为5类，流动单元A、流动单元B、流动单元C、流动单元D和流动单元E。在流动单元基础上建立的储层孔渗测井解释模型预测的储层物性与岩心分析的物性的相关系数大于0.7，能满足精细储层建模的要求。
(2)研究区裂缝类型按成因主要分为两类：构造成因和构造-沉积复合成因。复合成因的层状裂缝最发育，断层有关的裂缝次之。将近60％裂缝属于复合成因，40％属于构造成因。裂缝按倾角主要分为两种类型，低角度裂缝和高角度裂缝。高角度裂缝走向以东北、南西和东南方向为主，平均为156°。其倾角分布范围为15°-85°之间，大多数在40°-75°之间分布，平均为54.2°。背斜的不同部位倾角不同，西北部大于东南部。低角度裂缝的倾角在0°-20°之间不等，常呈“薄饼状”，一些裂缝的宽度小于5cm。
(3)三孔隙度方法能够很好的预测该地区的低角度裂缝的密度，临界值为0.3。利用SPSS软件建立了低角度裂缝密度预测的多元回归方程，利用该方法预测的裂缝密度和岩性测量的裂缝之间的相关系数为0.716。在模型建立的基础上，编写了VB代码自动计算了除部分没有补偿中子、声波和密度测井的所有井，并分析了低角度裂缝的密度平面分布。高角度裂缝的预测模型通过分析裂缝发育指数和岩心裂缝密度的关系来建立。
(4)流动单元组合下的岩心尺度数值模拟结果是不同的流动单元储层单独开采下的采收率最高，其次是流动单元A和流动单元B组合，其他组合均低于这两种组合。流动单元组合下的小尺度模拟表明，在相对注入速度相等的情况下，井距对综合采出程度的影响较小；而注入速度对综合采出程度影响较大。
(5)四种不同倾角的裂缝模拟结果显示，在含水量达到98％时，流动单元A、流动单元B和流动单元C中采收率的关系是135°＜45°＜0°＜90°＜无裂缝模型。
(6)和没有裂缝的模型比较，具有裂缝的岩心当含水饱和度达到98％时，流动单元A和流动单元B的含水率减少，流动单元C的含水率增加。由于串通裂缝的层间流动影响，流动单元A的采收率增加，流动单元B和流动单元C的减少，总的采收率降低。
(7)离散裂缝等效模拟结果显示，裂缝的等效处理网格大小对基质渗透率较低的模型模拟结果影响较大。当垂向网格大小等于模型厚度时，模拟的含水率比离散裂缝模型的结果低20％。合理的垂向网格大小和基质渗透率成正比关系。
(8)多界面和流动单元共同控制储层地质建模方法能很好的控制砂体形态和砂体属性参数分布。流动单元A、B、C、D、E和泥岩的最大主变程分别为635m、429m、520m、558m、556m和783m。孔隙度在流动单元B、C、D和E中的最大主变程分别为483m、676m、516m和628m。孔隙度在流动单元B、C、D和E中的最小次变程分别为248m、200m、264m和232m。流动单元A、B、C、D和E的渗透率控制范围分别为52.5～162.6×10-3μm2、18.3～97×10-3μm2、6.6～51.6×10-3μm2、2.5～17.8×10-3μm2、1.4～16.8×10-3μm2。最后，建立了各种储层属性分布模型。
(9)运用随机模拟方法建立了不同裂缝参数条件下的裂缝分布模型，探讨了裂缝密度参数与裂缝模型结果参数的关系。当裂缝面积密度为0.05m2/m3、0.1m2/m3、0.5m2/m3、1m2/m3、2m2/m3和10m2/m3时，裂缝模型的平均孔隙度分别为0.0000、0.0002、0.0019、0.0028、0.0087和0.0461；平均渗透率分别为1×10-3μm2、10×10-3μm2、79×10-3μm2、119×10-3μm2、369×10-3μm2、1260×10-3μm2和1941×10-3μm2；裂缝形状系数平均值分别为3984m-2、371m-2、1476m-2、3×106m-2、8×105m-2、和1×105m-2。粗化的孔隙度和渗透率值与裂缝面密度成直线关系。
(10)利用数值模拟方法对(9)的模型在基质结合不同流动单元情况下的渗流特征进行了对比。发现不同模型的采收率和累计采油量与裂缝面积密度没有直线关系。采收率在基质为流动单元A、流动单元B和流动单元C时的拐点分别为2m2/m3、2m2/m3和1m2/m3；累计产油量在裂缝面密度大于0.5 m2/m3时递增，小于0.5 m2/m时递减。注采井连线与裂缝走向平行时的采收率和累计产油量均小于相应的在注采井连线与裂缝走向垂直下的模拟结果。而含水饱和度达到98％的时间关系为：当基质为流动单元A时，垂直模型晚于平行模型，当基质为流动单元B和C时，垂直模型早于平行模型。
(11)运用随机模拟方法建立了Ⅲ1油组的高角度裂缝和低角度裂缝离散分布模型。并利用Oda方法对该离散裂缝网络模型进行了粗化，得出了裂缝孔隙度、渗透率和裂缝形状指数。模型的平均孔隙度为0.4％，I、J和K方向的平均渗透率分别为11×10-3μm2、22×10-3μm2和13×10-3μm2。I、J和K方向的平均裂缝间距分别为2.78m、4.19m和0.58m。整个模型的裂缝形状指数为18836/m2.在该模型的基础上，利用数值模拟方法对剩余油分布进行了研究。剩余油分布主要受注采井网不完善、次级断层、砂体分布和流动单元储层非均质性等因素影响控制和影响。