CN114154430A - 一种压裂油藏co2驱油流动模拟方法 - Google Patents

一种压裂油藏co2驱油流动模拟方法 Download PDF

Info

Publication number
CN114154430A
CN114154430A CN202010926532.3A CN202010926532A CN114154430A CN 114154430 A CN114154430 A CN 114154430A CN 202010926532 A CN202010926532 A CN 202010926532A CN 114154430 A CN114154430 A CN 114154430A
Authority
CN
China
Prior art keywords
model
oil
fracture
phase
fracturing
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202010926532.3A
Other languages
English (en)
Inventor
杨勇
张世明
张庆福
张传宝
李友全
郭祥
王晓敏
王杰
王丽娜
韩凤蕊
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Petroleum and Chemical Corp
Exploration and Development Research Institute of Sinopec Shengli Oilfield Co
Original Assignee
China Petroleum and Chemical Corp
Exploration and Development Research Institute of Sinopec Shengli Oilfield Co
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China Petroleum and Chemical Corp, Exploration and Development Research Institute of Sinopec Shengli Oilfield Co filed Critical China Petroleum and Chemical Corp
Priority to CN202010926532.3A priority Critical patent/CN114154430A/zh
Publication of CN114154430A publication Critical patent/CN114154430A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • E—FIXED CONSTRUCTIONS
    • E21—EARTH OR ROCK DRILLING; MINING
    • E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B43/00—Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
    • E21B43/16—Enhanced recovery methods for obtaining hydrocarbons
    • E21B43/164—Injecting CO2 or carbonated water
    • E—FIXED CONSTRUCTIONS
    • E21—EARTH OR ROCK DRILLING; MINING
    • E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B43/00—Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
    • E21B43/25—Methods for stimulating production
    • E21B43/26—Methods for stimulating production by forming crevices or fractures
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00—Details relating to the application field
    • G06F2113/08—Fluids
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Geology (AREA)
  • Mining & Mineral Resources (AREA)
  • Fluid Mechanics (AREA)
  • General Physics & Mathematics (AREA)
  • Environmental & Geological Engineering (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Analysis (AREA)
  • Computing Systems (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Physics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Algebra (AREA)
  • Chemical & Material Sciences (AREA)
  • Chemical Kinetics & Catalysis (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明涉及低渗透油藏CO2驱提高采收率数值模拟的领域，特别涉及一种压裂油藏CO2驱油流动模拟方法。所述方法包括以下步骤：确定压裂产生的人工大裂缝位置及微裂缝区域，建立压裂油藏几何模型；建立混合模型表征不同尺度裂缝，进行网格剖分，形成网格***；建立组分模型描述CO2驱油过程，对其中的组分予以考虑，描述相态变化及CO2和原油的混相机理；制定求解方案，求解复杂的压裂油藏组分模型，使得模拟过程稳定。本发明方法在模拟CO2驱油流动过程的同时，既能准确描述人工大裂缝，也能描述压裂或注气过程中产生的小裂缝，克服了以上所述现有的压裂油藏CO2驱油方法面临的困难。

Description

一种压裂油藏CO2驱油流动模拟方法
技术领域
本发明涉及低渗透油藏CO2驱提高采收率数值模拟的领域，特别涉及一种压裂油藏CO2驱油流动模拟方法。
背景技术
CO2驱油是最具优势和潜力的低渗透油藏提高石油采收率技术，CO2驱油技术的发展对有效开采低渗透油藏，确保我国石油产量的稳定增长，维护我国石油安全具有重要意义。
然而在实施CO2驱油的过程中会难以避免的遇到裂缝。首先，裂缝作为最小的地质构造在地壳中分布广泛，几乎所有的油藏中都存在不同规模的裂缝。据不完全统计，世界上裂缝性油藏约占已探明总储量的一半；在我国，已探明的裂缝性油藏地质储量占总探明储量的28％以上。而且，在中低渗透油藏的开采过程中，为提高采油井的采油能力或注入井的注水(注气)能力，往往会进行人工压裂，人工压裂会产生不同尺度的压裂裂缝。一方面，在压裂过程中，除了大尺度压裂裂缝，也会产生众多小尺度的裂缝；而且在注水注气的过程中，注入压力偏大时也会产生众多小尺度的裂缝。压裂产生的裂缝会显著影响CO2驱油过程中的压力分布，组分分布及饱和度分布等。因此在进行压裂时，需要考虑压裂规模，压裂裂缝的尺度，数量，方位等对CO2驱油效果的影响。
对于裂缝对CO2驱油影响的探索，油藏流动的模拟成为探索油藏的重要方式，目前，主要有基于双重介质、等效连续介质和离散裂缝三种流动数学模型，双重介质模型认为介质中存在两个平行的渗流***：裂缝***为主要流动通道，基岩***为储集空间，模型在一定程度上刻画出了裂缝优先流特点；等效连续介质模型将整个介质视为一个连续***，通过等效参数来表征其非均质性；离散裂缝模型则对每条裂缝予以显示表征和模拟，计算精度高、拟真性好。
但是基于双重介质模型假设裂缝***呈均匀网络状分布，不能真实地反映裂缝的随机性和多尺度性，不适用于裂缝性油藏的精细流动模拟，等效连续介质模型的理论基础为尺度升级理论，本质上是对数学模型进行降阶处理，升级后的大尺度模型人为地抹平了介质的强非均质性和多尺度性，过于宏观简化，不能捕捉小尺度流动特征。离散裂缝模型则对每条裂缝予以显示表征和模拟，计算精度高、拟真性好，但网格剖分复杂。
对于人工压裂产生的主要的大裂缝，裂缝尺度大，且裂缝相对是孤立，此时双重介质膜模型不能准确反映裂缝对CO2驱油的影响，因此采用离散裂缝模型予以表征。然而基于非结构化网格剖分的离散裂缝模型在裂缝较多时存在计算量大的问题，难以推广至实际应用。对于压裂或注气产生的众多微裂缝，如果采用离散裂缝模型表征，会导致表征困难，计算量大。双重介质模型对于这种微裂缝具有较好的模拟效果，但是双重介质模型只考虑相互裂缝中的流动及基岩块与周围裂缝之间流动，忽略了基岩内部的流动过程。对于CO2驱油这种多相多组分的流动过程，该模型的描述不够精确。因此，需要寻求一种能够同时准确高效的表征人工大裂缝和微裂缝，同时便于和组分模型耦合的模型，构建压裂油藏CO2驱油模拟方法。
发明内容
本发明主要目的是提供一种压裂油藏CO2驱油流动模拟方法，本发明方法在模拟CO2驱油流动过程的同时，既能准确描述人工大裂缝，也能描述压裂或注气过程中产生的小裂缝，克服了以上所述现有的压裂油藏CO2驱油方法面临的困难。
为实现上述目的，本发明采用以下技术方案：
本发明提供一种压裂油藏CO2驱油流动模拟方法，所述方法包括以下步骤：
确定压裂产生的人工大裂缝位置及微裂缝区域，建立压裂油藏几何模型；
建立混合模型表征不同尺度裂缝，进行网格剖分，形成网格***；
建立组分模型描述CO2驱油过程，对其中的组分予以考虑，描述相态变化及CO2和原油的混相机理；
制定求解方案，求解复杂的压裂油藏组分模型，使得模拟过程稳定。
进一步优选地，所述建立混合模型的具体步骤包括：
对人工压裂产生的大裂缝建立嵌入式离散裂缝模型；
对压裂或注气过程中产生的微尺度裂缝建立MINC模型；
基于守恒定律，将基岩***、嵌入式离散裂缝***和MINC模型耦合起来，建立混合模型，将该混合模型作为人工压裂油藏流动模型，以此为基础建立流动方程。
进一步优选地，对压裂或注气过程中产生的微尺度裂缝建立MINC模型的方法包括：
首先建立双重介质模型：
Figure BDA0002666135390000031
Figure BDA0002666135390000032
F＝αmKmτm(pm-pf)pm≥pf
F＝αmKmτf(pm-pf)pm＜pf
其中，n＝1或2，分别代表双重介质模型里面的基岩***和裂缝***，vn***内流体流动速度，F为两个***间的窜流量，qn为***内源汇项，τn代表***内流体总体流度；Kn代表***绝对渗透率张量，pn为***内流体压力，αm为基岩***的形状因子，Km为基岩***的绝对渗透率张量，τm代表基岩***内流体流度系数，pm代表基岩压力，pf代表裂缝压力，τf代表裂缝***内流体流度系数；
其次，在双重介质模型的基础上将基岩单元进一步划分为几个嵌套的单元，构成MINC模型。
进一步优选地，对人工压裂产生的大裂缝建立嵌入式离散裂缝模型的方法包括：
对基岩进行剖分，然后将裂缝***嵌入基岩网格当中，
嵌入式离散裂缝模型是直接将裂缝***嵌入油藏网格当中，在进行网格剖分后，将裂缝嵌入网格中，然后计算裂缝***与网格***的相交关系，嵌入式离散模型的质量守恒方程为：
Figure BDA0002666135390000041
其中，kF为大裂缝渗透率，pF为大裂缝压力，μ为流体粘度，qF为大裂缝源汇项，qFf为两个***间的窜流量，qFF表示相交裂缝单元之间的窜流量，δFF表示大裂缝单元是否与其他裂缝单元相交，若大裂缝单元与其他裂缝单元相交δFF＝1，否则δFF＝0，VF分别是裂缝单元的体积。
进一步优选地，网格剖分的步骤具体包括：首先对于油藏整体，若裂缝为嵌入式，不需要作为基岩网格内边界，可直接对油藏进行网格剖分，然后嵌入裂缝***；
MINC模型的网格剖分主要就是对于基岩块的剖分，将基岩块划分为嵌套网格，然后采用正交网格剖分基岩单元；
由油藏基岩的网格***、MINC模型的网格单元和裂缝***，构成整体的网格***。
进一步优选地，计算嵌入式离散裂缝模型的参数方法：对于人工压裂的大裂缝需要确定其裂缝空间坐标，条数，间距等，然后在此基础上计算每条裂缝所穿过的网格单元以及裂缝与网格单元的连接信息，包括裂缝与网格的接触面积，裂缝所在网格到裂缝面的平均距离等信息，为后续进行混合模型内的组分流动模拟做准备。
进一步优选地，计算MINC模型参数方法：根据油田实际资料确定MINC区域的大小，然后计算MINC模型中每个子单元的体积分数，相邻关系，各子单元空间坐标，子单元体积等参数。
优选地，建立组分模型描述CO2驱油过程：
(1)通过实验结果确定组分数量，针对每一个组分建立质量守恒方程
Fi+＝qi＝θNi/θt
其中i表示组分，i＝1,…,nc,nc表示组分数量。F为质量流量，q代表源汇项，N为质量变化；N可以通过相摩尔密度ρ，饱和度S，油相中组分i的摩尔分数xi，气相中组分i的摩尔分数yi表示为
Figure BDA0002666135390000052
F可以写为
Figure BDA0002666135390000053
其中vβ为相β的达西速度
Figure BDA0002666135390000054
其中β表示油相，气相；
(2)各组分要满足组分质量守恒条件：
Figure BDA0002666135390000055
其中
Figure BDA0002666135390000061
表示逸度系数；逸度系数可以通过压力、温度和体积等数据计算；而温度T，压力P，体积V参数可以通过PR状态方程计算；
在求解的过程通过使方程满足质量守恒条件和热动力平衡条件，更新平衡比，最终得到平衡条件下的组分含量等参数；在计算中将裂缝也视为单元，可将组分模型、闪蒸计算在计算中适用于裂缝***。
优选地，所述模拟CO2驱油混相机理的具体步骤包括：当达到一定压力时，原油与CO2将发生混相，此时油气***相渗曲线将会发生改变，此时采用界面张力矫正法修正相渗曲线
Figure BDA0002666135390000062
Figure BDA0002666135390000063
其中：
Figure BDA0002666135390000064
其中，
Figure BDA0002666135390000065
和
Figure BDA0002666135390000066
为非混相下的油气相对渗透率；是krg气相渗透率，krog是油相相对渗透率，kroiw是束缚水饱和度下油或气的相对渗透率，n，no，ng是无因次可调参数，σ*是非混相下的油气界面张力值；
当原油与CO2逐渐混相时，气相和油相相对渗透率发生改变，达到完全混相时，界面张力为变为0，此时相渗曲线会变为两条直线。
优选地，所述求解复杂的压裂油藏组分模型的具体步骤包括：耦合组分模型和裂缝流动模型的压裂油藏组分模型，为了使气驱模拟过程收敛，采用隐式求解方案；为降低矩阵维度，将CO2驱油数学模型中的方程分为主方程和辅助方程，将变量分为主变量和辅助变量，数值求解过程中，先通过求解主方程得到主变量，然后通过辅助方程获得辅助变量。
与现有技术相比，本发明具有以下有益效果：
发明公开了一种压裂油藏CO2驱油流动模拟方法，针对低渗透油藏实施CO2驱油过程中遇到人工压裂裂缝以及众多微裂缝的情况，采用采用嵌入式离散裂缝模型表征大裂缝，采用MINC模型表征微裂缝，同时用组分模型模拟CO2驱油过程中组分相态的变化以及混相等机理。将上述模型耦合起来，描述压裂油藏中CO2驱油过程及机理。模型考虑大裂缝、微裂缝、CO2与原油混相机理及低渗透油藏特征的影响，更加符合实际。本发明方法在模拟CO2驱油流动过程的同时，既能准确描述人工大裂缝，也能描述压裂或注气过程中产生的小裂缝，克服了以上所述现有的压裂油藏CO2驱油方法面临的困难。
本发明方法为实施油田实施CO2驱油提供了一种高效的模拟方法，能够研究不同措施下CO2驱油效果，为制定压裂方案提供建议，实现油气藏的开发动态模拟，为油气藏的高效开发提供技术支撑。
附图说明
构成本发明的一部分的说明书附图用来提供对本发明的进一步理解，本发明的示意性实施例及其说明用于解释本发明，并不构成对本发明的不当限定。
图1为本发明所述压裂油藏CO2驱流动模拟方法一具体实施例的流程图；
图2为本发明一具体实施例所述表征微裂缝的MINC示意图；
图3为本发明一具体实施例所述嵌入式离散裂缝模型网格示意图。
具体实施方式
应该指出，以下详细说明都是示例性的，旨在对本发明提供进一步的说明。除非另有指明，本文使用的所有技术和科学术语具有与本发明所属技术领域的普通技术人员通常理解的相同含义。
需要注意的是，这里所使用的术语仅是为了描述具体实施方式，而非意图限制根据本发明的示例性实施方式。如在这里所使用的，除非上下文另外明确指出，否则单数形式也意图包括复数形式，此外，还应当理解的是，当在本说明书中使用术语“包含”和/或“包括”时，其指明存在特征、步骤、操作和/或它们的组合。
为了使得本领域技术人员能够更加清楚地了解本发明的技术方案，以下将结合具体的实施例详细说明本发明的技术方案。
如图1所示，所述一种压裂油藏CO2驱油流动模拟方法，包括以下步骤：
S1.根据人工压裂的实际情况，确定压裂产生的人工大裂缝位置，尺度，数量以及分布等情况，确定微裂缝区域，建立压裂油藏几何模型；
S2.根据裂缝分布情况，建立混合模型表征不同尺度裂缝：
所述建立混合模型的具体步骤包括：
对压裂或注气过程中产生的微尺度裂缝建立MINC模型，如图2所示；
对人工压裂产生的大裂缝建立嵌入式离散裂缝模型，如图3所示；
基于守恒定律，将基岩***、嵌入式离散裂缝***和MINC模型耦合起来，建立混合模型，将该混合模型作为人工压裂油藏流动模型，以此为基础建立流动方程。
对压裂或注气过程中产生的微尺度裂缝建立MINC模型的方法包括：
首先建立双重介质模型：
Figure BDA0002666135390000081
Figure BDA0002666135390000082
F＝αmKmτm(pm-pf)pm≥pf
F＝αmKmτf(pm-pf)pm＜pf
其中，n＝1或2，分别代表双重介质模型里面的基岩***和裂缝***，vn***内流体流动速度，F为两个***间的窜流量，qn为***内源汇项，τn代表***内流体总体流度；Kn代表***绝对渗透率张量，pn为***内流体压力，αm为基岩***的形状因子，Km为基岩***的绝对渗透率张量，τm代表基岩***内流体流度系数，pm代表基岩压力，pf代表裂缝压力，τf代表裂缝***内流体流度系数；
其次，在双重介质模型的基础上将基岩单元进一步划分为几个嵌套的单元，构成MINC模型。
计算MINC模型参数方法：根据油田实际资料确定MINC区域的大小，然后计算MINC模型中每个子单元的体积分数，相邻关系，各子单元空间坐标，子单元体积等参数。
对人工压裂产生的大裂缝建立嵌入式离散裂缝模型的方法包括：
对基岩进行剖分，然后将裂缝***嵌入基岩网格当中，
嵌入式离散裂缝模型是直接将裂缝***嵌入油藏网格当中，在进行网格剖分后，将裂缝嵌入网格中，然后计算裂缝***与网格***的相交关系，嵌入式离散模型的质量守恒方程为：
Figure BDA0002666135390000091
其中，kF为大裂缝渗透率，pF为大裂缝压力，μ为流体粘度，qF为大裂缝源汇项，qFf为两个***间的窜流量，qFF表示相交裂缝单元之间的窜流量，δFF表示大裂缝单元是否与其他裂缝单元相交，若大裂缝单元与其他裂缝单元相交δFF＝1，否则δFF＝0，VF分别是裂缝单元的体积。
计算嵌入式离散裂缝模型的参数方法：对于人工压裂的大裂缝需要确定其裂缝空间坐标，条数，间距等，然后在此基础上计算每条裂缝所穿过的网格单元以及裂缝与网格单元的连接信息，包括裂缝与网格的接触面积，裂缝所在网格到裂缝面的平均距离等信息，为后续进行混合模型内的组分流动模拟做准备。
S3.建立组分模型描述CO2驱油过程，对其中的组分予以考虑，描述相态变化及CO2和原油的混相机理。建立组分模型描述CO2驱油过程：
(1)通过实验结果确定组分数量，针对每一个组分建立质量守恒方程
Fi+＝qi＝θNi/θt
其中i表示组分，i＝1,…,nc,nc表示组分数量。F为质量流量，q代表源汇项，N为质量变化；N可以通过相摩尔密度ρ，饱和度S，油相中组分i的摩尔分数xi，气相中组分i的摩尔分数yi表示为
Figure BDA0002666135390000102
F可以写为：
Figure BDA0002666135390000103
其中vβ为相β的达西速度
Figure BDA0002666135390000104
其中β表示油相，气相；
(2)各组分要满足组分质量守恒条件：
Figure BDA0002666135390000105
其中
Figure BDA0002666135390000106
表示逸度系数；逸度系数可以通过压力、温度和体积等数据计算；而温度T，压力P，体积V参数可以通过PR状态方程计算；
在求解的过程通过使方程满足质量守恒条件和热动力平衡条件，更新平衡比，最终得到平衡条件下的组分含量等参数；在计算中将裂缝也视为单元，可将组分模型、闪蒸计算在计算中适用于裂缝***。
所述模拟CO2驱油混相机理的具体步骤包括：当达到一定压力时，原油与CO2将发生混相，此时油气***相渗曲线将会发生改变，此时采用界面张力矫正法修正相渗曲线：
Figure BDA0002666135390000111
Figure BDA0002666135390000112
其中：
Figure BDA0002666135390000113
其中，
Figure BDA0002666135390000114
和
Figure BDA0002666135390000115
为非混相下的油气相对渗透率；是krg气相渗透率，krog是油相相对渗透率，kroiw是束缚水饱和度下油或气的相对渗透率，n，no，ng是无因次可调参数，σ*是非混相下的油气界面张力值；
当原油与CO2逐渐混相时，气相和油相相对渗透率发生改变，达到完全混相时，界面张力为变为0，此时相渗曲线会变为两条直线。
S4.进行网格剖分，形成网格***。
网格剖分的步骤具体包括：首先对于油藏整体，若裂缝为嵌入式，不需要作为基岩网格内边界，可直接对油藏进行网格剖分，然后嵌入裂缝***。
MINC模型的网格剖分主要就是对于基岩块的剖分，将基岩块划分为嵌套网格，然后采用正交网格剖分基岩单元。
由油藏基岩的网格***、MINC模型的网格单元和裂缝***，构成整体的网格***。
S5.模型求解：制定求解方案，求解复杂的压裂油藏组分模型，使得模拟过程稳定。
所述求解复杂的压裂油藏组分模型的具体步骤包括：耦合组分模型和裂缝流动模型的压裂油藏组分模型，为了使气驱模拟过程收敛，采用隐式求解方案；为降低矩阵维度，将CO2驱油数学模型中的方程分为主方程和辅助方程，将变量分为主变量和辅助变量，数值求解过程中，先通过求解主方程得到主变量，然后通过辅助方程获得辅助变量。
上述实施例为本发明较佳的实施方式，但本发明的实施方式并不受上述实施例的限制，其他的任何未背离本发明的精神实质与原理下所作的改变、修饰、替代、组合、简化，均应为等效的置换方式，都包含在本发明的保护范围之内。

Claims (10)

1.一种压裂油藏CO2驱油流动模拟方法，其特征在于，所述方法包括以下步骤：
确定压裂产生的人工大裂缝位置及微裂缝区域，建立压裂油藏几何模型；
建立混合模型表征不同尺度裂缝，进行网格剖分，形成网格***；
建立组分模型描述CO2驱油过程，对其中的组分予以考虑，描述相态变化及CO2和原油的混相机理；
制定求解方案，求解复杂的压裂油藏组分模型，使得模拟过程稳定。
2.根据权利要求1所述方法，其特征在于，所述建立混合模型的具体步骤包括：
对人工压裂产生的大裂缝建立嵌入式离散裂缝模型；
对压裂或注气过程中产生的微尺度裂缝建立MINC模型；
基于守恒定律，将基岩***、嵌入式离散裂缝***和MINC模型耦合起来，建立混合模型，将该混合模型作为人工压裂油藏流动模型，以此为基础建立流动方程。
3.根据权利要求2所述方法，其特征在于，对压裂或注气过程中产生的微尺度裂缝建立MINC模型的方法包括：
首先建立双重介质模型：
Figure FDA0002666135380000011
Figure FDA0002666135380000012
F＝αmKmτm(pm-pf)pm≥pf
F＝αmKmτf(pm-pf)pm＜pf
其中，n＝1或2，分别代表双重介质模型里面的基岩***和裂缝***，vn***内流体流动速度，F为两个***间的窜流量，qn为***内源汇项，τn代表***内流体总体流度；Kn代表***绝对渗透率张量，pn为***内流体压力，αm为基岩***的形状因子，Km为基岩***的绝对渗透率张量，τm代表基岩***内流体流度系数，pm代表基岩压力，pf代表裂缝压力，τf代表裂缝***内流体流度系数；
其次，在双重介质模型的基础上将基岩单元进一步划分为几个嵌套的单元，构成MINC模型。
4.根据权利要求2所述方法，其特征在于，对人工压裂产生的大裂缝建立嵌入式离散裂缝模型的方法包括：
对基岩进行剖分，然后将裂缝***嵌入基岩网格当中，
嵌入式离散裂缝模型是直接将裂缝***嵌入油藏网格当中，在进行网格剖分后，将裂缝嵌入网格中，然后计算裂缝***与网格***的相交关系，嵌入式离散模型的质量守恒方程为：
Figure FDA0002666135380000021
其中，kF为大裂缝渗透率，pF为大裂缝压力，μ为流体粘度，qF为大裂缝源汇项，qFf为两个***间的窜流量，qFF表示相交裂缝单元之间的窜流量，δFF表示大裂缝单元是否与其他裂缝单元相交，若大裂缝单元与其他裂缝单元相交δFF＝1，否则δFF＝0，VF分别是裂缝单元的体积。
5.根据权利要求2所述方法，其特征在于，网格剖分的步骤具体包括：首先对于油藏整体，若裂缝为嵌入式，不需要作为基岩网格内边界，可直接对油藏进行网格剖分，然后嵌入裂缝***；
MINC模型的网格剖分主要就是对于基岩块的剖分，将基岩块划分为嵌套网格，然后采用正交网格剖分基岩单元；
由油藏基岩的网格***、MINC模型的网格单元和裂缝***，构成整体的网格***。
6.根据权利要求2所述方法，其特征在于，计算嵌入式离散裂缝模型的参数方法：对于人工压裂的大裂缝需要确定其裂缝空间坐标，条数，间距等，然后在此基础上计算每条裂缝所穿过的网格单元以及裂缝与网格单元的连接信息，包括裂缝与网格的接触面积，裂缝所在网格到裂缝面的平均距离等信息，为后续进行混合模型内的组分流动模拟做准备。
7.根据权利要求2所述方法，其特征在于，计算MINC模型参数方法：根据油田实际资料确定MINC区域的大小，然后计算MINC模型中每个子单元的体积分数，相邻关系，各子单元空间坐标，子单元体积等参数。
8.根据权利要求1所述方法，其特征在于，建立组分模型描述CO2驱油过程：
(1)通过实验结果确定组分数量，针对每一个组分建立质量守恒方程
Figure FDA0002666135380000031
其中i表示组分，i＝1,…,nc,nc表示组分数量。F为质量流量，q代表源汇项，N为质量变化；N可以通过相摩尔密度ρ，饱和度S，油相中组分i的摩尔分数xi，气相中组分i的摩尔分数yi表示为
Figure FDA0002666135380000034
F可以写为
Figure FDA0002666135380000032
其中vβ为相β的达西速度
Figure FDA0002666135380000033
其中β表示油相，气相；
(2)各组分要满足组分质量守恒条件：
Figure FDA0002666135380000041
其中
Figure FDA0002666135380000047
表示逸度系数；逸度系数可以通过压力、温度和体积等数据计算；而温度T，压力P，体积V参数可以通过PR状态方程计算；
在求解的过程通过使方程满足质量守恒条件和热动力平衡条件，更新平衡比，最终得到平衡条件下的组分含量等参数；在计算中将裂缝也视为单元，可将组分模型、闪蒸计算在计算中适用于裂缝***。
9.根据权利要求1所述方法，其特征在于，所述模拟CO2驱油混相机理的具体步骤包括：当达到一定压力时，原油与CO2将发生混相，此时油气***相渗曲线将会发生改变，此时采用界面张力矫正法修正相渗曲线
Figure FDA0002666135380000042
Figure FDA0002666135380000043
其中：
Figure FDA0002666135380000044
其中，
Figure FDA0002666135380000045
和
Figure FDA0002666135380000046
为非混相下的油气相对渗透率；是krg气相渗透率，krog是油相相对渗透率，kroiw是束缚水饱和度下油或气的相对渗透率，n，no，ng是无因次可调参数，σ*是非混相下的油气界面张力值；
当原油与CO2逐渐混相时，气相和油相相对渗透率发生改变，达到完全混相时，界面张力为变为0，此时相渗曲线会变为两条直线。
10.根据权利要求1所述方法，其特征在于，所述求解复杂的压裂油藏组分模型的具体步骤包括：耦合组分模型和裂缝流动模型的压裂油藏组分模型，为了使气驱模拟过程收敛，采用隐式求解方案；为降低矩阵维度，将CO2驱油数学模型中的方程分为主方程和辅助方程，将变量分为主变量和辅助变量，数值求解过程中，先通过求解主方程得到主变量，然后通过辅助方程获得辅助变量。
CN202010926532.3A 2020-09-04 2020-09-04 一种压裂油藏co2驱油流动模拟方法 Pending CN114154430A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010926532.3A CN114154430A (zh) 2020-09-04 2020-09-04 一种压裂油藏co2驱油流动模拟方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010926532.3A CN114154430A (zh) 2020-09-04 2020-09-04 一种压裂油藏co2驱油流动模拟方法

Publications (1)

Publication Number Publication Date
CN114154430A true CN114154430A (zh) 2022-03-08

Family

ID=80460185

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010926532.3A Pending CN114154430A (zh) 2020-09-04 2020-09-04 一种压裂油藏co2驱油流动模拟方法

Country Status (1)

Country Link
CN (1) CN114154430A (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114611339A (zh) * 2022-05-12 2022-06-10 中国石油大学(华东) 一种裂缝连通性识别方法
CN115935674A (zh) * 2022-12-20 2023-04-07 中国石油大学(北京) 基于co2驱油藏流体时空变化特征的多相带判别方法
CN117610313A (zh) * 2024-01-18 2024-02-27 西安石油大学 页岩油***工裂缝的网格加密方法及装置

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114611339A (zh) * 2022-05-12 2022-06-10 中国石油大学(华东) 一种裂缝连通性识别方法
CN115935674A (zh) * 2022-12-20 2023-04-07 中国石油大学(北京) 基于co2驱油藏流体时空变化特征的多相带判别方法
CN115935674B (zh) * 2022-12-20 2024-03-12 中国石油大学(北京) 基于co2驱油藏流体时空变化特征的多相带判别方法
CN117610313A (zh) * 2024-01-18 2024-02-27 西安石油大学 页岩油***工裂缝的网格加密方法及装置
CN117610313B (zh) * 2024-01-18 2024-04-26 西安石油大学 页岩油***工裂缝的网格加密方法及装置

Similar Documents

Publication Publication Date Title
CN107060746B (zh) 一种复杂裂缝性油藏流动模拟的方法
CN114154430A (zh) 一种压裂油藏co2驱油流动模拟方法
CN103400020B (zh) 一种测算多条相交离散裂缝流动状况的油藏数值模拟方法
CN109446706B (zh) 一种确定脉冲纤维加砂压裂支撑剂团铺置形态的方法
CN111814364A (zh) 一种岩溶储层演化数值模拟方法
CN111104766B (zh) 基于离散裂缝模型的油水两相非达西渗流数值模拟方法
CN109372478B (zh) 一种确定非混相气驱油开采方式的实验方法及装置
CN103399970B (zh) 用离散裂缝线处理进行数模测算油藏流动状况的方法
CN112541287B (zh) 疏松砂岩压裂充填防砂增产调剖一体化设计方法
CN110984973A (zh) 缝洞型碳酸盐岩气藏单井控制储量的确定方法
CN112069654B (zh) 一种碳酸盐岩酸化数值模拟方法
CN204140039U (zh) 裂缝性油藏可视化水驱油物理模型和物理模拟实验装置
Dawei et al. Productivity simulation of hydraulically fractured wells based on hybrid local grid refinement and embedded discrete fracture model
Qasem et al. Recovery performance of partially fractured reservoirs by capillary imbibition
CN112528503A (zh) 一种废弃矿井瓦斯抽采数值模拟分析方法
Zhao et al. Physical simulation of waterflooding development in large-scale fractured-vuggy reservoir considering filling characteristics
CN109187321A (zh) 一种大尺寸型盐水平连通井组水溶造腔模拟实验方法
CN110991016B (zh) 不规则边界油藏两口体积压裂水平井渗流模型的建立方法
CN116167302B (zh) 一种天然气水合物增产模拟中人工复杂裂缝的描述方法
CN109558614A (zh) 页岩气藏多尺度裂缝内气体流动的模拟方法及***
CN107169227A (zh) 一种分段压裂水平井的粗网格模拟方法及***
CN107366529A (zh) 非均质水驱油藏无因次采液指数曲线表征方法
CN112901129B (zh) 一种页岩油藏基质裂缝***高温高压多测点吞吐物理模拟装置
Huang et al. Research advance on prediction and optimization for fracture propagation in stimulated unconventional reservoirs
Zhang et al. Staged fracturing of horizontal wells in continental tight sandstone oil reservoirs: a case study of Yanchang Formation in Western Ordos Basin, China

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination