CN105201479A - 一种页岩储层水平井分段压裂射孔簇参数优化设计方法 - Google Patents
一种页岩储层水平井分段压裂射孔簇参数优化设计方法 Download PDFInfo
- Publication number
- CN105201479A CN105201479A CN201510648540.5A CN201510648540A CN105201479A CN 105201479 A CN105201479 A CN 105201479A CN 201510648540 A CN201510648540 A CN 201510648540A CN 105201479 A CN105201479 A CN 105201479A
- Authority
- CN
- China
- Prior art keywords
- sigma
- stress
- fracture
- crack
- perforation
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 37
- 238000013461 design Methods 0.000 title abstract description 11
- 230000008859 change Effects 0.000 claims abstract description 13
- 230000007246 mechanism Effects 0.000 claims abstract description 10
- 238000013178 mathematical model Methods 0.000 claims abstract description 6
- 206010017076 Fracture Diseases 0.000 claims description 76
- 208000010392 Bone Fractures Diseases 0.000 claims description 73
- 238000000926 separation method Methods 0.000 claims description 48
- 239000011435 rock Substances 0.000 claims description 27
- 230000000977 initiatory effect Effects 0.000 claims description 23
- 229910052500 inorganic mineral Inorganic materials 0.000 claims description 20
- 239000011707 mineral Substances 0.000 claims description 20
- 230000006378 damage Effects 0.000 claims description 16
- 230000008569 process Effects 0.000 claims description 15
- 238000005336 cracking Methods 0.000 claims description 11
- 238000010168 coupling process Methods 0.000 claims description 9
- 230000008878 coupling Effects 0.000 claims description 8
- 238000005859 coupling reaction Methods 0.000 claims description 8
- 238000004088 simulation Methods 0.000 claims description 7
- 208000013201 Stress fracture Diseases 0.000 claims description 6
- 238000004458 analytical method Methods 0.000 claims description 6
- 230000015572 biosynthetic process Effects 0.000 claims description 4
- 238000006073 displacement reaction Methods 0.000 claims description 3
- 238000013210 evaluation model Methods 0.000 claims description 3
- 230000000694 effects Effects 0.000 abstract description 8
- 238000010276 construction Methods 0.000 abstract description 3
- 230000001131 transforming effect Effects 0.000 abstract 2
- 230000006698 induction Effects 0.000 abstract 1
- 230000002250 progressing effect Effects 0.000 abstract 1
- 239000012530 fluid Substances 0.000 description 24
- 239000011159 matrix material Substances 0.000 description 12
- 238000005516 engineering process Methods 0.000 description 8
- 230000009466 transformation Effects 0.000 description 7
- 238000005457 optimization Methods 0.000 description 6
- 230000001965 increasing effect Effects 0.000 description 5
- 238000004519 manufacturing process Methods 0.000 description 5
- 230000035699 permeability Effects 0.000 description 5
- 238000011161 development Methods 0.000 description 4
- 230000018109 developmental process Effects 0.000 description 4
- 239000000463 material Substances 0.000 description 4
- 229910021532 Calcite Inorganic materials 0.000 description 3
- 238000009826 distribution Methods 0.000 description 3
- 229910000514 dolomite Inorganic materials 0.000 description 3
- 239000010459 dolomite Substances 0.000 description 3
- 210000002569 neuron Anatomy 0.000 description 3
- 239000010453 quartz Substances 0.000 description 3
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N silicon dioxide Inorganic materials O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 description 3
- 210000004027 cell Anatomy 0.000 description 2
- 230000003993 interaction Effects 0.000 description 2
- 230000000704 physical effect Effects 0.000 description 2
- 239000011148 porous material Substances 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000007704 transition Effects 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
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 1
- 206010010149 Complicated fracture Diseases 0.000 description 1
- 101100001675 Emericella variicolor andJ gene Proteins 0.000 description 1
- 208000027418 Wounds and injury Diseases 0.000 description 1
- 230000001133 acceleration Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 239000002734 clay mineral Substances 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000006866 deterioration Effects 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 230000005489 elastic deformation Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 238000001730 gamma-ray spectroscopy Methods 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000008595 infiltration Effects 0.000 description 1
- 238000001764 infiltration Methods 0.000 description 1
- 208000014674 injury Diseases 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 230000035515 penetration Effects 0.000 description 1
- 239000003209 petroleum derivative Substances 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 239000007790 solid phase Substances 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
Abstract
本发明涉及一种页岩储层水平井分段压裂射孔簇参数优化设计方法,它包括以下步骤:S1、选出高可压性压裂井段;S2、建立水力压裂裂缝扩展的数学模型,分析裂缝延伸的诱导应力场模型,并分析裂缝转向机理;S3、分析主裂缝的延伸情况,选出能使各主裂缝均匀向前扩展的簇间距;S4、分析主裂缝周围应力场的变化情况,选出水平主应力比小于1.3的射孔簇间距;并结合步骤S3确定最优的射孔簇间距;S5、利用交替压裂原理,进一步增加裂缝的复杂性,将中间射孔簇孔密设为10-16孔/m,两侧射孔簇孔密设为16孔/m以上。本发明可以优化射孔簇间距,获得最大的裂缝改造体积,增加储层改造的实际效果,避免因盲目设置射孔簇间距引起的施工成本偏高和压裂效果差的情况。
Description
技术领域
本发明涉及石油天然气开发技术领域。尤其涉及一种页岩储层水平井分段压裂射孔簇参数优化设计方法。
背景技术
页岩储层在我国储量增长、能源安全方面正发挥越来越重要的作用。页岩储层具有渗透率低、强度高、储层分布不均匀等特点,水平井分段压裂增产改造技术是开发该类储层、增加单井控制储量和产能的一项重要技术。常规水平井分段压裂采用单段射孔、单段压裂方式,从减小缝间干扰方面出发,考虑尽可能增大射孔间距,形成多条垂直于井筒的横切缝。
目前在页岩储层开发中己经认识到由于基质向裂缝的供油气能力较差,仅靠单一的压裂主缝很难取得预期的增产效果,往往需要特殊的压裂技术形成复杂裂缝,增加储层基质向水力裂缝供油气能力,提高压裂增产效果。而近些年来开始推广应用的“缝网压裂”和“体积改造”技术就是利用裂缝间的干扰,形成复杂的应力场,使裂缝在延伸过程中形成分叉缝或者发生转向,沟通天然裂缝,增加裂缝复杂程度,以期获得更好的增产效果。
针对“缝网压裂”和“体积改造”技术,国外首先开展了裂缝诱导应力场的研究。Sneddon和Elliott较早的在理论上研究了裂缝周围应力场,推导出了在无限大弹性体中裂缝周围的应力场计算公式。而对于扁平型裂缝(penny-shapedcrack),Sneddon建立了新的方程来描述它们对周围应力场的影响。L.D.Palmer分析了煤层气藏中初次压裂裂缝产生的诱导应力,同时研究了裂缝形态和净压力对储层渗透率的影响。N.P.Roussel等研究了水力压裂裂缝诱导应力和孔隙压力诱导应力对水平井重复压裂新缝起裂和延伸的影响。同时,国内学者也对缝网压裂做过相关研究。ChanganM.Du等人提出采用水力压裂技术在双孔隙页岩气藏中形成人工缝网的方法,利用微地震技术记录“改造体积”并校正。雷群等人利用储层水平主应力差与裂缝延伸净压力的关系,在远井地带也实现了复杂缝网,改善了油气流动通道。但是地层因素、多裂缝及裂缝间距对水平井压裂结果有着重要影响,而国内外现有研究均未考虑以上因素。
发明内容
本发明的目的在于克服现有技术的不足,提供一种可以优化射孔簇间距,获得最大裂缝改造体积,增加储层改造实际效果的页岩储层水平井分段压裂射孔簇参数优化设计方法。
本发明的目的是通过以下技术方案实现的:一种页岩储层水平井分段压裂射孔簇参数优化设计方法,它包括以下步骤:
S1、使用储层岩石力学参数以及矿物组成参数来计算储层可压性指数,获取井段可压性指数剖面,选出高可压性压裂井段;
S2、建立水力压裂裂缝扩展的数学模型,分析裂缝延伸的诱导应力场模型,并分析裂缝转向机理;
S3、模拟不同射孔簇间距下,主裂缝的延伸情况,选出最优裂缝尺寸下的簇间距;
S4、分析不同射孔簇间距下,主裂缝周围应力场的变化情况,分析不同射孔簇间距下压裂形成“缝网”的可行性,选出能够形成“缝网”的簇间距;并结合步骤S3中形成最优裂缝尺寸的射孔簇间距范围,确定最优的射孔簇间距;
S5、利用交替压裂的应力干扰,形成更加复杂的裂缝,将中间射孔簇孔密设为10-16孔/m、两侧射孔簇孔密为16孔/m以上。
进一步地,所述步骤S1包括以下子步骤:
S11、计算储层岩石的弹性参数脆性指数EEn和矿物含量脆性指数BBn,求弹性参数脆性指数EEn和矿物含量脆性指数BBn的平均值,得到储层岩石的综合脆性指数BTot;
S12、根据储层岩石的可压裂性评价模型:
Frac=ωBTot/((1-ω)sinφKICKIIC)(5)
结合测井数据,绘制储层岩石沿井眼轴线的可压裂性指数剖面,选出高可压裂性区域进行压裂改造,高可压性指数区域是根据实际情况选择相对较高区域。其中Frac是可压性指数,无量纲;BTot为综合脆性指数,无量纲;是内摩擦角,单位为度;ω为数值范围0~1的权重系数,无量纲;KIC是I型断裂韧性,单位为MPa·m1/2;KIIC是II型断裂韧性,单位为MPa·m1/2。
进一步地,所述步骤S2包括以下子步骤:
S21、建立水力压裂裂缝扩展的数学模型,根据所述步骤S1选出的储层的地应力、岩石力学特征、物性参数等数据和地质分层特征,采用有限元、离散元、边界元或位移不连续等方法,考虑多条裂缝同时扩展的情形,建立页岩储层水平井分段压裂裂缝起裂与动态扩展的渗流-应力-断裂损伤耦合数学模型;
S22、分析裂缝延伸的诱导应力场模型,原始地应力由最大水平主应力、最小水平主应力和垂向应力组成,后续起裂裂缝周围的应力场由先起裂裂缝产生的诱导应力场与原地应力场叠加组成,根据叠加原理,分段压裂产生的第n条裂缝周围的复合应力场模型为:
式中,σ'H(n)、σ'h(n)、σ'v(n)为第n条裂缝周围的复合应力分量,单位为MPa;σH、σh、συ分别为储层水平最大应力、水平最小应力和垂向应力,单位为MPa;ν为泊松比,无量纲;σ1x(in)、σ1y(in)和σ1z(in)分别为第i条裂缝对第n条裂缝产生的诱导应力X轴向、Y轴向和Z轴向分量,单位均为MPa;
S23、分析转向机理,分段压裂产生的第n条裂缝发生转向的条件为:该条裂缝受到的最小水平主应力方向的诱导应力与最大水平主应力方向的诱导应力之差大于或等于原始最大、最小水平主应力之差,公式表示为:
进一步地,所述步骤S3的具体过程是:
根据所述步骤S2中所建立的多裂缝动态扩展数学模型,模拟不同射孔簇间距下,主裂缝在缝长和缝高方向的延伸情况,并比较每簇裂缝最终形成的几何尺寸,以此判断诱导应力对主裂缝延伸的影响程度,选出能使各主裂缝均匀向前扩展的最优簇间距。
进一步地,所述步骤S4的具体过程是:
根据所述步骤S2分析不同射孔簇间距下,主裂缝周围应力场的变化情况,判断相邻两主裂缝之间及裂缝附近的水平主应力比σ'H(n)/σ'h(n),选出水平主应力比小于1.3的射孔簇间距,从而得到从应力干扰角度形成“缝网”所需的射孔簇间距;结合步骤S3中形成最优裂缝尺寸的射孔簇间距范围,确定最优的射孔簇间距;
进一步地,所述步骤S5的具体过程是:
模拟同时压裂和交替压裂时,中间裂缝的应力干扰规律,揭示中间裂缝的应力转向机理,利用交替压裂的应力干扰,以形成更加复杂的裂缝,将中间射孔簇孔密设为10-16孔/m、两侧射孔簇孔密为16孔/m以上。
本发明具有以下优点:
1、选择高可压性压裂井段进行优化设计,确保最大化的压后缝网体积和压后页岩气产量。页岩储层地质构造复杂,储层内部物性和非均质性差异大,储层可压裂性对压裂形成裂缝网络体积有着重要影响,如何筛选出可压裂性高的储层,是体积压裂改造形成更多复杂裂缝网络体积的前提。对于1千多米甚至接近2千米的长水平段,若压裂位置选择在储层可压裂性较差的井段,即使采用大规模的体积压裂和最优化的射孔簇间距,也不能得到预期的裂缝网络体积,致使压后产量较低,前期大规模的压裂成本投入难以收回。因此,要优选出水平井分段压裂射孔簇间距,就必须首先确定长水平段内的高可压裂性井段,在该井段内,再进行射孔簇间距优化。
2、根据页岩储层水平井分段压裂的施工工艺特征,考虑多簇射孔形成的水力裂缝的动态弹塑性损伤扩展,此优化设计方法可以模拟多簇裂缝同时扩展时相互影响的裂缝三维几何形态;
3、此优化设计方法中采用的渗流-应力-损伤耦合的有限元方法,可以实现多簇裂缝扩展过程中的裂缝诱导应力场的相互干扰;
4、此优化设计方法进一步利用了交替压裂时,中间射孔簇因射孔孔密小而裂缝后起裂的原理,来增加中间裂缝的应力干扰和转向程度,从而形成更加复杂的裂缝形态,对储层进行更加彻底的改造;
5、运用此优化设计方法,结合多簇裂缝动态扩展相互干扰形成的裂缝形态以及应力场,可优选合理的射孔簇间距,从而获得最大的裂缝改造体积,增加储层改造的实际效果,避免因盲目设置射孔簇间距引起的施工成本偏高和压裂效果较差的情况;
6、应用范围广:本发明不仅适用于页岩储层水平井分段压裂射孔簇间距的优化设计,同样适用于致密砂岩储层的多段体积压裂。
附图说明
图1是本发明一种页岩储层水平井分段压裂射孔簇参数优化设计方法的流程图;
图2是本发明页岩储层水平井段的可压裂性指数剖面;
图3是本发明中有限元模型;
图4是本发明中模型使用的基本参数;
图5是本发明中模型示意图;
图6是本发明中簇间距10m情况下,压裂30分钟后的裂缝立体示意图;
图7是本发明中簇间距20m情况下,压裂30分钟后的裂缝立体示意图;
图8是本发明中簇间距30m情况下,压裂30分钟后的裂缝立体示意图;
图9是本发明中簇间距40m情况下,压裂30分钟后的裂缝立体示意图;
图10是本发明中簇间距50m情况下,压裂30分钟后的裂缝立体示意图;
图11是本发明中簇间距10m情况下,压裂30分钟后的Z=0平面内不同高度处的水平应力差异系数;
图12是本发明中簇间距20m情况下,压裂30分钟后的Z=0平面内不同高度处的水平应力差异系数;
图13是本发明中簇间距30m情况下,压裂30分钟后的Z=0平面内不同高度处的水平应力差异系数;
图14是本发明中簇间距40m情况下,压裂30分钟后的Z=0平面内不同高度处的水平应力差异系数;
图15是本发明中簇间距50m情况下,压裂30分钟后的Z=0平面内不同高度处的水平应力差异系数。
具体实施方式
下面结合附图对本发明做进一步的描述,但本发明的保护范围不局限于以下所述。
如图1所示,一种页岩储层水平井分段压裂射孔簇参数优化设计方法,它包括以下步骤:
S1、使用储层岩石力学参数以及矿物组成参数来计算储层可压性指数,获取井段可压性指数剖面,选出高可压性压裂井段。页岩储层地质构造复杂,储层内部物性和非均质性差异大,储层可压裂性对压裂形成裂缝网络体积有着重要影响,如何筛选出可压裂性高的储层,是体积压裂改造形成更多复杂裂缝网络体积的前提。对于1千多米甚至接近2千米的长水平段,若压裂位置选择在储层可压裂性较差的井段,即使采用大规模的体积压裂和最优化的射孔簇间距,也不能得到预期的裂缝网络体积,致使压后产量较低,前期大规模的压裂成本投入难以收回。因此,要优选出水平井分段压裂射孔簇间距,就必须首先确定长水平段内的高可压裂性井段,在该井段内,再进行射孔簇间距优化,以确保最大化的压后缝网体积和压后页岩气产量。
S2、建立水力压裂裂缝扩展的数学模型,分析裂缝延伸的诱导应力场模型,并分析裂缝转向机理;
S3、模拟不同射孔簇间距下,主裂缝的延伸情况,选出最优裂缝尺寸下的簇间距;
S4、分析不同射孔簇间距下,主裂缝周围应力场的变化情况,分析不同射孔簇间距下压裂形成“缝网”的可行性,选出能够形成“缝网”的簇间距;并结合步骤S3中形成最优裂缝尺寸的射孔簇间距范围,确定最优的射孔簇间距。
S5、利用交替压裂的应力干扰,形成更加复杂的裂缝,将中间射孔簇孔密设为10-16孔/m、两侧射孔簇孔密为16孔/m以上。
所述步骤S1包括以下子步骤:
S11、计算储层岩石的弹性参数脆性指数EEn和矿物含量脆性指数BBn,求弹性参数脆性指数EEn和矿物含量脆性指数BBn的平均值,得到储层岩石的综合脆性指数BTot;
本实例中,页岩水平井完井井深4049m,垂深2965m,水平段长907.5m,该区块的岩石力学脆性指数采用Rickman的矿物含量脆性指数评价经验公式(RickmanR,MullenM,PetreE,etal.APracticalUseofShalePetrophysicsforStimulationDesignOptimization:AllShalePlaysAreNotClonesoftheBarnettShale.SPE115258,SPEAnnualTechnicalConferenceandExhibition,21-24September,Denver,Colorado,USA,2008)计算,Rickman的矿物含量脆性指数评价经验公式为:
式中,EEn为弹性参数脆性指数,无量纲;
En为弹性模量脆性指数,无量纲;E、Emax和Emin分别是页岩储层的弹性模量、页岩储层内最大弹性模量和最小弹性模量,单位为GPa;
νn为泊松比脆性指数,无量纲;ν、vmax和vmin分别是页岩储层的泊松比、页岩储层的最大泊松比和最小泊松比,无量纲。
通过声波测井资料获得产能段杨氏模量E和泊松比ν,该井所处区块的储层杨氏模量分布范围为10~60GPa;泊松比分布范围为0.1~0.50,利用公式(1)计算出岩石力学脆性指数。
采用Wang和Gale的矿物含量脆性指数BBn评价经验公式(Wang,F.P.,andJ.F.W.Gale.Screeningcriteriaforshale-gassystems:GulfCoastAssociationofGeologicalSocietiesTransactions,v.59,p.779-793,2009)计算矿物含量脆性指数BBn,Wang和Gale的矿物含量脆性指数BBn评价经验公式为:
BBn=(W石英、方解石+W白云岩)/W总质量(2)
式中,BBn为矿物含量脆性指数,无量纲;W石英、方解石为石英和方解石的质量,单位为Kg;W白云岩为白云石质量,单位为Kg;W总质量为矿物的总质量,单位为Kg。
总脆性指数BTot(无量纲),可以表示为:
BTot=(EEn+BBn)/2(3)
采用常规测井、元素俘获能谱测井(ECS)自然伽马能谱测井(NGS)和岩性密度测井等资料,计算目标产能段粘土矿物含量,石英类矿物含量,碳酸盐岩类矿物含量,利用公式(2)计算矿物含量脆性指数BBn,并根据公式(3)求取综合脆性指数BTot;
S12、采用巴西圆盘实验测定页岩试样的抗拉强度、I型和II型断裂韧性KIC和KIIC,建立内摩擦角I型和II型断裂韧性与岩石力学特征参数之间的关系式:
式中,σt为页岩的抗拉强度,单位为MPa;σn为裂缝面法向围压,单位为MPa;ρ为页岩密度,单位为Kg/m3;为岩石内摩擦角,单位为度;KIC是I型断裂韧性,单位为MPa·m1/2,KIIC是II型断裂韧性,单位为MPa·m1/2。
储层岩石的可压裂性评价模型为:
Frac=ωBTot/((1-ω)sinφKICKIIC)(5)
式中,Frac是可压性指数,无量纲;BTot为综合脆性指数,无量纲;是内摩擦角,单位为度;ω为数值范围0~1的权重系数,无量纲;
根据储层岩石力学参数以及矿物组成参数,并利用公式(4)和(5)计算出该水平井水平段的可压性指数剖面,如图2所示,从图2可压性指数剖面可以看出,可在该区域上选取4个高可压性指数区域。
所述步骤S2包括以下子步骤:
S21、建立水力压裂裂缝扩展的数学模型,根据所述步骤S1选出的储层的地应力、岩石力学特征、物性参数等数据和地质分层特征,采用有限元、离散元、边界元或位移不连续等方法,考虑多条裂缝同时扩展的情形,建立页岩储层水平井分段压裂裂缝起裂与动态扩展的渗流-应力-断裂损伤耦合数学模型。
水力压裂过程中,泵压随排量增加不断增大,作用于裂缝面上流体渗流压力也不断增大,使得流体向地层的滤失增加,导致岩石孔隙中的应力状态发生改变。而岩石中应力的变化必然引起储层孔隙度、流体渗流速度等参数的改变,反过来又会影响到裂缝面上渗流场孔隙压力的变化,储层岩石中这种流体渗流与岩石变形的相互制约,相互作用关系即称为渗流-应力耦合。以均质,各向同性的二维平面应变模型为基础,建立渗流-应力-断裂损伤流固耦合模型,水力裂缝的起裂与扩展采用粘弹性损伤cohesive单元描述。本实施例中,建立页岩储层水平井分段压裂裂缝起裂与动态扩展的渗流-应力-断裂损伤耦合数学模型过程如下(其中公式(6)~公式(9)为渗流-应力耦合模型的系列方程,公式(10)~公式(14)为裂缝动态扩展损伤模型的系列方程,公式(6)~公式(14)综合起来即为裂缝起裂与动态扩展的渗流-应力-断裂损伤耦合数学模型):
渗流和应力耦合通过有效应力与总应力的关系进行耦合,有效应力与总应力的关系为:
式中,为有效应力向量矩阵;σ为总应力矩阵,应力单位为MPa;pw是压力的绝对值,单位为MPa;I为二阶单位张量,无量纲;
取控制体微元为V,单位为m3,其表面为S,单位为㎡,岩石骨架的应力平衡方程:
式中,σ和分别为应力矩阵和虚应变率矩阵,t、和δv分别为表面力向量、体积力向量和虚速度向量;
将应力平衡方程式(7)离散化得到固相材料的有限元网格,同时使流体可以流经这些网格,流体渗流需要满足连续性方程:
式中,So为初始饱和度,无量纲;
mT为单位矩阵,mT=[1,1,1,0,0,0]T;
Dep为弹塑性矩阵;
KS为岩石骨架颗粒的压缩模量,单位GPa;
ko为初始渗透率张量与流体密度的乘积矩阵;
kr为比例渗透率,单位为m/s2;
ρo为流体密度,单位为Kg/m3;
g为重力加速度,单位为m/s2;
ξ为饱和度变化率,无量纲;
φ为孔隙度,无量纲;
Ko为流体体积模量,单位为GPa;
p为流体压力,单位为MPa;
运用粘弹性损伤cohesive单元模拟裂缝动态扩展,cohesive单元能够模拟裂缝中流体流动和流体向基质中渗流。渗流系数与体积应变的关系为:
式中,k为流体渗透系数,单位m/s;
k0为初始渗透系数,单位m/s;
n0为初始孔隙度,无量纲;
εV为体积应变,无量纲;
计算过程中,岩石的渗透性演化方程通过有限元软件的子程序二次开发功能予以实现。
本实施例中,采用二次应力起裂准则作为水力裂缝是否起裂的判断依据,即当cohesive孔压单元的三向应力与其对应的强度极限的比值的平方和达到1时,Cohesive孔压单元起裂,公式如下:
式中,tn、ts、tt分别为cohesive单元沿着法向(垂直于cohesive单元上下表面的方向)和两个切向的应力分量(二维情况下tt不存在),单位为MPa;为cohesive单元的抗拉强度,单位为MPa;和(二维情况下不存在)为两个切向方向的抗剪强度,单位为MPa;
Cohesive单元采用刚度退化描述单元损伤演化过程,其损伤演化模型为:
式中,和分别为cohesive孔压单元三个方向按照未损伤阶段线弹性变形时计算得到的应力,单位为MPa;tn、ts、tt为三个对应方向实际承受的压力,单位为MPa;D为无量纲损伤因子,取值0-1,D=0时材料未损伤,D=1时材料完全损伤。当裂缝中积累的能量达到岩石的临界破裂能量时单元起裂。可以采用下式计算裂缝起裂临界能量:
式中,GIC为裂缝起裂临界能量,单位为J;KIC为断裂韧性,单位为MPa·m1/2;ν为泊松比,无量纲;E为杨氏模量,单位为GPa。
对复合型裂缝起裂后的扩展,应用B-K准则,即由Benzeggagh和Kenane提出的裂缝扩展临界能量释放率准则,即:
式中,Gc为复合型裂缝临界断裂能量释放率,单位为N/mm;
为法向断裂临界应变能释放率,单位为N/mm;
为切向断裂临界能量释放率,单位为N/mm;
Gn、Gs、Gt分别为法向、第一切向和第二切向断裂能释放率,单位为N/mm;
B-K准则认为
η为与材料本身特性有关的常数,无量纲;
当裂缝尖端节点处计算的能量释放率大于B-K临界能量释放率是,Cohesive单元当前裂尖节点对绑定部分将解开,裂缝向前扩展。
裂缝内流体流动模型如下:
本实施例中模拟压裂液为牛顿流体,裂缝中的一部分流体通过cohesive单元的两个表面渗透进入地层。流体通过cohesive单元两个表面的法向流动计算公式为:
式中,qt和qb分别为流体流入cohesive单元上、下表面的流速,单位m/s;ct和cb分别为两个表面的滤失系数,无量纲;pt、pb、pi分别为cohesive单元两个表面的流体压力以及cohesive单元中面的流体压力,单位为MPa。
S22、分析裂缝延伸的诱导应力场模型,原始地应力由最大水平主应力、最小水平主应力和垂向应力组成,后续起裂裂缝周围的应力场由先起裂裂缝产生的诱导应力场与原地应力场叠加组成,根据叠加原理,分段压裂产生的第n条裂缝周围的复合应力场模型为:
式中,σ'H(n)、σ'h(n)、σ'v(n)为第n条裂缝周围的复合应力分量,单位为MPa;σH、σh、συ分别为储层水平最大应力、水平最小应力和垂向应力,单位为MPa;ν为泊松比,无量纲;σ1x(in)、σ1y(in)和σ1z(in)分别为第i条裂缝对第n条裂缝产生的诱导应力X轴向、Y轴向和Z轴向分量,单位均为MPa;
S23、分析转向机理
根据弹性力学和岩石破裂准则,水力裂缝总是产生于强度最弱、阻力最小的方向,即裂缝破裂面垂直于最小主应力方向。因此,分段压裂产生的初始裂缝是垂直于井筒方向的横向裂缝,而初始裂缝产生的诱导应力会对后续起裂裂缝周围应力场产生影响,在原地应力和诱导应力作用下,后续起裂裂缝在原来最小水平主应力方向上受到的应力可能会大于原来最大水平主应力方向的应力,即最大、最小水平主应力方向发生了变化,导致裂缝在延伸过程中发生转向,向平行于井筒方向延伸。分段压裂产生的第n条裂缝发生转向的条件为:该条裂缝受到的最小水平主应力方向的诱导应力与最大水平主应力方向的诱导应力之差大于或等于原始最大、最小水平主应力之差,公式表示为:
所述步骤S3的具体过程是:
根据储层的分层特征,建立三维模型如图3所示。模型使用如图4所示参数,压裂液流体粘度:30mPa.s;排量:12m3/min。cohesive单元参数如表1。
表1cohesive单元参数
根据所述步骤S2中所建立的多裂缝动态扩展数学模型(式(6)~式(14)),建立裂缝的有限元模型,如图5所示,模型纵向方向高度100m,横向方向长度100m,设置不同的簇间距进行压裂模拟,以优化簇间距。分析不同射孔簇间距下裂缝的最终的空间几何形态,如图6~图10所示,从图中可以看出,射孔簇间距在10m和20m的时候,射孔簇间裂缝应力干扰较强,中间簇裂缝的延伸受到很大程度的限制。而当簇间距为30m-50m时,中间簇裂缝的延伸几乎未受到限制,所有簇裂缝都正常均匀延伸。故在该储层情况下射孔簇间距应控制在30m或30m以上。
所述步骤S4的具体过程是:
根据所述步骤S2分析不同射孔簇间距下,主裂缝周围应力场的变化情况。判断相邻两主裂缝之间及裂缝附近的水平主应力比σ'H(n)/σ'h(n),水平主应力比小于1.3时,水力压裂能够形成充分的裂缝网络;水平主应力比为1.3~1.5时,水力压裂在高的净压力时能够形成较为充分的裂缝网络,此时裂缝延伸净压力需要大于水平主应力差;水平主应力比大于1.5时,水力压裂不能形成裂缝网络。根据以上可知,应选择水平主应力比小于1.3的射孔簇间距,此时从应力干扰角度能形成“缝网”。
从图11~图15可以看出:10m射孔簇间距时,裂缝干扰非常严重,出现应力反转现象;30m和40m射孔簇间距时,裂缝干扰同样较为严重,最大主应力出现应力反转现象,此时水平主应力比的大部分区域都小于1.3,对于主裂缝周围产生局部改造体积是非常有利的;相对于40m射孔簇间距,30m射孔簇间距的应力场干扰更强;50m射孔簇间距时,裂缝干扰情况较30m和40m射孔簇间距时的小,应力干扰情况不明显。因此此时簇间距应小于40m。
故综合步骤S3和步骤S4可知:从分段多簇、多簇同时压裂时主裂缝的几何形状和主裂缝周围形成复杂“缝网”的可行性来看,优选30~40m的射孔簇间距。
所述步骤S5的具体过程是:
模拟同时压裂和交替压裂时,中间裂缝的应力干扰规律,揭示中间裂缝的应力转向机理,利用交替压裂的应力干扰,以形成更加复杂的裂缝,进而将中间射孔簇孔密设为10-16孔/m、两侧射孔簇孔密为16孔/m以上。
本实施例中所用软件为ABAQUS软件,但实际应用中亦可使用其他模拟软件。
Claims (6)
1.一种页岩储层水平井分段压裂射孔簇参数优化设计方法,其特征在于,它包括以下步骤:
S1、使用储层岩石力学参数以及矿物组成参数来计算储层可压性指数,获取井段可压性指数剖面,选出高可压性压裂井段;
S2、建立水力压裂裂缝扩展的数学模型,分析裂缝延伸的诱导应力场模型,并分析裂缝转向机理;
S3、模拟不同射孔簇间距下,主裂缝的延伸情况,选出最优裂缝尺寸下的簇间距;
S4、分析不同射孔簇间距下,主裂缝周围应力场的变化情况,分析不同射孔簇间距下压裂形成“缝网”的可行性,选出能够形成“缝网”的簇间距;并结合步骤S3中形成最优裂缝尺寸的射孔簇间距范围,确定最优的射孔簇间距;
S5、利用交替压裂的应力干扰,形成更加复杂的裂缝,将中间射孔簇孔密设为10-16孔/m、两侧射孔簇孔密为16孔/m以上。
2.根据权利要求1所述的一种页岩储层水平井分段压裂射孔簇参数优化设计方法,其特征在于,所述步骤S1包括以下子步骤:
S11、计算储层岩石的弹性参数脆性指数EEn和矿物含量脆性指数BBn,求弹性参数脆性指数EEn和矿物含量脆性指数BBn的平均值,得到储层岩石的综合脆性指数BTot;
S12、根据储层岩石的可压裂性评价模型:
Frac=ωBTot/((1-ω)sinφKICKIIC)
结合测井数据,绘制储层岩石沿井眼轴线的可压裂性指数剖面,选出高可压裂性区域进行压裂改造,其中Frac是可压性指数,无量纲;BTot为综合脆性指数,无量纲;是内摩擦角,单位为度;ω为数值范围0~1的权重系数,无量纲;KIC是I型断裂韧性,单位为MPa·m1/2;KIIC是II型断裂韧性,单位为MPa·m1/2。
3.根据权利要求1所述的一种页岩储层水平井分段压裂射孔簇参数优化设计方法,其特征在于,所述步骤S2包括以下子步骤:
S21、建立水力压裂裂缝扩展的数学模型,根据所述步骤S1选出的储层的地应力、岩石力学特征、物性参数数据和地质分层特征,采用有限元、离散元、边界元或位移不连续方法,考虑多条裂缝同时扩展的情形,建立页岩储层水平井分段压裂裂缝起裂与动态扩展的渗流-应力-断裂损伤耦合数学模型;
S22、分析裂缝延伸的诱导应力场模型,原始地应力由最大水平主应力、最小水平主应力和垂向应力组成,后续起裂裂缝周围的应力场由先起裂裂缝产生的诱导应力场与原地应力场叠加组成,根据叠加原理,分段压裂产生的第n条裂缝周围的复合应力场模型为:
式中,σ'H(n)、σ'h(n)、σ'v(n)为第n条裂缝周围的复合应力分量,单位为MPa;σH、σh、συ分别为储层水平最大应力、水平最小应力和垂向应力,单位为MPa;ν为泊松比,无量纲;σ1x(in)、σ1y(in)和σ1z(in)分别为第i条裂缝对第n条裂缝产生的诱导应力X轴向、Y轴向和Z轴向分量,单位均为MPa;
S23、分析转向机理,分段压裂产生的第n条裂缝发生转向的条件为:该条裂缝受到的最小水平主应力方向的诱导应力与最大水平主应力方向的诱导应力之差大于或等于原始最大、最小水平主应力之差,公式表示为:
4.根据权利要求1所述的一种页岩储层水平井分段压裂射孔簇参数优化设计方法,其特征在于,所述步骤S3的具体过程是:
根据所述步骤S2中所建立的多裂缝动态扩展数学模型,模拟不同射孔簇间距下,主裂缝在缝长和缝高方向的延伸情况,并比较每簇裂缝最终形成的几何尺寸,以此判断诱导应力对主裂缝延伸的影响程度,选出能使各主裂缝均匀向前扩展的最优簇间距。
5.根据权利要求1所述的一种页岩储层水平井分段压裂射孔簇参数优化设计方法,其特征在于,所述步骤S4的具体过程是:
根据所述步骤S2分析不同射孔簇间距下,主裂缝周围应力场的变化情况,判断相邻两主裂缝之间及裂缝附近的水平主应力比σ'H(n)/σ'h(n),选出水平主应力比小于1.3的射孔簇间距,从而得到从应力干扰角度形成“缝网”所需的射孔簇间距;结合步骤S3中形成最优裂缝尺寸的射孔簇间距范围,确定最优的射孔簇间距。
6.根据权利要求1所述的一种页岩储层水平井分段压裂射孔簇参数优化设计方法,其特征在于,所述步骤S5的具体过程是:
模拟同时压裂和交替压裂时,中间裂缝的应力干扰规律,揭示中间裂缝的应力转向机理,利用交替压裂的应力干扰,以形成更加复杂的裂缝,将中间射孔簇孔密设为10-16孔/m、两侧射孔簇孔密为16孔/m以上。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510648540.5A CN105201479B (zh) | 2015-10-09 | 2015-10-09 | 一种页岩储层水平井分段压裂射孔簇参数优化设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510648540.5A CN105201479B (zh) | 2015-10-09 | 2015-10-09 | 一种页岩储层水平井分段压裂射孔簇参数优化设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105201479A true CN105201479A (zh) | 2015-12-30 |
CN105201479B CN105201479B (zh) | 2017-10-24 |
Family
ID=54949392
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510648540.5A Active CN105201479B (zh) | 2015-10-09 | 2015-10-09 | 一种页岩储层水平井分段压裂射孔簇参数优化设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105201479B (zh) |
Cited By (58)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105672973A (zh) * | 2016-01-26 | 2016-06-15 | 西南石油大学 | 一种煤层气藏整体体积压裂优化设计方法 |
CN105735960A (zh) * | 2016-03-22 | 2016-07-06 | 西南石油大学 | 一种低渗透油气藏水平井分段多簇压裂簇间距优化方法 |
CN105735961A (zh) * | 2016-02-18 | 2016-07-06 | 西南石油大学 | 一种低渗致密储层分层压裂级数优选方法 |
CN106097125A (zh) * | 2016-06-17 | 2016-11-09 | 延长油田股份有限公司 | 一种致密砂岩储层可压裂性评价方法 |
CN106194141A (zh) * | 2016-08-19 | 2016-12-07 | 中国石油天然气股份有限公司 | 基于体积压裂下缝网形成的工艺方法 |
CN106326591A (zh) * | 2016-08-31 | 2017-01-11 | 西南石油大学 | 水力压裂过程中裂缝内压裂液的压力场获取方法及装置 |
CN106650100A (zh) * | 2016-12-23 | 2017-05-10 | 西南石油大学 | 一种实验页岩藏储层水平井交替体积压裂方法 |
CN106761647A (zh) * | 2017-01-13 | 2017-05-31 | 黄仲尧 | 一种估算页岩储层压后平面改造面积的方法 |
CN106869911A (zh) * | 2017-02-24 | 2017-06-20 | 中石化重庆涪陵页岩气勘探开发有限公司 | 一种描述页岩储层可压性的评价方法 |
CN106909758A (zh) * | 2017-04-10 | 2017-06-30 | 西南石油大学 | 一种致密油储层水平井多段分簇射孔位置优化设计的新方法 |
CN106930743A (zh) * | 2017-05-06 | 2017-07-07 | 东北石油大学 | 一种水平井分段多簇压裂射孔簇位置优化设计方法 |
CN106951659A (zh) * | 2017-04-05 | 2017-07-14 | 中国石油大学(北京) | 一种高效实现煤层及煤层顶部砂岩合压的射孔方案 |
CN106979000A (zh) * | 2017-05-06 | 2017-07-25 | 东北石油大学 | 水平井分段多簇压裂各射孔簇破裂压力计算方法 |
CN106979001A (zh) * | 2017-06-06 | 2017-07-25 | 华美孚泰油气增产技术服务有限责任公司 | 厚层砂砾岩水平井立体缝网压裂优化方法 |
CN107045582A (zh) * | 2017-05-06 | 2017-08-15 | 东北石油大学 | 水平井分段多簇压裂诱导应力计算方法 |
CN107044277A (zh) * | 2017-06-06 | 2017-08-15 | 西南石油大学 | 低渗透非均质油藏水平井重复压裂增产潜力评价方法 |
CN107203667A (zh) * | 2017-05-23 | 2017-09-26 | 西南石油大学 | 水平井段内多簇压裂优化方法及系统 |
CN107237616A (zh) * | 2017-07-21 | 2017-10-10 | 中国石油集团川庆钻探工程有限公司长庆井下技术作业公司 | 一种水平井体积压裂簇间距优化方法 |
CN107476791A (zh) * | 2016-06-07 | 2017-12-15 | 中国石油化工股份有限公司 | 一种页岩气水平井分段压裂变密度簇射孔方法及射孔枪 |
CN107506840A (zh) * | 2016-06-14 | 2017-12-22 | 中国石油化工股份有限公司 | 一种页岩气水平井射孔簇间距的优化方法 |
CN107524432A (zh) * | 2016-06-21 | 2017-12-29 | 中国石油化工股份有限公司 | 一种压裂施工中的实时调参方法 |
CN107869344A (zh) * | 2016-09-23 | 2018-04-03 | 中国石油化工股份有限公司 | 低渗稠油油藏的储层改造方法 |
CN107916917A (zh) * | 2016-10-10 | 2018-04-17 | 中国石油化工股份有限公司 | 提高水力压裂双翼裂缝缝长的方法 |
CN108629089A (zh) * | 2018-04-17 | 2018-10-09 | 西南交通大学 | 一种基于三维连续-离散单元耦合的强夯加固地基模拟方法 |
CN108868748A (zh) * | 2018-04-28 | 2018-11-23 | 中国石油化工股份有限公司江汉油田分公司石油工程技术研究院 | 一种页岩气水平井重复压裂裂缝开启压力的计算方法 |
CN108984877A (zh) * | 2018-07-03 | 2018-12-11 | 中国石油化工股份有限公司胜利油田分公司石油工程技术研究院 | 一种致密储层复杂压裂裂缝形成条件的评估方法 |
CN109359376A (zh) * | 2018-10-10 | 2019-02-19 | 北京科技大学 | 水力压裂裂缝在页岩储层天然裂缝界面扩展判识方法 |
CN109472037A (zh) * | 2017-09-08 | 2019-03-15 | 中国石油化工股份有限公司 | 页岩气储层人工裂缝参数优选方法及系统 |
CN109632459A (zh) * | 2018-11-14 | 2019-04-16 | 中石化重庆涪陵页岩气勘探开发有限公司 | 一种页岩可压裂性评价方法 |
CN109723435A (zh) * | 2018-12-12 | 2019-05-07 | 长江大学 | 确定页岩气藏水平井重复压裂前现地应力场的方法及设备 |
CN109992864A (zh) * | 2019-03-22 | 2019-07-09 | 成都理工大学 | 非常规双重介质储层体积压裂数值模拟及参数优化方法 |
CN110083885A (zh) * | 2019-04-04 | 2019-08-02 | 中国石油大学(华东) | 一种体积压裂水平井簇间距范围确定方法和装置 |
CN110529088A (zh) * | 2019-08-30 | 2019-12-03 | 西南石油大学 | 一种基于薄片鉴定的岩石可压裂性剖面建立方法 |
CN111472730A (zh) * | 2020-05-06 | 2020-07-31 | 中国石油天然气股份有限公司 | 大段多簇压裂的射孔方案确定方法 |
US10761241B2 (en) | 2019-02-25 | 2020-09-01 | University Of Science And Technology Beijing | Methods for optimizing bunch distance of fractured horizontal wells of shale gas |
WO2020173253A1 (zh) * | 2019-02-25 | 2020-09-03 | 北京科技大学 | 页岩气压裂水平井簇间距优化方法 |
CN111797575A (zh) * | 2020-07-23 | 2020-10-20 | 西南石油大学 | 一种深层页岩压裂工艺参数优化方法 |
CN111898065A (zh) * | 2020-06-18 | 2020-11-06 | 长江大学 | 一种页岩脆性指数计算方法、设备及可读存储介质 |
CN112096362A (zh) * | 2020-09-27 | 2020-12-18 | 西南石油大学 | 非常规储层多簇射孔竞争起裂与扩展模拟方法及装置 |
CN112177598A (zh) * | 2019-07-01 | 2021-01-05 | 中国石油化工股份有限公司 | 一种考虑压裂液压缩性的地层起裂压力预测方法 |
CN112324412A (zh) * | 2020-11-02 | 2021-02-05 | 中国石油化工股份有限公司 | 一种体积压裂形成复杂缝网的方法 |
CN112343589A (zh) * | 2020-11-06 | 2021-02-09 | 武汉大学 | 研究水力压裂复杂缝网微裂缝相互作用机理的试验方法 |
CN112377104A (zh) * | 2020-10-30 | 2021-02-19 | 中国石油天然气集团有限公司 | 动态调整致密砂岩气藏水平井实钻轨迹的方法及装置 |
CN112647935A (zh) * | 2019-10-12 | 2021-04-13 | 中国石油化工股份有限公司 | 压裂裂缝参数计算方法及系统 |
CN112765807A (zh) * | 2021-01-15 | 2021-05-07 | 西南石油大学 | 一种裂缝穿越海陆过渡相岩性界面偏转角计算方法 |
CN112836446A (zh) * | 2021-01-28 | 2021-05-25 | 重庆科技学院 | 一种致密储层水平井分簇限流布孔优化设计方法 |
CN112878977A (zh) * | 2021-03-26 | 2021-06-01 | 重庆科技学院 | 一种致密储层水平井多簇压裂射孔簇间距优化方法 |
CN113033021A (zh) * | 2021-04-21 | 2021-06-25 | 北京艾迪博科油气技术有限公司 | 致密储层的改造评估方法和装置 |
CN113107362A (zh) * | 2021-05-31 | 2021-07-13 | 中国矿业大学 | 一种促进裂隙错动提升自支撑能力的立井施工方法 |
CN113389534A (zh) * | 2021-07-21 | 2021-09-14 | 西南石油大学 | 一种水平井密切割压裂裂缝扩展预测与设计参数优化方法 |
CN113431542A (zh) * | 2020-07-27 | 2021-09-24 | 中国石油化工股份有限公司 | 水平井压裂裂缝干扰强度的计算方法 |
CN113685162A (zh) * | 2021-07-22 | 2021-11-23 | 中国石油大学(北京) | 压裂参数的确定方法、装置、设备及存储介质 |
CN113756794A (zh) * | 2020-06-05 | 2021-12-07 | 中国石油天然气股份有限公司 | 分簇射孔器分布方式的确定方法、装置、设备及介质 |
CN113927738A (zh) * | 2021-10-15 | 2022-01-14 | 中国重汽集团青岛重工有限公司 | 一种搅拌罐螺旋叶片的优化设计方法及成型方法 |
CN114508334A (zh) * | 2020-11-17 | 2022-05-17 | 中国石油化工股份有限公司 | 一种基于三维地应力场分布的溶洞循缝沟通技术确定方法 |
CN115584963A (zh) * | 2022-09-20 | 2023-01-10 | 西南石油大学 | 一种非常规储层可压裂性综合评价方法 |
CN115639083A (zh) * | 2022-12-06 | 2023-01-24 | 新疆斐德莱布能源科技有限公司 | 一种矿场级的真三轴水力压裂模拟实验方法及装置 |
CN116306330A (zh) * | 2022-11-23 | 2023-06-23 | 西南石油大学 | 一种深层断层附近页岩气压裂簇间距优化方法 |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110005380B (zh) * | 2019-04-11 | 2020-08-11 | 中国石油大学(北京) | 非均质性页岩非均匀分簇射孔优化方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130048282A1 (en) * | 2011-08-23 | 2013-02-28 | David M. Adams | Fracturing Process to Enhance Propping Agent Distribution to Maximize Connectivity Between the Formation and the Wellbore |
CN103256046A (zh) * | 2013-04-28 | 2013-08-21 | 北京大学 | 非常规油气藏水平井全缝长压裂参数模拟的方法及装置 |
CN103256035A (zh) * | 2013-05-22 | 2013-08-21 | 中国石化集团华北石油局 | 一种致密气田水平井压裂缝地质设计方法 |
CN103527163A (zh) * | 2013-09-24 | 2014-01-22 | 西南石油大学 | 一种致密储层水平井体积压裂工艺 |
CN103953323A (zh) * | 2014-05-08 | 2014-07-30 | 西南石油大学 | 一种水平井产生缝网的水力压裂工艺 |
-
2015
- 2015-10-09 CN CN201510648540.5A patent/CN105201479B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130048282A1 (en) * | 2011-08-23 | 2013-02-28 | David M. Adams | Fracturing Process to Enhance Propping Agent Distribution to Maximize Connectivity Between the Formation and the Wellbore |
CN103256046A (zh) * | 2013-04-28 | 2013-08-21 | 北京大学 | 非常规油气藏水平井全缝长压裂参数模拟的方法及装置 |
CN103256035A (zh) * | 2013-05-22 | 2013-08-21 | 中国石化集团华北石油局 | 一种致密气田水平井压裂缝地质设计方法 |
CN103527163A (zh) * | 2013-09-24 | 2014-01-22 | 西南石油大学 | 一种致密储层水平井体积压裂工艺 |
CN103953323A (zh) * | 2014-05-08 | 2014-07-30 | 西南石油大学 | 一种水平井产生缝网的水力压裂工艺 |
Non-Patent Citations (4)
Title |
---|
尹建等: "水平井分段压裂射孔间距优化方法", 《石油钻探技术》 * |
曾凡辉等: "致密砂岩气藏水平井分段压裂优化设计与应用", 《石油学报》 * |
蔡儒帅等: "页岩储层可压性剖面模型研究与应用", 《大庆石油地质与开发》 * |
郭建春等: "川西页岩-砂岩交互水平井压裂参数优化设计", 《石油学报》 * |
Cited By (85)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105672973A (zh) * | 2016-01-26 | 2016-06-15 | 西南石油大学 | 一种煤层气藏整体体积压裂优化设计方法 |
CN105735961A (zh) * | 2016-02-18 | 2016-07-06 | 西南石油大学 | 一种低渗致密储层分层压裂级数优选方法 |
CN105735961B (zh) * | 2016-02-18 | 2017-12-05 | 西南石油大学 | 一种低渗致密储层分层压裂级数优选方法 |
CN105735960A (zh) * | 2016-03-22 | 2016-07-06 | 西南石油大学 | 一种低渗透油气藏水平井分段多簇压裂簇间距优化方法 |
CN105735960B (zh) * | 2016-03-22 | 2017-05-17 | 西南石油大学 | 一种低渗透油气藏水平井分段多簇压裂簇间距优化方法 |
CN107476791B (zh) * | 2016-06-07 | 2019-05-17 | 中国石油化工股份有限公司 | 一种页岩气水平井分段压裂变密度簇射孔方法及射孔枪 |
CN107476791A (zh) * | 2016-06-07 | 2017-12-15 | 中国石油化工股份有限公司 | 一种页岩气水平井分段压裂变密度簇射孔方法及射孔枪 |
CN107506840A (zh) * | 2016-06-14 | 2017-12-22 | 中国石油化工股份有限公司 | 一种页岩气水平井射孔簇间距的优化方法 |
CN106097125A (zh) * | 2016-06-17 | 2016-11-09 | 延长油田股份有限公司 | 一种致密砂岩储层可压裂性评价方法 |
CN107524432A (zh) * | 2016-06-21 | 2017-12-29 | 中国石油化工股份有限公司 | 一种压裂施工中的实时调参方法 |
CN107524432B (zh) * | 2016-06-21 | 2019-07-02 | 中国石油化工股份有限公司 | 一种压裂施工中的实时调参方法 |
CN106194141A (zh) * | 2016-08-19 | 2016-12-07 | 中国石油天然气股份有限公司 | 基于体积压裂下缝网形成的工艺方法 |
CN106326591B (zh) * | 2016-08-31 | 2019-09-06 | 西南石油大学 | 水力压裂过程中裂缝内压裂液的压力场获取方法及装置 |
CN106326591A (zh) * | 2016-08-31 | 2017-01-11 | 西南石油大学 | 水力压裂过程中裂缝内压裂液的压力场获取方法及装置 |
CN107869344A (zh) * | 2016-09-23 | 2018-04-03 | 中国石油化工股份有限公司 | 低渗稠油油藏的储层改造方法 |
CN107869344B (zh) * | 2016-09-23 | 2020-06-02 | 中国石油化工股份有限公司 | 低渗稠油油藏的储层改造方法 |
CN107916917A (zh) * | 2016-10-10 | 2018-04-17 | 中国石油化工股份有限公司 | 提高水力压裂双翼裂缝缝长的方法 |
CN106650100B (zh) * | 2016-12-23 | 2020-01-10 | 西南石油大学 | 一种实验页岩藏储层水平井交替体积压裂方法 |
CN106650100A (zh) * | 2016-12-23 | 2017-05-10 | 西南石油大学 | 一种实验页岩藏储层水平井交替体积压裂方法 |
CN106761647B (zh) * | 2017-01-13 | 2020-08-14 | 中国石油化工股份有限公司 | 一种估算页岩储层压后平面改造面积的方法 |
CN106761647A (zh) * | 2017-01-13 | 2017-05-31 | 黄仲尧 | 一种估算页岩储层压后平面改造面积的方法 |
CN106869911B (zh) * | 2017-02-24 | 2020-08-04 | 中石化重庆涪陵页岩气勘探开发有限公司 | 一种描述页岩储层可压性的评价方法 |
CN106869911A (zh) * | 2017-02-24 | 2017-06-20 | 中石化重庆涪陵页岩气勘探开发有限公司 | 一种描述页岩储层可压性的评价方法 |
CN106951659A (zh) * | 2017-04-05 | 2017-07-14 | 中国石油大学(北京) | 一种高效实现煤层及煤层顶部砂岩合压的射孔方案 |
CN106909758B (zh) * | 2017-04-10 | 2018-12-04 | 西南石油大学 | 一种致密油储层水平井多段分簇射孔位置优化设计的方法 |
CN106909758A (zh) * | 2017-04-10 | 2017-06-30 | 西南石油大学 | 一种致密油储层水平井多段分簇射孔位置优化设计的新方法 |
CN106930743A (zh) * | 2017-05-06 | 2017-07-07 | 东北石油大学 | 一种水平井分段多簇压裂射孔簇位置优化设计方法 |
CN106979000A (zh) * | 2017-05-06 | 2017-07-25 | 东北石油大学 | 水平井分段多簇压裂各射孔簇破裂压力计算方法 |
CN107045582A (zh) * | 2017-05-06 | 2017-08-15 | 东北石油大学 | 水平井分段多簇压裂诱导应力计算方法 |
CN107045582B (zh) * | 2017-05-06 | 2019-10-25 | 东北石油大学 | 水平井分段多簇压裂诱导应力计算方法 |
CN107203667A (zh) * | 2017-05-23 | 2017-09-26 | 西南石油大学 | 水平井段内多簇压裂优化方法及系统 |
CN107203667B (zh) * | 2017-05-23 | 2020-09-25 | 西南石油大学 | 水平井段内多簇压裂优化方法及系统 |
CN107044277B (zh) * | 2017-06-06 | 2019-04-26 | 西南石油大学 | 低渗透非均质油藏水平井重复压裂增产潜力评价方法 |
CN107044277A (zh) * | 2017-06-06 | 2017-08-15 | 西南石油大学 | 低渗透非均质油藏水平井重复压裂增产潜力评价方法 |
CN106979001A (zh) * | 2017-06-06 | 2017-07-25 | 华美孚泰油气增产技术服务有限责任公司 | 厚层砂砾岩水平井立体缝网压裂优化方法 |
CN107237616A (zh) * | 2017-07-21 | 2017-10-10 | 中国石油集团川庆钻探工程有限公司长庆井下技术作业公司 | 一种水平井体积压裂簇间距优化方法 |
CN109472037A (zh) * | 2017-09-08 | 2019-03-15 | 中国石油化工股份有限公司 | 页岩气储层人工裂缝参数优选方法及系统 |
CN108629089B (zh) * | 2018-04-17 | 2021-08-03 | 西南交通大学 | 一种基于连续-离散单元耦合的强夯加固地基模拟方法 |
CN108629089A (zh) * | 2018-04-17 | 2018-10-09 | 西南交通大学 | 一种基于三维连续-离散单元耦合的强夯加固地基模拟方法 |
CN108868748A (zh) * | 2018-04-28 | 2018-11-23 | 中国石油化工股份有限公司江汉油田分公司石油工程技术研究院 | 一种页岩气水平井重复压裂裂缝开启压力的计算方法 |
CN108868748B (zh) * | 2018-04-28 | 2022-03-25 | 中国石油化工股份有限公司 | 一种页岩气水平井重复压裂裂缝开启压力的计算方法 |
CN108984877A (zh) * | 2018-07-03 | 2018-12-11 | 中国石油化工股份有限公司胜利油田分公司石油工程技术研究院 | 一种致密储层复杂压裂裂缝形成条件的评估方法 |
CN109359376A (zh) * | 2018-10-10 | 2019-02-19 | 北京科技大学 | 水力压裂裂缝在页岩储层天然裂缝界面扩展判识方法 |
CN109632459B (zh) * | 2018-11-14 | 2021-04-13 | 中国石油化工股份有限公司 | 一种页岩可压裂性评价方法 |
CN109632459A (zh) * | 2018-11-14 | 2019-04-16 | 中石化重庆涪陵页岩气勘探开发有限公司 | 一种页岩可压裂性评价方法 |
CN109723435B (zh) * | 2018-12-12 | 2022-07-22 | 长江大学 | 确定页岩气藏水平井重复压裂前现地应力场的方法及设备 |
CN109723435A (zh) * | 2018-12-12 | 2019-05-07 | 长江大学 | 确定页岩气藏水平井重复压裂前现地应力场的方法及设备 |
US10761241B2 (en) | 2019-02-25 | 2020-09-01 | University Of Science And Technology Beijing | Methods for optimizing bunch distance of fractured horizontal wells of shale gas |
WO2020173253A1 (zh) * | 2019-02-25 | 2020-09-03 | 北京科技大学 | 页岩气压裂水平井簇间距优化方法 |
CN109992864A (zh) * | 2019-03-22 | 2019-07-09 | 成都理工大学 | 非常规双重介质储层体积压裂数值模拟及参数优化方法 |
US11567233B2 (en) * | 2019-03-22 | 2023-01-31 | Chengdu University Of Technology | Numerical simulation and parameter optimization method for volumetric fracturing of unconventional dual medium reservoir |
CN110083885A (zh) * | 2019-04-04 | 2019-08-02 | 中国石油大学(华东) | 一种体积压裂水平井簇间距范围确定方法和装置 |
CN112177598B (zh) * | 2019-07-01 | 2024-03-08 | 中国石油化工股份有限公司 | 一种考虑压裂液压缩性的地层起裂压力预测方法 |
CN112177598A (zh) * | 2019-07-01 | 2021-01-05 | 中国石油化工股份有限公司 | 一种考虑压裂液压缩性的地层起裂压力预测方法 |
CN110529088A (zh) * | 2019-08-30 | 2019-12-03 | 西南石油大学 | 一种基于薄片鉴定的岩石可压裂性剖面建立方法 |
CN112647935A (zh) * | 2019-10-12 | 2021-04-13 | 中国石油化工股份有限公司 | 压裂裂缝参数计算方法及系统 |
CN111472730A (zh) * | 2020-05-06 | 2020-07-31 | 中国石油天然气股份有限公司 | 大段多簇压裂的射孔方案确定方法 |
CN111472730B (zh) * | 2020-05-06 | 2022-07-05 | 中国石油天然气股份有限公司 | 大段多簇压裂的射孔方案确定方法 |
CN113756794B (zh) * | 2020-06-05 | 2024-03-01 | 中国石油天然气股份有限公司 | 分簇射孔器分布方式的确定方法、装置、设备及介质 |
CN113756794A (zh) * | 2020-06-05 | 2021-12-07 | 中国石油天然气股份有限公司 | 分簇射孔器分布方式的确定方法、装置、设备及介质 |
CN111898065B (zh) * | 2020-06-18 | 2023-06-06 | 长江大学 | 一种页岩脆性指数计算方法、设备及可读存储介质 |
CN111898065A (zh) * | 2020-06-18 | 2020-11-06 | 长江大学 | 一种页岩脆性指数计算方法、设备及可读存储介质 |
CN111797575A (zh) * | 2020-07-23 | 2020-10-20 | 西南石油大学 | 一种深层页岩压裂工艺参数优化方法 |
CN113431542A (zh) * | 2020-07-27 | 2021-09-24 | 中国石油化工股份有限公司 | 水平井压裂裂缝干扰强度的计算方法 |
CN112096362A (zh) * | 2020-09-27 | 2020-12-18 | 西南石油大学 | 非常规储层多簇射孔竞争起裂与扩展模拟方法及装置 |
CN112096362B (zh) * | 2020-09-27 | 2021-10-22 | 西南石油大学 | 非常规储层多簇射孔竞争起裂与扩展模拟方法及装置 |
CN112377104A (zh) * | 2020-10-30 | 2021-02-19 | 中国石油天然气集团有限公司 | 动态调整致密砂岩气藏水平井实钻轨迹的方法及装置 |
CN112324412A (zh) * | 2020-11-02 | 2021-02-05 | 中国石油化工股份有限公司 | 一种体积压裂形成复杂缝网的方法 |
CN112343589A (zh) * | 2020-11-06 | 2021-02-09 | 武汉大学 | 研究水力压裂复杂缝网微裂缝相互作用机理的试验方法 |
CN114508334A (zh) * | 2020-11-17 | 2022-05-17 | 中国石油化工股份有限公司 | 一种基于三维地应力场分布的溶洞循缝沟通技术确定方法 |
CN112765807A (zh) * | 2021-01-15 | 2021-05-07 | 西南石油大学 | 一种裂缝穿越海陆过渡相岩性界面偏转角计算方法 |
CN112765807B (zh) * | 2021-01-15 | 2022-03-01 | 西南石油大学 | 一种裂缝穿越海陆过渡相岩性界面偏转角计算方法 |
CN112836446B (zh) * | 2021-01-28 | 2022-04-26 | 重庆科技学院 | 一种致密储层水平井分簇限流布孔优化设计方法 |
CN112836446A (zh) * | 2021-01-28 | 2021-05-25 | 重庆科技学院 | 一种致密储层水平井分簇限流布孔优化设计方法 |
CN112878977A (zh) * | 2021-03-26 | 2021-06-01 | 重庆科技学院 | 一种致密储层水平井多簇压裂射孔簇间距优化方法 |
CN113033021A (zh) * | 2021-04-21 | 2021-06-25 | 北京艾迪博科油气技术有限公司 | 致密储层的改造评估方法和装置 |
CN113107362A (zh) * | 2021-05-31 | 2021-07-13 | 中国矿业大学 | 一种促进裂隙错动提升自支撑能力的立井施工方法 |
CN113107362B (zh) * | 2021-05-31 | 2023-09-05 | 中国矿业大学 | 一种促进裂隙错动提升自支撑能力的立井施工方法 |
CN113389534A (zh) * | 2021-07-21 | 2021-09-14 | 西南石油大学 | 一种水平井密切割压裂裂缝扩展预测与设计参数优化方法 |
CN113685162A (zh) * | 2021-07-22 | 2021-11-23 | 中国石油大学(北京) | 压裂参数的确定方法、装置、设备及存储介质 |
CN113927738B (zh) * | 2021-10-15 | 2023-02-17 | 中国重汽集团青岛重工有限公司 | 一种搅拌罐螺旋叶片的优化设计方法及成型方法 |
CN113927738A (zh) * | 2021-10-15 | 2022-01-14 | 中国重汽集团青岛重工有限公司 | 一种搅拌罐螺旋叶片的优化设计方法及成型方法 |
CN115584963A (zh) * | 2022-09-20 | 2023-01-10 | 西南石油大学 | 一种非常规储层可压裂性综合评价方法 |
CN116306330A (zh) * | 2022-11-23 | 2023-06-23 | 西南石油大学 | 一种深层断层附近页岩气压裂簇间距优化方法 |
CN115639083A (zh) * | 2022-12-06 | 2023-01-24 | 新疆斐德莱布能源科技有限公司 | 一种矿场级的真三轴水力压裂模拟实验方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN105201479B (zh) | 2017-10-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105201479B (zh) | 一种页岩储层水平井分段压裂射孔簇参数优化设计方法 | |
CN108319756B (zh) | 一种致密储层体积压裂缝网扩展模拟及表征方法 | |
CN105201484A (zh) | 一种直井分层压裂层段优选及施工参数优化设计方法 | |
CN105156103B (zh) | 一种岩屑-岩心-井眼-储层多尺度的页岩储层三维可压裂性评价方法 | |
CN106547034B (zh) | 一种确定致密储层岩石脆性指数的方法 | |
CN105134156B (zh) | 一种用于致密砂岩储层三维可压裂性模型的建模方法 | |
CN104750896B (zh) | 一种缝洞型碳酸盐岩油藏数值模拟方法 | |
CN108009705A (zh) | 一种基于支持向量机技术的页岩储层可压性评价方法 | |
CN106909758B (zh) | 一种致密油储层水平井多段分簇射孔位置优化设计的方法 | |
CN108681635B (zh) | 一种致密储层体积压裂可压性评价方法 | |
Bakhshi et al. | Numerical modeling and lattice method for characterizing hydraulic fracture propagation: a review of the numerical, experimental, and field studies | |
CN107044277A (zh) | 低渗透非均质油藏水平井重复压裂增产潜力评价方法 | |
Li et al. | A fully coupled model for hydraulic-fracture growth during multiwell-fracturing treatments: enhancing fracture complexity | |
CN106437674A (zh) | 仿水平井注水开发井网适配方法 | |
CN110469303B (zh) | 一种基于四类改造体积的体积压裂参数优化设计方法 | |
CN105735961B (zh) | 一种低渗致密储层分层压裂级数优选方法 | |
CN110094196B (zh) | 一种碳酸盐岩裸眼水平井分段酸压效果评价方法 | |
CN108661616B (zh) | 一种适用于砂岩油气藏的分层压裂方法 | |
Song et al. | Hydraulic fracturing simulation of complex fractures growth in naturally fractured shale gas reservoir | |
Lin et al. | Numerical simulation of the influence of stimulated reservoir volume on in-situ stress field | |
Li et al. | Improving fracture initiation predictions of a horizontal wellbore in laminated anisotropy shales | |
Yang et al. | Influence of reservoirs/interlayers thickness on hydraulic fracture propagation laws in low-permeability layered rocks | |
Jizhou et al. | Optimization of integrated geological-engineering design of volume fracturing with fan-shaped well pattern | |
CN106761647A (zh) | 一种估算页岩储层压后平面改造面积的方法 | |
CN112012710A (zh) | 一种砂泥岩交互地层水平井压裂三维裂缝扩展模拟方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
OL01 | Intention to license declared | ||
OL01 | Intention to license declared |