CN114818551B - 一种页岩油藏压裂缝网自适应半无限导流表征方法及装置 - Google Patents
一种页岩油藏压裂缝网自适应半无限导流表征方法及装置 Download PDFInfo
- Publication number
- CN114818551B CN114818551B CN202210765504.7A CN202210765504A CN114818551B CN 114818551 B CN114818551 B CN 114818551B CN 202210765504 A CN202210765504 A CN 202210765504A CN 114818551 B CN114818551 B CN 114818551B
- Authority
- CN
- China
- Prior art keywords
- infinite
- network
- model
- representing
- error threshold
- 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.)
- Active
Links
- 238000012512 characterization method Methods 0.000 title claims abstract description 38
- 239000003079 shale oil Substances 0.000 title claims abstract description 16
- 238000000034 method Methods 0.000 claims abstract description 41
- 230000011218 segmentation Effects 0.000 claims abstract description 39
- 238000004364 calculation method Methods 0.000 claims abstract description 22
- 238000002372 labelling Methods 0.000 claims abstract description 10
- 238000003860 storage Methods 0.000 claims abstract description 7
- 238000005192 partition Methods 0.000 claims description 19
- 238000007789 sealing Methods 0.000 claims description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims 1
- 238000004088 simulation Methods 0.000 abstract description 13
- 238000012545 processing Methods 0.000 abstract description 8
- 238000010586 diagram Methods 0.000 description 15
- 230000003044 adaptive effect Effects 0.000 description 9
- 238000004519 manufacturing process Methods 0.000 description 6
- 238000004891 communication Methods 0.000 description 5
- 239000012530 fluid Substances 0.000 description 5
- 239000011159 matrix material Substances 0.000 description 5
- 230000008569 process Effects 0.000 description 4
- 238000012546 transfer Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 3
- 230000006870 function Effects 0.000 description 3
- 208000035126 Facies Diseases 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 238000003795 desorption Methods 0.000 description 2
- 230000009977 dual effect Effects 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000035699 permeability Effects 0.000 description 2
- 230000000750 progressive effect Effects 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- 238000001179 sorption measurement Methods 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 1
- 230000005465 channeling Effects 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000005213 imbibition Methods 0.000 description 1
- 238000010946 mechanistic model Methods 0.000 description 1
- 230000005055 memory storage Effects 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 238000010295 mobile communication Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 229920001296 polysiloxane Polymers 0.000 description 1
- 230000001172 regenerating effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
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]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Forestry; Mining
-
- 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
- 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)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Business, Economics & Management (AREA)
- Mathematical Optimization (AREA)
- Marine Sciences & Fisheries (AREA)
- Fluid Mechanics (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Computing Systems (AREA)
- Life Sciences & Earth Sciences (AREA)
- Agronomy & Crop Science (AREA)
- Animal Husbandry (AREA)
- Mathematical Analysis (AREA)
- Mining & Mineral Resources (AREA)
- Algebra (AREA)
- Health & Medical Sciences (AREA)
- Economics (AREA)
- General Health & Medical Sciences (AREA)
- Human Resources & Organizations (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本申请公开了一种页岩油藏压裂缝网自适应半无限导流表征方法、装置及存储介质,涉及油藏数值模拟技术领域,解决了现有技术中不能有效表征大规模改造后形成的多尺度裂缝的问题。该方法包括:将嵌入式离散裂缝嵌入多重介质模型中;对多重介质模型进行分割,确定背景网络和每一个分割单元的质量守恒方程;根据质量守恒方程和多重介质模型的边界条件,计算每一个分割单元的压力,并在多重介质模型上进行标注;根据压力误差阈值和无限导流误差阈值搜索多重介质模型,确定无限导流裂缝和连通子网络;计算连通子网络与背景网格之间的传导率,确定半无限导流表征模型。该方法实现了对大规模改造后形成的多尺度裂缝能够进行有效的表征,且计算处理速度快。
Description
技术领域
本申请涉及油藏数值模拟技术领域,尤其涉及一种页岩油藏压裂缝网自适应半无限导流表征方法及装置。
背景技术
由于陆相页岩油藏不同储集介质的尺度存在极大的差异,多相流体流动规律具有很大的差异,通常需要对不同尺度和流动特征的储集介质进行区分表征,建立多尺度介质模型。在多尺度介质模型基础上,如何有效表征大规模改造后形成的多尺度裂缝成为页岩油藏压裂水平数值模拟和产能评估的一个难题。
在现有技术中,数值模拟方法对压裂多尺度形成的多尺度裂缝网络的表征可以分为等效介质模型、离散裂缝模型和混合模型三大类。等效介质模型稳定性好、求解速度快,但是精度较低;离散裂缝模型精度高但是求解速度慢,无法处理大规模缝网;混合模型当离散裂缝网格数较多时,计算效率仍比较具有挑战性。
发明内容
本申请实施例通过提供一种页岩油藏压裂缝网自适应半无限导流表征方法,解决了现有技术中不能有效表征大规模改造后形成的多尺度裂缝的问题,实现了对大规模改造后形成的多尺度裂缝能够进行有效的表征,且计算处理速度快。
第一方面,本发明实施例提供了一种页岩油藏压裂缝网自适应半无限导流表征方法,该方法包括:
将嵌入式离散裂缝嵌入多重介质模型中;
对所述多重介质模型进行分割,确定背景网络和每一个分割单元的质量守恒方程;
根据所述质量守恒方程和所述多重介质模型的边界条件,计算每一个所述分割单元的压力,并在所述多重介质模型上进行标注;
根据压力误差阈值和无限导流误差阈值搜索所述多重介质模型,确定无限导流裂缝和连通子网络,所述无限导流误差阈值包括:
其中,CSN表示所述连通子网络,Tmf表示所述背景网络中所述分割单元与所述分割单元内的嵌入式离散裂缝之间的传导率,表示所述连通子网络的平均压力值与所述背景网络中对应所述分割单元的压力值的差值,εq表示无限导流误差阈值,e表示无限导流误差,ΔPmf表示所述连通子网络中所述分割单元的压力值与所述背景网络中对应的压力值的差值;
计算所述连通子网络与所述背景网格之间的传导率,确定半无限导流表征模型。
结合第一方面,在一种可能的实现方式中,所述多重介质模型的边界条件包括:封闭边界条件和井控条件;
所述封闭边界条件是:在某个方向上,所述多重介质模型与外部没有能量交换,表示为:u(x,y)|Ω·n=0;其中,n表示边界处的法向量,u(x,y)|Ω表示边界处的流量速度。
结合第一方面,在一种可能的实现方式中,所述压力误差阈值包括:
其中,Pi表示第i个所述分割单元的压力值,Pj表示第j个所述分割单元的压力值,ep表示第i个所述分割单元和第j个所述分割单元的压力误差值,ΔPij表示第i个所述分割单元和第j个所述分割单元的压力差值,εP表示所述压力误差阈值。
结合第一方面,在一种可能的实现方式中,通过以下步骤确定所述连通子网络:
根据所述压力误差阈值和所述无限导流误差阈值对所述多重介质模型进行搜索后,将符合条件的裂缝进行合并,确定所述连通子网络。
结合第一方面,在一种可能的实现方式中,所述计算所述连通子网络与所述背景网格之间的传导率,包括:采用网格等效距离方法计算所述传导率。
第二方面,本发明实施例提供了一种页岩油藏压裂缝网自适应半无限导流表征装置,该装置包括:
嵌入模块,用于将嵌入式离散裂缝嵌入多重介质模型中;
分割计算模块,用于对所述多重介质模型进行分割,确定背景网络和每一个分割单元的质量守恒方程;
压力值标注模块,用于根据所述质量守恒方程和所述多重介质模型的边界条件,计算所述分割单元的压力,并在所述多重介质模型上进行标注;
无限导流裂缝确定模块,用于根据压力误差阈值和无限导流误差阈值搜索所述多重介质模型,确定无限导流裂缝和连通子网络,所述无限导流误差阈值包括:
其中,CSN表示所述连通子网络,Tmf表示所述背景网络中所述分割单元与所述分割单元内的嵌入式离散裂缝之间的传导率,表示所述连通子网络的平均压力值与所述背景网络中对应所述分割单元的压力值的差值,εq表示无限导流误差阈值,e表示无限导流误差,ΔPmf表示所述连通子网络中所述分割单元的压力值与所述背景网络中对应的压力值的差值;
表征模块确定模块,用于计算所述连通子网络与所述背景网格之间的传导率,确定半无限导流表征模型。
结合第二方面,在一种可能的实现方式中,所述多重介质模型包括封闭边界条件和井控条件;
所述封闭边界条件是:在某个方向上,所述多重介质模型与外部没有能量交换,表示为:u(x,y)|Ω·n=0;其中,n表示边界处的法向量,u(x,y)|Ω表示为边界处的流量速度。
结合第二方面,在一种可能的实现方式中,所述无限导流裂缝确定模块包括根据以下公式来计算所述压力误差阈值;
其中,Pi表示第i个所述分割单元的压力值,Pj表示第j个所述分割单元的压力值,ep表示第i个所述分割单元和第j个所述分割单元的压力误差值,ΔPij表示第i个所述分割单元和第j个所述分割单元的压力差值,εP表示所述压力误差阈值。
结合第二方面,在一种可能的实现方式中,所述无限导流裂缝确定模块具体通过以下步骤确定所述连通子网络:
根据所述压力误差阈值和所述无限导流误差阈值对所述多重介质模型进行搜索后,将符合条件的裂缝进行合并,确定所述连通子网络。
结合第二方面,在一种可能的实现方式中,所述表征模块确定模块采用网格等效距离方法计算所述传导率。
第三方面,本发明实施例提供了一种页岩油藏压裂缝网自适应半无限导流表征服务器,该服务器包括存储器和处理器;
所述存储器用于存储计算机可执行指令;
所述处理器用于执行所述计算机可执行指令,以实现第一方面或第一方面任一种可能的实现方式所述的方法。
第四方面,本发明实施例提供了一种计算机可读存储介质,所述计算机可读存储介质存储有可执行指令,计算机执行所述可执行指令时能够实现第一方面或第一方面任一种可能的实现方式所述的方法。
本发明实施例中提供的一个或多个技术方案,至少具有如下技术效果或优点:
本发明实施例通过采用了一种页岩油藏裂缝自适应半无限导流的表征方法,该方法发包括:将嵌入式离散裂缝嵌入多重介质模型中;对多重介质模型进行分割,确定背景网络和每一个分割单元的质量守恒方程;根据质量守恒方程和多重介质模型的边界条件,计算每一个分割单元的压力,并在多重介质模型上进行标注;根据压力误差阈值和无限导流误差阈值搜索多重介质模型,确定无限导流裂缝和连通子网络;计算连通子网络与背景网格之间的传导率,确定半无限导流表征模,本申请提供的方法中根据压力误差阈值和无限导流误差阈值搜索多重介质模型,在这一步骤中的压力误差阈值和无限导流误差阈值是能够进行人为定义的,能够将离散裂缝进行合并。不同的阈值设定导致计算速度快慢不同,合并的离散裂缝也不同,即得到的连通子网络也不同。同时计算机的计算压力随人为设定阈值变化。该方法解决了现有技术中不能有效表征大规模改造后形成的多尺度裂缝的问题,实现了对大规模改造后形成的多尺度裂缝能够进行有效的表征,且计算处理速度快。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对本发明实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本申请实施例提供的页岩油藏压裂缝网自适应半无限导流表征方法的步骤流程图;
图2为本申请实施例提供的多重介质模型示意图;
图3为本申请实施例提供的多重介质模型连通关系示意图;
图4为本申请实施例提供的一个分割单元中嵌入裂缝的示意图;
图5为本申请实施例提供的井控条件和封闭边界条件示意图;
图6为本申请实施例提供的将两个裂缝合并为无限导流裂缝的示意图;
图7为本申请实施例提供的单井概念模型示意图;
图8为本申请实施例提供的嵌入式离散裂缝示意图;
图9为本申请实施例提供的将离散裂缝和压力场嵌入单井概念模型后的示意图;
图10A至图10D为本申请实施例提供的不同阈值下连通子网络示意图;
图11为本申请实施例提供的单井概念模型模拟生产曲线对比图;
图12A至图12D为本申请实施例提供的连通子网络数目不同时半无限导流模型示意图;
图13为本申请实施例提供的单井概念模型数值模拟时间对比图;
图14为本申请实施例提供的页岩油藏压裂缝网自适应半无限导流表征装置示意图;
图15为本申请实施例提供的页岩油藏压裂缝网自适应半无限导流表征服务器示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述。显然,所描述的实施例是本发明的一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
由于陆相页岩油藏不同储集介质的尺度存在极大的差异,多项流体流动规律具有很大的差异,通常需要对不同尺度和流动特征的储集介质进行区分表征,建立多尺度介质模型。在多尺度介质模型基础上,如何有效表征大规模改造后形成了多尺度裂缝成为页岩油藏压裂水平井数值模拟和产能评价的一个难题。
裂缝表征的关键是确定多尺度介质与裂缝之间的传质表征。现有数值模拟方法对压裂多尺度形成的多尺度裂缝网络的表征可等效的分为等效介质模型、离散裂缝模型和混合模型三大类。
等效介质模型是现阶段裂缝性油藏数值模拟应用最广泛的模型,其基本原理是将裂缝和基质分开建模,通过窜流因子表征两者之间的流体质量交换,适合处理连通性好的天然裂缝储层。在双重介质模型基础上,Pruess和Narasimhan,Wu和Pruess分别于1985年及1988年对双重介质模型进行扩展,提出了一个多重介质模型,该模型精细划分基质块描述裂缝和基质间的压力梯度,提高了油藏渗吸、驱替以及缝洞型介质等情形下的精度。等效介质模型的优点是稳定性好、求解速度快,但是精度较低。
离散裂缝模型。为了更精确表征裂缝的影响,人们开始使用离散裂缝模型来表征裂缝。Karimi-Fard等人于2004年提出基于有限体积方法的离散裂缝模型。该方法采用对裂缝结合体进行降维简化处理,在保证精度的同时,提高了模型适用性,成功应用到了二维和三维油藏中的单相和多相流模拟。Lee等人和Moinfar等人提出了嵌入式离散裂缝模型,该模型将裂缝网络直接嵌入基岩结构化网格系统中,避免了复杂的非结构化网格剖分过程。离散裂缝模型的优点是精度高,缺点是求解较慢,无法处理大规模缝网。
混合模型。为了处理多尺度裂缝,许多学者开始把不同裂缝表征模型进行融合。为了区分压裂裂缝与天然裂缝对流动的影响,Xu等提出了一个改进的预测模型,在基质和嵌入式离散裂缝之间加入了天然裂缝等效连续介质,形成DP-EDFM模型;Zhang等提出了一个四重介质模型,描述了页岩气体从有机介质-无机介质-天然裂缝-人工裂缝的流动过程,在他们的模型中对人工裂缝也进行了连续介质等效处理;Jiang和Younis提出了一个通用的多重相互作用连续统一体耦合面向制造的设计/嵌入式离散裂缝的模型,并给出了构建模型的流程,他们的模型忽略了基质与嵌入式裂缝之间的流动。然而,由于巨大的体积和导流能力差异,面向制造的设计/嵌入式离散裂缝网格对于模型求解收敛性影响较大,当离散裂缝网格数较多时,计算效率仍比较具有挑战性。基于现有技术中存在的问题,本申请实施例提供了一种页岩油藏压裂缝网自适应半无限导流表征方法,该方法如图1所示包括以下步骤S101至S105。
S101:将嵌入式离散裂缝嵌入多重介质模型中。图2所示为多重介质模型示意图。
S102:对多重介质模型进行分割,确定背景网络和每一个分割单元的质量守恒方程。
S103:根据质量守恒方程和多重介质模型的边界条件,计算每一个分割单元的压力,并在多重介质模型上进行标注。
S104:根据压力误差阈值和无限导流误差阈值搜索多重介质模型,确定无限导流裂缝和连通子网络,无限导流误差阈值包括:
其中,CSN表示连通子网络,Tmf表示背景网络中分割单元与分割单元内的嵌入式离散裂缝之间的传导率,表示连通子网络的平均压力值与背景网络中对应分割单元的压力值的差值,εq表示无限导流误差阈值,e表示无限导流误差,ΔPmf表示连通子网络中分割单元的压力值与背景网络中对应的压力值的差值;
S105:计算连通子网络与背景网格之间的传导率,确定半无限导流表征模型。
在步骤S102中,将多重介质模型分割成多个分割单元,假设有机介质只与其他介质发生吸附解吸传质,无机介质和等效裂缝介质之间发生渗流和吸附解吸传质。为了更精细描述压力等变量在裂缝附近范围内的强烈变化,还可以继续嵌套分割每一种介质网络。图3所示是多重介质模型连通关系图。多重介质模型通过传导率表征不同介质间以及相邻网格间的传质关系,其连通关系为Qij=Tijλ(Pi-Pj),其中Tij表示网格i和网格j之间的传导率,λ表示流体流度,Pi表示网格i的压力,Pj表示网格j的压力。
为了能够提高流动模拟的精度,对于压裂形成的多尺度复杂缝网进行嵌入离散裂缝表征。如图4所示是一个分割单元中嵌入裂缝的示意图。对已经嵌入式离散裂缝的网格传导率采用等效距离方法计算,计算公式为:
其中,Amf表示嵌入裂缝的面积;kf表示嵌入单位裂缝的渗透率;km表示背景单元的渗透率;<d>表示背景网格与嵌入式离散裂缝之间的等效距离;n表示嵌入裂缝的法向量;x表示背景网格内的一点到裂缝的距离向量,S表示背景网格的体积。
在步骤S102中背景网络和每一个分割单元的质量守恒方程为:
其中,WI表示井指数,上标w表示编号为w的井,pw表示编号为w的源汇井的压力,Pj表示第j个分割单元的压力值。
在步骤S103中,多重介质模型的边界条件包括:封闭边界条件和井控条件;井控条件和封闭边界条件如图5所示。其中封闭边界条件是:在某个方向上,多重介质模型与外部没有能量交换,表示为:
u(x,y)|Ω·n=0 (5),
其中,n表示边界处的法向量,u(x,y)|Ω表示为边界处的流量速度。
井控条件如上述公式(6)表示。
联合上述公式(4)(5)(6),根据线性方程组求解算法,可以求出多重介质模型的压力场P。以上述的方法计算每一个分割单元的压力,并将计算所得的压力值在多重介质模型上进行标注。
压力误差阈值包括:
其中,Pi表示第i个分割单元的压力值,Pj表示第j个分割单元的压力值,ep表示第i个分割单元和第j个分割单元的压力误差值,ΔPij表示第i个分割单元和第j个分割单元的压力差值,εP表示压力误差阈值。
无限导流误差阈值包括:
其中,CSN表示连通子网络,Tmf表示背景网络中分割单元与分割单元内的嵌入式离散裂缝之间的传导率,表示连通子网络的平均压力值与背景网络中对应分割单元的压力值的差值,εq表示无限导流误差阈值,e表示无限导流误差,ΔPmf表示连通子网络中分割单元的压力值与背景网络中对应的压力值的差值。
通过以下步骤确定连通子网络:根据压力误差阈值和无限导流误差阈值对多重介质模型进行搜索后,将符合条件的裂缝进行合并,确定连通子网络。根据公式(7)和公式(8)联合计算,将满足公式(8)的裂缝字网络设置为无限导流裂缝,重新生成嵌入裂缝分割单元。将多个无限导流裂缝组成的网络称为连通子网络。
计算连通子网络与背景网格之间的传导率,包括:采用网格等效距离方法计算传导率。对连通子网络之间重新形成的分割单元f*,并计算重新形成的分割单元的体积,删除原有的网格传导率Tmf,重新计算网格的传导率,其公式表示为,其中,CSN表示连通子网络。
在本申请一个具体的应用实例中,建立单井概念模型,如图7所示为单井概念模型示意图。
考虑如图7所示单井概念模型,模型水平段长为1600米,7段体积压裂,段间距为200米,详细参数如表1所示。
表1机理模型参数
对单井概念模型添加嵌入式离散裂缝,图8为添加了嵌入式离散缝隙的模型示意图。
采用本申请提供的方法建立半无限导流表征模型。首先将压裂缝网裂缝嵌入单井概念模型,共形成2026个嵌入缝计算单元。根据如图1的流程,基于质量守恒和传导率连通表建立压力场,如图9所示。
确定表示压力误差阈值εp和无限导流误差阈值εq,搜索并确定连通子网络。需要注意,可以根据确定不同的εp和εq来改变连通子网络的数量,以此改变模型精度和计算速度,如图10A至10D所示。图10A的参数为εp=0.05,εq=0.05,NCSN=10,图10B的参数为εp=0.1,εq=0.05,NCSN=15,图10C的参数为εp=0.05,εq=0.1,NCSN=22,图10D的参数为εp=0.1,εq=0.1,NCSN=30。εp和εq越趋于0,则模型越接近全嵌入离散裂缝模型,εp和εq越趋于1,则模型越接近全无限导流模型。
同时,开展直接基于全嵌入离散裂缝模型的流动模拟,作为参考对照模型,并设计全无限导流对比模型作为对比模型。模拟结果如图11至图13所示。模拟生产曲线结果图11所示,横坐标代表生产时间(天),左纵坐标表示日产油(方/天),右纵坐标表示累产油(方),结果显示了本专利提出模型在连通子网络数NCSN为30时,与全嵌入离散裂缝模型较为接近,具有较高精度,相对误差仅为1.2%。随着连通子网络数的减小,较高估计了裂缝导流能力,导致预测产能偏高。对于全无限导流模型,由于其将所有裂缝视为全连通无限导流,偏差最大,预测累产油累积误差达到17%。模拟压力场结果如图12A至图12D所示,同样可以看出,当连通子网络数NCSN为30时,与对照模型十分接近。即便是NCSN为10的情况,无论在产能预测还是压力场分布,都比全无限导流模型要更加接近参考模型。由图13可以看出,由于大幅减少了小体积网格数和高传导率项,模拟时间均能大幅减少,相比全嵌入模型,连通子网络数NCSN为30的半无限导流模型可以提速3.3倍;而相比全无限导流模型,耗时多了大约11.3%。单井概念模型显示了方法的可靠性、稳定性,同时对计算速度有了提速的效果。
基于本申请提供的方法,对于压力误差阈值和无限导流误差阈值都是可变的,根据在实际应用的场景中对其进行不同的设定,随之连通子网络也随之变化,更进一步可以采用层层递进的计算方式来进行不同裂缝进行合并,这样能在一定程度上更近一步的减小运算量,加快计算速度,且不同的压力误差阈值和无限导流误差阈值带来的计算速率和计算的结果也不相同。对于具体的应用实例来说,取决于对计算速率和计算结果的相对取舍。总体上来说,本申请提供的方法,在现有技术的基础上加快了计算速度,并能够对大规模改造后形成的多尺度裂缝能够进行有效的表征。
本申请实施例提供了一种页岩油藏压裂缝网自适应半无限导流表征装置200,如图14所示该装置包括:嵌入模块201,分割计算模块201,压力值标注模块203,无限导流裂缝确定模块204,表征模块确定模块205。
嵌入模块201用于将嵌入式离散裂缝嵌入多重介质模型中。
分割计算模块202用于对多重介质模型进行分割,确定背景网络和每一个分割单元的质量守恒方程。
压力值标注模块203用于根据质量守恒方程和多重介质模型的边界条件,计算分割单元的压力,并在多重介质模型上进行标注;压力值标注模块203还包括封闭边界条件和井控条件;封闭边界条件是:在某个方向上,多重介质模型与外部没有能量交换,表示为:u(x,y)|Ω·n=0;其中,n表示边界处的法向量,u(x,y)|Ω表示为边界处的流量速度。
无限导流裂缝确定模块204用于根据压力误差阈值和无限导流误差阈值搜索多重介质模型,确定无限导流裂缝和连通子网络;多重介质模型根据以下公式来计算压力误差阈值;其中,Pi表示第i个分割单元的压力值,Pj表示第j个分割单元的压力值,ep表示第i个分割单元和第j个分割单元的压力误差值,ΔPij表示第i个分割单元和第j个分割单元的压力差值,εP表示压力误差阈值。
无限导流裂缝确定模块204根据以下公式来计算无限导流误差阈值;其中,CSN表示连通子网络,Tmf表示背景网络中分割单元与分割单元内的嵌入式离散裂缝之间的传导率,表示连通子网络的平均压力值与背景网络中对应分割单元的压力值的差值,εq表示无限导流误差阈值,e表示无限导流误差,ΔPmf表示连通子网络中分割单元的压力值与背景网络中对应的压力值的差值。
无限导流裂缝确定模块204通过以下步骤确定连通子网络:根据压力误差阈值和无限导流误差阈值对多重介质模型进行搜索后,将符合条件的裂缝进行合并,确定连通子网络。
表征模块确定模块205用于计算连通子网络与背景网格之间的传导率,确定半无限导流表征模型。表征模块确定模块205具体采用网格等效距离方法计算传导率。
本发明实施例提供了一种页岩油藏压裂缝网自适应半无限导流表征服务器,如图15所示,该服务器包括存储器301和处理器302;存储器301用于存储计算机可执行指令;处理器302用于执行计算机可执行指令,以实现以上所述的方法。
本发明实施例提供了一种计算机可读存储介质,计算机可读存储介质存储有可执行指令,计算机执行可执行指令时能够实现以上所述的方法。
上述存储介质包括但不限于随机存取存储器(英文:Random Access Memory;简称:RAM)、只读存储器(英文:Read-Only Memory;简称:ROM)、缓存(英文:Cache)、硬盘(英文:Hard Disk Drive;简称:HDD)或者存储卡(英文:Memory Card)。所述存储器可以用于存储计算机程序指令。
虽然本申请提供了如实施例或流程图所述的方法操作步骤,但基于常规或者无创造性的劳动可以包括更多或者更少的操作步骤。本实施例中列举的步骤顺序仅仅为众多步骤执行顺序中的一种方式,不代表唯一的执行顺序。在实际中的装置或客户端产品执行时,可以按照本实施例或者附图所示的方法顺序执行或者并行执行(例如并行处理器或者多线程处理的环境)。
上述实施例阐明的装置或模块,具体可以由计算机芯片或实体实现,或者由具有某种功能的产品来实现。为了描述的方便,描述以上装置时以功能分为各种模块分别描述。在实施本申请时可以把各模块的功能在同一个或多个软件和/或硬件中实现。当然,也可以将实现某功能的模块由多个子模块或子单元组合实现。
本申请中所述的方法、装置或模块可以以计算机可读程序代码方式实现控制器按任何适当的方式实现,例如,控制器可以采取例如微处理器或处理器以及存储可由该(微)处理器执行的计算机可读程序代码(例如软件或固件)的计算机可读介质、逻辑门、开关、专用集成电路(英文:Application Specific Integrated Circuit;简称:ASIC)、可编程逻辑控制器和嵌入微控制器的形式,控制器的例子包括但不限于以下微控制器:ARC 625D、Atmel AT91SAM、Microchip PIC18F26K20以及Silicone Labs C8051F320,存储器控制器还可以被实现为存储器的控制逻辑的一部分。
本申请所述装置中的部分模块可以在由计算机执行的计算机可执行指令的一般上下文中描述,例如程序模块。一般地,程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构、类等。也可以在分布式计算环境中实践本申请,在这些分布式计算环境中,由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中,程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。
通过以上的实施方式的描述可知,本领域的技术人员可以清楚地了解到本申请可借助软件加必需的硬件的方式来实现。基于这样的理解,本申请的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来,也可以通过数据迁移的实施过程中体现出来。该计算机软件产品可以存储在存储介质中,如ROM/RAM、磁碟、光盘等,包括若干指令用以使得一台计算机设备(可以是个人计算机,移动终端,服务器,或者网络设备等)执行本申请各个实施例或者实施例的某些部分所述的方法。
本说明书中的各个实施方式采用递进的方式描述,各个实施方式之间相同或相似的部分互相参见即可,每个实施方式重点说明的都是与其他实施方式的不同之处。本申请的全部或者部分可用于众多通用或专用的计算机系统环境或配置中。例如:个人计算机、服务器计算机、手持设备或便携式设备、平板型设备、移动通信终端、多处理器系统、基于微处理器的系统、可编程的电子设备、网络PC、小型计算机、大型计算机、包括以上任何系统或设备的分布式计算环境等等。
以上实施例仅用以说明本申请的技术方案,而非对本申请限制;尽管参照前述实施例对本申请进行了详细的说明,本领域普通技术人员应当理解:其依然可以对前述实施例所记载的技术方案进行修改,或者对其中部分或者全部技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本申请技术方案的范围。
Claims (8)
1.一种页岩油藏压裂缝网自适应半无限导流表征方法,其特征在于,包括:
将嵌入式离散裂缝嵌入多重介质模型中;
对所述多重介质模型进行分割,确定背景网络和每一个分割单元的质量守恒方程;
根据所述质量守恒方程和所述多重介质模型的边界条件,计算每一个所述分割单元的压力,并在所述多重介质模型上进行标注;
根据压力误差阈值和无限导流误差阈值搜索所述多重介质模型,确定无限导流裂缝和连通子网络,所述无限导流误差阈值包括:
其中,CSN表示所述连通子网络,Tmf表示所述背景网络中所述分割单元与所述分割单元内的嵌入式离散裂缝之间的传导率,表示所述连通子网络的平均压力值与所述背景网络中对应所述分割单元的压力值的差值,εq表示无限导流误差阈值,e表示无限导流误差,ΔPmf表示所述连通子网络中所述分割单元的压力值与所述背景网络中对应的压力值的差值;
计算所述连通子网络与背景网格之间的传导率,确定半无限导流表征模型。
2.根据权利要求1所述的方法,其特征在于,所述多重介质模型的边界条件包括:封闭边界条件和井控条件;
所述封闭边界条件是:在某个方向上,所述多重介质模型与外部没有能量交换,表示为:u(x,y)|Ω·n=0;其中,n表示边界处的法向量,u(x,y)|Ω表示为边界处的流量速度。
4.根据权利要求1所述的方法,其特征在于,通过以下步骤确定所述连通子网络:
根据所述压力误差阈值和所述无限导流误差阈值对所述多重介质模型进行搜索后,将符合条件的裂缝进行合并,确定所述连通子网络。
5.根据权利要求1所述的方法,其特征在于,所述计算所述连通子网络与所述背景网格之间的传导率,包括:采用网格等效距离方法计算所述传导率。
6.一种页岩油藏压裂缝网自适应半无限导流表征装置,其特征在于,包括:
嵌入模块,用于将嵌入式离散裂缝嵌入多重介质模型中;
分割计算模块,用于对所述多重介质模型进行分割,确定背景网络和每一个分割单元的质量守恒方程;
压力值标注模块,用于根据所述质量守恒方程和所述多重介质模型的边界条件,计算所述分割单元的压力,并在所述多重介质模型上进行标注;
无限导流裂缝确定模块,用于根据压力误差阈值和无限导流误差阈值搜索所述多重介质模型,确定无限导流裂缝和连通子网络,所述无限导流误差阈值包括:
其中,CSN表示所述连通子网络,Tmf表示所述背景网络中所述分割单元与所述分割单元内的嵌入式离散裂缝之间的传导率,表示所述连通子网络的平均压力值与所述背景网络中对应所述分割单元的压力值的差值,εq表示无限导流误差阈值,e表示无限导流误差,ΔPmf表示所述连通子网络中所述分割单元的压力值与所述背景网络中对应的压力值的差值;
表征模块确定模块,用于计算所述连通子网络与背景网格之间的传导率,确定半无限导流表征模型。
7.一种页岩油藏压裂缝网自适应半无限导流表征服务器,其特征在于,包括存储器和处理器;
所述存储器用于存储计算机可执行指令;
所述处理器用于执行所述计算机可执行指令,以实现权利要求1-5任一项所述的方法。
8.一种计算机可读存储介质,其特征在于,所述计算机可读存储介质存储有可执行指令,计算机执行所述可执行指令时能够实现如权利要求1-5任一项所述的方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210765504.7A CN114818551B (zh) | 2022-07-01 | 2022-07-01 | 一种页岩油藏压裂缝网自适应半无限导流表征方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210765504.7A CN114818551B (zh) | 2022-07-01 | 2022-07-01 | 一种页岩油藏压裂缝网自适应半无限导流表征方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114818551A CN114818551A (zh) | 2022-07-29 |
CN114818551B true CN114818551B (zh) | 2022-11-04 |
Family
ID=82523160
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210765504.7A Active CN114818551B (zh) | 2022-07-01 | 2022-07-01 | 一种页岩油藏压裂缝网自适应半无限导流表征方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114818551B (zh) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107169227A (zh) * | 2017-06-09 | 2017-09-15 | 中国石油大学(华东) | 一种分段压裂水平井的粗网格模拟方法及系统 |
CN108319756A (zh) * | 2017-12-29 | 2018-07-24 | 西安石油大学 | 一种致密储层体积压裂缝网扩展模拟及表征方法 |
CN111062129A (zh) * | 2019-12-16 | 2020-04-24 | 中国石油大学(华东) | 页岩油复杂缝网离散裂缝连续介质混合数值模拟方法 |
CN111104766A (zh) * | 2019-12-16 | 2020-05-05 | 中国石油大学(华东) | 基于离散裂缝模型的油水两相非达西渗流数值模拟方法 |
CN112084718A (zh) * | 2020-09-16 | 2020-12-15 | 西南石油大学 | 基于渗流差异的页岩气藏单相气体三孔三渗模型构建方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP3108088B1 (en) * | 2014-02-19 | 2018-08-22 | Repsol, S.A. | Method implemented in a computer for the numerical simulation of a porous medium |
CN108595782B (zh) * | 2018-04-03 | 2020-09-04 | 中国石油大学(北京) | 一种离散裂缝中基质与裂缝间的传质计算方法 |
CN109408855B (zh) * | 2018-09-03 | 2023-05-09 | 中国石油大学(华东) | 一种待压裂水平井的裂缝导流能力值计算方法和装置 |
CN109992864B (zh) * | 2019-03-22 | 2020-02-04 | 成都理工大学 | 非常规双重介质储层体积压裂数值模拟及参数优化方法 |
-
2022
- 2022-07-01 CN CN202210765504.7A patent/CN114818551B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107169227A (zh) * | 2017-06-09 | 2017-09-15 | 中国石油大学(华东) | 一种分段压裂水平井的粗网格模拟方法及系统 |
CN108319756A (zh) * | 2017-12-29 | 2018-07-24 | 西安石油大学 | 一种致密储层体积压裂缝网扩展模拟及表征方法 |
CN111062129A (zh) * | 2019-12-16 | 2020-04-24 | 中国石油大学(华东) | 页岩油复杂缝网离散裂缝连续介质混合数值模拟方法 |
CN111104766A (zh) * | 2019-12-16 | 2020-05-05 | 中国石油大学(华东) | 基于离散裂缝模型的油水两相非达西渗流数值模拟方法 |
CN112084718A (zh) * | 2020-09-16 | 2020-12-15 | 西南石油大学 | 基于渗流差异的页岩气藏单相气体三孔三渗模型构建方法 |
Non-Patent Citations (5)
Title |
---|
Dynamic and static comprehensive prediction method of natural fractures in fractured oil reservoirs: A case study of Triassic Chang 63 reservoirs in Huaqing Oilfield, Ordos Basin, NW China;Su, H (Su Hao)等;《Petroleum Exploration and Development》;20171231;第972-982页 * |
储集层多尺度裂缝高效数值模拟模型;苏皓;《石油学报》;20190531;第587-593页 * |
多尺度嵌入式离散裂缝模型模拟方法;张庆福等;《计算力学学报》;20180815(第04期);第112-118页 * |
页岩气藏压裂缝网扩展流动一体化模拟技术;盛广龙等;《西南石油大学学报(自然科学版)》;20211031;第84-96页 * |
页岩气藏数值模拟模型研究进展;王强等;《中国地质》;20191215(第06期);第38-53页 * |
Also Published As
Publication number | Publication date |
---|---|
CN114818551A (zh) | 2022-07-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
AU2011289788B2 (en) | Reservoir upscaling method with preserved transmissibility | |
CN109829217B (zh) | 压裂性裂缝油藏产能模拟方法及装置 | |
Huang et al. | Well-balanced finite volume scheme for shallow water flooding and drying over arbitrary topography | |
CA2961562C (en) | Formation fracture flow monitoring | |
Garzon et al. | Improvements of explicit crack surface representation and update within the generalized finite element method with application to three‐dimensional crack coalescence | |
CN113836695B (zh) | 一种基于无网格连接元的油藏数值模拟方法 | |
Zhang | Finite element generation of arbitrary 3-D fracture networks for flow analysis in complicated discrete fracture networks | |
Liang | A simplified adaptive Cartesian grid system for solving the 2D shallow water equations | |
Xue et al. | A fast numerical method and optimization of 3D discrete fracture network considering fracture aperture heterogeneity | |
CA2964250A1 (en) | Junction models for simulating proppant transport in dynamic fracture networks | |
CN115577562A (zh) | 一种裂缝油藏井位优化方法 | |
CN116108617A (zh) | 一种页岩油藏局部网格加密方法、装置及存储介质 | |
CN115146446A (zh) | 基于近似梯度算法和嵌入式离散裂缝模型的油藏优化方法 | |
CA2963485A1 (en) | Formation fracture flow monitoring | |
CN114818551B (zh) | 一种页岩油藏压裂缝网自适应半无限导流表征方法及装置 | |
El-Ghandour et al. | Groundwater management using a new coupled model of flow analytical solution and particle swarm optimization | |
CN114718556A (zh) | 人工裂缝参数的获取方法、装置及设备 | |
CN116167302A (zh) | 一种天然气水合物增产模拟中人工复杂裂缝的描述方法 | |
CN114970392A (zh) | 一种裂缝分形表征方法、系统及可读存储介质 | |
CN114460639B (zh) | 页岩油储层渗透率的预测方法及装置 | |
Blonsky et al. | Development of acidizing simulator for sandstone reservoirs | |
Kurbonov | 3D model and numerical algorithm for gas filtration in porous media | |
CN111852461B (zh) | 一种模拟致密油藏基岩与裂缝间非稳态窜流的方法 | |
CN113919247A (zh) | 一种复杂裂缝网络性油藏流动模拟方法 | |
Peng et al. | Mixed numerical method for bed evolution |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |