[发明专利]一种用于预测缝洞型储集层资源量的方法

摘要

本发明提出了一种用于预测缝洞型储集层资源量的方法。该方法利用反映储层非均质性的地震属性资料，识别缝洞型储层储集空间分布范围；在此基础上，利用储集层有效厚度、孔隙度及含油气饱和度等关键参数平面分布图进行变差分析，并通过网格化将非连续储集层进一步划分为微小单元，计算每个单元内的资源量，最后汇集所有有效单元资源量获得缝洞型储集层总的资源量。该方法能够有效获取缝洞型储集层分布范围，保证了面积参数及其它参数的有效性和准确性，资源评价结果更为准确可信。

主权项

1.一种用于预测缝洞型储集层资源量的方法，包括以下步骤：步骤S10：获取待评价区域内的多个井位的测井资料，根据测井资料获得这些井位处对应的目的层的储层有效厚度、有效孔隙度、含油饱和度和含气饱和度；步骤S20：获取待评价区域的叠前地震数据体或叠后地震数据体，根据叠前地震数据体或叠后地震数据体获得目的层的沿层切片，对沿层切片进行均匀的网格划分，识别出储集体在沿层切片上的分布区域，在储集体分布区域内设置有井位的单元格为已知单元格，其他为未知单元格；步骤S30：将每个已知单元格中的一个井位处的储层有效厚度、有效孔隙度、含油饱和度和含气饱和度的值依次分别作为目的层在该已知单元格范围内有效厚度、有效孔隙度、含油饱和度和含气饱和度的均值；步骤S40：随机构建模拟所有未知单元格的地质属性的顺序，按顺序依次对每个未知单元格的地质属性进行模拟，模拟第一个未知单元格的地质属性时，依次进行步骤S41和步骤S42，模拟第二个未知单元格至最后一个未知单元格的地质属性时，每对一个未知单元格的地质属性进行模拟则每次进行步骤S43再进行步骤S42；步骤S41：将已知单元格中的储层有效厚度的均值作为第一实测值加入到的第一实测值集合中，将已知单元格中的有效孔隙度的均值作为第二实测值加入到第二实测值集合中，将已知单元格中的含油饱和度的均值作为第三实测值加入到第三实测值集合中，步骤S42：根据第一实测值集合中各个第一实测值、各个第一实测值所对应的单元格的位置和当前待计算的未知单元格的位置采用克里金插值法计算出该未知单元格中储层有效厚度插值以及第一估计方差，以该未知单元格的储层有效厚度插值为期望值、以第一估计方差为方差建立该未知单元格储层有效厚度的正态分布模型，根据第二实测值集合中各个第二实测值、各个第二实测值所对应的单元格的位置和当前待计算的未知单元格的位置采用克里金插值法计算出该未知单元格中有效孔隙度插值以及第二估计方差，以该未知单元格的有效孔隙度插值为期望值、以第二估计方差为方差建立该未知单元格有效孔隙度的正态分布模型，根据第三实测值集合中各个第三实测值、各个第三实测值所对应的单元格的位置和当前待计算的未知单元格的位置采用克里金插值法计算出该未知单元格中含油饱和度插值以及第三估计方差，以该未知单元格的含油饱和度插值为期望值、以第三估计方差为方差建立该未知单元格含油饱和度的正态分布模型；步骤S43：采用蒙托卡罗模拟方法从每个未知单元格储层有效厚度的正态分布模型中随机抽取一个数作为相应未知单元格的第一实测值并将已知单元格中的储层有效厚度的均值作为第一实测值，用该步骤获得的第一实测值建立新的第一实测值集合，采用蒙托卡罗模拟方法从各个未知单元格有效孔隙度的正态分布模型中随机抽取一个数作为相应未知单元格的第二实测值并将已知单元格中的有效孔隙度的均值作为第二实测值，用该步骤获得的第二实测值建立新的第二实测值集合，采用蒙托卡罗模拟方法从各个未知单元格含油饱和度的正态分布模型中随机抽取一个数作为相应未知单元格的第三实测值并将已知单元格中的含油饱和度的均值作为第三实测值，用该步骤获得的第三实测值建立新的第三实测值集合；步骤S50：根据每个已知单元格的已知单元格的面积、储层有效厚度的均值、有效孔隙度的均值、有效孔隙度的均值、含油饱和度的均值、平均地面原油密度以及原油体积系数计算出每个已知单元格对应区域内的原油资源量，根据每个未知单元格的未知单元格的面积、储层有效厚度插值、有效孔隙度插值、含油饱和度插值、平均地面原油密度以及原油体积系数计算出每个已知单元格对应区域内的原油资源量，将每个已知单元格对应区域内的原油资源量和每个未知单元格对应区域内的原油资源量累加以获得整个待评价区域的累积原油资源量；步骤S60：重复步骤S40、S50至少一千次，以获得至少一千个累积原油资源量；步骤S70：根据所有累积原油资源量，采用中点法获得待评价区域的累积原油资源量最终估值。