二氧化碳驱压裂数模一体化方法技术

本发明专利技术提供一种二氧化碳驱压裂数模一体化方法，该二氧化碳驱压裂数模一体化方法包括：步骤1，根据油藏实际情况建立地质模型并进行网格剖分；步骤2，进行水力裂缝扩展数值模拟，捕捉裂缝扩展轨迹；步骤3，根据裂缝扩展轨迹及天然裂缝分布建立嵌入式离散裂缝模型；步骤4，耦合嵌入式离散裂缝模型和组分模型，得到压裂油藏二氧化碳驱油模型并进行数值模拟。该二氧化碳驱压裂数模一体化方法将水力压裂模拟和CO2驱油模拟结合在一起，实现CO2驱压裂数模一体化，更加快捷、准确地进行油藏数值模拟。可以为CO2驱油数值模拟技术提供新的方法，并可以为其他油藏数值模拟提供思路，具有较大意义。义。义。

【技术实现步骤摘要】
二氧化碳驱压裂数模一体化方法


[0001]本专利技术涉及水力压裂及二氧化碳驱油模拟
，特别是涉及到一种二氧化碳驱压裂数模一体化方法。

技术介绍

[0002]我国已探明低渗透地质储量逾140亿吨。低渗透油藏水驱动用困难，已开发低滲油田水驱开均采收率仅20％左右，低渗透油藏水驱开发存在注入性差、开采效果较差等问题。目前，低渗透油藏可采储量采出程度过半，亟需开发接替技术。
[0003]低渗透油气资源往往需要水力压裂后才具有较好的采收效果。水力压裂后可以提高油井的控油半径，而且注水井在高速注水压裂后还能补充地层能量，对于提高原油产量大有益处。
[0004]与此同时，低渗透油藏二氧化碳驱油技术也引起了广泛关注，被应用于国内外实际油田的开采中，取得了很好的效果。二氧化碳是一种优越的驱油剂，注二氧化碳驱油已被证实是一种有效的提高采收率的方法，二氧化碳拥有粘度低、易注入的特点，而且注入地层和原油接触后能够使原油体积膨胀，降低原油粘度，并抽提轻质组分，降低油水界面张力，改善流度比等，可以有效提高原油采收率。二氧化碳与原油混相后，能消除界面张力，大幅提高采收率。二氧化碳捕集驱油与埋存及提高采收率技术(CCS
‑
EOR)，通过二氧化碳捕集，注入地层驱油，回收再利用等技术，将二氧化碳永久封存于地下，可以有效降低温室气体的含量，实现增加原油产量与埋存的“双赢”，是目前最现实的碳减排和利用途径。开展二氧化碳驱油技术应用则是大型石油企业开展碳减排工作、践行国家绿色低碳发展理念的重要抓手。所以，发展二氧化碳驱技术契合国家低碳发展战略。
[0005]因此，水力压裂和二氧化碳驱油技术的结合，是一种低渗、致密油藏开采的有效技术。但是，目前该技术的数值模拟存在困难。
[0006]二氧化碳数值模拟技术经过发展已经能较好的模拟二氧化碳驱油过程，但是在模拟压裂油藏中二氧化碳驱油过程时，由于裂缝开度较小，采用常规离散裂缝模型表征裂缝会导致网格剖分复杂、计算量大，且难以与水力压裂技术结合。因此目前二氧化碳驱中压裂裂缝的表征尚存在较大的改进空间。模拟水力压裂过程也是比较复杂的，需要考虑岩石力学变形、流体在压裂裂缝中的流动、断裂扩展准则、水力裂缝与天然裂缝的相互作用等方面。国内外专家也在此方面做了大量研究，包括岩石力学变形的弹性理论，裂缝中流体的滤失理论以及断裂扩展模型。在此基础上，近年来发展了一些模型模拟压裂裂缝的形成。
[0007]但是，目前两种技术是分开、割裂的。在水力压裂模拟结束后，需要将裂缝扩展信息收集，然后重新建模，在进行二氧化碳驱油数值模拟。一方面，导致了较大的工作量，另一方面，两种工作的割裂会导致计算时出现误差。
[0008]在申请号：CN202010733943.0的中国专利申请中，涉及到一种陆相致密油藏强化体积改造缝网压裂的方法，该陆相致密油藏强化体积改造缝网压裂的方法包括：步骤1、建立含单压裂段的致密油藏全三维水平井压裂地质力学模型，进行水平井密切割优化设计强
化改造；步骤2、进行油藏增能渗吸强化改造；步骤3、建立致密油藏区块全三维水平井压裂地质力学模型，进行体积缝网压裂强化改造；步骤4、进行主裂缝簇式支撑高导流压裂强化改造。该陆相致密油藏强化体积改造缝网压裂的方法在体积压裂技术的基础上发展了强化体积改造缝网压裂技术，即在压裂工艺的各个环节充分考虑提高油藏初始产量和最终采收率的工艺手段，实现有效储层改造体积最大化。
[0009]在申请号：CN201611093557.X的中国专利申请中，涉及到一种提高CO2封存量和原油采收率的方法，包括：根据长细管驱油实验得到压力水平与驱油效率关系曲线；利用CO2驱油与封存靶区盖层和断层稳定性评价方法，确定靶区压力安全界限；利用油藏数值模拟方法，优选合理的压力保持水平；根据储层参数、生产参数和工艺方法，实现CO2纵向均衡驱替；根据靶区储层非均质性和剩余储量丰度分布情况，采用变井网井距、差异储层改造，实现CO2平面均衡驱替；通过对注入方式、注入速度、注采调控和泡沫调堵这些方式优选，调整注采流线。该专利技术为提高CO2封存量和原油采收率提供了可行方法，在节能减排和CO2资源化利用、提高原油采收率方面有很大应用前景。
[0010]在申请号：CN202010853616.9的中国专利申请中，涉及到一种页岩油藏开发方案的优选方法，包括以下步骤：根据目标油藏的地质模型，选择甜点区建立开发机理模型；根据所述开发机理模型，设计多种提高采收率的开发方式，将各种开发方式进行组合，形成多套整体开发方案，并通过数值模拟对各套整体开发方案进行采收率指标预测；预测未来多年国际油价，根据所述油价与所述采收率计算各套整体开发方案的累计净现值，选择累计净现值最大的整体开发方案作为所述目标油藏的优选开发方案。本专利技术优选的开发方案其经济效益更高，能够降低页岩油气开发的经济损失风险。
[0011]以上现有技术均与本专利技术有较大区别，未能解决我们想要解决的技术问题，为此我们专利技术了一种新的二氧化碳驱压裂数模一体化方法。

技术实现思路

[0012]本专利技术的目的是提供一种避免了水力压裂与二氧化碳驱分开建模的问题，能够高效、准确地进行二氧化碳驱油数值模拟的二氧化碳驱压裂数模一体化方法。
[0013]本专利技术的目的可通过如下技术措施来实现：二氧化碳驱压裂数模一体化方法，该二氧化碳驱压裂数模一体化方法包括：
[0014]步骤1，根据油藏实际情况建立地质模型并进行网格剖分；
[0015]步骤2，进行水力裂缝扩展数值模拟，捕捉裂缝扩展轨迹；
[0016]步骤3，根据裂缝扩展轨迹及天然裂缝分布建立嵌入式离散裂缝模型；
[0017]步骤4，耦合嵌入式离散裂缝模型和组分模型，得到压裂油藏二氧化碳驱油模型并进行数值模拟。
[0018]本专利技术的目的还可通过如下技术措施来实现：
[0019]在步骤1中，根据油藏的实际情况，收集油藏地质资料，建立油藏模型；根据模拟需要划分网格系统，计算网格系统的几何信息。
[0020]在步骤1中，收集的油藏地质资料包括天然裂缝分布、油藏几何形态、渗透率、孔隙度以及流体资料。
[0021]在步骤1中，计算的网格系统的几何信息，包括网格体积、网格与天然裂缝的相交
信息。
[0022]在步骤2中，水力裂缝扩展模拟需要考虑三个要素：流体在裂缝内的流动模拟、明确裂缝扩展准则以及岩体在流体压力和地应力的共同作用下变形；综合三个要素，采用扩展有限元法模拟裂缝扩展，获得裂缝扩展的动态数据，并捕捉裂缝轨迹信息。
[0023]在步骤2中，流体在裂缝内的流动模拟时，由于裂缝宽度较小，因此忽略沿裂缝开度方向的流动，仅考虑沿裂缝方向的流动：
[0024][0025]其中，q为单位厚度截面流量；w为裂缝面宽宽度；μ为压裂流体黏度；x为沿裂缝裤扩展方向的位移；p表示流体压力。
[0026]边界条件包括裂缝注入速率及尖端流量为0，表示为
[0027][0028]其中Q0为裂缝的注入速度；L为本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.二氧化碳驱压裂数模一体化方法，其特征在于，该二氧化碳驱压裂数模一体化方法包括：步骤1，根据油藏实际情况建立地质模型并进行网格剖分；步骤2，进行水力裂缝扩展数值模拟，捕捉裂缝扩展轨迹；步骤3，根据裂缝扩展轨迹及天然裂缝分布建立嵌入式离散裂缝模型；步骤4，耦合嵌入式离散裂缝模型和组分模型，得到压裂油藏二氧化碳驱油模型并进行数值模拟。2.根据权利要求1所述的二氧化碳驱压裂数模一体化方法，其特征在于，在步骤1中，根据油藏的实际情况，收集油藏地质资料，建立油藏模型；根据模拟需要划分网格系统，计算网格系统的几何信息。3.根据权利要求2所述的二氧化碳驱压裂数模一体化方法，其特征在于，在步骤1中，收集的油藏地质资料包括天然裂缝分布、油藏几何形态、渗透率、孔隙度以及流体资料。4.根据权利要求2所述的二氧化碳驱压裂数模一体化方法，其特征在于，在步骤1中，计算的网格系统的几何信息，包括网格体积、网格与天然裂缝的相交信息。5.根据权利要求1所述的二氧化碳驱压裂数模一体化方法，其特征在于，在步骤2中，水力裂缝扩展模拟需要考虑三个要素：流体在裂缝内的流动模拟、明确裂缝扩展准则以及岩体在流体压力和地应力的共同作用下变形；综合三个要素，采用扩展有限元法模拟裂缝扩展，获得裂缝扩展的动态数据，并捕捉裂缝轨迹信息。6.根据权利要求5所述的二氧化碳驱压裂数模一体化方法，其特征在于，在步骤2中，流体在裂缝内的流动模拟时，由于裂缝宽度较小，因此忽略沿裂缝开度方向的流动，仅考虑沿裂缝方向的流动：其中，q为单位厚度截面流量；w为裂缝面宽宽度；μ为压裂流体黏度；x为沿裂缝裤扩展方向的位移；p表示流体压力；边界条件包括裂缝注入速率及尖端流量为0，表示为其中Q0为裂缝的注入速度；L为裂缝半长。7.根据权利要求5所述的二氧化碳驱压裂数模一体化方法，其特征在于，在步骤2中，明确裂缝扩展准则时，采用最大周向应力准则，当裂缝尖端的应力强度因子大于岩石韧度时，裂缝将进行扩展，并且沿着最大周向应力方向扩展，根据最大周向应力准则，必须满足：其中O
M
为最大周向应力方向角，即裂缝开裂方向与裂缝面的夹角；K
I
为第一类应力强度因子；K
II
为第一类应力强度因子；K
IC
表示岩石断裂韧度。8.根据权利要求5所述的二氧化碳驱压裂数模一体化方法，其特征在于，在步骤2中，岩体在流体压力和地应力的共同作用下变形时，岩体在裂缝内流体压力和地应力的综合作用
下发生变形，岩石变形控制方程组为σ＝C:ε(u)其中边界条件包括外边界和裂缝边界条件，即其中σ为应力张量；f为体积力；ε为应变张量；C为弹性张量矩阵；u为位移；g为定位移；Γ
c
为位移外边界；Γ
c
为裂缝边界；σ
s
为切应力；为定切应力；σ
n
为正应力；p为流体压力。9.根据权利要求5所述的二氧化碳驱压裂数模一体化方法，其特征在于，在步骤2中，采用扩展有限元法模拟裂缝扩展时，考虑应力场和流动场的耦合；对于应力场，水力压裂裂缝的扩展是个动态的过程，如果采用传统有限元方法求解，网格单元需要不断重构，导致网格划分及计算过程复杂且容易出现计算误差；扩展有限元法在传统有限元的基础上增加自由度，避免了计算过程中的网格加密和网格重构；将位移分为连续位移和不连续位移，位移近似公式为其中，I为单元内位移节点集；N
i
为常规有限元单元节点位移形函数；u
i
为单元节点位移矢量；I
enr
和I
tip
分别为裂缝段和裂...

【专利技术属性】
技术研发人员：杨勇，张庆福，张传宝，李友全，王杰，王丽娜，郭祥，韩凤蕊，张东，于伟杰，
申请(专利权)人：中国石油化工股份有限公司胜利油田分公司勘探开发研究院，
类型：发明
国别省市：