致密油藏CO2吞吐有效作用半径的实验预测方法及系统与流程

本发明涉及非常规致密油藏开发实验领域，更具体地，涉及一种致密油藏co2吞吐有效作用半径的实验预测方法及系统。

背景技术：

目前，大部分国家将透率为1×10^-3μm²或3×10^-3μm²作为致密油藏储层渗透率的上限，一般最大不超过5×10^-3μm²。美国、俄罗斯、加拿大等国家均已成功开发致密油。由于致密油藏渗透率低、基质致密，致密油藏一般无自然产能，需要进行压裂改造才具备经济开采价值。目前，致密油藏开发主要采取水平井分段压裂来进行，初期产量较高，稳产期短，递减率大，但致密油藏的采收率亟待提高。

由于水平井分段压裂后，基质与裂缝接触面积较大，co2吞吐成为致密油藏提高采收率的首要选择。目前，美国bakken致密油藏等进行了前期co2吞吐机理研究，表明致密油藏co2吞吐具有较好效果。我国致密油藏开发尚处于起步阶段，对于致密油藏开发过程的认识尚存在不清楚之处。中石化华北分公司在鄂尔多斯盆地进行了小井组的co2吞吐试验，取得了一定的效果。对于致密油藏co2吞吐效果的评价一般采用物理模拟和数值模拟手段来进行，对于吞吐参数的模拟，一般采用单因素分析或正交实验的方法来进行优化。通过上述研究发现co2吞吐有效半径的预测能为致密油藏水平井分段压裂co2吞吐设计提供基础参数，但是，对于co2吞吐有效半径的预测，尚没有形成相应的方法。

因此，有必要开发一种致密油藏co2吞吐有效作用半径的实验预测方法。

公开于本发明背景技术部分的信息仅仅旨在加深对本发明的一般背景技术的理解，而不应当被视为承认或以任何形式暗示该信息构成已为本领域技术人员所公知的现有技术。

技术实现要素：

本发明提出了一种致密油藏co2吞吐有效作用半径的实验预测方法及系统，研发一种致密裂缝性油藏co2吞吐有效半径预测及实验方法，为致密油藏co2吞吐优化提供参考依据。

根据本发明的一方面，提出了一种致密油藏co2吞吐有效作用半径的实验预测方法，包括：

获取co2在饱和原油致密油藏中的扩散系数d；

基于油藏模型实验，获取致密油藏的岩心100％co2的浓度半径rco2；

基于所述油藏模型实验，确定焖井时间t与co2吞吐有效作用半径re的关系；

基于已知的扩散系数d、100％co2的浓度半径rco2及所述焖井时间t与co2吞吐有效作用半径re的关系，并根据数值微分计算方法，获取所述co2有效作用半径re。

优选地，所述焖井时间t对应的co2吞吐有效作用半径的关系为：

式中，rco2-100％co2的浓度半径；

d-扩散系数；

re-co2有效作用半径。

优选地，所述扩散系数d通过式(1)获取：

式中，j-扩散通量；

d-扩散系数；

c-浓度；

x-一维方向上的距离。

优选地，扩散量q随co2有效作用半径re变化时的变化量方程为：

式中，q-扩散量；

h-油层厚度。

优选地，通过式(6)进一步得到：

其中，基于数值微分计算方法对式(9)两边分别积分获得式(10)。

优选地，单位时间内，所述co2有效作用半径re时，扩散量q在所述co2有效作用半径re范围内co2浓度的方程式为：

式中，cr-距离井筒中心位置对应的co2浓度；

c-积分常数；

通过式(5)对re进行求导，获取式(6)。

优选地，在所述co2有效作用半径re范围内co2浓度分布方程为：

式中，cco2-co2有效作用半径re内的co2浓度；

ce-co2有效作用半径re的边缘浓度；

r-距离井筒中心的位置。

优选地，当co2有效作用半径re的co2浓度为0时，扩散速度方程为：

优选地，在径向流过程中，所述扩散半径为re时，扩散速度方程为：

式中，qco2-co2的扩散速度；

ce-co2有效作用半径re的边缘浓度；

当ce为0时，式(2)转化为式(3)。

根据本发明的另一方面，提出了一种致密油藏co2吞吐有效作用半径的实验预测系统，包括：

存储器，存储有计算机可执行指令；

处理器，所述处理器运行所述存储器中的计算机可执行指令，执行以下步骤：

获取co2在饱和原油致密油藏中的扩散系数d；

基于油藏模型实验，获取致密油藏的岩心100％co2的浓度半径rco2；

基于所述油藏模型实验，确定焖井时间t与co2吞吐有效作用半径re的关系；

基于已知的扩散系数d、100％co2的浓度半径rco2及所述焖井时间t与co2吞吐有效作用半径re的关系，并根据数值微分计算方法，获取所述co2有效作用半径re。

根据本发明的一种致密油藏co2吞吐有效作用半径的实验预测方法及系统，其优点在于：研发一种致密裂缝性油藏co2吞吐有效半径re预测及实验方法，为致密油藏co2吞吐优化提供参考依据。

本发明的方法和系统具有其它的特性和优点，这些特性和优点从并入本文中的附图和随后的具体实施例中将是显而易见的，或者将在并入本文中的附图和随后的具体实施例中进行详细陈述，这些附图和具体实施例共同用于解释本发明的特定原理。

附图说明

通过结合附图对本发明示例性实施例进行更详细的描述，本发明的上述以及其它目的、特征和优势将变得更加明显，其中，在本发明示例性实施例中，相同的参考标号通常代表相同部件。

图1示出了根据本发明的一个示例性实施例的一种致密油藏co2吞吐有效作用半径的实验预测方法的步骤的流程图。

图2示出了根据本发明的一个示例性实施例的co2扩散空间的示意图。

图3示出了根据本发明的一个示例性实施例的致密平板模型的示意图。

图4示出了根据本发明的一个示例性实施例的co2吞入过程中压力分布示意图。

图5示出了根据本发明的一个示例性实施例的室内实验co2焖井时间与扩散半径关系的示意图。

图6示出了根据本发明的一个示例性型实施例的矿场试验过程中co2焖井时间与扩散半径关系的示意图。

具体实施方式

下面将参照附图更详细地描述本发明。虽然附图中显示了本发明的优选实施例，然而应该理解，可以以各种形式实现本发明而不应被这里阐述的实施例所限制。相反，提供这些实施例是为了使本发明更加透彻和完整，并且能够将本发明的范围完整地传达给本领域的技术人员。

本实施例的一种致密油藏co2吞吐有效作用半径的实验预测方法，包括：

获取co2在饱和原油致密油藏中的扩散系数d；

基于油藏模型实验，获取致密油藏的岩心100％co2的浓度半径rco2；

基于油藏模型实验，确定焖井时间t与co2吞吐有效作用半径re的关系；

基于已知的扩散系数d、100％co2的浓度半径rco2及焖井时间t与co2吞吐有效作用半径re的关系，并根据数值微分计算方法，获取co2有效作用半径re。

当我们把注入的co2区域和地层看作为圆形时，渗流空间可以表征为图2所示情况，如图2所示，co2区域简化为一个半径为rco2的圆形区域，在这个区域内，为100％的co2存在，且在扩散过程中，该区域不发生变化，原油不发生流动；在扩散过程中，扩散的边缘为(即吞吐有效作用半径)re，边缘对应的co2浓度为0，re在扩散过程中，随时间的变化而扩大。

作为优选方案，所述扩散系数d通过式(1)获取：

式中，j-扩散通量；

d-扩散系数；

c-浓度；

x-一维方向上的距离。

在径向流过程中，通过推导可以获得当有效作用半径为re时，扩散速度方程为：

式中，qco2-co2的扩散速度；

ce-co2有效作用半径re的边缘浓度。

作为优选方案，当扩散边界co2浓度为0时，可以得到扩散速度方程为：

在re范围内，co2浓度分布公式为：

式中，cco2-co2有效作用半径re内的co2浓度；

ce-co2有效作用半径re的边缘浓度；

r-距离井筒中心的位置。

作为优选方案，在单位时间内，有效作用半径为re时，扩散量q就是re范围内co2总量，其表达式为：

式中，cr-距离井筒中心位置对应的co2浓度；

c-积分常数。

作为优选方案，上式对re进行求导，表示扩散量q随动态边界发生变化时的变化量方程：

式中，q-扩散量；

h-油层厚度。

则：

作为优选方案，对于式(9)两边分别积分得到焖井时间t与co2有效作用半径re的关系为：

式中，rco2-100％co2的浓度半径；

d-扩散系数；

re-co2有效作用半径。

式(10)为焖井时间与co2有效作用半径re的关系式，利用该式，在确定co2扩散系数、100％co2浓度半径后，根据数值微分计算方法，可以根据计算确定吞吐有效半径。

本发明研发一种致密裂缝性油藏co2吞吐有效半径re预测及实验方法，为致密油藏co2吞吐优化提供参考依据。

实施例

图1示出了根据本发明的一个示例性实施例的一种致密油藏co2吞吐有效作用半径的实验预测方法的步骤的流程图。

如图1所示，本发明的一种致密油藏co2吞吐有效作用半径的实验预测方法，包括：

获取co2在饱和原油致密油藏中的扩散系数d；

基于油藏模型实验，获取致密油藏的岩心100％co2的浓度半径rco2；

基于油藏模型实验，确定焖井时间t与co2吞吐有效作用半径re的关系；

基于已知的扩散系数d、100％co2的浓度半径rco2及焖井时间t与co2吞吐有效作用半径re的关系，并根据数值微分计算方法，获取co2有效作用半径re。

本实施例中，利用压力降落法原理，测量并确定co2在饱和原油致密岩心中的扩散系数为10^-9m²/s。

图3示出了根据本发明的一个示例性实施例的致密平板模型的示意图。

如图3所示，为了确定co2吞吐过程中100％co2浓度半径，选用致密平板二维模型进行实验，该模型尺寸为40cm×40cm×2.7cm，渗透率为4.37md，孔隙度10.15％，根据现场井网形态，按照与现场裂缝距离100米、裂缝半缝长100米相似的水平井及裂缝的模式。

图4示出了根据本发明的一个示例性实施例的co2吞入过程中压力分布示意图。

如图4所示，根据现场温压条件，进行65℃、20mpa条件下的co2吞吐实验，图4表示co2吞入过程中的压力分布情况，根据实验测定，co2注入初期100％co2区域的宽度约为3cm。

根据焖井时间t与吞吐有效作用半径re关系式，运用runge-kutta法进行数值微分计算，得到室内实验过程中扩散半径与焖井时间之间的关系。

表1焖井时间与扩散半径之间的关系

图5示出了根据本发明的一个示例性实施例的室内实验co2焖井时间与扩散半径关系的示意图。

根据表1及实验过程绘制焖井时间与扩散半径关系，如图5所示，按照径向扩散模型计算co2吞吐过程中的有效扩散pv数。按照实验模型与现场模拟单元相似的原则，确定现场模拟单元为200m╳200m╳15m，裂缝长度100m，裂缝间距100m，油层厚度为15米。根据有效扩散pv相等，计算不同焖井时间对应的扩散半径。

图6示出了根据本发明的一个示例性型实施例的矿场试验过程中co2焖井时间与扩散半径关系的示意图。

从图6可知，随着焖井时间的延长，co2有效扩散半径逐步增大，当焖井时间超过400小时时，即16天左右时，co2扩散半径增加幅度降低。当焖井时间为720小时，即30天时，co2有效扩散半径为6.5米。

以上已经描述了本发明的实施例，上述说明是示例性的，并非穷尽性的，并且也不限于所披露的实施例。在不偏离所说明的实施例的范围和精神的情况下，对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。本文中所用术语的选择，旨在最好地解释实施例的原理、实际应用或对市场中的技术的改进，或者使本技术领域的其它普通技术人员能理解本文披露的实施例。