CN109992864B - 非常规双重介质储层体积压裂数值模拟及参数优化方法 - Google Patents
非常规双重介质储层体积压裂数值模拟及参数优化方法 Download PDFInfo
- Publication number
- CN109992864B CN109992864B CN201910223536.2A CN201910223536A CN109992864B CN 109992864 B CN109992864 B CN 109992864B CN 201910223536 A CN201910223536 A CN 201910223536A CN 109992864 B CN109992864 B CN 109992864B
- Authority
- CN
- China
- Prior art keywords
- fracture
- crack
- unconventional
- natural
- fractures
- 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
- 238000000034 method Methods 0.000 title claims abstract description 28
- 238000004088 simulation Methods 0.000 title claims abstract description 22
- 238000005457 optimization Methods 0.000 title claims abstract description 9
- 238000010276 construction Methods 0.000 claims abstract description 21
- 230000008569 process Effects 0.000 claims abstract description 18
- 239000011159 matrix material Substances 0.000 claims abstract description 17
- 230000000694 effects Effects 0.000 claims abstract description 10
- 230000000977 initiatory effect Effects 0.000 claims abstract description 8
- 238000013178 mathematical model Methods 0.000 claims abstract description 7
- 238000011160 research Methods 0.000 claims abstract description 7
- 239000004568 cement Substances 0.000 claims abstract description 6
- 230000002860 competitive effect Effects 0.000 claims abstract description 6
- 206010017076 Fracture Diseases 0.000 claims description 196
- 208000010392 Bone Fractures Diseases 0.000 claims description 129
- 239000012530 fluid Substances 0.000 claims description 38
- 230000008878 coupling Effects 0.000 claims description 21
- 238000010168 coupling process Methods 0.000 claims description 21
- 238000005859 coupling reaction Methods 0.000 claims description 21
- 238000006073 displacement reaction Methods 0.000 claims description 16
- 230000009977 dual effect Effects 0.000 claims description 12
- 230000006399 behavior Effects 0.000 claims description 11
- 239000007788 liquid Substances 0.000 claims description 11
- 239000007787 solid Substances 0.000 claims description 8
- 239000011435 rock Substances 0.000 claims description 6
- 208000006670 Multiple fractures Diseases 0.000 claims description 5
- 238000004364 calculation method Methods 0.000 claims description 5
- 239000000126 substance Substances 0.000 claims description 5
- 238000000342 Monte Carlo simulation Methods 0.000 claims description 3
- 238000004422 calculation algorithm Methods 0.000 claims description 3
- 238000000926 separation method Methods 0.000 claims description 3
- 238000001914 filtration Methods 0.000 claims description 2
- 239000011148 porous material Substances 0.000 abstract description 8
- 230000010354 integration Effects 0.000 abstract description 5
- 238000005094 computer simulation Methods 0.000 abstract 1
- 230000001360 synchronised effect Effects 0.000 abstract 1
- 239000007789 gas Substances 0.000 description 8
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 8
- 230000007246 mechanism Effects 0.000 description 6
- 238000010586 diagram Methods 0.000 description 4
- 239000003345 natural gas Substances 0.000 description 4
- 239000002253 acid Substances 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000035515 penetration Effects 0.000 description 2
- 239000003209 petroleum derivative Substances 0.000 description 2
- 239000000243 solution Substances 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 230000000903 blocking effect Effects 0.000 description 1
- 238000009933 burial Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000001808 coupling effect Effects 0.000 description 1
- 238000005336 cracking Methods 0.000 description 1
- 238000013016 damping Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 230000000149 penetrating effect Effects 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 238000005086 pumping Methods 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000010008 shearing Methods 0.000 description 1
Images
Classifications
-
- G01V20/00—
-
- 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
-
- 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]
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/64—Geostructures, e.g. in 3D data cubes
- G01V2210/646—Fractures
Abstract
本发明公开了一种非常规双重介质储层体积压裂数值模拟及参数优化方法,步骤如下:基于双重介质孔隙弹性理论,考虑裂缝间的摩擦效应,发展基于显式时间积分的水力裂缝粘弹塑性损伤模型;采用在计算模型网格内部边界嵌入零厚度裂缝单元的方法模拟水力裂缝遇天然裂缝的随机相交和分岔,建立非常规双重介质储层体积压裂水力裂缝扩展的数学模型;编制非常规储层体积压裂复杂多裂缝破裂与竞争扩展的有限元程序,建立含天然裂缝的套管‑水泥环‑簇内射孔孔眼‑储层基质的水力压裂有限元模型;开展簇内近井地带复杂裂缝起裂与扩展规律研究,开展分段多簇水力主裂缝同步扩展规律研究,优化体积压裂施工工艺。本发明模拟因素全面,优化了油气开采施工方案。
Description
技术领域
本发明涉及石油与天然气开发技术领域,特别是涉及一种非常规双重介质储层体积压裂数值模拟及参数优化方法。
背景技术
长期以来,我国能源供需矛盾突出,石油天然气对外依存度逐年攀升,2017年我国石油对外依存度达到67.4%,天然气达到39.4%,大大超过国际公认的警戒线,严重威胁我国能源安全。据美国能源信息署EIA的全球调查显示,我国页岩气储量为36.1万亿立方米,排名全球第一位。作为一种典型的非常规天然气资源,我国页岩气资源开发潜力巨大。我国页岩气藏普遍埋藏较深,特别是四川盆地的五峰-龙马溪组平均埋深超3000m,部分已达到5000m,呈现地质构造复杂、页理发育、高应力和高强度等典型特征,致使储层施工压力高达80-110MPa。根据涪陵页岩气田近三年在一期浅层(垂深2000-3000m)的开发实践发现:非常规双重介质储层水平井分段多簇射孔体积压裂过程中,储层破裂压力和裂缝延伸压力高,常遇到施工压力超限储层压不开、近井地带裂缝形态复杂而引起砂堵的恶性事故。
近年来,国内外学者针对非常规双重介质储层气体积压裂复杂裂缝的起裂和扩展行为开展了大量的室内实验和数值模拟研究。然而,这些研究均假设多簇射孔压裂后,在每一个射孔簇处形成一条主缝,各射孔簇主裂缝共同扩展,研究重点集中在各簇主裂缝间相互干扰作用下的应力场和裂缝形态,少有学者研究从射孔簇内射孔孔眼起裂的多条裂缝在近井地带的扩张、扩展、相交和分岔行为。同时,非常规双重介质储层通常含有大量的天然裂缝,是典型的双重介质系统,而现有水力裂缝扩展模型中普遍将非常规双重介质储层处理为常规单重介质,甚至是不可渗透介质,这会对压裂过程中主裂缝附近流体压力分布的计算准确度产生影响,进而影响水力裂缝和天然裂缝相遇后以及主裂缝周围尚未连通天然裂缝的起裂与扩展行为。
发明内容
本发明的目的在于克服现有技术的缺点,考虑到非常规双重介质储层气体积压裂施工压力高、多数射孔簇内水力裂缝未得到有效扩展的科学难题,考虑非常规双重介质储层应力场、渗流场和裂缝内流场的相互耦合,构建显式时间积分的裂缝粘弹塑性损伤模型及网格单元,开发离散天然裂缝生成程序以及裂缝单元嵌入程序,创建非常规双重介质储层双重介质渗流-应力-损伤耦合的复杂裂缝交错扩展模型,揭示天然裂缝展布及参数、射孔参数、液体类型和排量提升方式等对簇内多条水力裂缝和天然裂缝交错扩展的排斥、吸引、分岔和相交影响机理,找出有利于降低近井摩阻、以及多裂缝扩展并串联形成一条优势主裂缝的簇内射孔方式、射孔参数、施工参数和工程地质条件,进而提供一种非常规双重介质储层体积压裂数值模拟及参数优化方法,它包括以下步骤:
S1.建立水力压裂裂缝扩展的数学模型,分析裂缝延伸的诱导应力场模型,发展基于显式时间积分求解策略的双重介质渗流-应力耦合模型及算法;
S2.针对裂缝的变形,考虑裂缝面之间的摩擦效应,发展显式时间积分的裂缝粘弹塑性损伤模型及网格单元;
S3.在非常规储层基质网格边界嵌入零厚度裂缝单元模拟水力裂缝遇天然裂缝的随机相交和分岔;
S4.编制井眼-射孔孔眼-多裂缝压裂液竞争分配程序,实现多个孔缝(包括射孔孔眼和裂缝)间压裂液的动态分配过程;
S5.建立非常规储层体积压裂复杂多裂缝破裂与竞争扩展的数值模型,并编制复杂多裂缝交错扩展程序;
S6.基于实际储层工程地质条件,建立含天然裂缝的套管-水泥环-簇内射孔孔眼-储层基质的网格模型,研究不同射孔方式(包括螺旋射孔、定向射孔和定面射孔)、射孔参数(包括孔径、穿深、孔密和相位角)、井眼参数(包括井斜角和方位角)和储层岩石力学特征(包括脆性-延性和酸损伤)等条件下,簇内孔眼的破裂位置、天然裂缝的开启、近井地带及射孔孔眼压力分布等,优化射孔方式、射孔参数、压裂初始阶段施工参数(包括排量、液量和泵注程序等)和压裂液物理化学性能;
S7.基于实际储层工程地质条件,针对非常规储层基质-天然裂缝双重介质,耦合天然裂缝和多条水力主裂缝,开展簇内近井地带多条水力主裂缝遇天然裂缝交错延伸的数值模拟,研究近井筒水力裂缝和天然裂缝交错扩展的排斥、吸引、分岔和相交规律,优化射孔方式及参数、液体类型(即压裂液物理化学性质)和排量提升方式等体积压裂施工工艺参数;
S8.基于实际储层工程地质条件,针对非常规储层基质-天然裂缝双重介质,耦合天然裂缝和多条水力主裂缝,开展多簇段的多条水力主裂缝遇天然裂缝交错延伸的数值模拟,研究多条水力主裂缝之间及其和天然裂缝交错扩展的排斥、吸引、分岔和相交规律;计算不同施工参数情况下缝网的演化规律,以形成最大总裂缝面积为目标,对体积压裂施工工艺参数进行优化。
在步骤S2中,显式时间积分的裂缝粘弹塑性损伤模型及网格单元:针对裂缝的变形,考虑裂缝面之间的摩擦效应,发展新的内聚力本构模型描述裂缝单元的变形破坏行为,构建零厚度、考虑裂缝面摩擦效应的裂缝粘弹塑性损伤模型,裂缝的起裂和扩展服从Traction-separation损伤准则;针对裂缝交叉和分岔扩展时有限元计算不易收敛的难题,推导三角形裂缝单元的显式时间积分有限元方程列式,并编制裂缝单元的计算程序;裂缝网格单元结构上、下表面节点具有流体压力(包括基质孔隙压力和裂隙孔隙压力)自由度和位移自由度,用以计算裂缝面的变形以及压裂液的滤失,中间层节点只具有流体压力自由度,用来计算裂缝内流体的流动。
在步骤S3中,水力裂缝遇天然裂缝的相交和分岔扩展数值实现:根据储层天然裂缝分布规律,采用蒙特卡洛方法生成离散天然裂缝网络,建立含离散天然裂缝的有限元几何模型,对整个研究区域进行网格剖分,并开发在所有单元间嵌入流固耦合裂缝单元的程序,在所有网格的边界上插入裂缝单元,用粘弹塑性损伤裂缝单元统一表征水力裂缝和天然裂缝;因此,水力裂缝能够沿着任意单元的边界扩展,从而实现多裂缝的相交和分岔扩展行为,并且无需重新划分网格;表征水力裂缝和天然裂缝的裂缝单元本质上是相同的,只是力学参数和水力参数大小上存在区别,裂缝单元和非常规双重介质储层基质单元之间通过共用节点的方式耦合。
本发明具有以下优点:
1.对整个研究区域进行网格剖分,并开发在所有单元间嵌入流固耦合裂缝单元的程序,在所有网格的边界上插入裂缝单元,用粘弹塑性损伤裂缝单元统一表征水力裂缝和天然裂缝。因此,水力裂缝能够沿着任意单元的边界扩展,从而实现多裂缝的相交和分岔扩展行为,并且无需重新划分网格。
2.考虑射孔簇附近非常规双重介质储层天然裂缝的影响,建立非常规双重介质储层体积压裂射孔簇内孔眼渗流-应力-损伤耦合模型,揭示非常规双重介质储层射孔簇内多裂缝的起裂机理。
3.综合考虑天然裂缝和多条水力主裂缝的相互耦合作用,揭示射孔簇内多条水力主裂缝与天然裂缝的排斥、吸引、分岔和相交的竞争扩展机理。
附图说明
以下附图旨在对本发明做示意性说明和解释,并不限定发明的范围,其中:
图1是本发明所示数值模拟及施工法流程示意图;
图2是裂缝粘弹塑性损伤本构模型示意图;
图3是裂缝网格单元结构示意图;
图4是储层天然裂缝分布,有限元网格划分策略,基质、裂缝单元通过共用节点耦合示意图;
图5是井眼-射孔孔眼-多裂缝压裂液竞争分配模型示意图;
图6是分段多簇压裂复杂裂缝网络计算结果。
具体实施方式
下面结合附图对本发明做进一步的描述,但本发明的保护范围不局限于以下所述。
如图1-6所示,一种非常规双重介质储层体积压裂数值模拟及参数优化方法,包括以下步骤:
S1.建立水力压裂裂缝扩展的数学模型,分析裂缝延伸的诱导应力场模型,发展基于显式时间积分求解策略的双重介质渗流-应力耦合模型及算法,即:根据Biot孔隙弹性理论,推导非常规双重介质储层基质骨架的应力平衡方程、几何方程和Biot有效应力模型,建立非常规双重介质储层的基质骨架变形控制方程。假设非常规双重介质储层孔隙流体流动服从达西定律,结合质量守恒方程,推导出非常规双重介质储层的渗流控制方程。联立基质骨架变形控制方程、渗流控制方程和渗透率应力敏感性张量,结合水平井分段多簇压裂簇内射孔的流体和应力边界条件,建立基于显式时间积分求解策略的非常规双重介质储层流固耦合数学模型,对压裂储层几何模型进行网格离散,并在单元的节点上引入未知位移和压力自由度,推导非常规双重介质储层流固耦合数学模型的有限元方程列式,开发流固耦合单元的显式时间积分计算程序;
S2.针对裂缝的变形,考虑裂缝面之间的摩擦效应,发展显式时间积分的裂缝粘弹塑性损伤模型及网格单元,即:针对裂缝的变形,考虑裂缝面之间的摩擦效应,发展新的内聚力本构模型描述裂缝单元的变形破坏行为,构建零厚度、考虑裂缝面摩擦效应的裂缝粘弹塑性损伤模型,其力学本构模型(如图2所示,图中:是储层岩石的抗拉强度(剪切强度),MPa;τf是裂缝面的摩擦应力,MPa;Kn(Kt)是拉伸(剪切)状态的界面硬化系数;是初始损伤位移,m;是完全破坏位移,m)中裂缝的起裂和扩展服从Traction-separation损伤准则;针对裂缝交叉和分岔扩展时有限元计算不易收敛的难题,推导三角形裂缝单元的显式时间积分有限元方程列式,并编制裂缝单元的计算程序;裂缝网格单元(如图3所示)的上、下表面节点具有流体压力自由度(包括基质孔隙压力及裂隙孔隙压力)和位移自由度,用以计算裂缝面的变形以及压裂液的滤失,中间层节点只具有流体压力自由度,用来计算裂缝内流体的流动;
S3.在非常规储层基质网格边界嵌入零厚度裂缝单元模拟水力裂缝遇天然裂缝的随机相交和分岔,即:根据储层天然裂缝分布规律,采用蒙特卡洛方法生成离散天然裂缝网络,建立含离散天然裂缝的有限元几何模型(如图4a所示),采用三维四面体网格对整个研究区域进行剖分,并开发在所有单元间嵌入流固耦合裂缝单元的程序,在所有网格的边界上插入裂缝单元(如图4b所示),用粘弹塑性损伤裂缝单元统一表征水力裂缝和天然裂缝;因此,水力裂缝能够沿着任意单元的边界扩展,从而实现多裂缝的相交和分岔扩展行为,并且无需重新划分网格,表征水力裂缝和天然裂缝的裂缝单元本质上是相同的,只是力学参数和水力参数大小上存在区别;裂缝单元和非常规双重介质储层基质单元之间通过共用节点的方式耦合(如图4c所示,图中:{u,pm,pf}是表面节点;{p}是中间层节点);
S4.编制井眼-射孔孔眼-多裂缝压裂液竞争分配程序,实现多个孔缝(包括射孔孔眼和裂缝)间压裂液的动态分配过程,即在压裂施工过程中,压裂液从井眼通过射孔孔眼进入多条水力裂缝(如图5所示,图中:Qk是进入第k条裂缝的流量,m3/s;是第k条裂缝在套管里的流体压力,Pa;是第k条裂缝套管外面紧邻套管处流体压力,Pa;是第k条裂缝射孔簇压降,Pa),进入各射孔孔眼的压裂液流量和各射孔孔眼外裂缝的起裂与扩展状况动态相关;根据基尔霍夫定律,将井筒内流体的流动视为一维流动问题,建立井筒流量和压力平衡方程;根据质量守恒定律,注入裂缝内流体的总质量等于所有裂缝内流体的质量和法向滤失量之和,建立多条裂缝内流体的质量守恒方程;将井筒流量和压力平衡方程、多条裂缝内流体的质量守恒方程联立,建立井眼-射孔孔眼-多裂缝压裂液竞争分配模型,通过开发一维管单元刻画压裂液流经射孔孔眼时的摩阻,将多个一维管单元按照图5的方式在一端通过共用注入点进行组合,实现多个射孔孔眼和裂缝间压裂液的动态分配过程;
S5.建立非常规储层体积压裂复杂多裂缝破裂与竞争扩展的数值模型,并编制复杂多裂缝交错扩展程序,即:对建立的有限元方程组在空间域和时间域上进行离散编程求解,采用显式时间积分方案克服水力裂缝扩展(尤其是复杂裂缝扩展)过程中强烈非线性导致的收敛困难问题:固体变形动力学方程在时间域上采用中心时间差分的积分方案离散,流体动力学方程在时间域上采用向前时间差分的积分方案离散,同时,将人工阻尼技术(可提高时间积分步长)与多核并行计算技术相结合,提高模型的求解效率;
S6.基于实际储层工程地质条件,建立含天然裂缝的套管-水泥环-簇内射孔孔眼-储层基质的网格模型,研究不同射孔方式(包括螺旋射孔、定向射孔和定面射孔)、射孔参数(包括孔径、穿深、孔密和相位角)、井眼参数(包括井斜角和方位角)和储层岩石力学特征(包括脆性-延性和酸损伤)条件下,簇内孔眼的破裂位置、天然裂缝的开启、近井地带及射孔孔眼压力分布等,优化射孔方式、射孔参数、压裂初始阶段施工参数(包括排量、液量和泵注程序)和压裂液物理化学性能等,即:根据非常规双重介质储层水平井分段多簇压裂的工艺特征,基于实际井眼尺寸、单簇射孔的射孔参数、非常规双重介质储层岩石力学参数及其本构模型,考虑套管-水泥环-簇内射孔孔眼的结构特征,建立非常规双重介质储层水平井分段多簇压裂簇内射孔的几何模型,并对该模型进行离散;
S7.在“非常规双重介质储层体积压裂射孔簇内孔眼渗流-应力-损伤耦合的数值模拟模型”的基础上,仅考虑单条水力主裂缝自孔眼扩展的情形,开展水力主裂缝与非常规双重介质储层层理界面/天然裂缝交错延伸的数值模拟,揭示水力裂缝遇非常规双重介质储层层理和高角度天然裂缝时的停止、穿透和分岔等非平面扩展机理,找出形成这三种相交模式的施工工艺参数和工程地质条件;忽略非常规双重介质储层的天然裂缝,仅模拟簇内多条水力主裂缝的排斥、吸引和相交的竞争扩展过程,研究射孔参数、液体类型和排量提升方式等对非平面裂缝形貌、流量分配、缝内压降等参数的影响,揭示射孔参数和压裂施工参数对裂缝非均衡扩展形态,以及多裂缝诱导应力场的相互干扰机理;针对非常规双重介质储层基质-天然裂缝双重介质,耦合非常规双重介质储层天然裂缝和多条水力主裂缝,考虑不同天然裂缝展布及几何参数、射孔方式及参数、液体类型和排量提升方式等条件,开展簇内近井地带多条水力主裂缝遇天然裂缝交错延伸的数值模拟,研究多条水力裂缝在远离井眼过程中水力裂缝与天然裂缝相遇、交错延伸和分离的模式,研究多裂缝相互干扰的非平面裂缝形态、井眼-裂缝内部流量分配、缝内净压力以及缝间应力场的变化规律,优化射孔方式、参数、液体类型和排量提升方式等体积压裂施工工艺参数;
S8.基于实际储层工程地质条件,针对非常规储层基质-天然裂缝双重介质,耦合天然裂缝和多条水力主裂缝,开展多簇段间多条水力主裂缝遇天然裂缝交错延伸的数值模拟,研究射孔簇间多条水力主裂缝和天然裂缝交错扩展的排斥、吸引、分岔和相交规律,优化液体类型和排量提升方式等体积压裂施工工艺参数;即:根据储层的天然裂缝分形描述结果,建立含天然裂缝的套管-水泥环-簇内射孔孔眼-非常规双重介质储层基质的网格模型;利用编制的裂缝单元嵌入程序,在模型所有网格的边界上插入显式时间积分的裂缝单元;通过天然裂缝面上接触力和摩擦力的数学模型,计算天然裂缝面的位移,考虑天然裂缝的剪切滑移和扩张行为,采用最大拉应力准则判断簇内孔眼的破裂,采用塑性损伤模拟裂缝的扩展,建立双重介质体积压裂射孔簇内孔眼渗流-应力-损伤耦合的数值模拟模型,实现簇内多裂缝的起裂与扩展模拟(如图6所示,图中:横轴是缝宽,m;纵轴是累计缝长,m;Ltotal是裂缝总长度,m;LNF是天然裂缝长度,m;LHF是水力裂缝长度,m);应用所编制的“复杂裂缝交错扩展程序”,研究不同射孔方式、射孔参数、井眼参数和储层岩石力学特征条件下,簇内孔眼的破裂位置、天然裂缝的开启、近井地带及射孔孔眼压力分布等,优化射孔方式、射孔参数、压裂初始阶段施工参数和压裂液物理化学性能等。
Claims (3)
1.一种非常规双重介质储层体积压裂数值模拟及参数优化方法,其特征在于,它包括以下步骤:
S1.建立非常规双重介质储层体积压裂水力压裂裂缝扩展的数学模型,发展基于显式时间积分求解策略的渗流-应力耦合模型及算法;
S2.针对裂缝的变形,考虑裂缝面之间的摩擦效应,发展显式时间积分的裂缝粘弹塑性损伤模型及网格单元;
S3.在非常规储层基质网格边界嵌入零厚度裂缝单元模拟水力裂缝遇天然裂缝的随机相交和分岔;
S4.编制井眼-射孔孔眼-多裂缝压裂液竞争分配程序,实现多个孔缝间压裂液的动态分配过程;
S5.建立非常规储层体积压裂复杂多裂缝破裂与竞争扩展的数值模型,并编制复杂多裂缝交错扩展程序;
S6.基于实际储层工程地质条件,建立含天然裂缝的套管-水泥环-簇内射孔孔眼-储层基质的网格模型,研究不同射孔方式、射孔参数、井眼参数和储层岩石力学特征条件下,簇内孔眼的破裂位置、天然裂缝的开启、近井地带及射孔孔眼压力分布,优化射孔方式、射孔参数、压裂初始阶段施工参数和压裂液物理化学性能;
S7.基于实际储层工程地质条件,针对非常规储层基质-天然裂缝双重介质,耦合天然裂缝和多条水力主裂缝,开展簇内近井地带多条水力主裂缝遇天然裂缝交错延伸的数值模拟,研究射孔簇内多条水力主裂缝和天然裂缝交错扩展的排斥、吸引、分岔和相交规律,优化以下体积压裂施工的工艺参数:射孔方式及参数、液体类型和排量提升方式;
S8.基于实际储层工程地质条件,针对非常规储层基质-天然裂缝双重介质,耦合天然裂缝和多条水力主裂缝,开展多簇段的多条水力主裂缝遇天然裂缝交错延伸的数值模拟,研究多条水力主裂缝之间及其和天然裂缝交错扩展的排斥、吸引、分岔和相交规律,优化体积压裂施工工艺参数。
2.根据权利要求1所述的非常规双重介质储层体积压裂数值模拟及参数优化方法,其特征在于:在步骤S2中,针对裂缝的变形,考虑裂缝面之间的摩擦效应,发展新的内聚力本构模型描述裂缝单元的变形破坏行为,构建零厚度、考虑裂缝面摩擦效应的裂缝粘弹塑性损伤模型,裂缝的起裂和扩展服从Traction-separation损伤准则;针对裂缝交叉和分岔扩展时有限元计算不易收敛的难题,推导三角形裂缝单元的显式时间积分有限元方程列式,并编制裂缝单元的计算程序;裂缝网格单元结构上、下表面节点具有流体压力自由度和位移自由度,用以计算裂缝面的变形以及压裂液的滤失,中间层节点只具有流体压力自由度,用来计算裂缝内流体的流动。
3.根据权利要求1所述的非常规双重介质储层体积压裂数值模拟及参数优化方法,其特征在于:在步骤S3中,根据储层天然裂缝分布规律,采用蒙特卡洛方法生成离散天然裂缝网络,建立含离散天然裂缝的有限元几何模型,对整个研究区域进行网格剖分,并开发在所有单元间嵌入流固耦合裂缝单元的程序,在所有网格的边界上插入裂缝单元,用粘弹塑性损伤裂缝单元统一表征水力裂缝和天然裂缝;因此,水力裂缝能够沿着任意单元的边界扩展,从而实现多裂缝的相交和分岔扩展行为,并且无需重新划分网格,表征水力裂缝和天然裂缝的裂缝单元本质相同,仅在力学参数和水力参数大小上存在区别,裂缝单元和非常规双重介质储层基质单元之间通过共用节点的方式耦合。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910223536.2A CN109992864B (zh) | 2019-03-22 | 2019-03-22 | 非常规双重介质储层体积压裂数值模拟及参数优化方法 |
US16/516,261 US11567233B2 (en) | 2019-03-22 | 2019-07-19 | Numerical simulation and parameter optimization method for volumetric fracturing of unconventional dual medium reservoir |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910223536.2A CN109992864B (zh) | 2019-03-22 | 2019-03-22 | 非常规双重介质储层体积压裂数值模拟及参数优化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109992864A CN109992864A (zh) | 2019-07-09 |
CN109992864B true CN109992864B (zh) | 2020-02-04 |
Family
ID=67130886
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910223536.2A Active CN109992864B (zh) | 2019-03-22 | 2019-03-22 | 非常规双重介质储层体积压裂数值模拟及参数优化方法 |
Country Status (2)
Country | Link |
---|---|
US (1) | US11567233B2 (zh) |
CN (1) | CN109992864B (zh) |
Families Citing this family (29)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110865013B (zh) * | 2019-11-20 | 2021-05-11 | 成都理工大学 | 双重煤岩孔隙结构模型的构建方法 |
CN111444612B (zh) * | 2020-03-26 | 2021-04-16 | 北京科技大学 | 一种致密油藏水平井多级压裂流场形态模拟方法 |
CN112096362B (zh) * | 2020-09-27 | 2021-10-22 | 西南石油大学 | 非常规储层多簇射孔竞争起裂与扩展模拟方法及装置 |
CN112345429B (zh) * | 2020-11-02 | 2021-08-06 | 北京科技大学 | 非常规油气藏非均匀地质参数的获取方法及装置 |
CN114458272A (zh) * | 2020-11-02 | 2022-05-10 | 中国石油天然气股份有限公司 | 储层压裂参数的确定方法及装置 |
CN112211607A (zh) * | 2020-11-12 | 2021-01-12 | 中国石油天然气集团有限公司 | 一种页岩气水平井压裂垂向改造体积控制方法 |
CN114542038A (zh) * | 2020-11-26 | 2022-05-27 | 中国石油天然气股份有限公司 | 一种水平井天然裂缝发育段的压裂方法 |
CN112883661B (zh) * | 2021-01-30 | 2023-07-07 | 太原理工大学 | 一种碎软低渗油气储层的压裂模拟方法 |
CN113177248B (zh) * | 2021-04-21 | 2022-07-19 | 武汉大学 | 隧道围岩破裂碎胀大变形失稳灾变过程数值模拟方法 |
CN113431562A (zh) * | 2021-07-19 | 2021-09-24 | 西南石油大学 | 一种裂缝性致密砂岩储层缝网酸压施工排量优化方法 |
CN113919147A (zh) * | 2021-09-29 | 2022-01-11 | 长江大学 | 一种油藏压裂裂缝网络扩展路径计算方法 |
CN114252343B (zh) * | 2021-12-17 | 2023-06-09 | 中国石油大学(华东) | 一种深层断溶体裂缝闭合机理实验装置 |
CN114239308B (zh) * | 2021-12-24 | 2022-08-12 | 西南石油大学 | 一种多尺度高密度压裂参数优化方法 |
WO2023151072A1 (zh) * | 2022-02-14 | 2023-08-17 | 烟台杰瑞石油装备技术有限公司 | 砂堵预测方法、装置、设备及存储介质 |
CN114638137B (zh) * | 2022-03-31 | 2023-03-28 | 王永亮 | 一种基于热-流-固-损伤耦合的干热岩产热预测方法 |
CN114547952B (zh) * | 2022-04-25 | 2022-07-08 | 西安交通大学 | 一种基于深度学习-界面流耦合的压水堆中子输运方法 |
CN114547953B (zh) * | 2022-04-26 | 2022-07-19 | 西南石油大学 | 一种基于优化设计图版的压裂施工参数优化方法及系统 |
CN114542043B (zh) * | 2022-04-28 | 2022-08-12 | 太原理工大学 | 基于压裂液粘度优化改进岩层压裂增渗的方法及装置 |
CN114692472B (zh) * | 2022-06-01 | 2022-08-12 | 中国石油大学(华东) | 储层改造辅助降压开发天然气水合物藏的数值模拟方法 |
CN114818551B (zh) * | 2022-07-01 | 2022-11-04 | 西安石油大学 | 一种页岩油藏压裂缝网自适应半无限导流表征方法及装置 |
CN114880957B (zh) * | 2022-07-11 | 2022-09-16 | 西南石油大学 | 一种非常规储层压裂射孔段簇参数组合优化方法 |
CN115288650B (zh) * | 2022-08-04 | 2023-06-27 | 西安石油大学 | 孔隙弹性介质中并行计算与模拟水力压裂的方法 |
CN115455770B (zh) * | 2022-09-14 | 2023-04-25 | 东北石油大学 | 一种基于子结构的油气田井筒结构应力精细模拟方法 |
CN115935687B (zh) * | 2022-12-27 | 2023-09-08 | 哈尔滨工程大学 | 一种计算翻边轴承耦合润滑与动力学特性参数的方法 |
CN116227287B (zh) * | 2023-02-27 | 2023-11-14 | 西南石油大学 | 一种基于线性互补方法的裂缝流体流动流固耦合模拟方法 |
CN117217106A (zh) * | 2023-07-31 | 2023-12-12 | 中国石油大学(北京) | 一种缝内暂堵后地质模型节点位移的确定方法及装置 |
CN116738622B (zh) * | 2023-08-09 | 2023-10-20 | 成都理工大学 | 一种能够降低能耗的射孔器设计方法 |
CN117216903A (zh) * | 2023-09-13 | 2023-12-12 | 西南石油大学 | 一种耦合储层流动的射孔参数优化方法 |
CN117495147B (zh) * | 2023-12-22 | 2024-03-22 | 中国石油大学(华东) | 一种考虑段间压裂工艺差异的裂缝网络扩展智能预测方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103122762A (zh) * | 2013-01-30 | 2013-05-29 | 中国石油天然气股份有限公司 | 一种非常规泥页岩油气藏有效压裂层段的检测方法及装置 |
CN104237104A (zh) * | 2014-09-26 | 2014-12-24 | 中国石油天然气股份有限公司 | 确定体积压裂改造程度的方法 |
CN104570128A (zh) * | 2013-10-16 | 2015-04-29 | 中国石油化工股份有限公司 | 一种基于岩石脆性因子的非常规油气储层预测方法 |
CN105201484A (zh) * | 2015-10-29 | 2015-12-30 | 西南石油大学 | 一种直井分层压裂层段优选及施工参数优化设计方法 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105156103B (zh) * | 2015-09-29 | 2018-06-05 | 西南石油大学 | 一种岩屑-岩心-井眼-储层多尺度的页岩储层三维可压裂性评价方法 |
CN105201479B (zh) * | 2015-10-09 | 2017-10-24 | 西南石油大学 | 一种页岩储层水平井分段压裂射孔簇参数优化设计方法 |
CN205172548U (zh) * | 2015-12-04 | 2016-04-20 | 西安石油大学 | 一种体积压裂纵横向裂缝缝网延伸扩展模拟装置 |
CN107545113B (zh) * | 2017-09-08 | 2020-02-07 | 西南石油大学 | 非常规油气藏水力压裂复杂缝网形成过程模拟方法 |
CN107705215B (zh) * | 2017-09-25 | 2018-07-27 | 西南石油大学 | 一种页岩储层重复压裂选井选段方法 |
CN108442911A (zh) * | 2018-02-28 | 2018-08-24 | 西南石油大学 | 一种页岩气水平井重复压裂水力裂缝参数优化设计方法 |
-
2019
- 2019-03-22 CN CN201910223536.2A patent/CN109992864B/zh active Active
- 2019-07-19 US US16/516,261 patent/US11567233B2/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103122762A (zh) * | 2013-01-30 | 2013-05-29 | 中国石油天然气股份有限公司 | 一种非常规泥页岩油气藏有效压裂层段的检测方法及装置 |
CN104570128A (zh) * | 2013-10-16 | 2015-04-29 | 中国石油化工股份有限公司 | 一种基于岩石脆性因子的非常规油气储层预测方法 |
CN104237104A (zh) * | 2014-09-26 | 2014-12-24 | 中国石油天然气股份有限公司 | 确定体积压裂改造程度的方法 |
CN105201484A (zh) * | 2015-10-29 | 2015-12-30 | 西南石油大学 | 一种直井分层压裂层段优选及施工参数优化设计方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109992864A (zh) | 2019-07-09 |
US20200301043A1 (en) | 2020-09-24 |
US11567233B2 (en) | 2023-01-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109992864B (zh) | 非常规双重介质储层体积压裂数值模拟及参数优化方法 | |
CN111322050B (zh) | 一种页岩水平井段内密切割暂堵压裂施工优化方法 | |
Lian et al. | A study on casing deformation failure during multi-stage hydraulic fracturing for the stimulated reservoir volume of horizontal shale wells | |
CN109933939B (zh) | 非常规双重介质储层多裂缝起裂与扩展的数值模拟方法 | |
Li et al. | A fully coupled model for hydraulic-fracture growth during multiwell-fracturing treatments: enhancing fracture complexity | |
CN105201484A (zh) | 一种直井分层压裂层段优选及施工参数优化设计方法 | |
CN104695928A (zh) | 一种裂缝性致密油藏直井体积改造产能评价方法 | |
Furui et al. | A Comprehensive Model of High-Rate Matrix-Acid Stimulation for Long Horizontal Wells in Carbonate Reservoirs: Part II—Wellbore/Reservoir Coupled-Flow Modeling and Field Application | |
Lin et al. | Numerical simulation of the influence of stimulated reservoir volume on in-situ stress field | |
CN111472730B (zh) | 大段多簇压裂的射孔方案确定方法 | |
CN114737943A (zh) | 非常规储层立体井网改造方法、装置、设备及介质 | |
Xu et al. | Mechanism study of casing deformation in multistage hydraulic fracturing shale reservoir | |
Hu et al. | Three‐dimensional complex fracture propagation simulation: Implications for rapid decline of production capacity | |
Yang et al. | Coupled hydro-mechanical analysis of gas production in fractured shale reservoir by random fracture network modeling | |
Yuan et al. | Experimental investigation on fracture geometry in multi-stage fracturing under tri-axial stresses | |
Fu et al. | 3D DEM simulation on the interference of multiple hydraulic fractures in horizontal wells | |
He et al. | Phase Field Modeling of Multiple Fracture Growth in Natural Fractured Reservoirs | |
Liu et al. | Optimization of multi-cluster fracturing model under the action of induced stress in horizontal wells | |
de Pater et al. | Stimulation for geothermal wells in the Netherlands | |
CN114607334A (zh) | 陆相页岩气储层压裂方法 | |
Su et al. | Analysis of casing stress during multistage fracturing of shale gas horizontal wells considering thermo‐hydro‐mechanical coupling | |
Yang et al. | Visualization and characterization of experimental hydraulic fractures interacting with karst fracture-cavity distributions | |
Luo et al. | Optimization of the number of openhole packers in horizontal wells completed by inflow control devices for different types of reservoirs | |
Chen et al. | Study on Activation Behavior of Complex Fracture Network in Weak-Planar Rock | |
Haifeng et al. | Steady-state optimization model and practical design of the fracture network system in tight sand gas reservoirs |
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 | ||
TA01 | Transfer of patent application right | ||
TA01 | Transfer of patent application right |
Effective date of registration: 20191213 Address after: Three road 610000 Sichuan city of Chengdu Province, No. 1 East Applicant after: Chengdu University of Technology Applicant after: Southwest Petroleum University Address before: Three road 610059 Sichuan city of Chengdu province Chenghua District Erxian Qiaodong No. 1 Applicant before: Chengdu University of Technology |
|
GR01 | Patent grant | ||
GR01 | Patent grant |