CN110929413B - 一种模拟裂缝性储层中水力裂缝壁面不连续性行为的方法 - Google Patents
一种模拟裂缝性储层中水力裂缝壁面不连续性行为的方法 Download PDFInfo
- Publication number
- CN110929413B CN110929413B CN201911225898.1A CN201911225898A CN110929413B CN 110929413 B CN110929413 B CN 110929413B CN 201911225898 A CN201911225898 A CN 201911225898A CN 110929413 B CN110929413 B CN 110929413B
- Authority
- CN
- China
- Prior art keywords
- fracture
- natural
- hydraulic
- equation
- natural fracture
- 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 21
- 230000035699 permeability Effects 0.000 claims abstract description 33
- 238000004364 calculation method Methods 0.000 claims abstract description 22
- 238000006073 displacement reaction Methods 0.000 claims abstract description 18
- 238000009826 distribution Methods 0.000 claims abstract description 11
- 230000008878 coupling Effects 0.000 claims abstract description 7
- 238000010168 coupling process Methods 0.000 claims abstract description 7
- 238000005859 coupling reaction Methods 0.000 claims abstract description 7
- 230000008569 process Effects 0.000 claims abstract description 6
- 239000011148 porous material Substances 0.000 claims description 16
- 239000012530 fluid Substances 0.000 claims description 15
- 239000011159 matrix material Substances 0.000 claims description 10
- 239000000463 material Substances 0.000 claims description 8
- 238000001914 filtration Methods 0.000 claims description 4
- 238000013459 approach Methods 0.000 claims description 3
- 230000017105 transposition Effects 0.000 claims description 3
- 238000005065 mining Methods 0.000 claims 1
- 230000005012 migration Effects 0.000 abstract description 4
- 238000013508 migration Methods 0.000 abstract description 4
- 230000006399 behavior Effects 0.000 description 9
- 230000000694 effects Effects 0.000 description 5
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 4
- 238000011160 research Methods 0.000 description 4
- 230000015572 biosynthetic process Effects 0.000 description 3
- 238000010276 construction Methods 0.000 description 3
- 239000011435 rock Substances 0.000 description 3
- 238000004088 simulation Methods 0.000 description 3
- 238000007796 conventional method Methods 0.000 description 2
- 239000007789 gas Substances 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 239000003345 natural gas Substances 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 239000003795 chemical substances by application Substances 0.000 description 1
- 230000001808 coupling effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 229920001971 elastomer Polymers 0.000 description 1
- 239000000806 elastomer Substances 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000008439 repair process Effects 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Images
Classifications
-
- 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
- E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
- E21B49/003—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells by analysing drilling variables or conditions
-
- 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
-
- 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
- E21B47/00—Survey of boreholes or wells
- E21B47/003—Determining well or borehole volumes
-
- 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
- E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
- E21B49/005—Testing the nature of borehole walls or the formation by using drilling mud or cutting data
-
- G01V20/00—
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- 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
- E21B2200/00—Special features related to earth drilling for obtaining oil, gas or water
- E21B2200/20—Computer models or simulations, e.g. for reservoirs under production, drill bits
-
- 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
技术领域
本发明涉及水力裂缝技术领域,特别涉及一种模拟裂缝性储层中水力裂缝壁面不连续性行为的方法。
背景技术
针对于非常规储层的水力压裂技术研究,目前大多数研究者把原本为孔弹性的地层考虑为弹性地层,把流体的流动仅考虑为裂缝内的流动而忽略了孔隙中的流动。这样的考虑是建立在弹性力学的基础上,但是却与非常规储层的实际情况相差甚远。所以得到的水力裂缝壁面在大多数的研究中都是光滑的、连续的面。但是,通过很多实验可知(王跃.基于大型水力压裂实验系统的室内模拟研究[J].地球物理学进展,2017,32(1):408-413;张健,张国祥,李良等.页岩水力压裂模拟实验研究[J].合肥工业大学学报(自然科学版),2019,42(4):541-545),在对储层岩石实施水力压裂措施后形成的裂缝壁面并不是光滑、连续分布的,而是呈现出凹凸起伏、非连续的状态。这种现象在大多数研究中被忽略,但是水力裂缝壁面的不连续性对于支撑剂在其中的运移有很大的影响。所以常规的方法已经不能够准确地模拟出水力裂缝在压裂过程中的真实状态。
再者,目前大多数的研究关注的是水力裂缝在非常规储层中的延伸轨迹以及水力裂缝与天然裂缝的相交行为(李晓,赫建明,尹超等.页岩结构面特征及其对水力压裂的控制作用[J].天然气地球科学,2019,40(3):653-660;Chen X,Li Y,Zhao J,Xu W,FuD.Numerical investigation for simultaneous growth ofhydraulic fractures inmultiple horizontal wells[J].Journal ofNatural Gas Science and Engineering,2018,51:44-52)。其目的是研究水力裂缝与天然裂缝间的相交作用机制,然后通过模拟预测压裂过程中储层的成网能力以及形成的改造体积。但是受计算方法的限制,当天然裂缝条数成倍数增加时,计算时间将呈指数递增,对于计算机的性能以及配置要求也将大幅度提升(Roussel,N.P.and Sharma,M.M.Role of Stress Reorientation in the Successof Refracture Treatments in Tight Gas Sands[J].SPE Prod&Oper,2012,27(4):346–355)。所以原有的模拟方法就不再适用于现场工程实际。除此之外,由于储层被考虑为弹性体,而非孔弹性介质,在计算流体滤失以及储层应力时,得到的结果会存在较大的误差。这也是常规方法的弊端(Gao Q,Ghassemi A.Pore pressure and stress distributionsaround a hydraulic fracture in heterogeneous rock[J].Rock Mechanics and RockEngineering,2017,50(12),3157-3173)。
发明内容
针对上述问题,本发明旨在提供一种模拟裂缝性储层中水力裂缝壁面不连续性行为的方法,通过考虑天然裂缝、孔弹性效应、水力驱动作用以及压裂时流体与地质应力间的耦合作用,基于Biot理论、流固耦合原理、有限差分理论以及多相渗流原理等建立一个全新的可以模拟多簇水力裂缝延伸以及分析水力裂缝壁面不连续性的数学模型,利用该模型结合施工参数、初始条件以及边界条件预测在不同参数下压裂动态过程中水力裂缝壁面受非均匀应力影响产生的偏移以及不均匀扩展,从而通过分析裂缝壁面的不连续性来分析壁面的粗糙程度。最终可以为支撑剂在水力裂缝内的运移提供指导与参考。
本发明的技术方案如下:
一种模拟裂缝性储层中水力裂缝壁面不连续性行为的方法,包括以下步骤:
获取裂缝性储层的地质参数,结合天然裂缝长、宽、高的物理条件,基于裂缝连续体模型建立天然裂缝的物理模型;
基于水力裂缝内流动方程、滤失方程、宽度方程以及物质平衡方程建立水力裂缝扩展计算方程;
通过摩尔库伦定律建立天然裂缝破坏模型,根据所述天然裂缝破坏模型得到天然裂缝开度,再根据所述天然裂缝开度计算得到天然裂缝渗透率,最后通过所述裂缝连续体模型将所述天然裂缝渗透率转化为多孔介质的渗透率;
将所述水力裂缝扩展计算方程通过裂缝扩展准则以及流体滤失与所述多孔介质的渗透率进行耦合,得到考虑天然裂缝影响的耦合天然裂缝的孔弹性模型;
通过所述耦合天然裂缝的孔弹性模型,即可得到水力裂缝壁面应力、位移分布情况,通过位移大小,分析水力裂缝壁面的偏移以及不连续性行为。
进一步地,所述地质参数通过测井或压裂获得,具体包括天然裂缝初始开度、基质初始渗透率、初始孔隙度、弹性模量、泊松比。
进一步地,所述裂缝连续体模型为:
式中:
kij为基质渗透率张量,m2;
knf为天然裂缝渗透率,m2;
n1、n2、n3的计算方法为:
式中:
ξ为逼近角,°;
所述天然裂缝渗透率knf的计算方法为:
式中:
wNF为天然裂缝开度,m;
d为天然裂缝间距,m。
进一步地,通过将所述水力裂缝内流动方程、滤失方程、宽度方程代入所述物质平衡方程即可获得所述水力裂缝扩展计算方程,所述水力裂缝内流动方程、滤失方程、宽度方程代入以及物质平衡方程具体依次如下所示:
wHF=(u+-u-)·n (6)
式中:
pHF为缝内压力,Pa;
s为裂缝长度方向坐标,m;
μ为流体粘度,mPa·s;
hf为水力裂缝高度,m;
wHF为水力裂缝开度,m;
q为水力裂缝内流量,m3/s;
qL为压裂液滤失速度,m/s;
S为压裂液滤失通过的面积,m2;
f为滤失系数,1/m;
pp为储层孔隙压力,MPa;
u+、u-分别为水力裂缝左右两面上的位移,m;
n为水力裂缝面上的单位法向量,无量纲;
t为压裂的时间,min。
进一步地,所述天然裂缝破坏模型为:
pNF≥σn+Kt (9)
wNF=a0+aNFT+aNFS (10)
式中:
στ为天然裂缝壁面所受剪切应力,MPa;
τ0为天然裂缝内聚力,MPa;
σn为天然裂缝面上的正应力,MPa;
pNF为天然裂缝中流体压力,MPa;
Kt为天然裂缝抗张强度,MPa;
a0为天然裂缝的初始开度,m;
aNFT为天然裂缝的张性开度,m;
aNFS为天然裂缝的剪切开度,m。
进一步地,所述基本摩擦角在30°~40°范围内。
进一步地,所述耦合天然裂缝的孔弹性模型为:
▽·[μ▽u+μ▽uT+λItr(▽u)]+▽·σi-b▽pp+b▽pi=0 (11)
式中:
▽为拉普拉斯算子;
u为位移张量,m;
T为矩阵的转置;
λ为拉梅常数,MPa;
Itr为积分符号;
σi为应力张量,MPa;
b为Biot有效系数,无量纲;
pi为初始孔隙压力,MPa;
M为Biot模量,MPa;
Kdr为排采过程中的体积模量,MPa。
与现有技术相比,本发明具有如下优点:
针对非常规储层实施分段多簇压裂措施过程,本发明考虑了天然裂缝对储层流体流动、水力裂缝扩展以及应力分布的影响,还耦合考虑了储层的孔弹性效应,由此可以准确地计算得到不同裂缝簇数、不同压裂时间水力裂缝壁面位移的偏移与不连续性行为。且当天然裂缝条数成倍数增加,本发明通过所述裂缝连续体模型直接计算天然裂缝的基质渗透率张量,解决了现有技术中,需要单独计算每条天然裂缝,计算量大的缺点。综上所述,本发明对分析水力压裂时支撑剂的运移以及压后的裂缝导流能力具有重要作用,另外对于油田气田开发以及丰富水力压裂基础理论也具有借鉴意义。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本发明天然裂缝的物理模型示意图;
图2为实施例1压裂结束后的应力分布云图;
图3为实施例1压裂结束后的位移分布云图;
图4为实施例1水力裂缝壁面的不连续分布示意图。
具体实施方式
下面结合附图和实施例对本发明进一步说明。需要说明的是,在不冲突的情况下,本申请中的实施例及实施例中的技术特征可以相互结合。
一种模拟裂缝性储层中水力裂缝壁面不连续性行为的方法,包括以下步骤:
首先,获取裂缝性储层的地质参数,结合天然裂缝长、宽、高的物理条件,基于裂缝连续体模型建立天然裂缝的物理模型,所述地质参数通过测井或压裂获得,具体包括天然裂缝初始开度、基质初始渗透率、初始孔隙度、弹性模量、泊松比。所述裂缝连续体模型为:
式中:
kij为基质渗透率张量,m2;
knf为天然裂缝渗透率,m2;
n1、n2、n3的计算方法为:
式中:
ξ为逼近角,°;
所述天然裂缝渗透率knf的计算方法为:
式中:
wNF为天然裂缝开度,m;
d为天然裂缝间距,m。
其次,基于水力裂缝内流动方程、滤失方程、宽度方程以及物质平衡方程建立水力裂缝扩展计算方程,所述水力裂缝扩展计算方程通过将所述水力裂缝内流动方程、滤失方程、宽度方程代入所述物质平衡方程即可获得。所述水力裂缝内流动方程、滤失方程、宽度方程代入以及物质平衡方程具体依次如下所示:
wHF=(u+-u-)·n (6)
式中:
pHF为缝内压力,Pa;
s为裂缝长度方向坐标,m;
μ为流体粘度,mPa·s;
hf为水力裂缝高度,m;
wHF为水力裂缝开度,m;
q为水力裂缝内流量,m3/s;
qL为压裂液滤失速度,m/s;
S为压裂液滤失通过的面积,m2;
f为滤失系数,1/m;
pp为储层孔隙压力,MPa;
u+、u-分别为水力裂缝左右两面上的位移,m;
n为水力裂缝面上的单位法向量,无量纲;
t为压裂的时间,min。
再次,通过摩尔库伦定律建立天然裂缝破坏模型,根据所述天然裂缝破坏模型得到天然裂缝开度,再根据所述天然裂缝开度计算得到天然裂缝渗透率,最后通过所述裂缝连续体模型将所述天然裂缝渗透率转化为多孔介质的渗透率。所述天然裂缝破坏模型为:
pNF≥σn+Kt (9)
wNF=a0+aNFT+aNFS (10)
式中:
στ为天然裂缝壁面所受剪切应力,MPa;
τ0为天然裂缝内聚力,MPa;
σn为天然裂缝面上的正应力,MPa;
pNF为天然裂缝中流体压力,MPa;
Kt为天然裂缝抗张强度,MPa;
a0为天然裂缝的初始开度,m;
aNFT为天然裂缝的张性开度,m;
aNFS为天然裂缝的剪切开度,m。
可选地,所述基本摩擦角在30°~40°范围内。
然后,将所述水力裂缝扩展计算方程通过裂缝扩展准则以及流体滤失与所述多孔介质的渗透率进行耦合,得到考虑天然裂缝影响的耦合天然裂缝的孔弹性模型,所述耦合天然裂缝的孔弹性模型为:
▽·[μ▽u+μ▽uT+λItr(▽u)]+▽·σi-b▽pp+b▽pi=0 (11)
式中:
▽为拉普拉斯算子;
u为位移张量,m;
T为矩阵的转置;
λ为拉梅常数,MPa;
Itr为积分符号;
σi为应力张量,MPa;
b为Biot有效系数,无量纲;
pi为初始孔隙压力,MPa;
M为Biot模量,MPa;
Kdr为排采过程中的体积模量,MPa。
最后,通过所述耦合天然裂缝的孔弹性模型,即可得到水力裂缝壁面应力、位移分布情况,通过位移大小,分析水力裂缝壁面的偏移以及不连续性行为。
实施例1
以涪陵某区块典型的非常规储层作为待模拟裂缝性储层。通过现场测井、试井,获取该区块的地质参数和施工参数如表1所示:
表1待模拟裂缝性储层地质参数和施工参数
根据上述地质参数结合天然裂缝长、宽、高的物理条件,建立如图1所示的天然裂缝物理模型。
根据式(8)-式(10)的天然裂缝破坏模型得到天然裂缝开度,再通过式(3)计算得到天然裂缝渗透率,最后通过式(1)和式(2)将所述天然裂缝渗透率转化为多孔介质的渗透率。
将式(4)-式(6)代入式(7)得到水力裂缝扩展计算方程,将所述水力裂缝扩展计算方程通过裂缝扩展准则以及流体滤失与所述多孔介质的渗透率进行耦合,得到式(11)、式(12)的考虑天然裂缝影响的耦合天然裂缝的孔弹性模型;
通过所述耦合天然裂缝的孔弹性模型,即可得到如图2所示的水力裂缝壁面应力、如图3所示的位移分布情况,通过位移大小,即可分析水力裂缝壁面的偏移以及不连续性行为,在水力压裂25min、50min、75min时,水力裂缝壁面的偏移以及不连续性分布如图4所示。
以上所述,仅是本发明的较佳实施例而已,并非对本发明作任何形式上的限制,虽然本发明已以较佳实施例揭露如上,然而并非用以限定本发明,任何熟悉本专业的技术人员,在不脱离本发明技术方案范围内,当可利用上述揭示的技术内容作出些许更动或修饰为等同变化的等效实施例,但凡是未脱离本发明技术方案的内容,依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化与修饰,均仍属于本发明技术方案的范围内。
Claims (3)
1.一种模拟裂缝性储层中水力裂缝壁面不连续性行为的方法,其特征在于,包括以下步骤:
获取裂缝性储层的地质参数,结合天然裂缝长、宽、高的物理条件,基于裂缝连续体模型建立天然裂缝的物理模型,所述裂缝连续体模型为:
式中:kij为基质渗透率张量,m2;knf为天然裂缝渗透率,m2;
n1、n2、n3的计算方法为:
所述天然裂缝渗透率knf的计算方法为:
式中:wNF为天然裂缝开度,m;d为天然裂缝间距,m;
基于水力裂缝内流动方程、滤失方程、宽度方程以及物质平衡方程建立水力裂缝扩展计算方程,具体的通过将所述水力裂缝内流动方程、滤失方程、宽度方程代入所述物质平衡方程即可获得所述水力裂缝扩展计算方程,所述水力裂缝内流动方程、滤失方程、宽度方程代入以及物质平衡方程具体依次如下所示:
wHF=(u+-u-)·n (6)
式中:为偏微分符号;pHF为缝内压力,Pa;s为裂缝长度方向坐标,m;μ为流体粘度,mPa·s;hf为水力裂缝高度,m;wHF为水力裂缝开度,m;q为水力裂缝内流量,m3/s;qL为压裂液滤失速度,m/s;S为压裂液滤失通过的面积,m2;f为滤失系数,1/m;pp为储层孔隙压力,MPa;u+、u-分别为水力裂缝左右两面上的位移,m;n为水力裂缝面上的单位法向量,无量纲;t为压裂的时间,min;
通过摩尔库伦定律建立天然裂缝破坏模型,所述天然裂缝破坏模型为:
pNF≥σn+Kt (9)
wNF=a0+aNFT+aNFS (10)
式中:στ为天然裂缝壁面所受剪切应力,MPa;τ0为天然裂缝内聚力,MPa;为基本摩擦角,°;σn为天然裂缝面上的正应力,MPa;pNF为天然裂缝中流体压力,MPa;Kt为天然裂缝抗张强度,MPa;a0为天然裂缝的初始开度,m;aNFT为天然裂缝的张性开度,m;aNFS为天然裂缝的剪切开度,m;
根据所述天然裂缝破坏模型得到天然裂缝开度,再根据所述天然裂缝开度计算得到天然裂缝渗透率,最后通过所述裂缝连续体模型将所述天然裂缝渗透率转化为多孔介质的渗透率;
将所述水力裂缝扩展计算方程通过裂缝扩展准则以及流体滤失与所述多孔介质的渗透率进行耦合,得到考虑天然裂缝影响的耦合天然裂缝的孔弹性模型,所述耦合天然裂缝的孔弹性模型为:
式中:为拉普拉斯算子;u为位移张量,m;T为矩阵的转置;λ为拉梅常数,MPa;Itr为积分符号;σi为应力张量,MPa;b为Biot有效系数,无量纲;pi为初始孔隙压力,MPa;M为Biot模量,MPa;Kdr为排采过程中的体积模量,MPa;
通过所述耦合天然裂缝的孔弹性模型,即可得到水力裂缝壁面应力、位移分布情况,通过位移大小,分析水力裂缝壁面的偏移以及不连续性行为。
2.根据权利要求1所述的模拟裂缝性储层中水力裂缝壁面不连续性行为的方法,其特征在于,所述地质参数通过测井或压裂获得,具体包括天然裂缝初始开度、基质初始渗透率、初始孔隙度、弹性模量、泊松比。
3.根据权利要求1所述的模拟裂缝性储层中水力裂缝壁面不连续性行为的方法,其特征在于,所述基本摩擦角在30°~40°范围内。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911225898.1A CN110929413B (zh) | 2019-12-04 | 2019-12-04 | 一种模拟裂缝性储层中水力裂缝壁面不连续性行为的方法 |
US16/929,108 US11118450B2 (en) | 2019-12-04 | 2020-07-15 | Method for simulating the discontinuity of the hydraulic fracture wall in fractured reservoirs |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911225898.1A CN110929413B (zh) | 2019-12-04 | 2019-12-04 | 一种模拟裂缝性储层中水力裂缝壁面不连续性行为的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110929413A CN110929413A (zh) | 2020-03-27 |
CN110929413B true CN110929413B (zh) | 2020-11-10 |
Family
ID=69857760
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911225898.1A Active CN110929413B (zh) | 2019-12-04 | 2019-12-04 | 一种模拟裂缝性储层中水力裂缝壁面不连续性行为的方法 |
Country Status (2)
Country | Link |
---|---|
US (1) | US11118450B2 (zh) |
CN (1) | CN110929413B (zh) |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112241593B (zh) * | 2020-10-21 | 2021-11-30 | 西南石油大学 | 一种基于多重时间步的裂缝性储层滤失计算方法 |
CN112360420B (zh) * | 2020-11-20 | 2022-05-10 | 中国石油天然气股份有限公司 | 一种实现水平井分段多簇压裂多裂缝均衡扩展系统及方法 |
CN114592840B (zh) * | 2020-12-04 | 2023-10-27 | 中国石油天然气股份有限公司 | 暂堵压裂方法及其应用 |
CN112730198B (zh) * | 2021-01-15 | 2022-04-01 | 中国地质大学(武汉) | 一种页岩天然裂缝渗透率计算方法 |
CN113009571B (zh) * | 2021-02-18 | 2022-03-08 | 中国矿业大学(北京) | 双相介质中的水平裂缝的反射系数和透射系数确定方法 |
CN112761609B (zh) * | 2021-02-19 | 2022-02-01 | 西南石油大学 | 一种水力压裂作业中高效铺置支撑剂的优化方法 |
CN113591420B (zh) * | 2021-09-08 | 2022-05-13 | 中国石油大学(北京) | 真三轴水力压裂实验的仿真方法及处理器 |
CN114441337B (zh) * | 2022-01-14 | 2022-10-18 | 中国石油大学(北京) | 一种金属断裂韧性裂纹尖端张开位移的检测方法 |
CN114547998B (zh) * | 2022-02-28 | 2024-03-22 | 西南石油大学 | 一种通过耦合储层流动确定水平井压裂改造体积的方法 |
CN115288650B (zh) * | 2022-08-04 | 2023-06-27 | 西安石油大学 | 孔隙弹性介质中并行计算与模拟水力压裂的方法 |
CN115898375A (zh) * | 2022-12-20 | 2023-04-04 | 西南石油大学 | 一种模拟裂缝流固耦合变形的颗粒运移可视化实验装置及方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106869892A (zh) * | 2017-03-15 | 2017-06-20 | 西南石油大学 | 一种重复压裂井缝内暂堵起裂的判断方法 |
CN108319756A (zh) * | 2017-12-29 | 2018-07-24 | 西安石油大学 | 一种致密储层体积压裂缝网扩展模拟及表征方法 |
CN108468538A (zh) * | 2018-03-23 | 2018-08-31 | 西南石油大学 | 一种页岩水压裂缝扩展预测方法 |
CN108920824A (zh) * | 2018-07-02 | 2018-11-30 | 中国地质大学(北京) | 一种窄河道储层的产量递减分析方法 |
Family Cites Families (33)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8756016B2 (en) * | 2009-01-29 | 2014-06-17 | Schlumberger Technology Corporation | Method and system to estimate fracture aperture in horizontal wells |
US20100256964A1 (en) * | 2009-04-07 | 2010-10-07 | Schlumberger Technology Corporation | System and technique to quantify a fracture system |
US8494827B2 (en) * | 2009-09-25 | 2013-07-23 | Exxonmobil Upstream Research Company | Method of predicting natural fractures and damage in a subsurface region |
US8898044B2 (en) * | 2009-11-25 | 2014-11-25 | Halliburton Energy Services, Inc. | Simulating subterranean fracture propagation |
US8392165B2 (en) * | 2009-11-25 | 2013-03-05 | Halliburton Energy Services, Inc. | Probabilistic earth model for subterranean fracture simulation |
US8386226B2 (en) * | 2009-11-25 | 2013-02-26 | Halliburton Energy Services, Inc. | Probabilistic simulation of subterranean fracture propagation |
CA2820942C (en) * | 2010-12-10 | 2016-12-20 | Baker Hughes Incorporated | A method to improve reservoir simulation and recovery from fractured reservoirs |
WO2012083463A1 (en) * | 2010-12-22 | 2012-06-28 | Dusseault Maurice B | Multi-stage fracture injection process for enhanced resource production from shales |
WO2012097405A1 (en) * | 2011-01-20 | 2012-07-26 | Commonwealth Scientific And Industrial Research Organisation | Hydraulic fracturing |
DK177547B1 (da) * | 2011-03-04 | 2013-10-07 | Maersk Olie & Gas | Fremgangsmåde og system til brønd- og reservoir-management i udbygninger med åben zone såvel som fremgangsmåde og system til produktion af råolie |
US8762118B2 (en) * | 2011-03-07 | 2014-06-24 | Schlumberger Technology Corporation | Modeling hydraulic fractures |
US9618652B2 (en) * | 2011-11-04 | 2017-04-11 | Schlumberger Technology Corporation | Method of calibrating fracture geometry to microseismic events |
WO2015003028A1 (en) * | 2011-03-11 | 2015-01-08 | Schlumberger Canada Limited | Method of calibrating fracture geometry to microseismic events |
US9715026B2 (en) * | 2011-03-11 | 2017-07-25 | Schlumberger Technology Corporation | System and method for performing microseismic fracture operations |
MX2014000772A (es) * | 2011-07-28 | 2014-05-01 | Schlumberger Technology Bv | Sistema y metodo para llevar a cabo operaciones de fractura del hoyo. |
WO2013055930A1 (en) * | 2011-10-11 | 2013-04-18 | Schlumberger Canada Limited | System and method for performing stimulation operations |
US10422208B2 (en) * | 2011-11-04 | 2019-09-24 | Schlumberger Technology Corporation | Stacked height growth fracture modeling |
US20160265331A1 (en) * | 2011-11-04 | 2016-09-15 | Schlumberger Technology Corporation | Modeling of interaction of hydraulic fractures in complex fracture networks |
MX340326B (es) * | 2012-03-06 | 2016-07-05 | Ion Geophysical Corp | Modelo que predice la facturacion del esquisto. |
RU2591857C1 (ru) * | 2012-08-24 | 2016-07-20 | Шлюмбергер Текнолоджи Б.В. | Система и способ выполнения работ по стимуляции недр |
US9416642B2 (en) * | 2013-02-01 | 2016-08-16 | Halliburton Energy Services, Inc. | Modeling subterranean rock blocks in an injection treatment simulation |
EP3108088B1 (en) * | 2014-02-19 | 2018-08-22 | Repsol, S.A. | Method implemented in a computer for the numerical simulation of a porous medium |
US10788604B2 (en) * | 2014-06-25 | 2020-09-29 | Schlumberger Technology Corporation | Fracturing and reactivated fracture volumes |
US10519769B2 (en) * | 2014-09-10 | 2019-12-31 | Fracture ID, Inc. | Apparatus and method using measurements taken while drilling to generate and map mechanical boundaries and mechanical rock properties along a borehole |
WO2016122792A1 (en) * | 2015-01-28 | 2016-08-04 | Schlumberger Canada Limited | Method of performing wellsite fracture operations with statistical uncertainties |
WO2017027342A1 (en) * | 2015-08-07 | 2017-02-16 | Schlumberger Technology Corporation | Method of performing complex fracture operations at a wellsite having ledged fractures |
AU2015413610A1 (en) * | 2015-11-02 | 2018-04-19 | Halliburton Energy Services, Inc. | Three-dimensional geomechanical modeling of casing deformation for hydraulic fracturing treatment design |
US10571605B2 (en) * | 2016-03-28 | 2020-02-25 | Exxonmobil Upstream Research Company | Method and system for modeling and simulating a fractured reservoir |
US10572611B2 (en) * | 2016-04-29 | 2020-02-25 | Exxonmobil Upstream Research Company | Method and system for characterizing fractures in a subsurface region |
WO2017216594A1 (en) * | 2016-06-15 | 2017-12-21 | Schlumberger Technology Corporation | Induced seismicity |
WO2018069744A1 (en) * | 2016-10-14 | 2018-04-19 | Schlumberger Technology Corporation | Geologic structural model generation |
CN106869896B (zh) * | 2017-04-26 | 2018-05-08 | 西南石油大学 | 一种低渗透油气藏直井体积压裂储层改造体积预测方法 |
US20200095858A1 (en) * | 2017-09-14 | 2020-03-26 | Saudi Arabian Oil Company | Modeling reservoir permeability through estimating natural fracture distribution and properties |
-
2019
- 2019-12-04 CN CN201911225898.1A patent/CN110929413B/zh active Active
-
2020
- 2020-07-15 US US16/929,108 patent/US11118450B2/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106869892A (zh) * | 2017-03-15 | 2017-06-20 | 西南石油大学 | 一种重复压裂井缝内暂堵起裂的判断方法 |
CN108319756A (zh) * | 2017-12-29 | 2018-07-24 | 西安石油大学 | 一种致密储层体积压裂缝网扩展模拟及表征方法 |
CN108468538A (zh) * | 2018-03-23 | 2018-08-31 | 西南石油大学 | 一种页岩水压裂缝扩展预测方法 |
CN108920824A (zh) * | 2018-07-02 | 2018-11-30 | 中国地质大学(北京) | 一种窄河道储层的产量递减分析方法 |
Non-Patent Citations (2)
Title |
---|
Simulation of complex fracture networks influenced by natural fractures in shale gas reservoir;ZhaoJinzhou;《Natural Gas Industry B》;20141113;第1卷(第1期);第89-95页 * |
多场耦合作用对水力裂缝形态的影响;胡永全;《大庆石油地质与开发》;20180427;第37卷(第3期);第82-89页 * |
Also Published As
Publication number | Publication date |
---|---|
US11118450B2 (en) | 2021-09-14 |
US20200347721A1 (en) | 2020-11-05 |
CN110929413A (zh) | 2020-03-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110929413B (zh) | 一种模拟裂缝性储层中水力裂缝壁面不连续性行为的方法 | |
CN107545113B (zh) | 非常规油气藏水力压裂复杂缝网形成过程模拟方法 | |
Zhao et al. | Latest research progresses in network fracturing theories and technologies for shale gas reservoirs | |
CN110750930B (zh) | 一种基于裂缝连续体模型预测裂缝性储层应力演化的方法 | |
Lian et al. | A study on casing deformation failure during multi-stage hydraulic fracturing for the stimulated reservoir volume of horizontal shale wells | |
CN108868748B (zh) | 一种页岩气水平井重复压裂裂缝开启压力的计算方法 | |
CN109284571B (zh) | 一种二氧化碳置换页岩气多尺度多场耦合渗流数学建模方法 | |
Zhou et al. | Hydraulic fracture propagation direction during volume fracturing in unconventional reservoirs | |
CN104806233B (zh) | 一种预测弱面地层坍塌压力当量密度窗口的方法 | |
Zhang et al. | Employing a quad-porosity numerical model to analyze the productivity of shale gas reservoir | |
CN109209333B (zh) | 页岩气多井群高效开采间距优化方法 | |
CN112727534B (zh) | 基于真三向应力和渗透率动态变化的瓦斯钻孔布置方法 | |
CN106897544B (zh) | 对来自压裂的非常规地层的油气生产进行建模的方法 | |
RU2666573C1 (ru) | Способ разработки нефтяной залежи с проведением повторного гидроразрыва пласта с изменением направления трещины | |
US8788250B2 (en) | Method of improved reservoir simulation of fingering systems | |
Yu et al. | On how asymmetric stimulated rock volume in shales may impact casing integrity | |
Lin et al. | Numerical simulation of the influence of stimulated reservoir volume on in-situ stress field | |
CN105114065A (zh) | 水力压裂裂缝延伸过程的模拟方法 | |
Yu et al. | Impact of asymmetric stimulated rock volume on casing deformation in multi-stage fracturing; a case study | |
CN112012710A (zh) | 一种砂泥岩交互地层水平井压裂三维裂缝扩展模拟方法 | |
CN110781617A (zh) | 基于数模的缝内暂堵压裂裂缝复杂性识别方法 | |
JPWO2001063092A1 (ja) | 石油産出方法 | |
CN114580100B (zh) | 压裂水平井全井筒压力计算方法、设备和计算机可读储存介质 | |
Shi et al. | Perforation cluster spacing optimization with hydraulic fracturing-reservoir simulation modeling in shale gas reservoir | |
CN115510778A (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |