油藏地质模型的优化方法、装置、电子设备及存储介质与流程

1.本公开涉及石油勘探开发技术领域，特别地涉及一种油藏地质模型的优化方法、装置、电子设备及存储介质。


背景技术：

2.缝洞型油藏是一种高度离散型油藏。作为一种特殊类型的石油储层，缝洞型油藏有极强的非均质性，其储集空间类型主要为大型溶洞、溶蚀孔洞和不同尺度的裂缝。由于缝洞型油藏的强非均质性和储集空间的多样性，导致建立既符合储层发育的地质规律又满足油藏数值模拟需要的油藏地质模型仍是该类油藏开发的技术难题，其中一个突出的矛盾表现为油藏地质模型在数值模拟时与真实生产动态数据的不一致，从而影响历史拟合的效果。
3.目前应用动态数据约束建模的研究多针对储层类型简单的碎屑岩储层，对于复杂的缝洞型碳酸盐岩储层不适用，并且这些研究的目的多在于产量压力数据的历史拟合。所以在缝洞型碳酸盐岩储层的实际生产中，常出现油藏地质模型中个别井的井控储量与根据动态资料计算的动态储量相差较大的情况。因此，目前井控储量的模拟还有所欠缺，使用的油藏地质模型精度和准确性较差，导致数值模拟效果较差。


技术实现要素：

4.针对上述问题，本公开提供一种油藏地质模型的优化方法、装置、电子设备及存储介质，解决了现有技术中油藏地质模型的井控储量与动态储量相差较大，以及油藏地质模型精度和准确性较差导致数值模拟效果较差的问题。
5.第一方面，本公开提供了一种油藏地质模型的优化方法，所述方法包括：
6.步骤s1：确定油藏地质模型的井控边界；
7.步骤s2：对所述井控边界范围内的井控储量进行计算，并确定所述井控边界范围内的井控储量的影响因素；
8.步骤s3：根据所述井控边界范围内的井控储量的影响因素，选取相应的优化方法对所述油藏地质模型进行优化。
9.根据本公开的实施例，可选的，上述油藏地质模型的优化方法中，所述步骤s1包括以下步骤：
10.步骤s11：确定所述油藏地质模型的目标井附近第一预设范围内的邻近井，并依次连接所述目标井与各个所述邻近井连线的中点，以形成一个闭合的多边形并将其作为第一边界；
11.步骤s12：确定所述第一边界范围内的所有断层，通过将其中与所述目标井的距离大于预设距离的断层作为边界线的一部分来修正所述第一边界，从而得到第二边界；
12.步骤s13：判断所述目标井是否为酸化压裂的井，当所述目标井为酸化压裂的井且其压裂范围超过所述第二边界时，修正所述第二边界，以得到第三边界并使所述压裂范围
位于所述第三边界范围内；
13.步骤s14：确定所述第三边界与所述油藏地质模型内的孤立溶洞的位置关系，当所述第三边界穿过所述孤立溶洞且所述孤立溶洞属于油藏储集体时，修正所述第三边界，以得到第四边界并使所述孤立溶洞位于所述第四边界范围内；
14.步骤s15：根据所述目标井在生产过程中的压力数据，确定所述目标井的压力传播范围边界，当所述压力传播范围边界位于所述第四边界范围内时，将所述压力传播范围边界作为边界线的一部分来修正所述第四边界，从而得到最终的井控边界。
15.根据本公开的实施例，可选的，上述油藏地质模型的优化方法中，所述油藏地质模型包括三维正交网格体系，所述对所述井控边界范围内的井控储量进行计算的步骤，包括以下步骤：
16.步骤s21：根据所述三维正交网格体系，对所述井控边界范围内的所有溶洞进行网格划分，以得到若干方格；
17.步骤s22：对所有方格内的含油体积进行逐个计算并累加，以得到所述井控边界范围内的含油体积；
18.步骤s23：结合地面原油密度和地面原油体积系数将所述井控边界范围内的含油体积换算成所述井控边界范围内的井控储量。
19.根据本公开的实施例，可选的，上述油藏地质模型的优化方法中，根据下式计算所述井控边界范围内的井控储量：
[0020][0021]
其中，n0为所述井控边界范围内的井控储量；
[0022]
j为所述井控边界范围内的总方格数；
[0023]
v
b
(k)为第k个方格的体积；
[0024]
n
tg
(k)为第k个方格的净毛比；
[0025]
φ(k)为第k个方格的孔隙度；
[0026]
s0(k)为第k个方格的含油饱和度；
[0027]
b0为地面原油体积系数；
[0028]
ρ0为地面原油密度。
[0029]
根据本公开的实施例，可选的，上述油藏地质模型的优化方法中，所述步骤s3中，
[0030]
当所述井控边界范围内的井控储量的影响因素为体积参数时，所述优化方法为体积优化法；
[0031]
当所述井控边界范围内的井控储量的影响因素为属性参数时，所述优化方法为属性优化法；
[0032]
当所述井控边界范围内的井控储量的影响因素同时包括体积参数和属性参数时，所述优化方法为体积优化法和属性优化法相结合的方法。
[0033]
根据本公开的实施例，可选的，上述油藏地质模型的优化方法中，所述步骤s3中，当所述井控边界范围内的井控储量的影响因素为体积参数时，选取体积优化法对所述油藏地质模型进行优化，包括以下步骤：
[0034]
步骤s311：以位于所述井控边界范围内且储量受体积参数影响的溶洞为待优化的第一溶洞；
[0035]
步骤s312：对各个第一溶洞的体积分别进行扰动，在每次扰动后重新计算所述井控边界范围内的井控储量并计算其与井控储量目标值的差值，以得到第一差值，并将所述第一差值的绝对值与第一预设结束阈值进行比较；
[0036]
步骤s313：当至少一个第一溶洞的体积在经历扰动后的第一差值的绝对值小于所述第一预设结束阈值时，结束扰动并得到优化后的所述油藏地质模型；
[0037]
步骤s314：当所有的第一溶洞的体积在经历扰动后的第一差值的绝对值都大于或等于所述第一预设结束阈值时，返回步骤s312，以再一次对各个第一溶洞的体积分别进行扰动。
[0038]
根据本公开的实施例，可选的，上述油藏地质模型的优化方法中，所述步骤s312中，对各个第一溶洞的体积分别进行扰动，包括以下步骤：
[0039]
将各个第一溶洞的体积均匀地向外扩大或向内缩小一个方格大小的一层且保持第一溶洞的形状不变；
[0040]
确认第一溶洞的属性参数；其中，对于体积扩大的第一溶洞，扩大部分的方格的属性参数为其周围方格的相应属性参数的平均值，以使扩大后的第一溶洞的属性参数在空间上是逐渐过渡的；对于体积缩小的第一溶洞，各方格的属性参数不变。
[0041]
根据本公开的实施例，可选的，上述油藏地质模型的优化方法中，所述步骤s3中，当所述井控边界范围内的井控储量的影响因素为属性参数时，选取属性优化法对所述油藏地质模型进行优化，包括以下步骤：
[0042]
步骤s321：以位于所述井控边界范围内且储量受属性参数影响的溶洞为待优化的第二溶洞；
[0043]
步骤s322：对所有第二溶洞的属性参数同时进行扰动，在每次扰动后重新计算所述井控边界范围内的井控储量，并分析其与所述井控储量目标值的差值，以得到第二差值；
[0044]
步骤s323：根据当前一次扰动后的第二差值的绝对值与前一次扰动后的第二差值的绝对值的比较结果选择接受或拒绝当前一次扰动；其中，当当前一次扰动后的第二差值的绝对值小于前一次扰动后的第二差值的绝对值时，接受当前一次扰动，并将当前一次扰动过程中所使用的扰动因子和对应的第二差值存储为最新的扰动因子和最新的第二差值；
[0045]
步骤s324：确认当前扰动次数，并将其与预设扰动次数进行比较；
[0046]
步骤s325：当所述当前扰动次数大于或等于所述预设扰动次数时，结束扰动并将当前存储的最新的扰动因子对应的所述油藏地质模型作为优化后的所述油藏地质模型；
[0047]
步骤s326：当所述当前扰动次数小于所述预设扰动次数时，比较当前存储的最新的第二差值的绝对值与第二预设结束阈值的大小：
[0048]
当当前存储的最新的第二差值的绝对值小于所述第二预设结束阈值时，结束扰动并将当前存储的最新的扰动因子对应的所述油藏地质模型作为优化后的所述油藏地质模型；
[0049]
当当前存储的最新的第二差值的绝对值大于或等于所述第二预设结束阈值时，返回步骤s322，以对所有第二溶洞的属性参数进行下一次扰动。
[0050]
根据本公开的实施例，可选的，上述油藏地质模型的优化方法中，所述步骤s322中，当第二溶洞为放空的和未钻遇的溶洞时，对所有第二溶洞的属性参数同时进行扰动，包括以下步骤：
[0051]
将所有第二溶洞的所有方格的属性参数同时乘以对应的扰动因子，以成比例地改变所有第二溶洞的属性参数；
[0052]
其中，扰动因子为定值且为所述井控储量目标值与优化前的井控储量比值的0.8倍至1.2倍。
[0053]
根据本公开的实施例，可选的，上述油藏地质模型的优化方法中，所述步骤s322中，当第二溶洞为钻遇的非放空溶洞时，对所有第二溶洞的属性参数同时进行扰动，包括以下步骤：
[0054]
将所有第二溶洞的所有方格的属性参数同时乘以对应的扰动因子，以成比例地改变第二溶洞的属性参数；
[0055]
其中，位于所述目标井位置处的第二溶洞的扰动因子为1，位于所述目标井附近第二预设范围内的第二溶洞的扰动因子随第二溶洞与所述目标井的距离的增大而减小，位于所述目标井附近所述第二预设范围外的第二溶洞的扰动因子为定值且小于1。
[0056]
根据本公开的实施例，可选的，上述油藏地质模型的优化方法中，所述步骤s3中，当所述井控边界范围内的井控储量的影响因素包括体积参数和属性参数时，选取体积优化法和属性优化法相结合的方法对所述油藏地质模型进行优化，包括以下步骤：
[0057]
步骤s331：以位于所述井控边界范围内且储量受体积参数和属性参数共同影响的溶洞为待优化的第三溶洞；
[0058]
步骤s332：通过体积优化法对第三溶洞进行优化并计算体积优化前后所述井控边界范围内的井控储量变化值，以得到井控储量第一变化值；
[0059]
步骤s333：通过属性优化法对第三溶洞进行优化并计算属性优化前后所述井控边界范围内的井控储量变化值，以得到井控储量第二变化值；其中，所述井控储量第二变化值等于井控储量变化值的期望值与所述井控储量第一变化值的差值，所述井控储量变化值的期望值为所述井控储量目标值与体积优化前的井控储量的差值。
[0060]
第二方面，本公开提供了一种油藏地质模型的优化装置，所述装置包括：
[0061]
边界确定模块，用于确定油藏地质模型的井控边界；
[0062]
储量计算和影响因素确定模块，用于对所述井控边界范围内的井控储量进行计算，并确定所述井控边界范围内的井控储量的影响因素；
[0063]
模型优化模块，用于根据所述井控边界范围内的井控储量的影响因素，选取相应的优化方法对所述油藏地质模型进行优化。
[0064]
第三方面，本公开提供一种电子设备，包括存储器和处理器，所述存储器上存储有计算机程序，该计算机程序被所述处理器执行时，执行如第一方面中任意一项所述的油藏地质模型的优化方法。
[0065]
第四方面，本公开提供一种存储介质，该存储介质存储的计算机程序，可被一个或多个处理器执行，可用来实现如第一方面中任一项所述的油藏地质模型的优化方法。
[0066]
与现有技术相比，上述方案中的一个或多个实施例可以具有如下优点或有益效果：
[0067]
本公开提供的一种油藏地质模型的优化方法、装置、电子设备及存储介质，所述方法包括确定油藏地质模型的井控边界；对所述井控边界范围内的井控储量进行计算，并确定所述井控边界范围内的井控储量的影响因素；根据所述井控边界范围内的井控储量的影
响因素，选取相应的优化方法对所述油藏地质模型进行优化。该方法对油藏地质模型进行局部的优化，使优化后的油藏地质模型的井控储量更接近于根据动态资料计算的目标井控储量，提高了油藏地质模型精度和准确性，从而优化了模拟效果。而且该方法在保持原有储集体形态不变的情况下对不同类型的储集体的体积参数和/或属性参数进行优化，灵活度高，操作性强。
附图说明
[0068]
在下文中将基于实施例并参考附图来对本公开进行更详细的描述：
[0069]
图1为本公开实施例提供的一种油藏地质模型的优化方法的流程示意图；
[0070]
图2为本公开实施例提供的井控边界示意图；
[0071]
图3为本公开实施例提供的某一目标井ta1的pda双对数曲线示意图；
[0072]
图4为本公开实施例提供的某一目标井ta1的压力传播范围边界示意图；
[0073]
图5为本公开实施例提供的溶洞体积扰动示意图；
[0074]
图6为本公开实施例提供的选取体积优化法对油藏地质模型进行优化的流程示意图；
[0075]
图7为本公开实施例提供的选取属性优化法对油藏地质模型进行优化的流程示意图；
[0076]
图8为本公开实施例提供的放空的和未钻遇的溶洞的孔隙度扰动示意图；
[0077]
图9为本公开实施例提供的钻遇的非放空溶洞的孔隙度扰动示意图；
[0078]
图10为本公开实施例提供基于塔河地区某单元油藏地质模型示意图；
[0079]
图11为本公开实施例提供的基于塔河地区某单元油藏地质模型的井控边界示意图；
[0080]
图12为本公开实施例提供的基于塔河地区某单元油藏地质模型优化前井控边界范围内的储集体平面分布示意图；
[0081]
图13为本公开实施例提供的基于塔河地区某单元油藏地质模型优化前井控边界范围内的储集体空间分布示意图；
[0082]
图14为本公开实施例提供的基于塔河地区某单元油藏地质模型优化后井控边界范围内的储集体平面分布示意图；
[0083]
图15为本公开实施例提供的基于塔河地区某单元油藏地质模型优化后井控边界范围内的储集体空间分布示意图；
[0084]
图16为本公开实施例提供的一种油藏地质模型的优化装置的结构示意图；
[0085]
在附图中，相同的部件使用相同的附图标记，附图并未按照实际的比例绘制。
具体实施方式
[0086]
以下将结合附图及实施例来详细说明本公开的实施方式，借此对本公开如何应用技术手段来解决技术问题，并达到相应技术效果的实现过程能充分理解并据以实施。本公开实施例以及实施例中的各个特征，在不相冲突前提下可以相互结合，所形成的技术方案均在本公开的保护范围之内。
[0087]
同时，在以下说明中，出于解释的目的而阐述了许多具体细节，以提供对本发明实
施例的彻底理解。然而，对本领域的技术人员来说显而易见的是，本发明可以不用这里的具体细节或者所描述的特定方式来实施。
[0088]
实施例一
[0089]
图1为本公开实施例提供的一种油藏地质模型的优化方法的流程示意图，请参阅图1，本实施例提供一种油藏地质模型的优化方法，包括：
[0090]
步骤s1：确定油藏地质模型的井控边界。
[0091]
具体的，步骤s1包括以下步骤：
[0092]
步骤s11：确定所述油藏地质模型的目标井附近第一预设范围内的邻近井，并依次连接所述目标井与各个所述邻近井连线的中点，以形成一个闭合的多边形并将其作为第一边界；
[0093]
步骤s12：确定所述第一边界范围内的所有断层，通过将其中与所述目标井的距离大于预设距离的断层作为边界线的一部分来修正所述第一边界，从而得到第二边界；
[0094]
步骤s13：判断所述目标井是否为酸化压裂的井，当所述目标井为酸化压裂的井且其压裂范围超过所述第二边界时，修正所述第二边界，以得到第三边界并使所述压裂范围位于所述第三边界范围内；
[0095]
步骤s14：确定所述第三边界与所述油藏地质模型内的孤立溶洞的位置关系，当所述第三边界穿过所述孤立溶洞且所述孤立溶洞属于油藏储集体时，修正所述第三边界，以得到第四边界并使所述孤立溶洞位于所述第四边界范围内；
[0096]
步骤s15：根据所述目标井在生产过程中的压力数据，确定所述目标井的压力传播范围边界，当所述压力传播范围边界位于所述第四边界范围内时，将所述压力传播范围边界作为边界线的一部分来修正所述第四边界，从而得到最终的井控边界。
[0097]
需要说明的是，如图2所示，当所述压力传播范围边界位于外时，不考虑所述目标井的压力传播的问题，将所述第四边界作为最终的井控边界。而如果所述压力传播范围边界只有部分位于所述第四边界范围内，则只将位于所述第四边界范围内的那部分压力传播范围边界作为边界线的一部分来修正所述第四边界，以得到最终的井控边界。缝洞型碳酸盐岩储层非均质性严重，储集体的规模与展布方向的不确定性较大，在确定井控边界时，不仅要考虑单井(所述目标井)与所述邻近井的距离和连通关系，还要考虑所述断层的遮挡和所述孤立溶洞的连通等因素。
[0098]
一般情况下，第一预设范围通常取1.5倍平均井距。对于酸化压裂的井，认为压裂范围内的裂缝储集体都是连通的，所以最终的井控边界需包括压裂范围内的裂缝。
[0099]
步骤s15之前得到的都是静态边界，步骤s15得到的是动态边界。
[0100]
随着油井生产的进行，压力波不断往外传播，压力波的传播是有一定范围的，只有该压力传播范围内的油藏储量属于所述目标井的井控储量，该压力传播范围外的油藏都不属于所述目标井的井控储量。压力传播范围边界为封闭边界，本实施例中，以方形的压力传播范围边界为例，根据所述目标井在生产过程中的压力数据，可以得到压力传播到各边界时的特征及其传播时间，根据这一时间段可以确定所述目标井与压力传播范围某一条边界的距离，进而依次确定封闭边界的4条边界与所述目标井的距离。具体的压力传播范围边界的确定过程如下：通过所述目标井的pda(production data analysis，生产数据分析)双对数曲线进行探边测试方法，确定拟稳态条件下压力传播范围，当压力传播到传播范围边界
时，pda曲线上会出现边界特征，然后根据所述边界特征的时间段对应的与所述目标井的距离确定传播范围边界的位置(即与所述目标井的距离)。
[0101]
以某一目标井ta1井为例，如图3所示，首先通过的pda双对数曲线进行探边测试方法，pda双对数曲线上依次出现了压力传播范围的4条边界(东、西、南、北四个方向)的特征(特征1、2、3、4)，分别将这四个特征的时间段对应的距离作为压力传播范围的4条边界与目标井的距离，如图4所示，从而得到压力传播范围边界的示意图。然后结合静态边界(第四边界)，得到最终的井控边界。
[0102]
步骤s2：对所述井控边界范围内的井控储量进行计算，并确定所述井控边界范围内的井控储量的影响因素。
[0103]
所述油藏地质模型是基于三维正交网格体系建立的。
[0104]
具体的，所述步骤s2中，所述对所述井控边界范围内的井控储量进行计算的步骤，包括以下步骤：
[0105]
步骤s21：根据所述三维正交网格体系，对所述井控边界范围内的所有溶洞进行网格划分，以得到若干方格；
[0106]
步骤s22：对所有方格内的含油体积进行逐个计算并累加，以得到所述井控边界范围内的含油体积；
[0107]
步骤s23：结合地面原油密度和地面原油体积系数将所述井控边界范围内的含油体积换算成所述井控边界范围内的井控储量。
[0108]
需要说明的是，因为所述油藏地质模型是基于三维正交网格体系建立的，所以所述油藏地质模型的溶洞被自动网格化，得到若干方格。
[0109]
具体，所述井控边界范围内的井控储量的计算公式如下：
[0110][0111]
其中，n0为所述井控边界范围内的井控储量，单位为t；
[0112]
j为所述井控边界范围内的总方格数；
[0113]
v
b
(k)为第k个方格的体积，单位为m3；
[0114]
n
tg
(k)为第k个方格的净毛比；
[0115]
φ(k)为第k个方格的孔隙度；
[0116]
s0(k)为第k个方格的含油饱和度；
[0117]
b0为地面原油体积系数；
[0118]
ρ0为地面原油密度，单位为103kg/m3。
[0119]
步骤s3：根据所述井控边界范围内的井控储量的影响因素，选取相应的优化方法对所述油藏地质模型进行优化。
[0120]
根据步骤s2中计算井控储量的公式，可以发现该公式是对所述井控边界范围内所有溶洞的所有方格的储量的求和，从公式来看，所述井控边界范围内的井控储量的影响因素可能为方格的体积v
b
和总方格数j等体积参数，也可能为方格的净毛比n
tg
、孔隙度φ和含油饱和度s0等属性参数。各种参数的不确定性的来源不同，主要来源于储集体体积建模过程的不确定性以及储集体属性建模过程中的不确定性两个方面。其中，总方格数j和方格体积v
b
这两个参数的不确定性来源于储集体(溶洞)体积建模过程中的不确定性，净毛比n
tg
、孔隙度φ和含油饱和度s0这三个参数的不确定性来源于储集体属性建模过程中的不确定
性。
[0121]
当井控储量与动态储量出现矛盾时，应当结合实际生产具体分析矛盾的原因。例如，对于放空溶洞，其属性参数的赋值是使用动态类比法半定量确定的，当用此属性参数计算的井控储量和动态储量产生矛盾时，孔隙度、含油饱和度等属性参数的不确定性影响最大，应当选取属性优化方法对此类溶洞进行优化。例如，对于位于目标井处的溶洞，此类溶洞的属性参数是相对准确的，溶洞的体积(总方格数j和方格体积v
b
)的不确定性更大，应当选取体积优化方法对此类溶洞进行优化。井点(目标井)一定距离外的属性参数的不确定较大，此时，应当选取属性优化方法对此类溶洞进行优化。
[0122]
较为准确的属性模型是得到相对准确井控储量的前提，但在实际建模过程中诸多因素都会增加属性模型的不确定性，例如，有些井虽钻遇溶洞，由于钻头放空，导致没有相应的孔隙度测井解释，该溶洞孔隙度的赋值多是通过类比确定的，具有一定的不确定性。当溶洞的属性参数成为影响井控储量计算的主要因素时，需要调整所述井控边界范围内的属性模型以达到对井控储量优化的目的，即应当选取属性优化方法对此类溶洞进行优化。在缝洞型油藏中，绝大部分的储集空间为大型溶洞，因此，通过体积优化法和/或属性优化法对储集体的体积和属性进行优化时，优化的目标是大型溶洞。
[0123]
具体的，所述优化方法的选取原则如下：当所述井控边界范围内的井控储量的影响因素为体积参数时，所述优化方法为体积优化法；当所述井控边界范围内的井控储量的影响因素为属性参数时，所述优化方法为属性优化法；当所述井控边界范围内的井控储量的影响因素同时包括体积参数和属性参数时，所述优化方法为体积优化法和属性优化法相结合的方法。
[0124]
下面对这三种优化方式进行说明。
[0125]
第一种优化方式：所述步骤s3中，当所述井控边界范围内的井控储量的影响因素为属性参数时，选取属性优化法对所述油藏地质模型进行优化，包括以下步骤：
[0126]
步骤s311：以位于所述井控边界范围内且储量受体积参数影响的溶洞为待优化的第一溶洞；
[0127]
步骤s312：对各个第一溶洞的体积分别进行扰动，在每次扰动后重新计算所述井控边界范围内的井控储量并计算其与井控储量目标值的差值，以得到第一差值，并将所述第一差值的绝对值与第一预设结束阈值进行比较；
[0128]
步骤s313：当至少一个第一溶洞的体积在经历扰动后的第一差值的绝对值小于所述第一预设结束阈值时，结束扰动并得到优化后的所述油藏地质模型；
[0129]
步骤s314：当所有的第一溶洞的体积在经历扰动后的第一差值的绝对值都大于或等于所述第一预设结束阈值时，返回步骤s312，以再一次对各个第一溶洞的体积分别进行扰动。
[0130]
其中，对各个第一溶洞的体积分别进行扰动，包括以下步骤：将各个第一溶洞的体积均匀地向外扩大或向内缩小一个方格大小的一层且保持第一溶洞的形状不变；确认第一溶洞的属性参数；其中，对于体积扩大的第一溶洞，扩大部分的方格的属性参数为其周围方格的相应属性参数的平均值，以使扩大后的第一溶洞的属性参数在空间上是逐渐过渡的；对于体积缩小的第一溶洞，各方格的属性参数不变。
[0131]
具体的，扩大第一溶洞方法为，在优化前的油藏地质模型中，设所述井控边界范围
内代表某个第一溶洞的所有方格的集合为s，则对其中任意一个方格u(i，j，k)，有：
[0132]
u(i，j，k)∈s{所述井控边界范围内代表某个第一溶洞的方格}
[0133]
对所有u(i，j，k)∈s，令u(i
±
1，j
±
1，k
±
1)∈s，如图5所示，即将第一溶洞的体积均匀地向外扩大一个方格大小的一层且保持第一溶洞的形状不变。而且扩大部分的净毛比、孔隙度、含油饱和度等属性参数为其周围方格相应属性参数的平均值，以确保扩大后第一溶洞的属性参数在空间上是逐渐过渡的。所以，对于属性均匀的第一溶洞，扩大部分的属性参数与其它部分一致。
[0134]
而对于体积缩小的第一溶洞，只是第一溶洞的方格数发生了变化，各方格的属性参数并未发生变化。
[0135]
需要说明的是，体积扰动方法中可以对多个第一溶洞一起进行扰动，也可以逐个进行。
[0136]
如图6所示，以扩大第一溶洞的体积为例，本实施例中提供一种逐个扰动的方法，包括：
①
逐个将每个第一溶洞的体积进行扰动，在每个第一溶洞扰动后重新计算所述井控边界范围内的井控储量并计算其与井控储量目标值的差值，以得到第一差值并将其绝对值与所述第一预设结束阈值进行比较；
②
当当前的第一溶洞扰动后的第一差值的绝对值小于所述第一预设结束阈值时，结束扰动并得到优化后的所述油藏地质模型；
③
当所有的第一溶洞扰动后的第一差值的绝对值都大于或等于所述第一预设结束阈值时，返回步骤
①
，以再一次对每个第一溶洞的体积进行扰动。
[0137]
第二种优化方式：如图7所示，所述步骤s3中，当所述井控边界范围内的井控储量的影响因素为属性参数时，选取属性优化法对所述油藏地质模型进行优化，包括以下步骤：
[0138]
步骤s321：以位于所述井控边界范围内且储量受属性参数影响的溶洞为待优化的第二溶洞；
[0139]
步骤s322：对所有第二溶洞的属性参数同时进行扰动，在每次扰动后重新计算所述井控边界范围内的井控储量，并分析其与所述井控储量目标值的差值，以得到第二差值；
[0140]
步骤s323：根据当前一次扰动后的第二差值的绝对值与前一次扰动后的第二差值的绝对值的比较结果选择接受或拒绝当前一次扰动；其中，当当前一次扰动后的第二差值的绝对值小于前一次扰动后的第二差值的绝对值时，接受当前一次扰动，并将当前一次扰动过程中所使用的扰动因子和对应的第二差值存储为最新的扰动因子和最新的第二差值；
[0141]
步骤s324：确认当前扰动次数，并将其与预设扰动次数进行比较；
[0142]
步骤s325：当所述当前扰动次数大于或等于所述预设扰动次数时，结束扰动并将当前存储的最新的扰动因子对应的所述油藏地质模型作为优化后的所述油藏地质模型；
[0143]
步骤s326：当所述当前扰动次数小于所述预设扰动次数时，比较当前存储的最新的第二差值的绝对值与第二预设结束阈值的大小：
[0144]
当当前存储的最新的第二差值的绝对值小于所述第二预设结束阈值时，结束扰动并将当前存储的最新的扰动因子对应的所述油藏地质模型作为优化后的所述油藏地质模型；
[0145]
当当前存储的最新的第二差值的绝对值大于或等于所述第二预设结束阈值时，返回步骤s322，以对所有第二溶洞的属性参数进行下一次扰动。
[0146]
其中，上述属性优化方法对放空溶洞和非放空溶洞的优化过程又有所区别。
[0147]
对于所述井控边界范围内的放空的和未钻遇的溶洞储集体，对此类溶洞的属性参数进行扰动时，每次扰动将第二溶洞的所有方格的属性参数同时乘以对应的扰动因子c，以成比例地改变所有第二溶洞的属性参数；其中，扰动因子c为定值且为所述井控储量目标值与优化前的井控储量比值(原始井控储量)的0.8倍至1.2倍。以属性参数中的孔隙度的扰动为例，钻遇的非放空溶洞的孔隙度扰动示意图如图8所示。对孔隙度优化过程来说，称对溶洞孔隙度模型的一次变更为一次扰动。
[0148]
对于所述井控边界范围内钻遇的非放空溶洞储集体，对此类溶洞的属性参数进行扰动时，每次扰动都将第二溶洞的所有方格的属性参数同时乘以对应的扰动因子c，以成比例地改变所有第二溶洞的属性参数。由于所述目标井(井点)处有相对准确的属性参数解释数据，离井点的距离越近，属性参数的不确定性越低，所以井点(目标井)处的第二溶洞的属性参数保持不变，扰动因子c为1，而后向外第二预设范围内的第二溶洞的扰动因子c随第二溶洞与井点的距离的增大而减小(其减小趋势呈正弦曲线)，再往外的第二溶洞的扰动因子c变为常数即为定值，且可以得出该常数小于1。以属性参数中的孔隙度的扰动为例，定义孔隙度扰动模型，以及孔隙度扰动示意图如图9所示。
[0149]
第三种优化方式：当所述井控边界范围内的井控储量的影响因素包括体积参数和属性参数时，选取体积优化法和属性优化法相结合的方法对所述油藏地质模型进行优化，包括以下步骤：
[0150]
步骤s331：以位于所述井控边界范围内且储量受体积参数和属性参数共同影响的溶洞为待优化的第三溶洞；
[0151]
步骤s332：通过体积优化法对第三溶洞进行优化并计算体积优化前后所述井控边界范围内的井控储量变化值，以得到井控储量第一变化值；
[0152]
步骤s333：通过属性优化法对第三溶洞进行优化并计算属性优化前后所述井控边界范围内的井控储量变化值，以得到井控储量第二变化值；其中，所述井控储量第二变化值等于井控储量变化值的期望值与所述井控储量第一变化值的差值，所述井控储量变化值的期望值为所述井控储量目标值与体积优化前的井控储量的差值；
[0153]
步骤s334：计算所述井控储量第一变化值与所述井控储量第二变化值的比值以得到扰动系数α。
[0154]
也就是说，
△
b＝
△
b
v
+
△
b
p
，α＝(
△
b
v
)/(
△
b
p
)，其中，扰动系数α越大，则表示第三溶洞的储量受体积参数不确定性的影响越大，扰动系数α越小，则表示第三溶洞的储量受属性参数不确定性的影响越大。得到扰动系数α之后，可以根据扰动系数α的大小选择只通过体积优化法或属性优化法重新对第三溶洞进行优化，这样灵活度更高，操作性更强。此处，体积优化法和属性优化法的具体步骤和上述提供的体积优化法和属性优化法的具体步骤相同，此处不再赘述。
[0155]
需要说明的是，井控储量是油田开发技术方案制定和调整的重要依据之一，其确定方法主要有静态法和动态法。静态法是利用地质模型确定单井控制储量即本公开实施例中计算井控储量所使用的方法。动态法是依据地质储量同动态数据的关系，利用动态数据反推地质储量，方法主要包括物质平衡法、递减曲线法、水驱曲线法等。
[0156]
而本公开中就是根据动态资料通过动态法计算得到的动态储量为目标井控储量，对地质模型进行优化。使优化后的油藏地质模型的井控储量更接近于根据动态资料计算的
目标井控储量，提高了油藏地质模型精度和准确性，从而优化了模拟效果。
[0157]
具体的，所述井控储量目标值根据如下步骤计算得到：
[0158]
根据所述目标井的状态参数，确定所述目标井的孔隙介质模型；
[0159]
获取所述目标井的属性参数，基于所述目标井的孔隙介质模型，结合试井数据和/或生产动态数据，确定所述目标井的单井控制面积；
[0160]
根据所述目标井的属性参数和单井控制面积，确定所述目标井的单井控制储量。
[0161]
其中，孔隙介质模型的确定步骤具体包括：根据所述目标井的状态参数，对所述目标井进行分类，得到单井分类结果；根据所述单井分类结果，确定所述目标井的地质特征；根据所述地质特征，确定所述目标井的孔隙介质模型。孔隙介质模型为三孔双渗模型、三孔单渗模型、双孔双渗模型、双孔单渗模型或单孔介质模型中的一种。而其中，根据所述目标井的状态参数，对所述目标井进行分类聚类分析的方式，所利用的状态参数包括：钻井过程的放空情况、钻井过程中的漏失情况、过井地震剖面特征、酸压过程泵压状态、测井解释裂缝信息、储集体连通情况、岩心储集体发育情况、初期日产油量、累积产油量、油压变化幅度和产量递减情况等参数。
[0162]
具体地，对于某口井来说，若钻井过程中放空超0.5米，那么放空情况所对应的权值为1；若放空不超0.5米，那么放空情况所对应的权值为0；若无放空现象，那么放空情况所对应的权值为-1。
[0163]
若该口井在钻井过程中的漏失泥浆超千方，那么漏失情况所对应的权值为1；若漏失现象不超千方，那么漏失情况所对应的0；若无漏失现象，那么漏失情况所对应的权值为-1。
[0164]
若该口井的过井地震剖面特征呈串珠状反射特征，那么该状态参数(即过井地震剖面)所对应的权值1；若过井地震剖面特征呈其它类型反射特征，那么该状态参数所对应的权值为0；若过井地震剖面特征呈无反射特征，那么该状态参数所对应的权值为-1。
[0165]
若该口井在酸压过程中，泵压降低速度超5％，那么泵压所对应的权值为1；若泵压降低速度不超5％，那么泵压所对应的权值为0；若泵压为其它情况，那么泵压所对应的权值为-1。
[0166]
若该口井的测井解释裂缝信息中，存在10条以上裂缝发育，那么裂缝所对应的权值为1；若存在1～10条裂缝发育，那么裂缝所对应的权值为0；若不存在裂缝发育，那么裂缝所对应的权值为-1。
[0167]
对于该口井的储集体连通情况，则可以通过该口井的地质模型来判别。具体地，若钻遇区溶洞与1条以上大裂缝连通，那么则说明连通良好，储集体连通情况所对应的权值为1；若钻遇区溶孔仅与微裂缝连通，那么连通性一般，储集体连通情况所对应的权值为0；若钻遇区仅存在微裂缝，那么连通性较差，储集体连通情况所对应的权值为-1。
[0168]
对于该口井的岩心储集体发育情况来说，若孔、洞、缝均较发育，那么岩心储集体发育情况所对应的权值为1；若孔、缝较发育，那么岩心储集体发育情况所对应的权值为0；如果为其它情况，那么岩心储集体发育情况所对应的权值为-1。
[0169]
对于生产井的初期日产油量，本实施例中，以生产井的生产首月的平均值来进行分析。若初期日产油量大于100m3，那么初期日产油量所对应的权值为1；若初期日产油量在10～100m3，那么初期日产油量所对应的权值为0；若初期日产油量小于10m3，那么初期日产
油量所对应的权值为-1。
[0170]
对于生产井的累积产油量，本实施例中，以累积到当前时间的产油量来进行分析。若累积产油量大于10
×
104m3，那么累积产油量所对应的权值为1；若累积产油量在1
×
104～10
×
104m3之间，那么累积产油量所对应的权值为0；若累积产油量小于1
×
104m3，那么累积产油量所对应的权值为-1。
[0171]
对于生产井的油压变化幅度，本实施例中，以连续两天监测得到的油压来进行分析。若油压变化幅度小于5％，那么油压变化幅度所对应的权值为1；若油压变化幅度在0.5％～5％之间，那么油压变化幅度所对应的权值为0；若油压变化幅度小于0.5％，那么油压变化幅度所对应的权值为-1。
[0172]
对于生产井的产量递减情况，若产量递减小于1％，那么产量递减情况所对应的权值为1；若产量递减为1％～10％，那么产量递减情况所对应的权值为0；若产量递减大于10％，那么产量递减情况所对应的权值为-1。
[0173]
这样，便可以根据各口生产井的各个状态参数确定出各自所对应的权重。随后，本实施例中，基于各口生产井的各个状态参数所对应的权重，采用聚类分析的方式对各口生产井进行分类，从而得到单井分类结果。具体的，基于数据矩阵，利用聚类分析的方式便可以得到单井分类结果，即各口生产井的类型。
[0174]
其中，目标井可以分为第一类井(即ⅰ类井)、第二类井(即ⅱ类井)和第三类井(即ⅲ类井)。
[0175]
其中，第一类井为钻遇较好的缝洞储集体发育带的生产井，其地质特征包括：储集体以溶洞为主、周围裂缝发育较密、存在能够起导流作用的大裂缝和储集体之间的连通性较好。
[0176]
第二类井为钻遇裂缝较发育区的生产井，其地质特征包括：钻遇区存在小的溶孔和/或溶洞、存在较大的裂缝和储集体之间的连通性较好。
[0177]
第三类井为钻遇仅发育孤立溶洞和/或溶孔、小裂缝等储集体的生产井，其地质特性为储集体间的连通性很差，流体很难流动。
[0178]
而所述目标井的单井控制面积的确定步骤中，当所述目标井为已进行试井试验的生产井，则根据试井数据和/或生产动态数据确定所述目标井的单井控制面积；如果目标井为未进行试井试验的生产井，则根据生产动态数据确定所述目标井的单井控制面积。
[0179]
具体的所述目标井的单井控制面积的确定步骤具体包括：根据所述生产动态数据，确定所述目标井的井底压力；根据所述井底压力和孔隙介质模型，确定所述目标井的控制形状；基于所述控制形状，结合所述井控储量参数，确定所述目标井的单井控制面积。其中，所述井控储量参数包括：储层孔隙度、原油粘度和综合压缩系数。
[0180]
当所述控制形状为圆形或类圆形时，根据如下表达式计算所述目标井的单井控制面积：
[0181][0182]
当所述控制形状为矩形或类矩形时，根据如下表达式计算所述目标井的单井控制面积：
[0183][0184]
其中，a表示单井控制面积，k表示地层渗透率，t
s
表示拟稳态时间段的开始时间，t表示生产时间，φ表示储层孔隙度，μ表示原油粘度，c
t
表示综合压缩系数。
[0185]
而所述目标井的控制形状是根据所述井底压力，通过曲线拟合确定所述目标井的控制形状。
[0186]
当所述控制形状为圆形或类圆形时，在确定所述目标井的单井控制面积时，还根据拟合曲线确定所述目标井的拟稳态时间段，并根据所述拟稳态时间段中的井底压力来对所述井控储量参数进行修正，得到修正后的井控储量参数；从而根据修正后的井控储量参数，确定所述目标井的单井控制面积。而如果所述控制形状为圆形或类圆形但拟合曲线中并未含有拟稳态时间段时，则根据含水率小于预设含水率阈值的生产段来确定目标井的单井控制面积，从而计算出所述井控储量目标值。
[0187]
本公开实施例提供的一种油藏地质模型的优化方法，所述方法包括确定油藏地质模型的井控边界；对所述井控边界范围内的井控储量进行计算，并确定所述井控边界范围内的井控储量的影响因素；根据所述井控边界范围内的井控储量的影响因素，选取相应的优化方法对所述油藏地质模型进行优化。该方法对油藏地质模型进行局部的优化，使优化后的油藏地质模型的井控储量更接近于根据动态资料计算的目标井控储量，提高了油藏地质模型精度和准确性，从而优化了模拟效果。而且该方法在保持原有溶洞储集体形态不变的情况下对不同类型的储集体的体积参数和/或属性参数进行优化，灵活度高，操作性强。
[0188]
实施例二
[0189]
在实施例一的基础上，本实施例通过具体实施案例来对实施例一中所述的方法进行说明。
[0190]
本实施例是基于塔河地区某单元油藏地质模型和动态数据，选取该单元的一口生产井w井，用上述方法对油藏地质模型进行优化。
[0191]
如图10所示，从所建立的油藏地质模型中可以看出，w井钻遇一溶洞，该溶洞的储量为17万吨，远远低于动态法所预测(动态储量)的65万吨。由于w井进行了酸化压裂措施，在井筒附近产生了一系列的裂缝，连通了附近一定范围内的溶洞储集体。因此动态法所预测的动态储量应是w井连通的溶洞储集体的储量之和。
[0192]
在对w井的井油藏地质模型优化前，应用上述方法确定w井的井控边界。首先，确定w井附近一定距离(第一预设范围)不同方位的邻近井(tk1至tk5)，依次连接w井与各邻近井连线的中点，形成一个闭合的多边形作为第一边界(即初始井控边界)，然后考虑的遮挡作用。观察到w井和tk3、tk4、tk5三口井之间有断层，且w井到这三条断层的距离大于断层遮挡距离阈值d
f
(即预设距离)，因此将这三条断层作为边界线的一部分，以得到更新后的井控边界即第二边界。因为酸化压裂范围位于第二边界内，且第二边界未穿过孤立溶洞，而且本实施例中，不考虑目标井w井的压力传播范围边界的影响，所以，第二边界即为最终的井控边界，如图11所示。
[0193]
然后，在确定了井控边界后，要进一步确定井控边界内的溶洞储集体，根据w井射孔段深度和该井与各溶洞储集体的连通性，筛选掉那些和w井射孔段深度相差较大，和w井不连通的溶洞储集体，得到w井井控范围内的溶洞储集体，如图12和13所示，并将这些溶洞
储集体作为后续优化的对象。可以看到，w井的井控边界范围内共有4个溶洞储集体，其中，w井钻遇一个溶洞，其余3个溶洞分布在w井一定范围之内。利用井控储量计算公式计算得到这四个溶洞的储量。从计算结果可以看出井控边界范围内的井控储量为49.2万吨，与动态储量(65万吨)还有一定差距，考虑到该井钻遇的溶洞为非放空洞且位于井点处，具有相对准确的属性参数值，其储量受体积参数不确定性影响更大，所以选取体积优化法对油藏地质模型进行优化，优化过程中，在井控边界范围内的溶洞储集体在保持原有形态不变的前提下，对溶洞的体积进行了优化，如图14和15所示。表1为优化前后各溶洞的储量统计表。
[0194]
表1
[0195][0196]
如表1所示，经过以上优化，w井的井控储量和动态储量(65万吨)基本一致，既保证了溶洞储集体的形态和分布受各类静态资料的约束，其储量又与生产动态数据吻合。
[0197]
实施例三
[0198]
图16为本公开实施例提供的一种油藏地质模型的优化装置的结构示意图，请参阅图16，本实施例提供一种油藏地质模型的优化装置100，包括边界确定模块101、储量计算和影响因素确定模块102，以及模型优化模块103。
[0199]
边界确定模块101，用于确定油藏地质模型的井控边界；
[0200]
储量计算和影响因素确定模块102，用于对所述井控边界范围内的井控储量进行计算，并确定所述井控边界范围内的井控储量的影响因素；
[0201]
模型优化模块103，用于根据所述井控边界范围内的井控储量的影响因素，选取相应的优化方法对所述油藏地质模型进行优化。
[0202]
边界确定模块101确定油藏地质模型的井控边界，储量计算和影响因素确定模块102对所述井控边界范围内的井控储量进行计算，并确定所述井控边界范围内的井控储量的影响因素，最后模型优化模块103根据所述井控边界范围内的井控储量的影响因素，选取相应的优化方法对所述油藏地质模型进行优化。
[0203]
基于上述各模块执行油藏地质模型的优化方法的具体实施例已在实施例一中详述，此处不再赘述。
[0204]
实施例四
[0205]
本技术实施例提供了一种电子设备，该电子设备可以是手机、电脑或平板电脑等，包括存储器和处理器，所述存储器上存储有计算器程序，该计算机程序被处理器执行时实现如实施例一中所述的油藏地质模型的优化方法。可以理解，电子设备还可以包括，多媒体组件，输入/输出(i/o)接口，以及通信组件。
[0206]
其中，处理器用于执行如实施例一中的油藏地质模型的优化方法中的全部或部分
步骤。存储器用于存储各种类型的数据，这些数据例如可以包括电子设备中的任何应用程序或方法的指令，以及应用程序相关的数据。
[0207]
所述处理器可以是专用集成电路(application specific integrated circuit，简称asic)、数字信号处理器(digital signal processor，简称dsp)、数字信号处理设备(digital signal processing device，简称dspd)、可编程逻辑器件(programmable logic device，简称pld)、现场可编程门阵列(field programmable gate array，简称fpga)、控制器、微控制器、微处理器或其他电子元件实现，用于执行上述实施例一中的油藏地质模型的优化方法。
[0208]
所述存储器可以由任何类型的易失性或非易失性存储设备或者它们的组合实现，例如静态随机存取存储器(static random access memory，简称sram)，电可擦除可编程只读存储器(electrically erasable programmable read-only memory，简称eeprom)，可擦除可编程只读存储器(erasable programmable read-only memory，简称eprom)，可编程只读存储器(programmable read-only memory，简称prom)，只读存储器(read-only memory，简称rom)，磁存储器，快闪存储器，磁盘或光盘。
[0209]
实施例五
[0210]
本实施例还提供一种计算机可读存储介质，如闪存、硬盘、多媒体卡、卡型存储器(例如，sd或dx存储器等)、随机访问存储器(ram)、静态随机访问存储器(sram)、只读存储器(rom)、电可擦除可编程只读存储器(eeprom)、可编程只读存储器(prom)、磁性存储器、磁盘、光盘、服务器、app应用商城等等，其上存储有计算机程序，所述计算机程序被处理器执行时可以实现如下方法步骤：
[0211]
步骤s1：确定油藏地质模型的井控边界；
[0212]
步骤s2：对所述井控边界范围内的井控储量进行计算，并确定所述井控边界范围内的井控储量的影响因素；
[0213]
步骤s3：根据所述井控边界范围内的井控储量的影响因素，选取相应的优化方法对所述油藏地质模型进行优化。
[0214]
上述方法步骤的具体实施例过程可参见实施例一，本实施例在此不再重复赘述。
[0215]
综上，本公开提供的一种油藏地质模型的优化方法、装置、电子设备及存储介质，所述方法包括确定油藏地质模型的井控边界；对所述井控边界范围内的井控储量进行计算，并确定所述井控边界范围内的井控储量的影响因素；根据所述井控边界范围内的井控储量的影响因素，选取相应的优化方法对所述油藏地质模型进行优化。该方法对油藏地质模型进行局部的优化，使优化后的油藏地质模型的井控储量更接近于根据动态资料计算的目标井控储量，提高了油藏地质模型精度和准确性，从而优化了模拟效果。而且该方法在保持原有储集体形态不变的情况下对不同类型的储集体的体积参数和/或属性参数进行优化，灵活度高，操作性强。
[0216]
在本公开实施例所提供的几个实施例中，应该理解到，所揭露的方法，也可以通过其它的方式实现。以上所描述的方法实施例仅仅是示意性的。
[0217]
需要说明的是，在本文中，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、物品或者设备所固有
的要素。在没有更多限制的情况下，由语句“包括一个
……”
限定的要素，并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。
[0218]
虽然本公开所揭露的实施方式如上，但所述的内容只是为了便于理解本公开而采用的实施方式，并非用以限定本公开。任何本公开所属技术领域内的技术人员，在不脱离本公开所揭露的精神和范围的前提下，可以在实施的形式上及细节上作任何的修改与变化，但本公开的专利保护范围，仍须以所附的权利要求书所界定的范围为准。