一种基于凝胶泡沫的数值模拟方法
一种基于凝胶泡沫的数值模拟方法
【技术领域】
[0001] 本发明涉及油藏数值模拟技术应用领域，尤其涉及一种基于凝胶泡沫的数值模拟 方法。
【背景技术】
[0002] 泡沫（例如空气泡沫、氮气泡沫、二氧化碳泡沫等）调驱作为一种提高采收率的方 法在众多油田都得到了广泛的应用。但是在实际应用过程中，相关人员发现：普通泡沫封堵 能力有限，当遇到渗透率较大的窜流通道时，普通泡沫并不能有效封堵。另外，可以在封堵 渗透率较大的窜流通道上使用主要堵剂类型为冻胶类，冻胶虽然具备较好的封堵能力，但 是注入能力较差且经济成本较高。凝胶泡沫作为介于普通泡沫与冻胶之间的堵剂类型，兼 具了普通泡沫注入性能优异经济成本较低和冻胶封堵能力强的双重特点，但现有技术中没 有基于凝胶泡沫的数值模拟方式，无法为凝胶泡沫作为堵剂的实际应用提供支持。
【发明内容】
[0003] 本发明的目的是提供一种基于凝胶泡沫的数值模拟方法，以解决现有技术中没有 基于凝胶泡沫的数值模拟方式，以致于无法为凝胶泡沫作为堵剂的实际应用提供支持的问 题。
[0004] 本发有中一种基于凝胶泡沫的数值模拟方法包括：
[0005] 步骤1，构建油田或模拟油田实验装置的数学模型，将所述数学模型设置为3相9 组分形式，所述3相包括水相、油相和气相；所述9组分包括水、油、氮气、起泡剂、普通泡沫、 聚合物、交联剂、凝胶、凝胶泡沫；为各组分浓度设置初始值，为所述数学模型中的压力设置 初始值，为各组分粘度设置初始值，设置凝胶泡沫生成反应的反应速率和凝胶泡沫破灭反 应的反应速率，设置普通泡沫生成反应的反应速率和普通泡沫破灭反应的反应速率，设置 油水相对渗透率曲线以及气油相对渗透率曲线；设置时间步长和时间步数，以所述时间步 长和所述时间步数循环执行步骤2至步骤5 :
[0006] 步骤2,求解水相压力方程，油相压力方程，气相压力方程；
[0007] 步骤3,计算气相流速和凝胶泡沫液膜张力；其中，凝胶泡沫液膜张力为凝胶泡沫 内的气相压力与凝胶泡沫外壁的水相压力的差；
[0008] 步骤4,比较所述气相流速与临界流速的大小，在比较结果为所述气相流速大于 所述临界流速时执行泡沫生成反应；并且，比较凝胶泡沫液膜张力与临界压力的大小，并在 比较结果为所述凝胶泡沫液膜张力大于所述临界压力时执行泡沫破灭反应；其中，所述泡 沫生成反应包括普通泡沫生成反应和凝胶泡沫生成反应；将所述普通泡沫生成反应设置为 以下方式：由水、氮气、起泡剂生成普通泡沫；将所述凝胶泡沫生成反应设置为以下方式： 由水、氮气、起泡剂、交联剂、聚合物生成普通泡沫；所述泡沫破灭反应包括普通泡沫破灭反 应和凝胶泡沫破灭反应；将所述普通破灭生成反应设置为以下方式：由普通泡沫分解为氮 气、水、起泡剂；将所述凝胶泡沫生成反应设置为以下方式：由凝胶泡沫分解为水、氮气、起 泡剂、凝胶；
[0009] 步骤5,求解各组分浓度和各组分的粘度；修正凝胶泡沫的粘度和普通泡沫的粘 度，修正所述气油相对渗透率曲线。
[0010] 进一步地，本方法还可以具有以下特点：
[0011] 所述凝胶泡沫生成反应的反应速率的设置方法为：根据不同时间步的凝胶泡沫的 粘度的变化，设置凝胶泡沫的反应速率，具体为：根据实验数据获得相邻时间间隔相同的N 个时间点上的N个凝胶泡沫粘度，根据所述N个时间点与所述N个凝胶泡沫粘度，确定所述 凝胶泡沫生成反应的反应速率的值，使所述反应速率作为拟合函数的一参数使所述拟合函 数在时间点与凝胶泡沫粘度此两个坐标轴上能够拟合所述N个凝胶泡沫粘度，所述N为大 于1的整数。
[0012] 进一步地，本方法还可以具有以下特点：
[0013] 所述步骤1中设置气油相对渗透率曲线的过程中包括：将气相相对渗透率设置为 原始气相相对渗透率；
[0014] 所述步骤5中还包括：通过气相相对渗透率修正系数R对气相相对渗透率数值进 行修正，所述步骤5中还包括：判断普通泡沫或凝胶泡沫破灭时，将气相相对渗透率设置为 原始气相相对渗透率。
[0015] 进一步地，本方法还可以具有以下特点：
[0016] 获取所述气相相对渗透率修正系数R的方法包括：通过室内实验测试得到不同渗 透率的凝胶泡沫残余阻力系数，根据所述凝胶泡沫残余阻力系数通过历史拟合方法确定所 述气相相对渗透率修正系数R。
[0017] 进一步地，本方法还可以具有以下特点：
[0018] 根据下述方法对所述气相相对渗透率进行修正：
[0019] 满足
或cws>cwslim时，根据下述公式修正：
[0020] 满足& >(C+4或Cws彡cwslim时，根据下述公式进行修正：
[0022] 式中，为修正前的气相相对渗透率，为修正后的气相相对渗透率，< 为临界 含水饱和度，Cwslim为存在泡膜的最低起泡剂浓度，e为饱和度小范围变化量，SWS含水饱 和度也称为水相饱和度，Cws为起泡性剂浓度。
[0023] 进一步地，本方法还可以具有以下特点：
[0024] 所述气相相对渗透率修正系数R是80至106范围内的一个值。
[0025] 进一步地，本方法还可以具有以下特点：
[0026] 所述气相相对渗透率修正系数R的值为86。
[0027] 进一步地，本方法还可以具有以下特点：
[0028] 通过实验方法获得所述普通泡沫破灭反应的反应速率和所述凝胶泡沫破灭反应 的反应速率，具体包括：
[0029] 打开循环水浴锅并加热到80至100摄氏度之间，连接循环水浴锅和起泡仪；
[0030] 配制起泡剂溶液10-30毫升；
[0031] 向起泡仪注入所述起泡液；
[0032] 按照10至30毫升每分钟的速度注入氮气，气液比为1 : 1至2:1之间，观察普通 泡沫和凝胶泡沫的起泡情况，记录普通泡沫和凝胶泡沫的起泡体积及半衰期，重复测量多 次，取测量到的多个半衰期取平均值，得到普通泡沫的半衰期和凝胶泡沫的半衰期，根据所 述普通泡沫的半衰期计算得到所述普通泡沫破灭反应的反应速度，根据所述凝胶泡沫的半 衰期计算得到所述泡沫破灭反应的反应速度。
[0033] 进一步地，本方法还可以具有以下特点：
[0034] 设置所述普通泡沫破灭反应的反应速率的值位于0. 01至0. 012范围内，
[0035] 设置所述凝胶泡沫破灭反应的反应速率的值位于0. 001至0. 003范围内。
[0036] 进一步地，本方法还可以具有以下特点：
[0037] 设置所述普通泡沫破灭反应的反应速率的值为0. 011，
[0038] 设置所述凝胶泡沫破灭反应的反应速率的值为0. 002。
[0039] 本发明的基于凝胶泡沫的数值模拟方法，在引入凝胶泡沫的基础上建立了用于调 驱分析的数学模型，提出了有效表征凝胶泡沫调驱的数值模拟方法，为油藏进一步提高采 收率的数值模拟提供了有效的方法和途径。
【附图说明】
[0040] 图1是本发明中基于凝胶泡沫的数值模拟方法的流程图。 具体实施例
[0041] 为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例 中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是 本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员 在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。需要 说明的是，在不冲突的情况下，本申请中的实施例及实施例中的特征可以相互任意组合。
[0042] 本方法包括：
[0043] 步骤1， 构建油田或模拟油田实验装置的数学模型，将此数学模型设置为3相9组 分形式，其中3相包括水相、油相和气相；9组分包括水、油、氮气、起泡剂、普通泡沫、聚合 物、交联剂、凝胶、凝胶泡沫。为各组分浓度设置初始值，为所述数学模型中的压力设置初始 值，为各组分粘度设置初始值，设置凝胶泡沫生成反应的反应速率和凝胶泡沫破灭反应的 反应速率，设置普通泡沫生成反应的反应速率和普通泡沫破灭反应的反应速率；设置油水 相对渗透率曲线以及气油相对渗透率曲线；设置时间步长和时间步数，以上述时间步长循 环执行步骤2至步骤5,执行次数与上述时间步数相同。例如，可以设置时间步长为1小时， 时间步数为100步。
[0044] 步骤2,求解水相压力方程，油相压力方程，气相压力方程。
[0045] 步骤3,计算气相流速和凝胶泡沫液膜张力；其中，凝胶泡沫液膜张力为凝胶泡沫 内的气相压力与凝胶泡沫外壁的水相压力的差。
[0046] 步骤4,比较气相流速与临界流速的大小，在比较结果为气相流速大于临界流速时 执行泡沫生成反应；并且，比较凝胶泡沫液膜张力与临界压力的大小，并在比较结果为凝胶 泡沫液膜张力大于临界压力时执行泡沫破灭反应。
[0047] 本步骤4中，泡沫生成反应包括普通泡沫生成反应和凝胶泡沫生成反应；将普通 泡沫生成反应设置为以下方式：由水、氮气、起泡剂生成普通泡沫，即水+氮气+起泡剂=普 通泡沫。
[0048] 将凝胶泡沫生成反应设置为以下方式：由水、氮气、起泡剂、交联剂、聚合物生成普 通泡沫，即水+氮气+起泡剂+交联剂+聚合物=凝胶泡沫。
[0049] 泡沫破灭反应包括普通泡沫破灭反应和凝胶泡沫破灭反应。凝胶泡沫和普通泡沫 的破灭速度受控于泡沫的动态半衰期。
[0050] 将普通破灭生成反应设置为以下方式：由普通泡沫分解为氮气、水、起泡剂，即普 通泡沫=氮气+水+起泡剂。
[0051] 将凝胶泡沫生成反应设置为以下方式：由凝胶泡沫分解为水、氮气、起泡剂、凝胶， 即凝胶泡沫=水+氮气+起泡剂+凝胶。
[0052] 本步骤4中，根据步骤1中设置的普通泡沫生成反应的反应速率执行普通泡沫生 成反应，根据步骤1中设置的凝胶泡沫生成反应的反应速率执行凝胶泡沫生成反应，根据 步骤1中设置的普通泡沫破灭反应的反应速率执行普通泡沫破灭反应，根据步骤1中设置 的凝胶泡沫破灭反应的反应速率执行凝胶泡沫破灭反应。
[0053] 步骤5,求解各组分浓度和各组分的粘度；修正凝胶泡沫的粘度和普通泡沫的粘 度，修正气油相对渗透率曲线。
[0054] 下面详细说明本方法。
[0055] 步骤1中构建的数学模型满足以下条件：
[0056] (1)模型设置为3相9组分，具体组分设置与相平衡见表1
[0057] 表1模型组分设置
[0059] (2)驱替过程为等温过程，化学反应无热量消耗和释放。
[0060] (3)注入流体和油藏流体满足达西定律。
[0061] (4)瞬间达到相间平衡。
[0062] (5)组分弥散规律遵循Fick定律。
[0063] (6)满足理想状态的混合规则，即混合时体积变化为零。
[0064] 步骤2中求解水相压力方程，油相压力方程，气相压力方程的方法为已知技术，此 处不再赘述。
[0065] 普通泡沫生成后，普通泡沫壁面的聚合物与交联剂发生化学反应，生成凝胶泡沫， 不同的化学反应时间生成的凝胶泡沫粘度不同，凝胶泡沫生成后，在不破灭的情况下，其粘 度不变。本发明的步骤1中根据不同时间步的凝胶泡沫的粘度的变化设置凝胶泡沫的反应 速率，使模型能够准确模拟凝胶泡沫的粘度变化。具体的，凝胶泡沫生成反应的反应速率 的设置方法为：根据不同时间步的凝胶泡沫的粘度的变化，设置凝胶泡沫的反应速率，具体 为：根据实验数据获得相邻时间间隔相同的N个时间点上的N个凝胶泡沫粘度，根据N个时 间点与N个凝胶泡沫粘度，确定凝胶泡沫生成反应的反应速率的值，使反应速率作为拟合 函数的一参数使拟合函数在时间点与凝胶泡沫粘度此两个坐标轴上能够拟合N个凝胶泡 沫粘度，N为大于1的整数。
[0066]例如：根据室内实验，得到实验参数，如表2所示：
[0067]表2凝胶泡沫粘度设置参数
[0069]根据实验数据获得相邻时间间隔相同的3个时间点上的3个凝胶泡沫粘度，根据 此3个时间点与3个凝胶泡沫粘度，确定凝胶泡沫生成反应的反应速率的值，使反应速率作 为拟合函数的一参数使拟合函数在时间点与凝胶泡沫粘度此两个坐标轴上能够拟合上述3 个凝胶泡沫粘度。
[0070]步骤1中设置气油相对渗透率曲线的过程中包括：将气相相对渗透率设置为原始 气相相对渗透率。步骤5中还包括：判断普通泡沫或凝胶泡沫破灭时，将气相相对渗透率设 置为原始气相相对渗透率。
[0071]凝胶泡沫的封堵机理主要有两个，第一是增加气相粘度，第二是降低气相相对渗 透率，在泡沫气相相对渗透率的修正中同样需要考虑泡沫生成、聚井、破灭对泡沫体系的影 响。泡沫的破灭与凝胶泡沫液膜张力有关，进而与含水饱和度有关。当有泡沫存在时，只有 一部分气相在流动，泡沫的有效渗透率受到体系渗透率降低程度的影响。所以步骤5中还 包括：通过气相相对渗透率修正系数R对气相相对渗透率数值进行修正。具体的，根据下述 方法对所述气相相对渗透率进行修正：
[0074]式中，为修正前的气相相对渗透率，为修正后的气相相对渗透率为临界 含水饱和度，Cwslim为存在泡膜的最低起泡剂浓度，e为饱和度小范围变化量，SWS含水饱 和度也称为水相饱和度，Cws为起泡性剂浓度。
[0075] 在此修正方法中，获取气相相对渗透率修正系数R的方法包括两种：
[0076] 第一种：通过室内实验测试得到不同渗透率的凝胶泡沫残余阻力系数，举例如 下：
[0077] 通过室内实验测试得到不同渗透率的凝胶泡沫残余阻力系数，如表3所示：
[0078] 表3不同渗透率对应的凝胶泡沫残余阻力系数
[0081] 然后，根据凝胶泡沫残余阻力系数通过历史拟合方法确定气相相对渗透率修正系 数R。
[0082] 第二种：采用优选范围（此优先范围通过多次实验确定出），气相相对渗透率修正 系数R的优选范围为80至106,优选的采用86。
[0083] 步骤1中获得普通泡沫破灭反应的反应速率和凝胶泡沫破灭反应的反应速率的 方法有两种：
[0084] 第一种，通过实验方法获得。
[0085] 泡沫破灭的反应速率取决于泡沫的半衰期。实验方法具体方法包括：
[0086] 打开循环水浴锅并加热到80至100摄氏度之间，连接循环水浴锅和起泡仪；
[0087] 配制起泡剂溶液10-30毫升；
[0088] 向起泡仪注入起泡液；
[0089] 按照10至30毫升每分钟的速度注入氮气，气液比为1 : 1至1.5:1之间，观察普 通泡沫和凝胶泡沫的起泡情况，记录普通泡沫和凝胶泡沫的起泡体积及半衰期，重复测量 多次，取测量到的多个半衰期取平均值，得到普通泡沫的半衰期和凝胶泡沫的半衰期，根据 普通泡沫的半衰期计算得到普通泡沫破灭反应的反应速度，根据凝胶泡沫的半衰期计算得 到泡沫破灭反应的反应速度。
[0090] 使用上述方法产生的泡沫比较细小致密，这种结构有利于泡沫的稳定性，大而松 散的气泡排液速度快，液膜变薄处来不及修复，容易破灭。在实验中观察到，在0. 05%~ 〇. 4%的较宽浓度范围内，起泡剂的起泡体积随浓度增加呈现一个上升再到略有下降的平 缓趋势，这是由于在低于临界胶束浓度时，随着浓度的增加，注入一定氮气后，表面活性剂 分子在氮气/水界面上的吸附量增大，界面张力降低，产生的泡沫体系自由能降低， 能形成 更多的气液界面，产生更多的气泡，使泡沫体积增加，在临界胶束浓度时达到最高值，但是 在达到临界胶束浓度后，起泡性能不再提高，这主要是因为溶液表面张力不再减小，体系的 表面能趋于不变，另外形成的胶束在一定程度会增加液相的粘度，会使起泡体积受到一定 影响。
[0091] 第二种，采用优选范围（此优先范围通过多次实验确定出），普通泡沫破灭反应的 反应速率的值的优选范围为0. 01至0. 012,优选采用0. 011。凝胶泡沫破灭反应的反应速 率的值的优选范围为〇. 001至〇. 003,优先采用0. 002。
[0092] 本发明的基于凝胶泡沫的数值模拟方法，在引入凝胶泡沫的基础上建立了用于调 驱分析的数学模型，提出了有效表征凝胶泡沫调驱的数值模拟方法，对油藏进一步提高采 收率的数值模拟提供了有效的方法和途径。
[0093] 上面描述的内容可以单独地或者以各种方式组合起来实施，而这些变型方式都在 本发明的保护范围之内。
[0094] 需要说明的是，在本文中，术语"包括"、"包含"或者其任何其他变体意在涵盖非排 他性的包含，从而使得包括一系列要素的物品或者设备不仅包括那些要素，而且还包括没 有明确列出的其他要素，或者是还包括为这种物品或者设备所固有的要素。在没有更多限 制的情况下，由语句"包括……"限定的要素，并不排除在包括所述要素的物品或者设备中 还存在另外的相同要素。
[0095] 以上实施例仅用以说明本发明的技术方案而非限制，仅仅参照较佳实施例对本发 明进行了详细说明。本领域的普通技术人员应当理解，可以对本发明的技术方案进行修改 或者等同替换，而不脱离本发明技术方案的精神和范围，均应涵盖在本发明的权利要求范 围当中。
【主权项】
1. 一种基于凝胶泡沫的数值模拟方法，其特征在于，所述方法包括： 步骤1，构建油田或模拟油田实验装置的数学模型，将所述数学模型设置为3相9组分 形式，所述3相包括水相、油相和气相；所述9组分包括水、油、氮气、起泡剂、普通泡沫、聚合 物、交联剂、凝胶、凝胶泡沫；为各组分浓度设置初始值，为所述数学模型中的压力设置初始 值，为各组分粘度设置初始值，设置凝胶泡沫生成反应的反应速率和凝胶泡沫破灭反应的 反应速率，设置普通泡沫生成反应的反应速率和普通泡沫破灭反应的反应速率，设置油水 相对渗透率曲线以及气油相对渗透率曲线；设置时间步长和时间步数，以所述时间步长和 所述时间步数循环执行步骤2至步骤5 : 步骤2,求解水相压力方程，油相压力方程，气相压力方程； 步骤3,计算气相流速和凝胶泡沫液膜张力；其中，凝胶泡沫液膜张力为凝胶泡沫内的 气相压力与凝胶泡沫外壁的水相压力的差； 步骤4,比较所述气相流速与临界流速的大小，在比较结果为所述气相流速大于所述临 界流速时执行泡沫生成反应；并且，比较凝胶泡沫液膜张力与临界压力的大小，并在比较结 果为所述凝胶泡沫液膜张力大于所述临界压力时执行泡沫破灭反应；其中，所述泡沫生成 反应包括普通泡沫生成反应和凝胶泡沫生成反应；将所述普通泡沫生成反应设置为以下方 式：由水、氮气、起泡剂生成普通泡沫；将所述凝胶泡沫生成反应设置为以下方式：由水、氮 气、起泡剂、交联剂、聚合物生成普通泡沫；所述泡沫破灭反应包括普通泡沫破灭反应和凝 胶泡沫破灭反应；将所述普通破灭生成反应设置为以下方式：由普通泡沫分解为氮气、水、 起泡剂；将所述凝胶泡沫生成反应设置为以下方式：由凝胶泡沫分解为水、氮气、起泡剂、 凝胶； 步骤5,求解各组分浓度和各组分的粘度；修正凝胶泡沫的粘度和普通泡沫的粘度，修 正所述气油相对渗透率曲线。2. 如权利要求1所述的方法，其特征在于， 所述凝胶泡沫生成反应的反应速率的设置方法为：根据不同时间步的凝胶泡沫的粘度 的变化，设置凝胶泡沫的反应速率，具体为：根据实验数据获得相邻时间间隔相同的N个时 间点上的N个凝胶泡沫粘度，根据所述N个时间点与所述N个凝胶泡沫粘度，确定所述凝胶 泡沫生成反应的反应速率的值，使所述反应速率作为拟合函数的一参数使所述拟合函数在 时间点与凝胶泡沫粘度此两个坐标轴上能够拟合所述N个凝胶泡沫粘度，所述N为大于1 的整数。3. 如权利要求1所述的方法，其特征在于， 所述步骤1中设置气油相对渗透率曲线的过程中包括：将气相相对渗透率设置为原始 气相相对渗透率； 所述步骤5中还包括：通过气相相对渗透率修正系数R对气相相对渗透率数值进行修 正，所述步骤5中还包括：判断普通泡沫或凝胶泡沫破灭时，将气相相对渗透率设置为原始 气相相对渗透率。4. 如权利要求3所述的方法，其特征在于， 获取所述气相相对渗透率修正系数R的方法包括：通过室内实验测试得到不同渗透率 的凝胶泡沫残余阻力系数，根据所述凝胶泡沫残余阻力系数通过历史拟合方法确定所述气 相相对渗透率修正系数R。5. 如权利要求4所述的方法，其特征在于， 根据下述方法对所述气相相对渗透率进行修正： 满足pi. - 4 <丨< (S: + U或Cws彡C wslim时，根据下述公式修正：满足心 > 认:+ 4或C11J Cii5hm时，根据下述公式进行修正： ? ' I式中K为修正前的气相相对渗透率，为修正后的气相相对渗透率，< 为临界含水 饱和度，Cwslim为存在泡膜的最低起泡剂浓度，ε为饱和度小范围变化量，Sw为含水饱和度 也称为水相饱和度，Cws为起泡性剂浓度。6. 如权利要求3所述的方法，其特征在于， 所述气相相对渗透率修正系数R是80至106范围内的一个值。7. 如权利要求6所述的方法，其特征在于， 所述气相相对渗透率修正系数R的值为86。8. 如权利要求1所述的方法，其特征在于， 通过实验方法获得所述普通泡沫破灭反应的反应速率和所述凝胶泡沫破灭反应的反 应速率，具体包括： 打开循环水浴锅并加热到80至100摄氏度之间，连接循环水浴锅和起泡仪； 配制起泡剂溶液10-30毫升； 向起泡仪注入所述起泡液； 按照10至30毫升每分钟的速度注入氮气，气液比为I : 1至2:1之间，观察普通泡沫 和凝胶泡沫的起泡情况，记录普通泡沫和凝胶泡沫的起泡体积及半衰期，重复测量多次，取 测量到的多个半衰期取平均值，得到普通泡沫的半衰期和凝胶泡沫的半衰期，根据所述普 通泡沫的半衰期计算得到所述普通泡沫破灭反应的反应速度，根据所述凝胶泡沫的半衰期 计算得到所述泡沫破灭反应的反应速度。9. 如权利要求1所述的方法，其特征在于， 设置所述普通泡沫破灭反应的反应速率的值位于〇. 01至〇. 012范围内， 设置所述凝胶泡沫破灭反应的反应速率的值位于〇. 001至〇. 003范围内。10. 如权利要求9所述的方法，其特征在于， 设置所述普通泡沫破灭反应的反应速率的值为〇. 011， 设置所述凝胶泡沫破灭反应的反应速率的值为〇. 002。
【专利摘要】本发明公开了一种基于凝胶泡沫的数值模拟方法，本方法包括：构建3相9组分形式的数学模型；求解水相、油相、气相压力方程；计算气相流速和凝胶泡沫液膜张力；比较气相流速与临界流速的大小，在比较结果为气相流速大于临界流速时执行泡沫生成反应；比较凝胶泡沫液膜张力与临界压力的大小，并在比较结果为凝胶泡沫液膜张力大于临界压力时执行泡沫破灭反应；求解各组分浓度和各组分的粘度；修正凝胶泡沫的粘度和普通泡沫的粘度，修正气油相对渗透率曲线。本方法在引入凝胶泡沫的基础上建立了用于调驱分析的数学模型，提出了有效表征凝胶泡沫调驱的数值模拟方法，为油藏进一步提高采收率的数值模拟提供了有效手段。
【IPC分类】G06F19/00
【公开号】CN104899438
【申请号】CN201510295337
【发明人】王硕亮, 张媛, 于春磊, 付强, 李俊键, 刘同敬
【申请人】中国地质大学（北京）
【公开日】2015年9月9日
【申请日】2015年6月2日
【技术领域】
[0001] 本发明涉及油藏数值模拟技术应用领域，尤其涉及一种基于凝胶泡沫的数值模拟 方法。
【背景技术】
[0002] 泡沫（例如空气泡沫、氮气泡沫、二氧化碳泡沫等）调驱作为一种提高采收率的方 法在众多油田都得到了广泛的应用。但是在实际应用过程中，相关人员发现：普通泡沫封堵 能力有限，当遇到渗透率较大的窜流通道时，普通泡沫并不能有效封堵。另外，可以在封堵 渗透率较大的窜流通道上使用主要堵剂类型为冻胶类，冻胶虽然具备较好的封堵能力，但 是注入能力较差且经济成本较高。凝胶泡沫作为介于普通泡沫与冻胶之间的堵剂类型，兼 具了普通泡沫注入性能优异经济成本较低和冻胶封堵能力强的双重特点，但现有技术中没 有基于凝胶泡沫的数值模拟方式，无法为凝胶泡沫作为堵剂的实际应用提供支持。
【发明内容】
[0003] 本发明的目的是提供一种基于凝胶泡沫的数值模拟方法，以解决现有技术中没有 基于凝胶泡沫的数值模拟方式，以致于无法为凝胶泡沫作为堵剂的实际应用提供支持的问 题。
[0004] 本发有中一种基于凝胶泡沫的数值模拟方法包括：
[0005] 步骤1，构建油田或模拟油田实验装置的数学模型，将所述数学模型设置为3相9 组分形式，所述3相包括水相、油相和气相；所述9组分包括水、油、氮气、起泡剂、普通泡沫、 聚合物、交联剂、凝胶、凝胶泡沫；为各组分浓度设置初始值，为所述数学模型中的压力设置 初始值，为各组分粘度设置初始值，设置凝胶泡沫生成反应的反应速率和凝胶泡沫破灭反 应的反应速率，设置普通泡沫生成反应的反应速率和普通泡沫破灭反应的反应速率，设置 油水相对渗透率曲线以及气油相对渗透率曲线；设置时间步长和时间步数，以所述时间步 长和所述时间步数循环执行步骤2至步骤5 :
[0006] 步骤2,求解水相压力方程，油相压力方程，气相压力方程；
[0007] 步骤3,计算气相流速和凝胶泡沫液膜张力；其中，凝胶泡沫液膜张力为凝胶泡沫 内的气相压力与凝胶泡沫外壁的水相压力的差；
[0008] 步骤4,比较所述气相流速与临界流速的大小，在比较结果为所述气相流速大于 所述临界流速时执行泡沫生成反应；并且，比较凝胶泡沫液膜张力与临界压力的大小，并在 比较结果为所述凝胶泡沫液膜张力大于所述临界压力时执行泡沫破灭反应；其中，所述泡 沫生成反应包括普通泡沫生成反应和凝胶泡沫生成反应；将所述普通泡沫生成反应设置为 以下方式：由水、氮气、起泡剂生成普通泡沫；将所述凝胶泡沫生成反应设置为以下方式： 由水、氮气、起泡剂、交联剂、聚合物生成普通泡沫；所述泡沫破灭反应包括普通泡沫破灭反 应和凝胶泡沫破灭反应；将所述普通破灭生成反应设置为以下方式：由普通泡沫分解为氮 气、水、起泡剂；将所述凝胶泡沫生成反应设置为以下方式：由凝胶泡沫分解为水、氮气、起 泡剂、凝胶；
[0009] 步骤5,求解各组分浓度和各组分的粘度；修正凝胶泡沫的粘度和普通泡沫的粘 度，修正所述气油相对渗透率曲线。
[0010] 进一步地，本方法还可以具有以下特点：
[0011] 所述凝胶泡沫生成反应的反应速率的设置方法为：根据不同时间步的凝胶泡沫的 粘度的变化，设置凝胶泡沫的反应速率，具体为：根据实验数据获得相邻时间间隔相同的N 个时间点上的N个凝胶泡沫粘度，根据所述N个时间点与所述N个凝胶泡沫粘度，确定所述 凝胶泡沫生成反应的反应速率的值，使所述反应速率作为拟合函数的一参数使所述拟合函 数在时间点与凝胶泡沫粘度此两个坐标轴上能够拟合所述N个凝胶泡沫粘度，所述N为大 于1的整数。
[0012] 进一步地，本方法还可以具有以下特点：
[0013] 所述步骤1中设置气油相对渗透率曲线的过程中包括：将气相相对渗透率设置为 原始气相相对渗透率；
[0014] 所述步骤5中还包括：通过气相相对渗透率修正系数R对气相相对渗透率数值进 行修正，所述步骤5中还包括：判断普通泡沫或凝胶泡沫破灭时，将气相相对渗透率设置为 原始气相相对渗透率。
[0015] 进一步地，本方法还可以具有以下特点：
[0016] 获取所述气相相对渗透率修正系数R的方法包括：通过室内实验测试得到不同渗 透率的凝胶泡沫残余阻力系数，根据所述凝胶泡沫残余阻力系数通过历史拟合方法确定所 述气相相对渗透率修正系数R。
[0017] 进一步地，本方法还可以具有以下特点：
[0018] 根据下述方法对所述气相相对渗透率进行修正：
[0019] 满足
或cws>cwslim时，根据下述公式修正：
[0020] 满足& >(C+4或Cws彡cwslim时，根据下述公式进行修正：
[0022] 式中，为修正前的气相相对渗透率，为修正后的气相相对渗透率，< 为临界 含水饱和度，Cwslim为存在泡膜的最低起泡剂浓度，e为饱和度小范围变化量，SWS含水饱 和度也称为水相饱和度，Cws为起泡性剂浓度。
[0023] 进一步地，本方法还可以具有以下特点：
[0024] 所述气相相对渗透率修正系数R是80至106范围内的一个值。
[0025] 进一步地，本方法还可以具有以下特点：
[0026] 所述气相相对渗透率修正系数R的值为86。
[0027] 进一步地，本方法还可以具有以下特点：
[0028] 通过实验方法获得所述普通泡沫破灭反应的反应速率和所述凝胶泡沫破灭反应 的反应速率，具体包括：
[0029] 打开循环水浴锅并加热到80至100摄氏度之间，连接循环水浴锅和起泡仪；
[0030] 配制起泡剂溶液10-30毫升；
[0031] 向起泡仪注入所述起泡液；
[0032] 按照10至30毫升每分钟的速度注入氮气，气液比为1 : 1至2:1之间，观察普通 泡沫和凝胶泡沫的起泡情况，记录普通泡沫和凝胶泡沫的起泡体积及半衰期，重复测量多 次，取测量到的多个半衰期取平均值，得到普通泡沫的半衰期和凝胶泡沫的半衰期，根据所 述普通泡沫的半衰期计算得到所述普通泡沫破灭反应的反应速度，根据所述凝胶泡沫的半 衰期计算得到所述泡沫破灭反应的反应速度。
[0033] 进一步地，本方法还可以具有以下特点：
[0034] 设置所述普通泡沫破灭反应的反应速率的值位于0. 01至0. 012范围内，
[0035] 设置所述凝胶泡沫破灭反应的反应速率的值位于0. 001至0. 003范围内。
[0036] 进一步地，本方法还可以具有以下特点：
[0037] 设置所述普通泡沫破灭反应的反应速率的值为0. 011，
[0038] 设置所述凝胶泡沫破灭反应的反应速率的值为0. 002。
[0039] 本发明的基于凝胶泡沫的数值模拟方法，在引入凝胶泡沫的基础上建立了用于调 驱分析的数学模型，提出了有效表征凝胶泡沫调驱的数值模拟方法，为油藏进一步提高采 收率的数值模拟提供了有效的方法和途径。
【附图说明】
[0040] 图1是本发明中基于凝胶泡沫的数值模拟方法的流程图。 具体实施例
[0041] 为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例 中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是 本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员 在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。需要 说明的是，在不冲突的情况下，本申请中的实施例及实施例中的特征可以相互任意组合。
[0042] 本方法包括：
[0043] 步骤1， 构建油田或模拟油田实验装置的数学模型，将此数学模型设置为3相9组 分形式，其中3相包括水相、油相和气相；9组分包括水、油、氮气、起泡剂、普通泡沫、聚合 物、交联剂、凝胶、凝胶泡沫。为各组分浓度设置初始值，为所述数学模型中的压力设置初始 值，为各组分粘度设置初始值，设置凝胶泡沫生成反应的反应速率和凝胶泡沫破灭反应的 反应速率，设置普通泡沫生成反应的反应速率和普通泡沫破灭反应的反应速率；设置油水 相对渗透率曲线以及气油相对渗透率曲线；设置时间步长和时间步数，以上述时间步长循 环执行步骤2至步骤5,执行次数与上述时间步数相同。例如，可以设置时间步长为1小时， 时间步数为100步。
[0044] 步骤2,求解水相压力方程，油相压力方程，气相压力方程。
[0045] 步骤3,计算气相流速和凝胶泡沫液膜张力；其中，凝胶泡沫液膜张力为凝胶泡沫 内的气相压力与凝胶泡沫外壁的水相压力的差。
[0046] 步骤4,比较气相流速与临界流速的大小，在比较结果为气相流速大于临界流速时 执行泡沫生成反应；并且，比较凝胶泡沫液膜张力与临界压力的大小，并在比较结果为凝胶 泡沫液膜张力大于临界压力时执行泡沫破灭反应。
[0047] 本步骤4中，泡沫生成反应包括普通泡沫生成反应和凝胶泡沫生成反应；将普通 泡沫生成反应设置为以下方式：由水、氮气、起泡剂生成普通泡沫，即水+氮气+起泡剂=普 通泡沫。
[0048] 将凝胶泡沫生成反应设置为以下方式：由水、氮气、起泡剂、交联剂、聚合物生成普 通泡沫，即水+氮气+起泡剂+交联剂+聚合物=凝胶泡沫。
[0049] 泡沫破灭反应包括普通泡沫破灭反应和凝胶泡沫破灭反应。凝胶泡沫和普通泡沫 的破灭速度受控于泡沫的动态半衰期。
[0050] 将普通破灭生成反应设置为以下方式：由普通泡沫分解为氮气、水、起泡剂，即普 通泡沫=氮气+水+起泡剂。
[0051] 将凝胶泡沫生成反应设置为以下方式：由凝胶泡沫分解为水、氮气、起泡剂、凝胶， 即凝胶泡沫=水+氮气+起泡剂+凝胶。
[0052] 本步骤4中，根据步骤1中设置的普通泡沫生成反应的反应速率执行普通泡沫生 成反应，根据步骤1中设置的凝胶泡沫生成反应的反应速率执行凝胶泡沫生成反应，根据 步骤1中设置的普通泡沫破灭反应的反应速率执行普通泡沫破灭反应，根据步骤1中设置 的凝胶泡沫破灭反应的反应速率执行凝胶泡沫破灭反应。
[0053] 步骤5,求解各组分浓度和各组分的粘度；修正凝胶泡沫的粘度和普通泡沫的粘 度，修正气油相对渗透率曲线。
[0054] 下面详细说明本方法。
[0055] 步骤1中构建的数学模型满足以下条件：
[0056] (1)模型设置为3相9组分，具体组分设置与相平衡见表1
[0057] 表1模型组分设置
[0059] (2)驱替过程为等温过程，化学反应无热量消耗和释放。
[0060] (3)注入流体和油藏流体满足达西定律。
[0061] (4)瞬间达到相间平衡。
[0062] (5)组分弥散规律遵循Fick定律。
[0063] (6)满足理想状态的混合规则，即混合时体积变化为零。
[0064] 步骤2中求解水相压力方程，油相压力方程，气相压力方程的方法为已知技术，此 处不再赘述。
[0065] 普通泡沫生成后，普通泡沫壁面的聚合物与交联剂发生化学反应，生成凝胶泡沫， 不同的化学反应时间生成的凝胶泡沫粘度不同，凝胶泡沫生成后，在不破灭的情况下，其粘 度不变。本发明的步骤1中根据不同时间步的凝胶泡沫的粘度的变化设置凝胶泡沫的反应 速率，使模型能够准确模拟凝胶泡沫的粘度变化。具体的，凝胶泡沫生成反应的反应速率 的设置方法为：根据不同时间步的凝胶泡沫的粘度的变化，设置凝胶泡沫的反应速率，具体 为：根据实验数据获得相邻时间间隔相同的N个时间点上的N个凝胶泡沫粘度，根据N个时 间点与N个凝胶泡沫粘度，确定凝胶泡沫生成反应的反应速率的值，使反应速率作为拟合 函数的一参数使拟合函数在时间点与凝胶泡沫粘度此两个坐标轴上能够拟合N个凝胶泡 沫粘度，N为大于1的整数。
[0066]例如：根据室内实验，得到实验参数，如表2所示：
[0067]表2凝胶泡沫粘度设置参数
[0069]根据实验数据获得相邻时间间隔相同的3个时间点上的3个凝胶泡沫粘度，根据 此3个时间点与3个凝胶泡沫粘度，确定凝胶泡沫生成反应的反应速率的值，使反应速率作 为拟合函数的一参数使拟合函数在时间点与凝胶泡沫粘度此两个坐标轴上能够拟合上述3 个凝胶泡沫粘度。
[0070]步骤1中设置气油相对渗透率曲线的过程中包括：将气相相对渗透率设置为原始 气相相对渗透率。步骤5中还包括：判断普通泡沫或凝胶泡沫破灭时，将气相相对渗透率设 置为原始气相相对渗透率。
[0071]凝胶泡沫的封堵机理主要有两个，第一是增加气相粘度，第二是降低气相相对渗 透率，在泡沫气相相对渗透率的修正中同样需要考虑泡沫生成、聚井、破灭对泡沫体系的影 响。泡沫的破灭与凝胶泡沫液膜张力有关，进而与含水饱和度有关。当有泡沫存在时，只有 一部分气相在流动，泡沫的有效渗透率受到体系渗透率降低程度的影响。所以步骤5中还 包括：通过气相相对渗透率修正系数R对气相相对渗透率数值进行修正。具体的，根据下述 方法对所述气相相对渗透率进行修正：
[0074]式中，为修正前的气相相对渗透率，为修正后的气相相对渗透率为临界 含水饱和度，Cwslim为存在泡膜的最低起泡剂浓度，e为饱和度小范围变化量，SWS含水饱 和度也称为水相饱和度，Cws为起泡性剂浓度。
[0075] 在此修正方法中，获取气相相对渗透率修正系数R的方法包括两种：
[0076] 第一种：通过室内实验测试得到不同渗透率的凝胶泡沫残余阻力系数，举例如 下：
[0077] 通过室内实验测试得到不同渗透率的凝胶泡沫残余阻力系数，如表3所示：
[0078] 表3不同渗透率对应的凝胶泡沫残余阻力系数
[0081] 然后，根据凝胶泡沫残余阻力系数通过历史拟合方法确定气相相对渗透率修正系 数R。
[0082] 第二种：采用优选范围（此优先范围通过多次实验确定出），气相相对渗透率修正 系数R的优选范围为80至106,优选的采用86。
[0083] 步骤1中获得普通泡沫破灭反应的反应速率和凝胶泡沫破灭反应的反应速率的 方法有两种：
[0084] 第一种，通过实验方法获得。
[0085] 泡沫破灭的反应速率取决于泡沫的半衰期。实验方法具体方法包括：
[0086] 打开循环水浴锅并加热到80至100摄氏度之间，连接循环水浴锅和起泡仪；
[0087] 配制起泡剂溶液10-30毫升；
[0088] 向起泡仪注入起泡液；
[0089] 按照10至30毫升每分钟的速度注入氮气，气液比为1 : 1至1.5:1之间，观察普 通泡沫和凝胶泡沫的起泡情况，记录普通泡沫和凝胶泡沫的起泡体积及半衰期，重复测量 多次，取测量到的多个半衰期取平均值，得到普通泡沫的半衰期和凝胶泡沫的半衰期，根据 普通泡沫的半衰期计算得到普通泡沫破灭反应的反应速度，根据凝胶泡沫的半衰期计算得 到泡沫破灭反应的反应速度。
[0090] 使用上述方法产生的泡沫比较细小致密，这种结构有利于泡沫的稳定性，大而松 散的气泡排液速度快，液膜变薄处来不及修复，容易破灭。在实验中观察到，在0. 05%~ 〇. 4%的较宽浓度范围内，起泡剂的起泡体积随浓度增加呈现一个上升再到略有下降的平 缓趋势，这是由于在低于临界胶束浓度时，随着浓度的增加，注入一定氮气后，表面活性剂 分子在氮气/水界面上的吸附量增大，界面张力降低，产生的泡沫体系自由能降低， 能形成 更多的气液界面，产生更多的气泡，使泡沫体积增加，在临界胶束浓度时达到最高值，但是 在达到临界胶束浓度后，起泡性能不再提高，这主要是因为溶液表面张力不再减小，体系的 表面能趋于不变，另外形成的胶束在一定程度会增加液相的粘度，会使起泡体积受到一定 影响。
[0091] 第二种，采用优选范围（此优先范围通过多次实验确定出），普通泡沫破灭反应的 反应速率的值的优选范围为0. 01至0. 012,优选采用0. 011。凝胶泡沫破灭反应的反应速 率的值的优选范围为〇. 001至〇. 003,优先采用0. 002。
[0092] 本发明的基于凝胶泡沫的数值模拟方法，在引入凝胶泡沫的基础上建立了用于调 驱分析的数学模型，提出了有效表征凝胶泡沫调驱的数值模拟方法，对油藏进一步提高采 收率的数值模拟提供了有效的方法和途径。
[0093] 上面描述的内容可以单独地或者以各种方式组合起来实施，而这些变型方式都在 本发明的保护范围之内。
[0094] 需要说明的是，在本文中，术语"包括"、"包含"或者其任何其他变体意在涵盖非排 他性的包含，从而使得包括一系列要素的物品或者设备不仅包括那些要素，而且还包括没 有明确列出的其他要素，或者是还包括为这种物品或者设备所固有的要素。在没有更多限 制的情况下，由语句"包括……"限定的要素，并不排除在包括所述要素的物品或者设备中 还存在另外的相同要素。
[0095] 以上实施例仅用以说明本发明的技术方案而非限制，仅仅参照较佳实施例对本发 明进行了详细说明。本领域的普通技术人员应当理解，可以对本发明的技术方案进行修改 或者等同替换，而不脱离本发明技术方案的精神和范围，均应涵盖在本发明的权利要求范 围当中。
【主权项】
1. 一种基于凝胶泡沫的数值模拟方法，其特征在于，所述方法包括： 步骤1，构建油田或模拟油田实验装置的数学模型，将所述数学模型设置为3相9组分 形式，所述3相包括水相、油相和气相；所述9组分包括水、油、氮气、起泡剂、普通泡沫、聚合 物、交联剂、凝胶、凝胶泡沫；为各组分浓度设置初始值，为所述数学模型中的压力设置初始 值，为各组分粘度设置初始值，设置凝胶泡沫生成反应的反应速率和凝胶泡沫破灭反应的 反应速率，设置普通泡沫生成反应的反应速率和普通泡沫破灭反应的反应速率，设置油水 相对渗透率曲线以及气油相对渗透率曲线；设置时间步长和时间步数，以所述时间步长和 所述时间步数循环执行步骤2至步骤5 : 步骤2,求解水相压力方程，油相压力方程，气相压力方程； 步骤3,计算气相流速和凝胶泡沫液膜张力；其中，凝胶泡沫液膜张力为凝胶泡沫内的 气相压力与凝胶泡沫外壁的水相压力的差； 步骤4,比较所述气相流速与临界流速的大小，在比较结果为所述气相流速大于所述临 界流速时执行泡沫生成反应；并且，比较凝胶泡沫液膜张力与临界压力的大小，并在比较结 果为所述凝胶泡沫液膜张力大于所述临界压力时执行泡沫破灭反应；其中，所述泡沫生成 反应包括普通泡沫生成反应和凝胶泡沫生成反应；将所述普通泡沫生成反应设置为以下方 式：由水、氮气、起泡剂生成普通泡沫；将所述凝胶泡沫生成反应设置为以下方式：由水、氮 气、起泡剂、交联剂、聚合物生成普通泡沫；所述泡沫破灭反应包括普通泡沫破灭反应和凝 胶泡沫破灭反应；将所述普通破灭生成反应设置为以下方式：由普通泡沫分解为氮气、水、 起泡剂；将所述凝胶泡沫生成反应设置为以下方式：由凝胶泡沫分解为水、氮气、起泡剂、 凝胶； 步骤5,求解各组分浓度和各组分的粘度；修正凝胶泡沫的粘度和普通泡沫的粘度，修 正所述气油相对渗透率曲线。2. 如权利要求1所述的方法，其特征在于， 所述凝胶泡沫生成反应的反应速率的设置方法为：根据不同时间步的凝胶泡沫的粘度 的变化，设置凝胶泡沫的反应速率，具体为：根据实验数据获得相邻时间间隔相同的N个时 间点上的N个凝胶泡沫粘度，根据所述N个时间点与所述N个凝胶泡沫粘度，确定所述凝胶 泡沫生成反应的反应速率的值，使所述反应速率作为拟合函数的一参数使所述拟合函数在 时间点与凝胶泡沫粘度此两个坐标轴上能够拟合所述N个凝胶泡沫粘度，所述N为大于1 的整数。3. 如权利要求1所述的方法，其特征在于， 所述步骤1中设置气油相对渗透率曲线的过程中包括：将气相相对渗透率设置为原始 气相相对渗透率； 所述步骤5中还包括：通过气相相对渗透率修正系数R对气相相对渗透率数值进行修 正，所述步骤5中还包括：判断普通泡沫或凝胶泡沫破灭时，将气相相对渗透率设置为原始 气相相对渗透率。4. 如权利要求3所述的方法，其特征在于， 获取所述气相相对渗透率修正系数R的方法包括：通过室内实验测试得到不同渗透率 的凝胶泡沫残余阻力系数，根据所述凝胶泡沫残余阻力系数通过历史拟合方法确定所述气 相相对渗透率修正系数R。5. 如权利要求4所述的方法，其特征在于， 根据下述方法对所述气相相对渗透率进行修正： 满足pi. - 4 <丨< (S: + U或Cws彡C wslim时，根据下述公式修正：满足心 > 认:+ 4或C11J Cii5hm时，根据下述公式进行修正： ? ' I式中K为修正前的气相相对渗透率，为修正后的气相相对渗透率，< 为临界含水 饱和度，Cwslim为存在泡膜的最低起泡剂浓度，ε为饱和度小范围变化量，Sw为含水饱和度 也称为水相饱和度，Cws为起泡性剂浓度。6. 如权利要求3所述的方法，其特征在于， 所述气相相对渗透率修正系数R是80至106范围内的一个值。7. 如权利要求6所述的方法，其特征在于， 所述气相相对渗透率修正系数R的值为86。8. 如权利要求1所述的方法，其特征在于， 通过实验方法获得所述普通泡沫破灭反应的反应速率和所述凝胶泡沫破灭反应的反 应速率，具体包括： 打开循环水浴锅并加热到80至100摄氏度之间，连接循环水浴锅和起泡仪； 配制起泡剂溶液10-30毫升； 向起泡仪注入所述起泡液； 按照10至30毫升每分钟的速度注入氮气，气液比为I : 1至2:1之间，观察普通泡沫 和凝胶泡沫的起泡情况，记录普通泡沫和凝胶泡沫的起泡体积及半衰期，重复测量多次，取 测量到的多个半衰期取平均值，得到普通泡沫的半衰期和凝胶泡沫的半衰期，根据所述普 通泡沫的半衰期计算得到所述普通泡沫破灭反应的反应速度，根据所述凝胶泡沫的半衰期 计算得到所述泡沫破灭反应的反应速度。9. 如权利要求1所述的方法，其特征在于， 设置所述普通泡沫破灭反应的反应速率的值位于〇. 01至〇. 012范围内， 设置所述凝胶泡沫破灭反应的反应速率的值位于〇. 001至〇. 003范围内。10. 如权利要求9所述的方法，其特征在于， 设置所述普通泡沫破灭反应的反应速率的值为〇. 011， 设置所述凝胶泡沫破灭反应的反应速率的值为〇. 002。
【专利摘要】本发明公开了一种基于凝胶泡沫的数值模拟方法，本方法包括：构建3相9组分形式的数学模型；求解水相、油相、气相压力方程；计算气相流速和凝胶泡沫液膜张力；比较气相流速与临界流速的大小，在比较结果为气相流速大于临界流速时执行泡沫生成反应；比较凝胶泡沫液膜张力与临界压力的大小，并在比较结果为凝胶泡沫液膜张力大于临界压力时执行泡沫破灭反应；求解各组分浓度和各组分的粘度；修正凝胶泡沫的粘度和普通泡沫的粘度，修正气油相对渗透率曲线。本方法在引入凝胶泡沫的基础上建立了用于调驱分析的数学模型，提出了有效表征凝胶泡沫调驱的数值模拟方法，为油藏进一步提高采收率的数值模拟提供了有效手段。
【IPC分类】G06F19/00
【公开号】CN104899438
【申请号】CN201510295337
【发明人】王硕亮, 张媛, 于春磊, 付强, 李俊键, 刘同敬
【申请人】中国地质大学（北京）
【公开日】2015年9月9日
【申请日】2015年6月2日
最新回复(0)