CN114510854A - 一种循缝找洞的酸压数值模拟结果准确性评价方法 - Google Patents

一种循缝找洞的酸压数值模拟结果准确性评价方法 Download PDF

Info

Publication number
CN114510854A
CN114510854A CN202011286785.5A CN202011286785A CN114510854A CN 114510854 A CN114510854 A CN 114510854A CN 202011286785 A CN202011286785 A CN 202011286785A CN 114510854 A CN114510854 A CN 114510854A
Authority
CN
China
Prior art keywords
acid
numerical simulation
fracture
model
acid 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
CN202011286785.5A
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
Sinopec Northwest Oil Field Co
Original Assignee
China Petroleum and Chemical Corp
Sinopec Northwest Oil Field 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, Sinopec Northwest Oil Field Co filed Critical China Petroleum and Chemical Corp
Priority to CN202011286785.5A priority Critical patent/CN114510854A/zh
Publication of CN114510854A publication Critical patent/CN114510854A/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/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • E—FIXED CONSTRUCTIONS
    • E21—EARTH DRILLING; MINING
    • E21B—EARTH DRILLING, e.g. DEEP 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
    • 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]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
    • G06T17/05—Geographic models
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/10—Numerical modelling
    • 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

Abstract

本发明所循缝找洞的酸压数值模拟结果准确性评价方法，通过建立能够进行酸压数值模拟的缝洞型三维应力场地质模型，并在该模型中进行循缝找洞酸压数值模拟计算，明确沟通溶洞的数量，根据沟通溶洞的储集体体积计算指定产量下的理论稳定供液时间；并将理论稳定供液时间与实际生产过程中稳定供液时间进行对比，将循缝找洞酸压数值模拟的准确性从定性上升到了定量，进而提高循缝找洞酸压数值模拟的准确性。

Description

一种循缝找洞的酸压数值模拟结果准确性评价方法
技术领域
本发明涉及缝洞型油藏的开发技术领域，尤其涉及一种循缝找洞的酸压数值模拟结果准确性评价方法。
背景技术
在我国油气资源较为丰富，且大部分都储集在碳酸盐岩储层中。碳酸盐岩储层比较突出的特征就是发育有大量的天然裂缝、溶洞，勘探开发实践表明溶洞是油气资源的储集空间，而天然裂缝则是油气的主要渗流通道。对于缝洞型油气藏，储层改造成败的关键就在于裂缝能否与溶洞储集体沟通，而数值模拟是一种研究缝洞沟通的有效手段。
目前已经有大量关于缝洞型碳酸盐岩储层酸压数值模拟的研究，但一般都是针对单缝单洞进行机理研究，不适用于实际工程情况。同时，酸压数值模拟结果准确性一般用微地震的方法判断，即实际施工后通过微地震信号判断裂缝是否与溶洞沟通，微地震信号受缝、洞等储层非均质性结构的影响非常大，往往不能准确表征实际的沟通情况，所以此种准确性评价方法的评价结果严重受限于微地震信号的接收与处理精度，且此种评价方法成本较高。
因此，有必要针对缝洞型碳酸盐岩储层循缝找洞提出一种全新的数值模拟结果评价方法。
发明内容
本发明提供一种循缝找洞的酸压数值模拟结果准确性评价方法，用于碳酸盐岩储层的酸压数值模拟结果评价，通过将理论稳定供液时间与实际稳定供液时间进行对比，使得循缝找洞酸压数值模拟的准确性从定性上升到了定量，进而提高循缝找洞酸压数值模拟的准确性。
本发明技术方案如下：
一种循缝找洞的酸压数值模拟结果准确性评价方法，包括以下步骤：
S1，在有限元处理软件中，建立目标缝洞型油气藏储层的三维地质模型；
S2，根据收集的目标缝洞型油气藏储层的地应力场数据，在有限元处理软件中建立目标缝洞型油气藏储层基于所述三维地质模型的三维应力场分布模型；
S3，在所述三维应力场分布模型中，标记目标缝洞型油气藏储层的溶洞位置和油气井位置，形成含天然裂缝及溶洞的三维应力场分布模型，并将包括达西尺度模型和孔隙尺度模型的多尺度数学模型赋予所述含天然裂缝及溶洞的三维应力场分布模型，形成能够进行酸压数值模拟的缝洞三维应力场地质模型；
S4，计算或获取指定产量下的实际稳定供液时间；
S5，在所述能够进行酸压数值模拟的缝洞三维应力场地质模型中，进行在S4中指定产量下的实际稳定供液时间的生产施工条件相同的条件进行循缝找洞酸压数值模拟计算，并确定指定产量下的理论稳定供液时间；
S6,判断循缝找洞的酸压数值模拟结果准确性；若理论稳定供液时间和实际稳定供液时间的相对误差绝对值在5％之内，则数值模拟结果准确；若该相对误差绝对值大于5％，则数值模拟结果不准确。
作为优选，所述S4中，进行不同生产施工条件下的循缝找洞酸压数值模拟计算时，输入实际生产施工条件下的对应数值后，所述能够进行酸压数值模拟的缝洞三维应力场地质模型形成主裂缝并开启天然裂缝，且主裂缝和开启的天然裂缝沟通溶洞，循缝找洞酸压数值模拟计算结束后，最终确定能够沟通的溶洞。
作为优选，所述理论稳定供液时间为酸压数值模拟计算中沟通的所有溶洞的总体积与指定产量的比值。
作为优选，所述指定产量下的实际稳定供液时间，是指在与所述酸压数值模拟计算条件相同的实际生产施工条件下，按照指定产量进行实际生产时，保持该指定产量恒定的时间。
作为优选，所述理论稳定供液时间和实际稳定供液时间的相对误差为，所述理论稳定供液时间和所述实际稳定供液时间的差值的绝对值与所述实际稳定供液时间的比值。
作为优选，所述S3中，将包括达西尺度模型和孔隙尺度模型的多尺度数学模型赋予所述三维应力场分布模型，形成能够进行酸压数值模拟的缝洞三维应力场地质模型，包括如下步骤：S3.1建立达西尺度模型：达西尺度模型用于描述厘米级至微米级多孔介质的模型，酸液在厘米级至微米级多孔介质内为达西渗流；
(1)基质中的酸液流动
酸液以一定速度注入地层，酸液中的氢离子在流体中的宏观运动(对流作用)和浓度梯度(扩散作用)下，由孔隙介质流体中传质到碳酸盐岩表面发生反应，进而改变地层孔隙度和渗透率。酸液在基质地层中的流动由达西定律控制：
Figure RE-GDA0002965964320000021
流体压力分布由不可压缩流体的连续性方程控制：
Figure RE-GDA0002965964320000022
氢离子在流体中的浓度分布由对流扩散方程控制，其中对流扩散方程分两种情况，情况一为酸液未将岩石完全溶蚀(ε<1)，情况二为酸液将岩石完全溶蚀(ε＝1)；
酸液未将岩石完全溶蚀(ε<1)情况下的对流扩散方程，需考虑氢离子在岩石表面的消耗和孔隙度的改变：
Figure RE-GDA0002965964320000031
酸液将岩石完全溶蚀(ε＝1)情况下的对流扩散方程，不需考虑氢离子在岩石表面的消耗和孔隙度的改变：
Figure RE-GDA0002965964320000032
式中，
Figure RE-GDA0002965964320000033
为达西速度矢量，m/s；k为地层渗透率，m2；μ为酸液粘度，Pa·s；P为酸液压力，Pa；ε为地层孔隙度；t为反应时间，s；Cf为岩石孔隙中的酸液浓度，mol/m3；De为酸液扩散张量，m2/s；kc为酸液局部传质系数，m/s；av为单位体积岩石具有的孔隙面积，m2/m3； Cs为岩石表面的酸液浓度，mol/m3；α为单位摩尔的酸液所能溶蚀的岩石质量，kg/mol；ρs为岩石密度，kg/m3；
(2)裂缝中的酸液流动
酸液在天然裂缝中的流动与在基质中的流动对应着不同的机理，天然裂缝中流动为自由流，基质中的流动为达西定律控制的多孔介质渗流；酸液在天然裂缝中的流动根据等效渗透率的概念将天然裂缝看作渗透率较大的区域，利用上述数学模型，研究酸液在裂缝性地层中的溶蚀现象，并将对压力场影响较大的天然裂缝进行网格加密，加快计算速度并保证计算的收敛性；
S3.2建立孔隙尺度模型，用于为所述达西尺度模型提供包括渗透率、孔隙半径、比表面积、孔隙度、扩散张量、传质系数的参数支持；
S3.3，将构建的含天然裂缝及溶洞的三维应力场分布模型进行网格化，将所述多尺度数学模型计算得到的参数赋予所述含天然裂缝及溶洞的三维应力场分布模型的网格单元、节点上；得到网格化后及赋值属性参数后的能够进行酸化数值模拟的缝洞三维应力场地质模型。
作为优选，在S3.2中孔隙尺度模型为所述达西尺度模型提供包括渗透率、孔隙半径、比表面积、孔隙度、扩散张量、传质系数的参数支持中，渗透率、孔隙半径、比表面积与孔隙度的关系包括：
地层渗透率、孔隙半径、比表面积大小与孔隙度直接相关，采用经验公式描述地层渗透率、孔隙半径、比表面积与孔隙度之间的关系为：
Figure RE-GDA0002965964320000041
式中，ε(0<ε<1)、k、rp、av分别是孔隙度、渗透率、孔隙半径、比表面积；ε0、k0、 r0、a0分别是初始的孔隙度、渗透率、孔隙半径、比表面积；β是与孔隙结构有关的常数，取β＝1；
一般情况下，酸液在地层中的流动速度很小，可以看作是层流运动，可将裂缝看作较细的圆管，利用圆管层流公式来考察流量；
Figure RE-GDA0002965964320000042
式中，Q为流量，cm3/s；ΔP为驱动压差，0.1MPa；D为裂缝直径，cm；A为截面积，cm2；μ为流体粘度，Pa·s；L为裂缝长度，cm；
利用达西公式计算得到的流量如下：
Figure RE-GDA0002965964320000043
取缝宽为0.2cm时，基于公式(6)、公式(7)计算得到天然裂缝对应的等效渗透率为k＝D2/32＝125×103μm2；令天然裂缝孔隙度εmax＝0.999(ε0＝0.05、k0＝0.32)，代入公式(5)得渗透率K：
Figure RE-GDA0002965964320000051
k与K是同一个数量级，即从流动阻力的角度考虑，将宽度为0.2cm的天然裂缝等效为孔隙度为0.999的基质是合理的，对于缝宽超过0.2cm的天然裂缝，由于天然裂缝相对于基质其渗透率已相当高，此时缝宽不是限制裂缝导流能力(渗透率与缝宽的乘积)的因素，所以将天然裂缝等效为缝宽0.2cm的裂缝处理即可。
作为优选，在S3.2中孔隙尺度模型为所述达西尺度模型提供包括渗透率、孔隙半径、比表面积、孔隙度、扩散张量、传质系数的参数支持中，由于地层的沉积压实作用使孔隙结构在横、纵方向上不同，所述扩散张量包括水平扩散张量DeX和垂向扩散张量DeT；所述扩散张量与分子扩散系数Dm、孔隙结构、流动速度等有关，其中佩克莱特数是表示对流与扩散比值的无量纲数；
Figure RE-GDA0002965964320000052
扩散张量分为水平扩散张量DeX和垂向扩散张量DeT的公式为：
DeX＝(αos+λx·PeP)Dm (10)
DeT＝(αos+λT·PeP)Dm (11)
式中，αos是与孔隙结构相关的常数，取αos＝0.5；由岩心扩散实验得到λx＝0.5，λT＝0.1； Dm是分子扩散系数；Pep是佩克莱特数。
作为优选，在S3.2中孔隙尺度模型为所述达西尺度模型提供包括渗透率、孔隙半径、比表面积、孔隙度、扩散张量、传质系数的参数支持中，
传质系数kc的计算方法为：
Figure RE-GDA0002965964320000053
式中，Sh∞为渐进Sherwood数；Rep为孔隙尺度雷诺数，
Figure RE-GDA0002965964320000054
Sc为施密特数，
Figure RE-GDA0002965964320000061
本发明相对于现有技术优势在于：本发明所述的循缝找洞的酸压数值模拟结果准确性评价方法，用于碳酸盐岩储层的循缝找洞酸压数值模拟结果的评估，通过将理论稳定供液时间与实际稳定供液时间进行对比，使得循缝找洞酸压数值模拟的准确性从定性上升到了定量，进而提高循缝找洞酸压数值模拟结果评判的准确性同时，本发明所述循缝找洞的酸压数值模拟结果准确性评价方法，基于已有的实际生产施工条件下，按照指定产量进行实际生产的条件数据赋值给所述能够进行酸压数值模拟的缝洞三维应力场地质模型，进行循缝找洞酸压数值模拟计算，极大地降低了成本，具有较大的应用前景。所述能够进行酸压数值模拟的缝洞三维应力场地质模型通过给含随机天然裂缝和溶洞的地质模型赋予能够描述厘米级至微米级尺度的多孔介质双尺度数学模型，构建能够进行酸化数值模拟的缝洞型碳酸盐岩地质模型，模拟出酸液激活天然裂缝，沟通溶洞的过程，计算出酸蚀裂缝的延伸方向和条数，以及酸液循缝沟通溶洞的方位和数量，便于判断酸蚀裂缝沟通井周潜在储集体的能力，为现场施工设计提供依据
附图说明
图1是本发明循缝找洞的酸压数值模拟结果准确性评价方法的流程图；
图2是本发明循缝找洞的酸压数值模拟结果准确性评价方法模拟出的缝洞沟通情况示意图；
图3是本发明循缝找洞的酸压数值模拟结果准确性评价方法中实际酸压施工后的稳定供液示意图。
具体实施方式
为了便于理解本发明，下面结合具体实施例和对比例，对本发明进行更详细的说明。
S1，在有限元处理软件中，建立目标缝洞型碳酸盐岩油藏储层的三维地质模型；
S2，根据收集的目标缝洞型碳酸盐岩油气藏储层的地应力场数据，在有限元处理软件中建立目标缝洞型碳酸盐岩油气藏储层基于所述三维地质模型的三维应力场分布模型；
S3，在所述三维应力场分布模型中，标记目标缝洞型碳酸盐岩油气藏储层的溶洞位置和油气井位置，形成含天然裂缝及溶洞的三维应力场分布模型，并将包括达西尺度模型和孔隙尺度模型的多尺度数学模型赋予所述含天然裂缝及溶洞的三维应力场分布模型，形成能够进行酸压数值模拟的缝洞三维应力场地质模型；根据测井等数据确定其能够进行酸压数值模拟的缝洞三维应力场地质模型如图2所示，目标油气藏储层为200m×200m的碳酸盐岩缝洞型油藏，且该区块内共发育有6个溶洞及多条天然裂缝。各溶洞半径如下：溶洞1(7.3m)、溶洞2(4.1m)、溶洞3(6.9m)、溶洞4(3.2m)、溶洞5(3.1m)、溶洞6(6.5m)。
S4，计算或获取指定产量下的实际稳定供液时间；所述指定产量下的实际稳定供液时间，是指在与所述酸压数值模拟计算条件相同的实际生产施工条件下，按照指定产量进行实际生产时，保持该指定产量恒定的时间。根据与所述酸压数值模拟计算条件相同的实际施工下获得的已有的数据，确定指定产量为30m3/d时的实际稳定供液时间如图2所示，为93d；
S5，在如图2所示所述能够进行酸压数值模拟的缝洞三维应力场地质模型中，根据合适的循缝找洞酸压工艺，进行在S4中指定产量为30m3/d时的实际稳定供液时间的生产施工条件相同的条件进行循缝找洞酸压数值模拟计算，通过酸压数值模拟计算形成主裂缝并开启天然裂缝，使主裂缝和开启的天然裂缝沟通溶洞，最终确定能够沟通的溶洞为溶洞1、溶洞5、溶洞6；
通过沟通到的溶洞1、溶洞5、溶洞6，的体积，明确指定产量为30m3/d时的理论稳定供液时间，由于溶洞1、溶洞5、溶洞6的总体积约为2905m3，故而理论稳定供液时间＝2905m3/30m3/d，约为97d；
S6，判断循缝找洞的酸压数值模拟结果准确性；理论稳定供液时间和实际稳定供液时间的相对误差＝(97d-93d)/93d≈4.3％。由于相对误差小于5％，故循缝找洞的酸压数值模拟结果较为准确。
应当指出，以上所述具体实施方式可以使本领域的技术人员更全面地理解本发明，但不以任何方式限制本发明。因此，尽管本说明书参照附图和实施例对本发明已进行了详细的说明，但是，本领域技术人员应当理解，仍然可以对本发明进行修改或者等同替换，总之，一切不脱离本发明的精神和范围的技术方案及其改变，其均应涵盖在本发明专利的保护范围当中。

