Processing math: 100%

滑脱和应力敏感效应对页岩气开发动态影响的数值模拟研究

杨金辉, 李立, 李钟洋, 鞠斌山

杨金辉, 李立, 李钟洋, 鞠斌山. 滑脱和应力敏感效应对页岩气开发动态影响的数值模拟研究[J]. 石油钻探技术, 2017, 45(1): 83-90. DOI: 10.11911/syztjs.201701015
引用本文: 杨金辉, 李立, 李钟洋, 鞠斌山. 滑脱和应力敏感效应对页岩气开发动态影响的数值模拟研究[J]. 石油钻探技术, 2017, 45(1): 83-90. DOI: 10.11911/syztjs.201701015
YANG Jinhui, LI Li, LI Zhongyang, JU Binshan. Numerical Simulation on the Effects of Slippage and Stress Sensibility on the Performance of Shale Gas Development[J]. Petroleum Drilling Techniques, 2017, 45(1): 83-90. DOI: 10.11911/syztjs.201701015
Citation: YANG Jinhui, LI Li, LI Zhongyang, JU Binshan. Numerical Simulation on the Effects of Slippage and Stress Sensibility on the Performance of Shale Gas Development[J]. Petroleum Drilling Techniques, 2017, 45(1): 83-90. DOI: 10.11911/syztjs.201701015

滑脱和应力敏感效应对页岩气开发动态影响的数值模拟研究

基金项目: 

国家重大科技专项“大型油气田及煤层气开发”之课题“新一代油藏数值模拟软件”(编号:2011ZX05009-006)和中央高校基本科研业务费专项资金(编号:2652015142)联合资助。

详细信息
    作者简介:

    杨金辉(1993-),男,山东东营人,2015年毕业于中国地质大学(武汉)石油工程专业,油气田开发工程专业在读硕士研究生,主要从事油气田开发工程技术研究工作。通信作者鞠斌山,

    杨金辉(1993-),男,山东东营人,2015年毕业于中国地质大学(武汉)石油工程专业,油气田开发工程专业在读硕士研究生,主要从事油气田开发工程技术研究工作。通信作者鞠斌山,

  • 中图分类号: TE33+2

Numerical Simulation on the Effects of Slippage and Stress Sensibility on the Performance of Shale Gas Development

  • 摘要: 滑脱和应力敏感效应是影响页岩气藏开发动态的2个重要因素,为定量研究其对产气量的影响,运用等效孔隙介质模型,建立了同时考虑滑脱和应力敏感效应的页岩气藏两相渗流数学模型,采用有限差分方法对该模型进行了数值求解,用Fortran语言开发了页岩气藏数值模拟器,在定产量和定压衰竭式开采条件下应用该模拟器研究了滑脱与应力敏感效应对页岩气井产能的影响规律,并将计算结果与美国Haynesville页岩气藏W3井的实际生产数据进行了对比。结果表明:衰竭式开发导致页岩气藏的压力和渗透率均降低,生产井附近储层压力和渗透率急剧降低;在投产第1年内,产气量虽迅速下降,但滑脱和应力敏感效应对产气量影响较小;生产500 d后,滑脱和应力敏感效应对产气量产生明显影响,滑脱效应使产气量增加,应力敏感导致产气量降低;模拟结果与生产数据的衰减趋势吻合度较高,证明该页岩气藏数值模拟器计算结果具有较高的可靠性。页岩气藏开发数值模拟器的成功研制,为页岩气井产能预测提供了一种新的有效手段。
    Abstract: Gas slippage and stress sensibility are two important factors which influence the development performances in shale gas reservoirs.To highlight their quantitative effects on gas production rates,a mathematical model with consideration of both slippage and stress sensibility in a two-phase flow in shale gas reservoirs was set up by using an equivalent porous media model, in which a finite difference method was used to solve the model. In addition, a numerical simulator of shale gas reservoirs was developed in Fortran(computer language).The effect of stress sensibility and gas slippage on gas well productivity was quantitatively studied by using the innovative simulator at a constant production rate and a specified bottom flow pressure in declining production.Numerical simulation results were then compared with field production data of Well W3 in a Haynesville shale gas reservoir.The numerical results showed that both formation pressure and permeability declined during depletion development of the gas reservoir and dramatical decline was observed near the wellbore region.In the first production year,shale gas production rate declined sharply;but slippage and stress sensibility displayed minor effect on the production rate.The effect of stress sensibility and gas slippage on production became more significant after 500 days of production.The decline curve closely tracked the simulation and history matching,and it thus validated the result calculated from a shale gas reservoir simulator with relatively high reliability.The successful development of a numerical simulator for shale gas reservoirs migth provide an innovative and effective approach for predicting gas production rates and decline curves.
  • 非常规油气储层采用水力压裂技术进行改造时,支撑剂随压裂液进入地层并在裂缝中运移、沉降,其在裂缝中的铺置形态极大地影响裂缝的导流能力进而决定储层改造效果[1-4]。在页岩储层压裂中,大规模采用低黏滑溜水容易形成复杂缝网,由于储层两向应力差较大,导致分支裂缝宽度较小,支撑剂难以进入(放射性示踪结果证实,泵入地层的支撑剂仍主要分布在与井筒直连的主缝中)。为了改善主裂缝内的铺置效果,需要研究支撑剂在主缝内的运移与铺置规律以及影响砂堤最终堆积形态的主控因素 [5-10]

    截至目前,针对支撑剂的运移与铺置问题,国内外已进行了大量室内物模试验与数值模拟研究。不过,研究人员仍以视频记录支撑剂铺置过程为主,考虑裂缝形态、压裂液黏度、支撑剂类型、加砂浓度、注入口位置、泵注排量、泵注工艺等因素,通过形态参数提取、砂堤曲线拟合演化等手段,宏观描述支撑剂的铺置过程,定性抽提出不同因素作用下的砂堤堆积模式[11-15],并未从支撑剂运动的角度解释砂堤生长。数字粒子图像测速(particle image velocimetry,PIV)技术因其具有无接触测量的特点,在固液两相流测速中得到大量应用[16]。因此,Li Genghong等人[17-20]将PIV技术用于支撑剂运动程测试,但发现存在测试效果差、不能同时兼顾固相与液相运动和单一颗粒速度难以准确测定等问题。

    基于以上情况,笔者在传统平板可视化支撑剂输送试验中采用PIV技术,捕获支撑剂颗粒与示踪粒子的瞬态运动状态,通过基于曝光差异的分离算法将支撑剂颗粒图像与示踪粒子图像分离,运用粒子跟踪测速(particle track velocimetry,PTV)技术实现对压裂液–支撑剂两相运动场的解析,并分析了泵注排量、压裂液黏度、注入位置、支撑剂类型对砂堤堆积特征区域形态的影响,以期为压裂工艺参数设计优化提供参考依据。

    在平板裂缝内开展不同试验参数下的支撑剂宏观铺置试验,分析砂堤生长过程。在砂堤形态不再发生变化后,采用PIV技术测定裂缝内不同区域两相粒子的运动速度。采集的两相粒子图像经过增强–分离处理得到仅含支撑剂颗粒的单相粒子图像,并基于PTV算法对支撑剂颗粒运动进行解析以获得支撑剂颗粒运动场。结合铺置试验结果,分析不同因素对砂堤形态的影响,总结2种典型的砂堤生长模式。

    基于缝内流动相似准则,将现场泵注排量换算为试验装置缝内流速:

    ve=500hewehfwfnvf (1)

    式中:ve为试验装置缝内流速,L/min;vf为现场施工排量,m3/min;he为裂缝装置缝高,m;we为裂缝装置缝宽,m;hf为真实压裂裂缝缝高,m;wf为真实压裂裂缝缝宽,m;n为裂缝数量。

    参照川南页岩气现场施工参数,排量12~16 m3/min,主体施工液体黏度3~5 mPa·s,进入地层后实际黏度2~3 mPa·s,页岩储层水力裂缝缝高10~12 m[5],平均缝宽6 mm,分6~12簇射孔[8],假设每簇均匀进液,根据式(1),换算出室内试验主体泵注排量为10.00~26.67 L/min。由于试验中支撑剂微观运动测试需要采用片状光源从裂缝顶部照亮缝内流场,砂比过大时,会造成过曝光现象,继而影响后续速度解释,确定本试验最大砂比为4%,且所有试验均在此砂比条件下进行。支撑剂运移与铺置试验方案见表1(试验组1,2和3研究泵注排量的影响,试验组4,5和6研究黏度的影响,试验组7,8和9研究注入位置的影响,试验组1和6研究支撑剂类型的影响)。

    表  1  支撑剂运移与铺置试验方案
    Table  1.  Scheme of proppant transport and placement test
    试验组泵注排量/(L∙min–1支撑剂压裂液黏度/(mPa·s)注入位置PIV拍摄区域(砂堤位置)
    112.040/70目陶粒2.5上、中、下孔前缘、中部、后缘
    218.040/70目陶粒2.5上、中、下孔前缘、中部
    324.040/70目陶粒2.5上、中、下孔前缘、中部
    412.070/140目石英砂6.0上、中、下孔中部
    512.070/140目石英砂1.0上、中、下孔中部
    612.070/140目石英砂2.5上、中、下孔中部、后缘
    712.070/140目石英砂2.5上孔前缘
    812.070/140目石英砂2.5中孔前缘
    912.070/140目石英砂2.5下孔前缘
    下载: 导出CSV 
    | 显示表格

    试验材料主要包括支撑剂、减阻剂(成分为聚丙烯酰胺)、滑溜水,如图1所示。支撑剂选用40/70目陶粒(粒径212~420 μm,视密度2 660 kg/m3)与70/140目石英砂(粒径104~212 μm,视密度2 630 kg/m3),试验前通过调整减阻剂的加量配制不同黏度的液体。

    图  1  支撑剂运移与铺置试验材料
    Figure  1.  Materials used in proppant transport and placement test

    试验采用自主研发的多尺度支撑剂输送物理模拟系统[21]进行。该系统主要由控制模块、循环泵注模块、平板裂缝模块、基于数字PIV的定量化测试模块、数据采集及辅助测试模块等5部分组成,如图2所示。

    图  2  多尺度支撑剂输送物理模拟系统示意
    Figure  2.  Physical simulation system of multi-scale proppant transport

    参考传统水力压裂形成对称双翼单缝设计了平板裂缝模块。该模块井筒装置与4块100 cm×30 cm裂缝单元拼接而成,裂缝总长400 cm,高30 cm,井筒设有3个等间距孔眼,孔眼间距10 cm,孔眼直径10 mm,采用球阀控制孔眼的开闭。平板裂缝由高强度铝框架、PC透明板、高透光亚克力条、橡胶密封条组成,通过加密螺栓紧固以保证良好的密封性,采用不同宽度的高透光亚克力条控制裂缝宽度,其高透光性保证了激光能从裂缝顶部照透整个缝内流场,如图3所示。

    图  3  裂缝模块紧固与顶部透光的实现
    Figure  3.  Tightening of fracture module and realization of top light transmission

    循环泵注模块由2个300 L带搅拌器的配液罐、1个300 L混砂搅拌罐、1个沉降循环罐、1个定量加砂装置、1台离心泵和2台螺杆泵组成,最大泵注排量可达150 L/min。测试与控制模块采用中控柜系统集成控制,能够在线调控泵注排量、配液罐和混砂搅拌罐的搅拌速率及配液罐补液速率,同时实时监测泵注流量和缝内压力。

    试验内容主要包括支撑剂宏观铺置试验、支撑剂微观运动测速和砂堤宏观演化–颗粒微观运动联合分析。

    在配液罐中配制300 L试验液体,并加入一定量的中空玻璃珠作为示踪粒子(其密度与液体相同,其运动场可等效为流体流动场),根据泵注排量与砂比计算加砂速率。打开系统循环阀门并采用小排量泵注试验液体,使平板裂缝内充满液体,并通过裂缝入口、出口处的排空气管排出循环中的气体。循环5 min后开始加砂,通过调整配液罐补液速率,保证配液罐内液量不变以维持恒定砂比。调整视频采集系统,保证每块平板裂缝都能清晰成像。

    铺置过程中,采用PIV系统对缝内两相流动进行运动测速,并基于支撑剂颗粒与示踪粒子曝光特性差异,运用图像处理技术与PTV算法,进行支撑剂单相图片分离与运动速度场解析,解析流程如图4所示(由于本文只针对支撑剂运动开展研究,不再对分离后的示踪粒子图片进行后续速度解析)。

    图  4  压裂液–支撑剂两相流场解析流程
    Figure  4.  Analytical process of fracturing fluid–proppant two-phase flow field

    对视频采集结果进行分析,选取相同时刻(t1t2t3)进行拼接得到完整的砂堤形态图,并统计砂堤平衡高度、前缘倾角与沉降剖面角,结合砂堤表面支撑剂运动速度场,分析不同条件下砂堤前缘、中部和后缘形态差异及其主控因素,总结砂堤生长模式。

    泵注过程中,选取相同时间节点对砂堤剖面进行记录,以展示支撑剂的铺置过程。泵注排量为12 L/min时,支撑剂在注入口附近堆积,并在第1块板内先增长至平衡高度;当前缘砂堤达到平衡高度后,支撑剂从砂堤顶部过流通道流出并向裂缝深处运移沉降,此时后第2块板内底部仍为凹形层状铺置,如图5(a)所示。泵注排量为18 L/min时,砂堤首次出现平衡高度的位置较12 L/min靠后,如图5(b)所示。泵注排量为24 L/min时,砂堤生长速度较慢,支撑剂沉降趋势减弱,部分支撑剂未能在裂缝内沉降,直接从裂缝尾端流出,如图5(c)所示。

    图  5  不同泵注排量下的支撑剂铺置过程
    Figure  5.  Proppant placement process under different displacements

    泵注清水时,支撑剂铺置过程与试验组1类似(排量12 L/min,压裂液黏度2.5 mPa·s,注入40/70目陶粒),如图6(a)所示。泵注2.5 mPa·s滑溜水时,平板裂缝内砂堤均匀铺置;泵注中前期,靠近入射口位置砂堤增长至中孔位置时下孔完全堵塞,中、上孔注入速度增大,使之后的支撑剂速度增大向后端运移;泵注后期,裂缝前端砂堤增长至平衡高度后,中孔注入的的部分支撑剂越过砂堤前缘坡面与上孔入射的支撑剂汇流向后端运移,一小部分回卷到入口,如图6(b)所示。泵注6.0 mPa·s低黏滑溜水时,由于裂缝尺寸较小,支撑剂难以在裂缝内沉降堆积,如图6(c)所示。

    图  6  不同黏度压裂液下的支撑剂铺置过程
    Figure  6.  Proppant placement process under different viscosity

    采用上孔注入时,在第1块板中部产生强烈的涡流,一部分回卷至注入位置附近堆积,一部分往裂缝远端输送。当砂堤生长至裂缝高度一半时,中孔、下孔堵塞导致上孔注入速率增加,到砂堤达一定高度后,受涡流主导的支撑剂无法平衡重力效应继续向砂堤顶部运移,支撑剂会在凹陷处开始缓慢堆积直至与后面均匀铺置的砂堤高度一致。泵注中后期,支撑剂难以堆积在入射口位置,大部分将被携带至远端形成层状铺置(见图7(a))。采用中孔注入时,支撑剂运移与铺置过程与上孔注入时类似(见图7(b)),由于注入位置高度降低,涡流的位置更靠近入口,回流支撑剂的占比减少,注入口下端砂堤更平缓,前缘砂堤形成速率更慢、高度更低。采用下孔注入时,砂堤呈现整体均匀增长,且入口位置无砂堤堆积,如图7(c)所示。

    图  7  不同注入位置下的支撑剂铺置过程
    Figure  7.  Proppant placement process at different inject positions

    在泵注排量12 L/min、压裂液黏度2.5 mPa·s条件下,40/70目陶粒与70/140目石英砂的支撑剂铺置过程如图5(a)、图6(b)所示。由于支撑剂粒径、密度差异,相同滑溜水携带条件下,支撑剂沉降速度不同,进而影响铺置过程。

    砂堤前缘形态很大程度上决定了裂缝与井筒的连通性,如何保证并提高近井筒区域支撑剂有效堆积一直是难题,下面选取砂堤前缘倾角作为砂堤前缘形态的评价参数进行试验。砂堤前缘倾角统计结果如图8所示。

    图  8  砂堤前缘倾角统计结果
    Figure  8.  Statistics of leading edge inclination angle of sand embankment

    图8可知,注入位置与泵注排量对砂堤前缘形态影响最大,其次是压裂液黏度、支撑剂粒径。泵注排量越大,砂堤前缘倾角越小。泵注排量从12 L/min增至18 L/min时,砂堤前缘倾角减小31.43%;排量从18 L/min增至24 L/min时,砂堤前缘倾角减小11.48%;泵注排量较大时,继续提高排量对砂堤前缘形态影响较小,降低排量更容易使支撑剂在裂缝前端堆积形成具有更高倾角的砂堤前缘。

    针对注入位置与泵注排量,开展了砂堤前缘支撑剂颗粒速度测试,结果表明:

    1)采用上孔注入时,支撑剂在裂缝顶部受强涡流作用,一部分支撑剂会回卷至裂缝底部,在注入孔下方形成有效堆积,使井筒与裂缝有最佳连通性;而采用下孔注入时,即使有支撑剂回卷至入口处,但裂缝底部的高速流体会将该部分支撑剂重新卷走往裂缝深部运移,支撑剂无法在注入孔上方形成有效堆积,进而使裂缝与井筒完全失去连接。不同注入位置下砂堤前缘近表面支撑剂颗粒的运动场如图9所示。

    图  9  不同注入位置下砂堤前缘近表面支撑剂颗粒运动场
    Figure  9.  Near-surface proppant particle motion field at leading edge of sand embankment at different injection positions

    2)泵注排量越大,裂缝前端受涡流的影响导致无砂区面积越大,支撑剂难以在入口位置形成堆积。泵注排量小时,无砂区呈现“缩颈”,携砂液平滑流进砂堤上方过流通道内,未出现大尺度涡流;随泵注排量增大,砂堤表面形状逐渐从“下凹”过渡至“上凸”,这是由于泵注排量增大导致出现强涡流,对砂堤前缘进行冲刷侵蚀。不同泵注排量下砂堤前缘近表面支撑剂颗粒的运动场如图10所示。

    图  10  不同排量下砂堤前缘近表面支撑剂颗粒运动场
    Figure  10.  Near-surface proppant particle motion field at leading edge of sand embankment under different displacements

    砂堤平衡高度决定了裂缝纵向有效支撑高度,笔者将平衡高度作为砂堤中部的形态主要评价参数。砂堤平衡高度统计结果如图11 所示。由图11可知:压裂液黏度与泵注排量对砂堤中部形态影响最大,其次是注入位置和支撑剂粒径;提高压裂液黏度,砂堤平衡高度呈先增大后减小的趋势。

    图  11  砂堤平衡高度统计结果
    Figure  11.  Statistics of equilibrium height of sand embankment

    测试不同黏度压裂液下砂堤表面支撑剂的运动速度(结果见图12),并统计支撑剂运动角度(见表2),以水平方向上下两侧颗粒运动角度来衡量支撑剂向前运动的趋势。角度分布结果表明:泵注清水时,支撑剂颗粒速度分布范围更大,支撑剂在砂堤表面出现了明显的迁跃运动,这是由于清水具有较大的摩阻,更容易卷起砂堤表面的支撑剂;而泵注滑溜水时,支撑剂颗粒速度更集中,速度方向具有较强的一致性,具有更强的向前运移的能力,对砂堤表面的冲蚀作用较弱。

    图  12  不同黏度压裂液下砂堤近表面支撑剂颗粒速度场
    Figure  12.  Velocity field of proppant particles near the surface of sand embankment under different viscosity
    表  2  不同黏度压裂液下砂堤表面支撑剂颗粒运动角度统计
    Table  2.  Statistics of movement angle of proppant particles near sand embankment surface under different viscosity
    压裂液黏度/
    (mPa·s)
    角度占比,%
    0°~10°和350°~360°0°~20°和340°~360°
    1.057.3690.68
    2.582.2296.13
    6.089.5697.02
    下载: 导出CSV 
    | 显示表格

    沉降坡面角反映砂堤后缘在裂缝内铺置的延伸距离,沉降坡面角越小,支撑剂运移距离越远,砂堤铺置长度越长。以沉降坡面角作为砂堤后缘形态特征描述参数,进行铺置试验,结果见图13图14

    图  13  砂堤沉降剖面角统计结果
    Figure  13.  Statistics of a sand embankment settlement profile angle
    图  14  2种砂堤后缘形态下支撑剂颗粒速度场
    Figure  14.  Velocity field of proppant particles under two kinds of sand embankment trailing shape

    试验结果表明,砂堤后缘存在2种典型形态:一种是砂堤后缘具有明显沉降剖面角;一种是砂堤后缘高度与平衡高度接近,没有明显沉降剖面角。分析认为,影响砂堤后缘差异的因素是:1)平衡高度与后缘高度是否具有明显差异;2)流体从砂堤上方过流通道内流出时,是否具有较强的携带支撑剂向远端运移的能力。图13所示结果表明,压裂液黏度对砂堤后缘形态影响最大,其次是泵注排量、支撑剂类型及注入位置。图14所示结果表明,形态1模式下支撑剂砂堤后缘的堆积主要受到涡流的控制,而形态2模式下支撑剂直接沉降并堆积。

    砂堤形貌演化过程表明,缝内砂堤生长经历3个阶段:1)前端砂堤形成阶段,泵注初期,大量支撑剂在注入口附近沉降形成初始砂堤;2)中部砂堤平衡高度建立阶段,缝内砂堤高度逐渐增长,达到平衡高度后,支撑剂通过裂缝顶部与砂堤表面之间的过流通道向裂缝深处运移;3)后缘砂堤延伸阶段,从过流通道内流出的支撑剂在裂缝后端运移沉降,砂堤整体向后延伸。

    从整体铺置过程来看,砂堤生长(支撑剂铺置)呈现2种模式:模式1,为“裂缝前端先堆积至平衡高度,再稳定向后端铺置”模式(见图15(a);图中T1T4表示支撑剂在铺置过程中的4个阶段,下同);模式2,为“砂堤整体纵向增长,稳定向后端铺置”模式(见图15(b))。其中,模式1在低泵注排量、低黏度压裂液、大粒径支撑剂时发生(试验组1和4),模式2在高泵注排量、中高黏度压裂液、小粒径支撑剂时发生(试验组2,3,5,6,7,8和9)。但这2种模式并非完全相互对立,而是可以在铺置过程中某个阶段发生转换。

    图  15  砂堤生长模式
    Figure  15.  Growth mode of sand embankment

    砂堤达到平衡高度后,过流通道内的流动速度大幅度提高,有利于支撑剂远距离输送。现阶段页岩储层压裂大规模采用“先加入70/140目石英砂打磨孔眼与铺砂+后加入40/70目陶粒尾注”的加砂模式,这会导致具有高导流能力的陶粒被携带至裂缝远端,无法在近井地带形成有效堆积。建议采取“先沉降架桥+后长距离运移+尾近井高导流”的思路,前期加入70/140目石英砂,打磨孔眼与适当加砂;中期加入40/70目陶粒,使砂堤快速达到平衡高度;后期提高压裂液黏度,加入70/140目石英砂,实现裂缝内长距离铺置;施工结束前降低泵注排量,并尾追40/70目陶粒,以获得高导流近井区域。

    1)基于可视化平板裂缝的多尺度支撑剂输送物理模拟系统,采用支撑剂宏观铺置与颗粒微观运动相结合的方式开展支撑剂运移试验,从定性描述与定量解释的角度分析砂堤生长过程与最终形态差异,为支撑剂输送研究提供了新的思路。

    2)缝内砂堤生长经历了前端砂堤形成阶段、中部砂堤平衡高度建立阶段和后缘砂堤延伸阶段等3个阶段。缝内砂堤生长存在2种模式:“裂缝前端先堆积至平衡高度,再稳定向后端铺置”模式,在低泵注排量、低黏度压裂液、大粒径支撑剂时发生;“砂堤整体纵向增长,稳定向后端铺置”模式,在高泵注排量、中高黏度压裂液、小粒径支撑剂时发生。2种模式在泵注的不同阶段出现并转换,可以从砂堤后缘形态上进行区分。

    3)砂堤最终形态与多个因素有关,注入位置与泵注排量对砂堤前缘形态影响最大,压裂液黏度对砂堤中部与砂堤后缘形态影响最大。

    4)现场施工时可考虑采用定向射孔模式,泵注前期采用高排量、中高黏度压裂液泵注主体支撑剂(70/140目石英砂),中期切换为高导流支撑剂(40/70目陶粒)架桥,后期提高压裂液黏度继续泵送70/140目石英砂,实现裂缝内长距离输送,最后降低排量尾追40/70目陶粒,以保证近井地带裂缝与井筒的高连通性。

    5)本文只考虑了支撑剂颗粒运动,并未对流体运动进行深入分析,后续研究建议考虑流体运动特征,进一步分析流体与颗粒之间的作用关系,从而更好地评价缝内携砂运动过程。

  • [1] 路保平.中国石化页岩气工程技术进步及展望[J].石油钻探技术,2013,41(5):1-8. LU Baoping.Sinopec engineering technical advance and its developing tendency in shale gas[J].Petroleum Drilling Techniques,2013,41(5):1-8.
    [2]

    PADILLA A.Social responsibility management systems:elevating performance for shale gas development[R].SPE 156728,2012.

    [3] 薛承瑾.国内页岩气有效开采值得关注的几个问题[J].石油钻探技术,2012,40(4):1-6. XUE Chengjin.Noteworthy issues on effective production of shale gas resource in China[J].Petroleum Drilling Techniques,2012,40(4):1-6.
    [4] 李亚洲,李勇明,罗攀,等.页岩气渗流机理与产能研究[J].断块油气田,2013,20(2):186-190. LI Yazhou,LI Yongming,LUO Pan,et al.Study on seepage mechanism and productivity of shale gas[J].Fault-Block Oil Gas Field,2013,20(2):186-190.
    [5]

    OZKAN E,RAGHAVAN R S,APAYDIN O G.Modeling of fluid transfer from shale matrix to fracture network[R].SPE 134830,2010.

    [6]

    WATSON A T,GATENS J M III S A,LEE W J,et al.An analytical model for history matching naturally fractured reservoir production data[R].SPE 18856,1990.

    [7]

    WARREN J E,ROOT P J.The behavior of naturally fracture reservoirs[J].SPE Journal,1963,3(3):245-255.

    [8] 肖晓春,潘一山.滑脱效应影响的低渗煤层气运移实验研究[J].岩土工程学报,2009,31(10):1554-1558. XIAO Xiaochun,PAN Yishan.Experimental study on gas slippage effects in hypotonic coal reservoir[J].Chinese Journal of Geotechnical Engineering,2009,31(10):1554-1558.
    [9] 刘萱,许浩,汤达祯,等.低渗煤层气藏中的滑脱效应及影响因素[J].科技创新导报,2014,11(20):32-33. LIU Xuan,XU Hao,TANG Dazhen,et al.Factors on gas slippage effects in low permeability coal reservoir[J].Science and Technology Innovation Herald,2014,11(20):32-33.
    [10]

    SHABRO V,TORRES-VERDIN C,JAVADPOUR F.Numerical simulation of shale-gas production:from pore-scale modeling of slip-flow,Knudsen diffusion,and Langmuir desorption to reservoir modeling of compressible fluid[R].SPE 144355,2001.

    [11]

    SWAMI V.Shale gas reservoir modeling:from nanopores to laboratory[R].SPE 163065,2012.

    [12]

    AZOM P N,JAVADPOUR F.Dual-continuum modeling of shale and tight gas reservoirs[R].SPE 159584,2012.

    [13]

    JAVADPOUR F,FISHER D,UNSWORTH M.Nanoscale gas flow in shale gas sediments[J].Journal of Canadian Petroleum Technology,2007,46(10):55-61.

    [14]

    JAVADPOUR F.Nanopores and apparent permeability of gas flow in mudrocks(shales and siltstone)[J].Journal of Canadian Petroleum Technology,2009,48(8):16-21.

    [15] 夏阳,金衍,陈勉,等.页岩气渗流数学模型[J].科学通报,2015,60(24):2259-2271. XIA Yang,JIN Yan,CHEN Mian,et al.Gas flow in shale reservoirs[J].Chin.Sci.Bull.,2015,60(24):2259-2271.
    [16] 邓佳,朱维耀,刘锦霞,等.考虑应力敏感性的页岩气产能预测模型[J].天然气地球科学,2013,24(3):456-460. DENG Jia,ZHU Weiyao,LIU Jinxia,et al.Productivity prediction model of shale gas considering stress sensitivity[J].Natural Gas Geoscience,2013,24(3):456-460.
    [17] 田冷,肖聪,顾岱鸿.考虑应力敏感与非达西效应的页岩气产能模型[J].天然气工业,2014,34(12):70-75. TIAN Leng,XIAO Cong,GU Daihong.A shale gas reservoir productivity model considering stress sensitivity and non-Darcy flow[J].Natural Gas Industry,2014,34(12):70-75.
    [18] 郭肖,任影,吴红琴.考虑应力敏感和吸附的页岩表观渗透率模型[J].岩性油气藏,2015,27(4):109-112. GUO Xiao,REN Ying,WU Hongqin.Apparent permeability model of shale gas considering stress sensitivity and adsorption[J].Lithologic Reservoirs,2015,27(4):109-112.
    [19] 郭小哲,周长沙.基于滑脱的页岩气藏压裂水平井渗透模型及产能预测[J].石油钻采工艺,2015,37(3):61-65. Guo Xiaozhe, Zhou Changsha, Seepage and productivity forecast based on slippage of fratured horizontal well in shale gas pool[J].Oil Drilling Production Technology, 2015, 37(3):61-65.
    [20] 孙海,姚军,孙致学,等.页岩气数值模拟技术进展及展望[J].油气地质与采收率,2012,19(1):46-49. SUN Hai,YAO Jun,SUN Zhixue.Recent development and prospect on numerical simulation of shale gas reservoirs[J].Petroleum Geology and Recovery Efficiency,2012,19(1):46-49.
    [21] 李玉梅,吕炜,宋杰,等.层理性页岩气储层复杂网络裂缝数值模拟研究[J].石油钻探技术,2016,44(4):108-113. LI Yumei,LYU Wei,SONG Jie,et al.Numerical simulation study on the complex network fractures of stratified shale gas reservoirs[J].Petroleum Drilling Techniques,2016,44(4):108-113.
    [22]

    CARLSON E S,MERCER J C.Devonian shale gas production:mechanisms and simple models[J].Journal of Petroleum Technology,1991,43(4):476-482.

    [23]

    MORIDIS G J,BLASINGAME T A,FREEMAN C M.Analysis of mechanisms of flow in fractured tight-gas and shale-gas reservoirs[R].SPE 139250,2010.

    [24]

    KLINKENBERG L J.The permeability of porous media to liquids and gases[J].Socar Proceedings,1941,2(2):200-213.

    [25] 张睿,宁正福,杨峰,等.微观孔隙结构对页岩应力敏感影响的实验研究[J].天然气地球科学,2014,25(8):1284-1289. ZHANG Rui,NING Zhengfu,YANG Feng,et al.Experimental study on microscopic pore structure controls on shale permeability under compaction process[J].Natural Gas Geoscience,2014,25(8):1284-1289.
    [26]

    JU Binshan,WU Yushu,FAN Tailiang.Study on fluid flow in nonlinear elastic porous media:experimental and modeling approaches[J].Journal of Petroleum Science Engineering,2011,76(3/4):205-211.

    [27] 任飞,王新海,谢玉银,等.考虑滑脱效应的页岩气井底压力特征[J].石油天然气学报,2013,35(3):124-126. REN Fei,WANG Xinhai,XIE Yuyin,et al.Characteristics of bottom pressure of shale-gas well considering slippage effect[J].Journal of Oil and Gas Technology,2013,35(3):124-126.
    [28] 韩大匡,陈钦雷,闫存章.油藏数值模拟基础[M].北京:石油工业出版社,1993. HAN Dakuang,CHEN Qinlei,YAN Cunzhang.Fundmentals of oil reservoir numerical simulation[M].Beijing:Petroleum Industry Press,1993.
    [29]

    JU Binshan,FAN Tailiang,WANG Xiaodong,et al.A new simulation framework for predicting the onset and effects of fines mobilization[J].Transport in Porous Media,2007,68(2):265-283.

    [30] 葛洪魁,申颍浩,宋岩,等.页岩纳米孔隙气体流动的滑脱效应[J].天然气工业,2014,34(7):46-54. GE Hongkui,SHEN Yinghao,SONG Yan,et al.Slippage effect of shale gas flow in nanoscale pores[J].Natural Gas Industry,2014,34(7):46-53.
    [31] 白玉湖,杨皓,陈桂华,等.页岩气产量递减典型曲线的不确定性分析方法[J].石油钻探技术,2013,41(4):97-100. BAI Yuhu,YANG Hao,CHEN Guihua,et al.An uncertainty analysis method on typical production decline curve for shale gas reservoirs[J].Petroleum Drilling Techniques,2013,41(4):97-100.
    [32] 张荻萩,李治平,苏皓.页岩气产量递减规律研究[J].岩性油气藏,2015,27(6):138-144. ZHANG Diqiu,LI Zhiping,SU Hao.Production decline trend of shale gas[J].Lithologic Reservoirs,2015,27(6):138-144.
  • 期刊类型引用(3)

    1. 刘拯君,张衍君,廖婉蓉,杨兵,杨柳,周德胜. 页岩储层盐离子扩散规律及其在压裂缝网诊断中的应用. 新疆石油天然气. 2025(01): 41-49 . 百度学术
    2. 端祥刚,胡志明,常进,石雨昕,吴振凯,许莹莹. 页岩储层无支撑缝网区流动能力影响因素研究与进展. 特种油气藏. 2025(01): 22-31 . 百度学术
    3. 雷刚,任建飞,刘佰衢,熊健,刘向君,程万,梁利喜,黎子阳. 超临界CO_2—水—岩相互作用下龙马溪组页岩物理特性演化及其对井壁稳定性影响. 天然气工业. 2025(05): 136-149 . 百度学术

    其他类型引用(0)

计量
  • 文章访问数:  19913
  • HTML全文浏览量:  107
  • PDF下载量:  14891
  • 被引次数: 3
出版历程
  • 收稿日期:  2016-10-09
  • 修回日期:  2016-12-26
  • 刊出日期:  1899-12-31

目录

/

返回文章
返回