一种碳酸盐岩油藏古暗河型储层建模方法与流程

本发明涉及油气田开发技术领域，具体涉及一种碳酸盐岩油藏古暗河型储层建模方法。

背景技术：

碳酸盐岩油藏喀斯特古暗河储层是重要的油藏类型，古暗河系统是油藏岩溶体系的核心，其形成控制因素多样，受地下水供给、断裂、潜水面和岩性差异的影响，溶蚀与改造并存，造成了暗河形态展布、结构复杂，规律性差。平面上为单支管道、多支管道网络系统，连续或断续，充填特征多样；剖面上呈多层或单层分布，单条地下河溶洞呈近圆形、椭圆形等特征。

目前，对此类油藏的建模多以离散型孔、缝、洞为主，相对连续地下暗河建模的研究较少。实际钻遇古岩溶管道的井较少，硬数据缺乏，变差函数分析难度大，使得传统两点地质统计学方法很难表征复杂几何形态和空间结构。

相关研究主要包括：

公开号为107219553a的专利公开了一种利用地震数据基于gr分频反演的暗河充填程度预测方法，该方法利用marr小波分频技术对原始地震数据进行分频处理，得到不同频段的分频数据体并从分频数据体分别提取不同频段的分频属性，根据储层不同厚度下振幅与频率关系构建核函数，利用支持向量机进行多次学习，建立分频属性与测井gr曲线的非线性映射关系，将不同频段分频属性与测井曲线的非线性映射关系合成到一起得到gr分频反演体，根据测井解释结果确定不同充填程度暗河样本对应的gr峰值分布概率，确定gr分频反演体的暗河充填程度。

鲁新便等人(“塔河油田奥陶系油藏喀斯特古河道发育特征描述”，鲁新便等，石油实验地质，第36卷第3期，第268-274页，2014)划分了古河道类型，建立了一套刻画与识别古河道的方法技术，并实际识别6条主干河道。

张娟等人(“塔河油田西部古暗河缝洞结构特征及控制因素”，张娟等，油气地质与采收率，第25卷第4期，第33-39页，2018)对塔河油田西部古暗河进行系统刻画与识别，平面上缝洞结构分为连续发育的条带状和不连续发育的团块状2大类，纵向上自上而下可分为未充填、砂泥质充填和角砾充填3种缝洞充填类型。

高翔等人(“地下暗河溶洞系统特征描述及地址建模-以塔里木盆地哈拉哈塘地区为例”，高翔等，断块油气田，第23卷第6期，第782-787页，2016)根据录井、测井、地震和生产资料等，针对该储层建立了地下暗河溶洞系统的地质概念模型，总结出地下暗河溶洞识别的3种地震反射特征：片状反射、杂乱反射和弱反射。

公开号为103077558a的专利公开了一种碳酸盐岩缝洞型油藏大型溶洞储集体分布模型建模方法，考虑了岩溶发育模式对大型溶洞储集体展布的控制作用，加强了地质规律的约束。

公开号103077548a的专利公开了一种碳酸盐岩缝洞型油藏溶蚀孔洞储集体分布模型的建模方法，考虑了溶蚀孔洞储集体在储集体空间尺度及规模上与大型溶洞等其它类型储集体的巨大差异，单独建立溶蚀孔洞储集体模型，同时建立了距大型溶洞的距离与溶蚀孔洞储集体发育的定量概率关系，客观反映了溶蚀孔洞储集体的分布规律。

胡向阳等人(“多尺度岩溶相控碳酸盐岩缝洞型油藏储集体建模方法”，胡向阳等，石油学报，第35页第2期，第340-346页，2014)提出了多尺度岩溶相控缝洞储集体建模方法，即在古岩溶发育模式控制下，根据缝洞尺度的差异，采用两步法建模：第1步建立4个单一类型储集体离散分布模型，即利用地震识别大型溶洞和大尺度裂缝，通过确定性建模方法，建立离散大型溶洞模型和离散大尺度裂缝模型；在岩溶相控约束下，基于溶洞发育概率体和井间裂缝发育概率体，采用随机模拟多属性协同模拟方法，建立溶蚀孔洞模型和小尺度离散裂缝模型。第二步采用同位条件赋值算法，将个单一类型储集体模型融合成多尺度离散缝洞储集体三维分布模型。

侯加根等人(“缝洞型碳酸盐岩储层多类多尺度建模方法研究-以塔河油田四区奥陶系油藏为例”，侯加根等，地学前缘，第19卷第2期，第59-66页，2012)提出碳酸盐岩缝洞型储层应该按照大型洞穴、溶蚀孔洞、大尺度裂缝、小尺度裂缝的“多类多尺度建模”的基本思路。

刘钰铭等人(“塔河油田古岩溶储集体三维建模”，刘钰铭等，中国石油大学学报(自然科学版)，第36卷第2期，第34-38页，2012)提出“垂向划带、平面分区、断裂优先”的古岩溶储集体三维展布建模方法。

吕心瑞等人(“碳酸盐岩储层多尺度缝洞体分类表征-以塔河油田s80单元奥陶系油藏为例”，吕心瑞等，石油与天然气地质，第38卷第4期，第813-821页，2017)以井点钻遇的河道位置做为硬数据，结合河道描述成果，以地震分频能量预测属性作为井间暗河形状模拟约束条件，采用确定性建模方法建立分支河道储集体分布模型。

以上研究均集中于暗河描述和规律认识，建模过程均采用地震属性进行约束，存在预测放大效应，地震预测属性均反映为异常体，预测结果不连续，预测河道宽、厚多在百米以上，与野外露头暗河规模(多在20m以内)认识差别大。

此外，现有技术还存在诸多的问题，包括：暗河形态的建模均采用地球物理预测成果进行约束，地震预测异常体存在明显放大效应，预测河道宽、厚与野外露头规模认识差别大；建模过程中缺少地质成因的约束，暗河模型多为外部轮廓，缺乏对暗河结构特征的表征；暗河属性模拟均基于井点硬数据，井间进行插值模拟，缺乏可靠的约束依据；这些均造成了暗河地质模型可靠性低，难以为此类油藏开发方案设计与调整提供可靠地质依据。

技术实现要素：

为了解决上述技术问题，本发明提供了一种多元约束的古暗河地质建模方法，整合井、地震、野外露头、生产动态、油藏综合描述成果等多元信息约束模型的构建，综合反映当前资料状况下对暗河的认识，表征暗河发育形态、规模大小及结构模式。

本发明的碳酸盐岩油藏古暗河型储层建模方法，具体涉及古暗河特征参数及分布规律的确定、训练图像制作、古暗河建模方法的确定、多元约束方法的设计、动态资料优化模型、物性参数赋值与模拟，最终实现碳酸盐岩油藏古暗河储层的三维地质建模。本发明的建模方法包括以下步骤：

步骤a.优选与研究区储层具有相似性的古岩溶暗河储层及现代岩溶暗河露头，观察形态、结构、测量暗河特征参数，统计相关分布规律，确定古暗河真实特征参数的及规律分析；

步骤b.将暗河进行分类，分别制作不同类型暗河的训练图像，用于表征不同情况的地下暗河特征；

步骤c.根据研究区的暗河类型选取合适的训练图像，采用多点地质统计学算法进行模拟；

步骤d.在多点地质统计学模拟过程中，以概率体的形式整合岩溶成因规律、地球物理预测结果等其它先验认识，实现古暗河储层建模过程的多元约束；

步骤e.确定单井物性参数，基于古暗河空间分布模型，采用序贯高斯模拟方法建立地下暗河物性参数模型；

步骤f.基于动态数据优化地质模型。

进一步的，所述步骤a主要包括：

步骤a1.分析露头区储层与研究油藏的一致性，确定具有相似性的露头区；

步骤a2.确定表征暗河特征的参数，形成暗河特征参数库；

步骤a3.分析沿暗河的充填、垮塌分布，绘制相应河道展布及结构平面图；

步骤a4.基于建立的暗河特征参数库，进行相关几何特征参数分布规律的分析，用于井间河道特征模拟的约束条件。

进一步的，所述步骤b，包括：

步骤b1.从多个角度对暗河分类，以反映不同类型暗河的主要特征；

步骤b2.根据暗河分类，分别制作不同类型的训练图像；

步骤b3.校正暗河宽度；

步骤b4.利用地质知识统计信息，结合现代岩溶考察和研究区地下暗河储层研究，从成因主控因素、垂向结构和平面形态等方面考虑，构建多类不同主控要素地下河训练图像模式库。

进一步的，在所述步骤b2中，制作训练图像包括：

a.综合现代地下河形态结合塔河露头古河道几何参数，同时考虑暗河垮塌埋藏效应，制作三维训练图像；或

b.通过密井网区地震属性雕刻，结合步骤a暗河模式修正的方式构建训练图像。

进一步的，在所述步骤b3中，对古暗河宽度的校正方法采用累积概率曲线方法，针对相似露头区暗河管道任意一点宽度，在其宽度累积频率分布曲线上找到对应的累积频率；在古暗河溶洞宽度累积频率曲线上找具有相同累积频率的古岩溶暗河宽度，即为校正后古暗河核部对应的宽度；

对古暗河除核部之外的影响范围可由暗河管道物理垮塌试验或野外统计规律获得，根据地下暗河的核部和影响范围的关系校正训练图像。

进一步的，在所述步骤c中，多点地质建模算法包括迭代方法和非迭代方法，所述迭代方法是基于系统优化的思想，用于相对连续目标的模拟；非迭代方法是基于图像构型的方法，用于对目标结构等特征的模拟。

进一步的，在所述步骤c中，还包括，采用模拟目标匹配训练图像的方式模拟几何形态，采用模拟目标匹配训练图像库的方法模拟结构特征，降低随机模拟的不确定性。

进一步的，所述步骤d中的建模过程的多元约束，包括：

步骤d1.岩溶成因的约束；

步骤d2.其它地震属性的约束；

步骤d3.岩溶成因规律和地震属性信息概率体融合，融合方法如下：

其中：p{}，代表概率；c，表示出现溶洞；k，表示岩溶成因信息；s，表示地震属性信息。

进一步的，所述步骤e，主要包括：

(1)对于单井有测井曲线部分，根据测井曲线解释得到单井孔隙度、渗透率参数；

(2)对于单井无曲线或曲线杂乱部分，结合生产动态特征进行参数赋值；

(3)以单井物性参数为硬数据，利用古暗河空间分布模型作为约束条件，采用序贯高斯模拟方法建立地下暗河物性参数模型。

进一步的，所述步骤e中的(2)中，对于单井无曲线或曲线杂乱部分，结合生产动态特征进行参数赋值包括：

根据单井累产、初产信息确定单井孔隙度信息；

根据生产试井或生产动态分析反演获取单井渗透率数据。

进一步的，所述步骤f包括：

(1)基于动态判别井间连通性优化地质模型；

对于动态判别的通过暗河连通的井组，通过基于目标的算法，局部优化地下河溶洞不连续之处，确保地下暗河的连通性；或

(2)基于动态储量优化地质模型；

应用井控动态储量数据优化井控静态地质储量，以单井动态地质储量为条件数据，在缝洞单元总体积不变的前提下，优化模型中个别井控溶洞孔隙度或体积，使模型中井控储量与动态生产数据一致。

进一步的，所述步骤f中的基于动态储量优化地质模型中，地质模型的优化方法为退火模拟方法，包括：

(1)构建目标函数设定合理误差范围；

(2)基于初始模型开展数值模拟判断生产动态信息，判定目标函数值的误差大小；

(3)如果误差范围小于或等于合理误差范围，则判定为合理模型；如误差范围大于合理误差范围，则继续优化，直到误差达到设定的合理范围；

其中，-代表目标误差，当其大于5％时，继续优化，小于或等于5％时，终止优化；od代表动态方法计算的井控储量；os代表基于地质模型计算的井控储量。

与现有技术相比，本发明的碳酸盐岩油藏古暗河型储层建模方法，针对传统方法难以有效表征古暗河结构特征、缺乏可靠约束依据、模型可靠性低的问题，综合井、地震、野外露头、生产动态、油藏综合描述成果等多元信息约束模型的构建，通过露头暗河特征参数测量与规律统计解决暗河几何尺寸建模缺乏依据的问题；通过多点地质统计学模拟精细表征暗河复杂的几何形态和结构特征，通过生产动态资料约束优化建模，提高了模型与生产的符合程度，进一步提高了模型精度。

上述技术特征可以各种技术上可行的方式组合以产生新的实施方案，只要能够实现本发明的目的。

附图说明

在下文中将基于仅为非限定性的实施例并参考附图来对本发明进行更详细的描述。其中：

图1显示了本发明的一种碳酸盐岩油藏古暗河型储层建模的方法的步骤框图；

图2显示了本发明实施例1中露头暗河展布及结构特征示意图；

图3显示了本发明实施例1中单支管道型暗河示意图，箭头代表水流方向；

图4显示了本发明实施例1中管道网络型暗河示意图，箭头代表水流方向；

图5显示了本发明实施例1中断裂廊道型地下暗河示意图，箭头代表水流方向；

图6显示了本发明实施例1中典型单支管道型暗河训练图像，表征了暗河形态、几何尺寸、水流方向、充填特征、结构特征及其组合关系等；

图7显示了本发明实施例1中典型古暗河管道网络训练图像，表征了管道网络的形态、几何尺寸、水流方向、及结构等；

图8显示了本发明实施例1中古暗河核部影响范围宽度校正示意图；

图9显示了本发明实施例1中古暗河除核部之外的影响范围校正示意图；

图10显示了本发明实施例1中不同岩溶古地貌(高地、斜坡、洼地)、不同岩溶带(表层、渗流、径流)下暗河发育概率；

图11显示了本发明实施例2中r1油藏上层河建模结果；

图12显示了本发明实施例2中r1油藏下层河建模结果；

图13显示了本发明实施例2中r1油藏河道整体建模结果；

图14显示了本发明实施例2中r1油藏日产油量模拟结果与实际产量对比图；

图15显示了本发明实施例2中r1油藏日含水率模拟结果与实际含水率对比图；

图16显示了本发明实施例3中校正前河道建模结果的宽度分布图；

图17显示了本发明实施例3中校正后河道建模结果的宽度分布图；

图18显示了本发明实施例3中基于地震雕刻建立的暗河模型，为单一的河道异常体；

图19显示了本发明实施例3中基于本方法建立的暗河模型，表征了河道的结构特征；

图20显示了本发明实施例4中优化前的古暗河孔隙度模型；

图21显示了本发明实施例4中优化后古暗河孔隙度模型；

图22显示了本发明实施例4中优化前后三口井的井控储量大小与动态储量的对比。

具体实施方式

以下将结合说明书附图和具体实施例对本发明做进一步详细说明。需要说明的是，只要不构成冲突，本发明中的各个实施例以及各实施例中的各个特征可以相互结合，所形成的技术方案均在本发明的保护范围之内。

实施例1

如图1所示，本发明的碳酸盐岩油藏古暗河型储层建模方法，包括以下步骤：

步骤a.优选与研究区储层具有相似性的古岩溶暗河储层及现代岩溶暗河露头，观察形态、结构、测量暗河特征参数，统计相关分布规律，确定古暗河真实特征参数的及规律分析。具体包括：

步骤a1.确定具有相似性的露头区：分析现代岩溶暗河露头区储层与研究油藏的古气候、古地理、古水系环境、可溶岩组分、暗河形态、结构模式等方面的一致性，将具有80％以上相似性的区域作为相似露头区。

步骤a2.确定表征暗河特征的参数，表征暗河特征的参数主要包括3大类：暗河结构类型(主干河、支流河、厅堂洞等)、几何参数(长、宽、厚、长宽比、宽厚比)及内部特征(充填物、充填程度、物性)等。

统计如表1所示的典型露头暗河结构及参数列表，综合多条露头及实钻井暗河参数，形成暗河特征参数库，为河道真实几何尺寸提供依据。

表1.典型露头暗河结构及参数列表

步骤a3.分析沿暗河的充填、垮塌分布，绘制如图2所示的相应河道展布及结构平面图，河道结构特征分析及模拟提供基础。

步骤a4.基于建立的暗河特征参数库，进行相关几何特征参数分布规律的分析，用于井间河道特征模拟的约束条件。

步骤b.将暗河进行分类，分别制作不同类型暗河的训练图像，用于表征不同情况的地下暗河特征。具体包括：

步骤b1.暗河分类；

可以从多个角度对其进行分类，目的是反映不同类型暗河的主要特征，进而通过训练图像进行表达。例如，根据暗河系统几何形态特征的不同，可将暗河类型分为单支管道、管道网络和断裂廊道管道3类，

①如图3所示为单支管道型暗河示意图，管道洞主体近圆形，地下暗河呈线性展布，主洞体直径2～10m，长度1～30km，仅发育有一个主进水口和一个主排泄口的地下河管道，由管道主体和洞体周边影响带构成断面结构。

②如图4所示为管道网络型暗河示意图，暗河管道洞主体近圆形或椭圆形，地下河呈网状分布，主洞体直径2～10m，长度10～数百千米；由主管道及其周边的溶蚀影响带组成，呈枝状分布，由主管道和至少1条支管道组成，中部发育落水洞、竖井等岩溶形态。

③如图5所示为断裂廊道型暗河示意图，地下河管道高宽比大于5，平面上一般沿断裂构造带曲折延伸，垂向空间很大，断面宽1～10m，高5～10m，长度0.5～5km。

步骤b2.训练图像制作；

a.综合现代地下河形态结合塔河露头古河道几何参数，同时考虑暗河垮塌埋藏效应，制作三维训练图像。或

b.通过密井网区地震属性雕刻，结合步骤a中暗河模式修正的方式构建训练图像。

训练图像可以反映暗河的几何形态、尺寸参数、展布规律和结构特征在三维上的组合关系。

如图6所示为典型单支管道型暗河训练图像，如图7所示为典型管道网络的训练图像，分别表征了单支管道型暗河、管道网络的几何形态、尺寸、水流方向、充填特征、结构特征及其组合关系等。

步骤b3.暗河宽度校正；

由于现代河宽度及其概率分布与塔河露头古暗河溶洞宽度有较大出入，需要对现代溶洞宽度数据进行校正；如图8所示，暗河宽度校正可以采用累积概率曲线方法，针对现代溶洞任意一点宽度(如：82米)，在其宽度累积频率分布曲线上找到对应的累积频率(55％)；在古暗河溶洞宽度累积频率曲线上找具有相同累积频率(55％)的古岩溶暗河宽度(3.7米)，即为校正后古暗河核部对应的宽度，将所有暗河宽度进行校正形成暗河宽度校正表，根据宽度校正表，对宽度进行校正。校正后溶洞宽度不仅规模与古溶洞一致，数学分布也完全一致。

另外，针对上述步骤b2中a获得的训练图像，同样可以采取该方法进行校正，校正后的结果与野外相似露头区观察到的暗河特征一致。

步骤b4.训练图像库构建；

利用地质知识统计信息，结合现代岩溶考察和研究区地下暗河储层研究，从成因主控因素、垂向结构和平面形态等方面考虑，构建了多类不同主控要素地下河训练图像模式库。

其中成因主控因素分为断层主控、潜水面主控、水动力-裂缝主控等；垂向结构因素分为单层地下河溶洞、多层地下河溶洞等；平面形态因素分为单支河道、网状河道等。各类河道的展布规律、几何形态、河道结构等各不相同。

根据建模对象的特点选取合适的训练图像，可以达到精确建模的目的。如图9所示为校正后的地下暗河储集体图，其中可以很清楚的看出古暗河除核部之外的影响范围。

步骤c.根据研究区的暗河类型选取合适的训练图像，采用多点地质统计学算法进行模拟。

(1)多点地质建模算法可选择迭代方法和非迭代方法两种。迭代方法是基于系统优化的思想，用于相对连续目标的模拟；非迭代方法是基于图像构型的方法，用于对目标结构等特征的模拟。

(2)采用模拟目标匹配训练图像的方式模拟几何形态，采用模拟目标匹配训练图像库的方法模拟结构特征，降低随机模拟的不确定性。

(3)根据研究对象可以选择模拟储层中单层暗河、多层地下暗河等。

步骤d.在多点地质统计学模拟过程中，为了提高模型与实际认识的符合率，以概率体的形式整合岩溶成因规律、地球物理预测结果等其它先验认识，实现多元约束方法。其具体包括：

步骤d1.岩溶成因规律的约束；

如统计如图10所示的不同岩溶古地貌(高地、斜坡、洼地)、不同岩溶带(表层、渗流、径流)下暗河发育概率，作为模拟过程中的协同约束条件，体现暗河成因规律的控制。

步骤d2.其它地震属性的约束；

对于前述步骤中未参与建模的地震属性，可统计单井地震属性与地下河溶洞频率响应关系，将地震属性数据转化为地下暗河溶洞发育概率体，作为建模的协同约束条件，体现不同地震属性的控制作用。

步骤d3.岩溶成因规律和地震属性信息概率体之间有不同也有重叠，故地下暗河溶洞发育综合概率体由二者形成的概率体融合得出。融合方法如下：

其中：p{}，代表概率；c，表示出现溶洞；k，表示岩溶成因信息；s，表示地震属性信息。

步骤e.利用钻井、测井和生产动态反演解释数据确定单井物性参数，基于古暗河空间分布模型，采用序贯高斯模拟方法建立地下暗河物性参数模型。

(1)对于单井有测井曲线部分，根据测井曲线解释得到单井孔隙度、渗透率参数。

(2)对于单井无曲线或曲线杂乱部分，结合生产动态特征进行参数赋值，具体地：如根据单井累产、初产信息确定单井孔隙度信息；根据生产试井或生产动态分析反演获取单井渗透率数据。

(3)以单井物性参数为硬数据，利用古暗河空间分布模型作为约束条件，采用序贯高斯模拟方法建立地下暗河物性参数模型。