Claims (9)

1.一种循缝找洞的酸压数值模拟结果准确性评价方法，其特征在于，包括以下步骤：
S1，在有限元处理软件中，建立目标缝洞型油气藏储层的三维地质模型；
S2，根据收集的目标缝洞型油气藏储层的地应力场数据，在有限元处理软件中建立目标缝洞型油气藏储层基于所述三维地质模型的三维应力场分布模型；
S3，在所述三维应力场分布模型中，标记目标缝洞型油气藏储层的溶洞位置和油气井位置，形成含天然裂缝及溶洞的三维应力场分布模型，并将包括达西尺度模型和孔隙尺度模型的多尺度数学模型赋予所述含天然裂缝及溶洞的三维应力场分布模型，形成能够进行酸压数值模拟的缝洞三维应力场地质模型；
S4，计算或获取指定产量下的实际稳定供液时间；
S5，在所述能够进行酸压数值模拟的缝洞三维应力场地质模型中，进行在S4中指定产量下的实际稳定供液时间的生产施工条件相同的条件进行循缝找洞酸压数值模拟计算，并确定指定产量下的理论稳定供液时间；
S6,判断循缝找洞的酸压数值模拟结果准确性；若理论稳定供液时间和实际稳定供液时间的相对误差绝对值在5％之内，则数值模拟结果准确；若该相对误差绝对值大于5％，则数值模拟结果不准确。
2.根据权利要求1所述循缝找洞的酸压数值模拟结果准确性评价方法，其特征在于，所述S5中，进行不同生产施工条件下的循缝找洞酸压数值模拟计算时，输入实际生产施工条件下的对应数值后，所述能够进行酸压数值模拟的缝洞三维应力场地质模型形成主裂缝并开启天然裂缝，且主裂缝和开启的天然裂缝沟通溶洞，循缝找洞酸压数值模拟计算结束后，最终确定能够沟通的溶洞。
3.根据权利要求2所述循缝找洞的酸压数值模拟结果准确性评价方法，其特征在于，所述理论稳定供液时间为酸压数值模拟计算中沟通的所有溶洞的总体积与指定产量的比值。
4.根据权利要求3所述循缝找洞的酸压数值模拟结果准确性评价方法，其特征在于，所述指定产量下的实际稳定供液时间，是指在与所述酸压数值模拟计算条件相同的实际生产施工条件下，按照指定产量进行实际生产时，保持该指定产量恒定的时间。
5.根据权利要求4所述循缝找洞的酸压数值模拟结果准确性评价方法，其特征在于，所述理论稳定供液时间和实际稳定供液时间的相对误差为，所述理论稳定供液时间和所述实际稳定供液时间的差值的绝对值与所述实际稳定供液时间的比值。
6.根据权利要求1所述循缝找洞的酸压数值模拟结果准确性评价方法，其特征在于，所述S3中，将包括达西尺度模型和孔隙尺度模型的多尺度数学模型赋予所述三维应力场分布模型，形成能够进行酸压数值模拟的缝洞三维应力场地质模型，包括如下步骤：
S3.1建立达西尺度模型：达西尺度模型用于描述厘米级至微米级多孔介质的模型，酸液在厘米级至微米级多孔介质内为达西渗流；
(1)基质中的酸液流动
酸液以一定速度注入地层，酸液中的氢离子在流体中的宏观运动(对流作用)和浓度梯度(扩散作用)下，由孔隙介质流体中传质到碳酸盐岩表面发生反应，进而改变地层孔隙度和渗透率，酸液在基质地层中的流动由达西定律控制：
Figure RE-FDA0002965964310000021
流体压力分布由不可压缩流体的连续性方程控制：
Figure RE-FDA0002965964310000022
氢离子在流体中的浓度分布由对流扩散方程控制，其中对流扩散方程分两种情况，情况一为酸液未将岩石完全溶蚀(ε<1)，情况二为酸液将岩石完全溶蚀(ε＝1)；
酸液未将岩石完全溶蚀(ε<1)情况下的对流扩散方程，需考虑氢离子在岩石表面的消耗和孔隙度的改变：
Figure RE-FDA0002965964310000023
酸液将岩石完全溶蚀(ε＝1)情况下的对流扩散方程，不需考虑氢离子在岩石表面的消耗和孔隙度的改变：
Figure RE-FDA0002965964310000024
式中，
Figure RE-FDA0002965964310000025
为达西速度矢量，m/s；k为地层渗透率，m2；μ为酸液粘度，Pa·s；P为酸液压力，Pa；ε为地层孔隙度；t为反应时间，s；Cf为岩石孔隙中的酸液浓度，mol/m3；De为酸液扩散张量，m2/s；kc为酸液局部传质系数，m/s；av为单位体积岩石具有的孔隙面积，m2/m3；Cs为岩石表面的酸液浓度，mol/m3；α为单位摩尔的酸液所能溶蚀的岩石质量，kg/mol；ρs为岩石密度，kg/m3；
(2)裂缝中的酸液流动
酸液在天然裂缝中的流动与在基质中的流动对应着不同的机理，天然裂缝中流动为自由流，基质中的流动为达西定律控制的多孔介质渗流；酸液在天然裂缝中的流动根据等效渗透率的概念将天然裂缝看作渗透率较大的区域，利用上述数学模型，研究酸液在裂缝性地层中的溶蚀现象，并将对压力场影响较大的天然裂缝进行网格加密，加快计算速度并保证计算的收敛性；
S3.2建立孔隙尺度模型，用于为所述达西尺度模型提供包括渗透率、孔隙半径、比表面积、孔隙度、扩散张量、传质系数的参数支持；
S3.3，将构建的含天然裂缝及溶洞的三维应力场分布模型进行网格化，将所述多尺度数学模型计算得到的参数赋予所述含天然裂缝及溶洞的三维应力场分布模型的网格单元、节点上；得到网格化后及赋值属性参数后的能够进行酸化数值模拟的缝洞三维应力场地质模型。
7.根据权利要求6所述循缝找洞的酸压数值模拟结果准确性评价方法，其特征在于，在S3.2中孔隙尺度模型为所述达西尺度模型提供包括渗透率、孔隙半径、比表面积、孔隙度、扩散张量、传质系数的参数支持中，渗透率、孔隙半径、比表面积与孔隙度的关系包括：
地层渗透率、孔隙半径、比表面积大小与孔隙度直接相关，采用经验公式描述地层渗透率、孔隙半径、比表面积与孔隙度之间的关系为：
Figure RE-FDA0002965964310000031
式中，ε(0<ε<1)、k、rp、av分别是孔隙度、渗透率、孔隙半径、比表面积；ε0、k0、r0、a0分别是初始的孔隙度、渗透率、孔隙半径、比表面积；β是与孔隙结构有关的常数，取β＝1；
一般情况下，酸液在地层中的流动速度很小，可以看作是层流运动，可将裂缝看作较细的圆管，利用圆管层流公式来考察流量；
Figure RE-FDA0002965964310000041
式中，Q为流量，cm3/s；ΔP为驱动压差，0.1MPa；D为裂缝直径，cm；A为截面积，cm2；μ为流体粘度，Pa·s；L为裂缝长度，cm；
利用达西公式计算得到的流量如下：
Figure RE-FDA0002965964310000042
取缝宽为0.2cm时，基于公式(6)、公式(7)计算得到天然裂缝对应的等效渗透率为k＝D2/32＝125×103μm2；令天然裂缝孔隙度εmax＝0.999(ε0＝0.05、k0＝0.32)，代入公式(5)得渗透率K：
Figure RE-FDA0002965964310000043
k与K是同一个数量级，即从流动阻力的角度考虑，将宽度为0.2cm的天然裂缝等效为孔隙度为0.999的基质是合理的，对于缝宽超过0.2cm的天然裂缝，由于天然裂缝相对于基质其渗透率已相当高，此时缝宽不是限制裂缝导流能力(渗透率与缝宽的乘积)的因素，所以将天然裂缝等效为缝宽0.2cm的裂缝处理即可。
8.根据权利要求7所述循缝找洞的酸压数值模拟结果准确性评价方法，其特征在于，在S3.2中孔隙尺度模型为所述达西尺度模型提供包括渗透率、孔隙半径、比表面积、孔隙度、扩散张量、传质系数的参数支持中，由于地层的沉积压实作用使孔隙结构在横、纵方向上不同，所述扩散张量包括水平扩散张量DeX和垂向扩散张量DeT；所述扩散张量与分子扩散系数Dm、孔隙结构、流动速度等有关，其中佩克莱特数是表示对流与扩散比值的无量纲数；
Figure RE-FDA0002965964310000044
扩散张量分为水平扩散张量DeX和垂向扩散张量DeT的公式为：
DeX＝(αos+λx·PeP)Dm (10)
DeT＝(αos+λT·PeP)Dm (11)
式中，αos是与孔隙结构相关的常数，取αos＝0.5；由岩心扩散实验得到λx＝0.5，λT＝0.1；Dm是分子扩散系数；Pep是佩克莱特数。
9.根据权利要求4所述循缝找洞的酸压数值模拟结果准确性评价方法，其特征在于，在S3.2中孔隙尺度模型为所述达西尺度模型提供包括渗透率、孔隙半径、比表面积、孔隙度、扩散张量、传质系数的参数支持中，
传质系数kc的计算方法为：
Figure RE-FDA0002965964310000051
式中，Sh∞为渐进Sherwood数；Rep为孔隙尺度雷诺数，
Figure RE-FDA0002965964310000052
Sc为施密特数，
Figure RE-FDA0002965964310000053
CN202011286785.5A 2020-11-17 2020-11-17 一种循缝找洞的酸压数值模拟结果准确性评价方法 Pending CN114510854A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011286785.5A CN114510854A (zh) 2020-11-17 2020-11-17 一种循缝找洞的酸压数值模拟结果准确性评价方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011286785.5A CN114510854A (zh) 2020-11-17 2020-11-17 一种循缝找洞的酸压数值模拟结果准确性评价方法

Publications (1)

Publication Number Publication Date
CN114510854A true CN114510854A (zh) 2022-05-17

Family

ID=81546431

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011286785.5A Pending CN114510854A (zh) 2020-11-17 2020-11-17 一种循缝找洞的酸压数值模拟结果准确性评价方法

Country Status (1)

Country Link
CN (1) CN114510854A (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114724641A (zh) * 2022-06-07 2022-07-08 中国石油大学(华东) 一种双矿物组分碳酸盐岩储层酸化数值模拟方法
CN116525021A (zh) * 2023-03-24 2023-08-01 中国石油大学(华东) 一种基于改进双尺度模型的碳酸盐岩酸化模拟方法
CN116579263A (zh) * 2023-05-17 2023-08-11 中国石油大学(北京) 基于油气井排采动态数据的综合分析方法
CN117494601A (zh) * 2023-11-06 2024-02-02 西南石油大学 一种基于嵌入式离散裂缝的缝洞型储层酸压效果评价方法

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114724641A (zh) * 2022-06-07 2022-07-08 中国石油大学(华东) 一种双矿物组分碳酸盐岩储层酸化数值模拟方法
CN116525021A (zh) * 2023-03-24 2023-08-01 中国石油大学(华东) 一种基于改进双尺度模型的碳酸盐岩酸化模拟方法
CN116525021B (zh) * 2023-03-24 2024-03-22 中国石油大学(华东) 一种基于改进双尺度模型的碳酸盐岩酸化模拟方法
CN116579263A (zh) * 2023-05-17 2023-08-11 中国石油大学(北京) 基于油气井排采动态数据的综合分析方法
CN116579263B (zh) * 2023-05-17 2024-01-30 中国石油大学(北京) 基于油气井排采动态数据的综合分析方法
CN117494601A (zh) * 2023-11-06 2024-02-02 西南石油大学 一种基于嵌入式离散裂缝的缝洞型储层酸压效果评价方法

Similar Documents

Publication Publication Date Title
CN114510854A (zh) 一种循缝找洞的酸压数值模拟结果准确性评价方法
CN109441422B (zh) 一种页岩气井间距优化开采方法
Aliyu et al. Sensitivity analysis of deep geothermal reservoir: Effect of reservoir parameters on production temperature
CN103413030B (zh) 一种缝洞型碳酸盐岩气藏动态分析方法及系统
CN104533370A (zh) 压裂水平井油藏、裂缝、井筒全耦合模拟方法
CN109241588B (zh) 一种基于拟连续地质力学模型的单裂缝扩展的模拟方法
CN110147561B (zh) 一种含天然裂缝致密油气储层体积压裂缝网预测方法
CN104992468A (zh) 缝洞型碳酸盐岩油气藏三维地质建模方法
CN110096718B (zh) 一种获取碳酸盐岩油藏中溶洞的体积的方法
CN111794740B (zh) 一种适用在缝洞型碳酸盐岩油藏动态储量计算的方法
CN114510882A (zh) 一种循缝找洞酸化数值模拟方法
Qun et al. A novel approach of tight oil reservoirs stimulation based on fracture controlling optimization and design
CN110984973A (zh) 缝洞型碳酸盐岩气藏单井控制储量的确定方法
WO2017031615A1 (zh) 一种构建双溶腔盐穴储库地面沉降预测模型的方法
CN105260543A (zh) 基于双孔模型的多重介质油气流动模拟方法及装置
CN107066718A (zh) 一种四维动态地应力模拟方法
Tian et al. Shale gas production from reservoirs with hierarchical multiscale structural heterogeneities
CN113255123A (zh) 煤层顶板水平井分段多簇压裂适用性地质条件的评价方法
Wang et al. Application of fuzzy analytic hierarchy process in sandstone aquifer water yield property evaluation
Yang et al. Influence of reservoirs/interlayers thickness on hydraulic fracture propagation laws in low-permeability layered rocks
CN111950112A (zh) 一种适用于底部封闭的碳酸盐岩储层动态分析方法
CN110096669B (zh) 一种获取碳酸盐岩油藏中裂缝体积的方法
CN107832482B (zh) 致密储层多尺度裂缝网络建模及模拟方法
Liu et al. A fractal model for characterizing hydraulic properties of fractured rock mass under mining influence
CN111950111A (zh) 一种适用于底部开放的碳酸盐岩储层动态分析方法

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