步骤f.基于动态数据优化地质模型。

建模过程中，古暗河与单井、地震数据以及地质规律吻合较好，但是模型与动态数据不吻合的情况，比如连通性、井控储量等。包括：

(1)基于动态判别井间连通性优化地质模型。对于动态判别的通过暗河连通的井组，通过基于目标的算法，局部优化地下河溶洞不连续之处，确保地下暗河的连通性。

(2)基于动态储量优化地质模型。应用井控动态储量数据优化井控静态地质储量，以单井动态地质储量为条件数据，在缝洞单元总体积不变的前提下，优化模型中个别井控溶洞孔隙度或体积，来使得模型中井控储量与动态生产数据一致。采用的优化方法可以采用退火模拟方法，构建目标函数

其中：-代表目标误差，当其大于5％时，继续优化，小于或等于5％时，终止优化；od代表动态方法计算的井控储量；os代表基于地质模型计算的井控储量。

按照本发明提供的碳酸盐岩油藏古暗河型储层建模方法，以r1油藏、r2油藏和r3油藏为例，进行了现场实施。

实施例2

r1油藏是一典型碳酸盐岩古暗河型储层，具有上、下两层河道，其中上层河道为典型的单支管道型暗河，下层河道为典型管道网络暗河。

依据本方法优选了相似露头，并测量如表1所示的相关参数，综合多学科资料构建训练图像，图6为r1油藏典型单支管道型暗河训练图像，表征了暗河形态、几何尺寸、水流方向、充填特征、结构特征及其组合关系等。图7为r1油藏典型古暗河管道网络训练图像，表征了管道网络的形态、几何尺寸、水流方向、及结构特征等。图10为r1油藏不同岩溶古地貌(高地、斜坡、洼地)、不同岩溶带(表层、渗流、径流)下暗河发育概率。图11为r1油藏上层河按照本方法的建模结果。图12为r1油藏下层河按照本方法的建模结果。图13为r1油藏河道整体建模结果。图14为r1油藏日产油量模拟结果与实际日产油量对比图，从图中可以看出，两者的吻合率达95％以上。图15为r1油藏日含水率模拟结果与实际含水率对比图，两者的吻合率达90％以上。将基于模型的模拟结果与实际生产状况进行对比，吻合率大幅提高，剩余油分布更加清晰，体现了模型的准确性和高效性。

实施例3

r2油藏为某油田典型碳酸盐岩古暗河储层，共钻井56口，暗河流向由北部流向南部，全长约19公里，整体发育在不整合面以下10-250米。

采用本方法构建的地质模型，第一层暗河发育在不整合面以下10-100m，北部为单支地下河，南部为网状河道；第二层地下河发育在不整合面以下100-250m，分支河道较少，覆盖范围局限。通过对10口抽稀井进行验证，与实钻井钻遇符合率由67％提高到93％。河道宽度由原地震雕刻的30-170米校正到10-50米，与野外相似露头区观察到的暗河特征更为相符，图16为校正前河道建模结果的宽度分布图，图17为校正后河道宽度分布图。模型的表征精度由单一河道异常体(图18为基于地震雕刻建立的暗河模型，为单一的河道异常体)到河道结构的建模(图19为基于本方法建立的暗河模型，表征了河道的结构特征)，新模型包含2层主河道，6条分支河道，5个溶蚀带和落水洞，表征了复杂结构，提高了模型精度。

实施例4

r3油藏为某油田典型碳酸盐岩古暗河储层，基于本方法建立了古暗河储层三维地质模型，图17为古暗河孔隙度模型，图18为基于动态储量优化后的地质模型，图19为三口典型井井控储量优化前后与动态储量的对比，可以看出优化后的井控储量与动态储量符合率显著增大，提高了模型精度。

而通过以上实例验证了本发明方法的正确性与优越性，与常规方法相比，通过本方法构建的地质模型有效表征了古暗河发育形态、规模大小及结构模式，综合反映了当前资料状况下对暗河的认识，进一步提高了模型精度，进而为油藏数值模拟和开发方案调整提供更为可靠的依据。

虽然本发明所揭露的实施方式如上，但所述的内容只是为了便于理解本发明而采用的实施方式，并非用以限定本发明。任何本发明所属技术领域内的技术人员，在不脱离本发明所揭露的精神和范围的前提下，可以在实施的形式上及细节上作任何的修改与变化，但本发明的专利保护范围，仍须以所附的权利要求书所界定的范围为准。