GEOFLUIDSgydF4y2Ba GeofluidsgydF4y2Ba 1468 - 8123gydF4y2Ba 1468 - 8115gydF4y2Ba HindawigydF4y2Ba 10.1155 / 2020/8886988gydF4y2Ba 8886988gydF4y2Ba 研究文章gydF4y2Ba 全面模型估计刺激储集层体积返排基于页岩气储层中的数据gydF4y2Ba 程ydF4y2Ba 气gydF4y2Ba 1gydF4y2Ba 王gydF4y2Ba ShaojungydF4y2Ba 2gydF4y2Ba 朱gydF4y2Ba 丹gydF4y2Ba 3gydF4y2Ba 任gydF4y2Ba GuoxuangydF4y2Ba 4gydF4y2Ba 张gydF4y2Ba 元gydF4y2Ba 1gydF4y2Ba https://orcid.org/0000 - 0002 - 4102 - 1954gydF4y2Ba 胡gydF4y2Ba 龚京红gydF4y2Ba 1gydF4y2Ba 徐gydF4y2Ba 金泽国际gydF4y2Ba 1gydF4y2Ba 北京重点实验室的非常规天然气地质评价和开发工程gydF4y2Ba 中国地质大学gydF4y2Ba 北京100083年gydF4y2Ba 中国gydF4y2Ba cug.edu.cngydF4y2Ba 2gydF4y2Ba 石油勘探与开发研究所gydF4y2Ba 北京100083年gydF4y2Ba 中国gydF4y2Ba cnpc.com.cngydF4y2Ba 3gydF4y2Ba 第五届石油生产工厂gydF4y2Ba 长庆油田gydF4y2Ba 西安710021年gydF4y2Ba 中国gydF4y2Ba 4gydF4y2Ba Leewen-Cobra国际能源(北京)技术gydF4y2Ba 有限公司gydF4y2Ba 有限公司gydF4y2Ba 100084年gydF4y2Ba 中国gydF4y2Ba 2020年gydF4y2Ba 9gydF4y2Ba 9gydF4y2Ba 2020年gydF4y2Ba 2020年gydF4y2Ba 12gydF4y2Ba 5gydF4y2Ba 2020年gydF4y2Ba 24gydF4y2Ba 7gydF4y2Ba 2020年gydF4y2Ba 28gydF4y2Ba 8gydF4y2Ba 2020年gydF4y2Ba 9gydF4y2Ba 9gydF4y2Ba 2020年gydF4y2Ba 2020年gydF4y2Ba 版权©2020陈Qi et al。gydF4y2Ba 这是一个开放的文章在知识共享归属许可下发布的,它允许无限制的使用,分布和繁殖在任何媒介,提供最初的工作是正确的引用。gydF4y2Ba

刺激水库卷(SRV)与多级水力压裂产生的水平钻孔控制页岩气储层的生产。尽管微震的数据被用来估计SRV,高价,有时会高估了。此外,没有考虑应力敏感性对电脑的影响异常超压的区域。因此,这项工作的目的是描述地下裂缝网络的渗透率应力敏感性通过页岩气井生产数据返排的早期阶段。排液的返排地区首次发现两个中国南方页岩气井的数据。然后,我们测量了渗透率应力敏感性的核心压裂后,耦合动态相对渗透率(DRP)计算获得一个准确和简单复制曲线。之后,一个全面的模型是建立在考虑动态两相相对渗透率和压力灵敏度函数。最后,我们比较了计算结果与微震的数据。结果表明,该模型可以合理地预测SRV使用返排压裂后的数据。此外,与微震的数据相比,应力敏感性应该包括,特别是异常超压块。 It is believed that this mathematical model is accurate and useful. The work provides an efficient approach to estimate stimulated reservoir volume in the shale gas reservoirs.

 中国国家自然科学基金gydF4y2Ba 51804282gydF4y2Ba 中央大学基础研究基金gydF4y2Ba 2652018209gydF4y2Ba 中国国家科技重大项目gydF4y2Ba 2017 zx05009 - 005gydF4y2Ba 1。介绍gydF4y2Ba

页岩气的发展得到了越来越多的关注与生产从传统的水库的衰落gydF4y2Ba 1gydF4y2Ba,gydF4y2Ba 2gydF4y2Ba]。由于超低渗透率页岩(从gydF4y2Ba 10gydF4y2Ba −gydF4y2Ba 23gydF4y2Ba 来gydF4y2Ba 10gydF4y2Ba −gydF4y2Ba 17gydF4y2Ba 米gydF4y2Ba 2gydF4y2Ba )[gydF4y2Ba 3gydF4y2Ba),水平钻孔和多级水力压裂已广泛用于深水救生艇,从而增加有效接触显著提高天然气产量从个体。因此,刺激的表征储层体积似乎是至关重要的管理和预测页岩储层性能(gydF4y2Ba 4gydF4y2Ba,gydF4y2Ba 5gydF4y2Ba]。gydF4y2Ba

而微震的数据被广泛用于验证储层体积,刺激SRV获得直接从微震的数据通常是高估了。例如,电脑,包括大型凸性假设(如果空白区域的gydF4y2Ba 6gydF4y2Ba- - - - - -gydF4y2Ba 9gydF4y2Ba]。的同时,实现微震的技术领域有时并不划算,这限制了领域广泛应用。因此,早期使用生产数据经过水力压裂返排一起早期数据的液压机液体正迅速成为一个关键仪器在评估SRV [gydF4y2Ba 10gydF4y2Ba- - - - - -gydF4y2Ba 13gydF4y2Ba返排),即瞬态分析的数据。gydF4y2Ba

两相流似乎是一个典型的、有代表性的流动行为后返排过程中水力压裂(gydF4y2Ba 2gydF4y2Ba,gydF4y2Ba 14gydF4y2Ba]。例如,大量的生产数据来自井河页岩角的观察表明,两相流返排过程中对每一个水力压裂后(gydF4y2Ba 15gydF4y2Ba,gydF4y2Ba 16gydF4y2Ba]。克拉克森(gydF4y2Ba 17gydF4y2Ba)和Williams-Kovacs克拉克森(gydF4y2Ba 18gydF4y2Ba]还发现在48小时内两相流的开放和生产在巴内特和马塞勒斯页岩。同样,两相流阶段发生在志留纪Longmaxi形成页岩气水井开井后在中国南方[gydF4y2Ba 19gydF4y2Ba]。和一个更大的推进生产数据研究的模式是,在处理生产数据从页岩气井的喇叭河流域,Adefidipe et al。gydF4y2Ba 20.gydF4y2Ba,gydF4y2Ba 21gydF4y2Ba和徐et al。gydF4y2Ba 22gydF4y2Ba,gydF4y2Ba 23gydF4y2Ba)生产数据划分为两个阶段,早期天然气生产出路)和晚期天然气生产(LGP)的“V”形气水比之间的关系(吉尼斯世界纪录)和累计天然气生产。他们认为出路阶段是特别的,值得研究,因为LGP阶段已经有了一个矩阵的天然气供应。gydF4y2Ba

利用两相流的特征在天然气生产出路)早期阶段,研究研究努力估计SRV耦合与各种类型的储层岩石物理模型。例如,Ezulike et al。gydF4y2Ba 24gydF4y2Ba- - - - - -gydF4y2Ba 26gydF4y2Ba集成一个明确确定动态相对渗透率(DRP)函数在线性双孔隙度模型来研究两相返排的数据。功能扩展了单相水流动方程气水两相流计算使用累积水和碳氢化合物的生产数据。克拉克森et al。gydF4y2Ba 17gydF4y2Ba,gydF4y2Ba 27gydF4y2Ba,gydF4y2Ba 28gydF4y2Ba]最初使用一种改进的方法对煤层气(CBM)分析多相流流动。他们认为,气体和水的流动通过生成的裂缝网络类似于同步流的气体和水从破碎的煤储层在长期生产。后,他们提出了大量的基于这种假设的两相流模型定量计算多流的数据,包括致密油、致密气,帮助分析页岩气井的流动特性,类似于分析两相(气、水)煤层气生产。在后续的研究工作(gydF4y2Ba 29日gydF4y2Ba- - - - - -gydF4y2Ba 32gydF4y2Ba),他们从上面的方法验证了分析结果通过整合各种技术,包括改性物质平衡方程(MBE)应用于before-breakthrough(《)(气体突破之前)单相rate-transient分析(等)和after-breakthrough (ABT)(后气体突破)多相rate-transient分析(等)。gydF4y2Ba

虽然上述模型在定量计算SRV通过多流的数据提供见解,将这些模型应用到领域的挑战仍然因为朗缪尔体积计算过程中需要增加,这将导致一定程度上的不确定性和不准确性。为了克服这些模型应用领域的复杂性,克拉克森et al。gydF4y2Ba 17gydF4y2Ba,gydF4y2Ba 28gydF4y2Ba和傅et al。gydF4y2Ba 33gydF4y2Ba)建立了一个线性匹配曲线关系来描述之间的直线关系rate-normalized压力(RNP)和物质平衡时间(MBT)。这是因为裂缝渗透率和裂缝体积的基础上可以得到拟合曲线的特征。此外,Alkouh et al。gydF4y2Ba 34gydF4y2Ba)修改的重对数坐标图RNP和RNP导数MBT对石油和天然气的情况下和提供了一个简单的方法来计算主要包括裂缝渗透率和裂缝半身的断裂参数。然而,在返排过程中压缩性骨折的影响没有考虑模型,这可能导致总压缩系数的估计误差。gydF4y2Ba

然而,由于一些缺陷的假设组成函数,非线性现象发生在模型应用到一些领域案件。和非零截距的现象导致研究人员质疑这些问题的物理性质。这可能是由于错误引起的应力敏感性在某些领域。特别是,有异常超压在中国南方页岩气储层gydF4y2Ba 35gydF4y2Ba),因此应力敏感性的影响,尤其是对两相相对渗透率,必须考虑(gydF4y2Ba 36gydF4y2Ba]。此外,所有研究人员并没有使生产数据与微震的验证和比较数据。gydF4y2Ba

我们因此旨在开发一个stress-dependent两相相对渗透率,可用于扩展现有的单相模型。同时,我们提出了一种新的多相模型来估计SRV和断裂特征。此外,来验证我们的模型,我们比较结果与微震的数据。gydF4y2Ba

在这项工作中,我们首先处理页岩气井的生产数据和获得气水比的特点。然后,我们计算了简单实用的两阶段动态相对渗透率(DRP)基于生产数据收集在这个领域属于页岩气井在四川南部。随后,渗透率的应力-应变曲线是通过实验获得。最后,物质平衡方程和扩散方程是用来处理生产数据推导出模型,从而计算裂缝参数。我们比较了计算结果与微震的数据来验证计算结果的合理性和准确性。gydF4y2Ba

 2。数据准备gydF4y2Ba 2.1。字段数据从一个典型的中国南方页岩气藏gydF4y2Ba

在这项工作中,我们选择了返排数据从两个生产井位于中国南方:Silurianare越低,页岩气形成的一个典型的海洋页岩气藏。开发大量micron-nanoscale毛孔在矿物颗粒和有机物(gydF4y2Ba 37gydF4y2Ba]。这些毛孔页岩气储存和运移的主要通道,它发挥着重要的作用在提高页岩气的存储性能和提供良好的存储空间海洋页岩气积累和浓缩gydF4y2Ba 38gydF4y2Ba]。页岩气储层的深度超过2000米在地下地层压力系数2.0。整个页岩储层的厚度约为100 m和孔隙率从2%到4%不等的矩阵渗透率的0.001。在早期阶段的生产gydF4y2Ba 10gydF4y2Ba ~gydF4y2Ba 20.gydF4y2Ba ×gydF4y2Ba 10gydF4y2Ba 4gydF4y2Ba 米gydF4y2Ba 3gydF4y2Ba /gydF4y2Ba dgydF4y2Ba 与在一个稳定的生产gydF4y2Ba 5gydF4y2Ba ~gydF4y2Ba 8gydF4y2Ba ×gydF4y2Ba 10gydF4y2Ba 4gydF4y2Ba 米gydF4y2Ba 3gydF4y2Ba /gydF4y2Ba dgydF4y2Ba 。gydF4y2Ba

一些错误的东西背斜构造带是安排在左边梯形行是这个地区的主要地质构造(gydF4y2Ba 39gydF4y2Ba]。这个钻井平台位于背斜构造。尽管存在两个压剪断裂在垫,整体上外侧密封性能导致页岩气的有效积累和保存(gydF4y2Ba 37gydF4y2Ba,gydF4y2Ba 38gydF4y2Ba]。图gydF4y2Ba 1gydF4y2Ba显示了水平井的布局。由于较低的石英含量优质页岩段底部的Longmaxi形成,自然骨折不发达整个水库gydF4y2Ba 40gydF4y2Ba]。整个水平井组1,4,5位于一形成,而2、3、6和7属于另一个形成。gydF4y2Ba

垫示意图。四川盆地南部七MFHWs钻。gydF4y2Ba

在这项工作中,我们使用3和4的数据来测试我们的模型对以下两个原因。首先,两井钻探页岩气储层在不同的方向和位于两个不同的阵型。其次,微地震监测是在两井一起进行综合钻井、完井数据和测试数据。图gydF4y2Ba 2gydF4y2Ba显示了两口井的天然气和水的生产速度。两个经验丰富的天然气井突破后立即打开和生产,这意味着没有单相流。不同之处在于,3生产一度经历了急剧下降(80小时),而4经历了一个更温和的波动,可能是由于不同的形成条件。很抱歉,没有进一步研究这一现象的原因本文的高频生产数据不足或完全相邻的测试数据。gydF4y2Ba

诊断情节两井页岩气地层属于四川南部。(一)产量图3。(b)生产速度图4。gydF4y2Ba

 2.2。生产气-水比例(吉尼斯世界纪录)gydF4y2Ba

为了探索吉尼斯世界纪录的生产数据的趋势,这一节简单的流程和分析生产数据后立即观察是否有气体突破页岩气井的生产(gydF4y2Ba 17gydF4y2Ba,gydF4y2Ba 27gydF4y2Ba),这意味着如果可以通过建立相关模型计算基于生产数据。gydF4y2Ba

这主要是因为有效裂缝网络系统是用天然气和水饱和阶段后两个自闭。这里的天然气来源被认为是来自三个方面:(1)最初现有初始气体活性天然裂缝,(2)气体压裂液的影响下而流离失所的强逆流水自吸到页岩矩阵由于巨大的压力差第一关井期间,和(3)气体积累了从自发的自吸的压裂液在第二关井期间(gydF4y2Ba 41gydF4y2Ba,gydF4y2Ba 42gydF4y2Ba]。gydF4y2Ba

同样,在建立模型之前,我们还需要使用气水产量比例增益v型趋势。因此,我们最初的两阶段生产数据处理3和4,最后获得的气水比的趋势(吉尼斯世界纪录)。图gydF4y2Ba 3gydF4y2Ba显示了一个v型趋势吉尼斯世界纪录随着时间的推移早期两相返排的特性数据。虽然通过曲线的拐点的两个井发生在不同的时间不同的形成条件和生产系统,同样的v型趋势仍然很明显。Ghanbari et al。gydF4y2Ba 15gydF4y2Ba]和Abbasi [gydF4y2Ba 16gydF4y2Ba]还发现v型吉尼斯世界纪录在喇叭河页岩井钻在加拿大。gydF4y2Ba

诊断情节两井页岩气地层属于四川南部。(一)吉尼斯世界纪录3情节。4 (b)吉尼斯世界纪录的阴谋。gydF4y2Ba

气-水比例增加然后减少,我们可以从最基本的理论研究。假定气水两相流满足达西定律:gydF4y2Ba (1)gydF4y2Ba 问gydF4y2Ba =gydF4y2Ba −gydF4y2Ba kgydF4y2Ba 一个gydF4y2Ba μgydF4y2Ba dgydF4y2Ba PgydF4y2Ba dgydF4y2Ba lgydF4y2Ba ,gydF4y2Ba 在哪里gydF4y2Ba 问gydF4y2Ba 是在压力下流体流量吗gydF4y2Ba PgydF4y2Ba ,米gydF4y2Ba3gydF4y2Ba/ s;gydF4y2Ba kgydF4y2Ba 液体渗透,gydF4y2Ba μgydF4y2Ba米gydF4y2Ba2gydF4y2Ba;gydF4y2Ba 一个gydF4y2Ba 流体流动区域,米gydF4y2Ba2gydF4y2Ba;gydF4y2Ba μgydF4y2Ba 流体的粘滞性,mP·s;gydF4y2Ba PgydF4y2Ba MPa的压力;和gydF4y2Ba lgydF4y2Ba 是流体流动长度,m。gydF4y2Ba

忽略毛细力后裂缝网络系统,气水比gydF4y2Ba (2)gydF4y2Ba 问gydF4y2Ba ggydF4y2Ba 问gydF4y2Ba wgydF4y2Ba =gydF4y2Ba μgydF4y2Ba wgydF4y2Ba μgydF4y2Ba ggydF4y2Ba PgydF4y2Ba kgydF4y2Ba rggydF4y2Ba 年代gydF4y2Ba ggydF4y2Ba kgydF4y2Ba rwgydF4y2Ba 年代gydF4y2Ba wgydF4y2Ba ∂gydF4y2Ba PgydF4y2Ba ggydF4y2Ba ∂gydF4y2Ba PgydF4y2Ba wgydF4y2Ba ≈gydF4y2Ba μgydF4y2Ba wgydF4y2Ba μgydF4y2Ba ggydF4y2Ba PgydF4y2Ba kgydF4y2Ba rggydF4y2Ba 年代gydF4y2Ba ggydF4y2Ba kgydF4y2Ba rwgydF4y2Ba 年代gydF4y2Ba wgydF4y2Ba ,gydF4y2Ba 在哪里gydF4y2Ba 问gydF4y2Ba ggydF4y2Ba 和gydF4y2Ba 问gydF4y2Ba wgydF4y2Ba 气相的流速和水阶段,分别mgydF4y2Ba3gydF4y2Ba/ d;gydF4y2Ba μgydF4y2Ba wgydF4y2Ba 和gydF4y2Ba μgydF4y2Ba ggydF4y2Ba PgydF4y2Ba 是气相和水相的粘度下吗gydF4y2Ba PgydF4y2Ba 压力,分别mP·s;gydF4y2Ba kgydF4y2Ba rggydF4y2Ba 年代gydF4y2Ba ggydF4y2Ba 和gydF4y2Ba kgydF4y2Ba rwgydF4y2Ba 年代gydF4y2Ba wgydF4y2Ba 是气相和水相的相对渗透率下各自的饱和,gydF4y2Ba μgydF4y2Ba米gydF4y2Ba2gydF4y2Ba;gydF4y2Ba 年代gydF4y2Ba ggydF4y2Ba 和gydF4y2Ba 年代gydF4y2Ba wgydF4y2Ba 气相的饱和度和水阶段,无量纲;和gydF4y2Ba PgydF4y2Ba ggydF4y2Ba 和gydF4y2Ba PgydF4y2Ba wgydF4y2Ba 是气相的压力和水阶段,分别MPa。gydF4y2Ba

开井生产后,气体粘度随压力的降低,而水粘度保持相对不变。根据气水比的变化随着时间的推移,气水相对渗透率降低的比率。因此,渗透率和饱和度之间的关系的基础上,可以推断,饱和度都有对应的变化趋势,即初始气体饱和度(gydF4y2Ba 年代gydF4y2Ba 胃肠道gydF4y2Ba )不为零,这也证明了最初的存在游离气在断裂之前开放。井筒储存效应消失时,气体流从矩阵裂缝网络,导致气体在断裂系统的补充。因此,气体饱和度增加,气水相对渗透率增加,气水比减少。这就是为什么吉尼斯世界纪录的斜率曲线先下降,然后上升。gydF4y2Ba

一般来说,吉尼斯世界纪录的诊断上的负斜率阴谋被称为早期天然气生产(出路和吉尼斯世界纪录的崛起诊断图被称为晚期天然气生产(LGP)。人们普遍认为,出路地区井筒储存效应的阶段,气体和水制作来自于有效的断裂与水平井。同时,水相对渗透率的增加大于气体粘度的降低。LGP阶段矩阵的结果气体转移到井筒效应后的裂缝网络变得可以忽略不计(gydF4y2Ba 10gydF4y2Ba,gydF4y2Ba 11gydF4y2Ba,gydF4y2Ba 17gydF4y2Ba,gydF4y2Ba 25gydF4y2Ba- - - - - -gydF4y2Ba 27gydF4y2Ba]。根据上面的分析,我们发现在出路阶段,地上的返排流只是从有效的裂缝系统,而LGP阶段开始后,矩阵气体提升者参与流动。因此,我们将建立一个分析模型的出路阶段计算SRV页岩气多级压裂后油井。gydF4y2Ba

 3所示。方法gydF4y2Ba 3.1。物质平衡方程(MBE)和扩散方程的断裂系统gydF4y2Ba

为了方便出路的数学模型的建立阶段,页岩气压裂井周围的裂缝网络简化成SRV地区组成的矩阵系统和断裂系统。如图gydF4y2Ba 4gydF4y2Ba,有效裂缝系统是由一个骨折部分和一个矩阵部分。断裂部分由人工液压骨折、压裂所产生的二次骨折,活跃自然骨折。矩阵部分页岩矩阵连接断裂部分。值得注意的是,解吸和吸附过程中忽略了物质平衡方程(MBE) [gydF4y2Ba 17gydF4y2Ba,gydF4y2Ba 43gydF4y2Ba]。gydF4y2Ba

原理图与多级水力压裂水平井的发展的物质平衡方程。虚线箭头显示流体流动方向,按顺序从矩阵骨折和骨折井筒。gydF4y2Ba

这个模型简化了复杂、活跃的自然和二次骨折,以及人工液压骨折,到一个简单的断裂系统。人工水力裂缝的长度作为整个储层体积,刺激的宽度和长度的水平井筒的长度。在整个库区刺激,主要断裂的高度,包括矩阵连接部分断裂系统,参与流动。在等效断裂系统,裂缝和压裂液饱和(水相)和天然气(气相)。假设不涉及矩阵气体流中的早期天然气生产(出路。gydF4y2Ba

假设裂缝系统可以近似为均匀/关闭/罐系统。和流体从骨折到水平井被认为是线性的。推动气水流动机制包括两个方面:(1)裂缝闭合和(2)扩张的流体(气水阶段)。Kuchuk et al。gydF4y2Ba 44gydF4y2Ba)也被认为是假设裂缝闭合和扩张的水和游离气在断裂系统驱动天然气和水在表面的积累表明pseudo-steady-state流入水库骨折的发生。gydF4y2Ba

一般来说,为了便于模型的建立和解决方案,我们做了以下假设:(1)在断裂系统忽略毛细管压力;(2)气体矩阵是可以忽略不计的出路阶段;(3)断裂系统近似为均匀/封闭系统;(4)达西定律适用于流体流动;(5)有效的断裂系统充斥着压裂液(水相)和天然气(气相)。gydF4y2Ba

有效压缩方面使用简化定义气体物质平衡方程(MBEs) [gydF4y2Ba 45gydF4y2Ba- - - - - -gydF4y2Ba 48gydF4y2Ba]。在这部作品中,物质平衡方程(MBE)的气相断裂网络是简化。气体和水的影响(压裂液)扩张和裂缝闭合被归类为“有效压缩”项由gydF4y2Ba (3)gydF4y2Ba CgydF4y2Ba ~gydF4y2Ba tgydF4y2Ba =gydF4y2Ba 1gydF4y2Ba −gydF4y2Ba GgydF4y2Ba pgydF4y2Ba GgydF4y2Ba figydF4y2Ba BgydF4y2Ba ggydF4y2Ba BgydF4y2Ba 胃肠道gydF4y2Ba 年代gydF4y2Ba 胃肠道gydF4y2Ba CgydF4y2Ba ggydF4y2Ba +gydF4y2Ba 1gydF4y2Ba −gydF4y2Ba WgydF4y2Ba pgydF4y2Ba WgydF4y2Ba figydF4y2Ba 年代gydF4y2Ba wigydF4y2Ba CgydF4y2Ba wgydF4y2Ba +gydF4y2Ba 1gydF4y2Ba VgydF4y2Ba figydF4y2Ba ∂gydF4y2Ba VgydF4y2Ba fgydF4y2Ba ∂gydF4y2Ba PgydF4y2Ba fgydF4y2Ba ,gydF4y2Ba 在哪里gydF4y2Ba GgydF4y2Ba pgydF4y2Ba 和gydF4y2Ba GgydF4y2Ba figydF4y2Ba 累计天然气生产和初始体积的气体在骨折,分别mgydF4y2Ba3gydF4y2Ba;gydF4y2Ba BgydF4y2Ba ggydF4y2Ba 和gydF4y2Ba BgydF4y2Ba 胃肠道gydF4y2Ba 是天然气的地层体积系数和天然气的地层体积系数与初始条件,分别mgydF4y2Ba3gydF4y2Ba/ mgydF4y2Ba3gydF4y2Ba;gydF4y2Ba 年代gydF4y2Ba 胃肠道gydF4y2Ba 和gydF4y2Ba 年代gydF4y2Ba wigydF4y2Ba 最初的气相饱和度和初始水相饱和度,分别,无量纲;gydF4y2Ba CgydF4y2Ba ggydF4y2Ba 和gydF4y2Ba CgydF4y2Ba wgydF4y2Ba 分别是压缩气体和水压缩性,爸爸gydF4y2Ba1gydF4y2Ba;gydF4y2Ba WgydF4y2Ba pgydF4y2Ba 和gydF4y2Ba WgydF4y2Ba figydF4y2Ba 累积水生产和初始体积的水在骨折,分别mgydF4y2Ba3gydF4y2Ba;gydF4y2Ba VgydF4y2Ba fgydF4y2Ba 和gydF4y2Ba VgydF4y2Ba figydF4y2Ba 是有效的体积有效的骨折和体积骨折在初始条件,mgydF4y2Ba3gydF4y2Ba;和gydF4y2Ba PgydF4y2Ba fgydF4y2Ba 破裂压力,Pa。gydF4y2Ba

CgydF4y2Ba ~gydF4y2Ba tgydF4y2Ba 类似于总压缩系数术语提出了传统多相试井。这是表示为一个函数的返排可测量的参数,包括累计天然气生产和累积水生产。每个术语的表达了不同的断裂系统的驱动机制:(1)扩张的气相,液相的扩张(2),(3)裂缝关闭。其中,gydF4y2Ba 1gydF4y2Ba /gydF4y2Ba VgydF4y2Ba figydF4y2Ba ∂gydF4y2Ba VgydF4y2Ba fgydF4y2Ba /gydF4y2Ba ∂gydF4y2Ba PgydF4y2Ba fgydF4y2Ba 类似于地层压缩系数在传统常规物质平衡分析(gydF4y2Ba 49gydF4y2Ba]。在处理骨折时,它指的逆断裂刚度取决于所需的弹性和应变能保持液压骨折开放(gydF4y2Ba 50gydF4y2Ba]。gydF4y2Ba

从矩阵系统忽略了气体流速断裂,气相物质平衡方程gydF4y2Ba (4)gydF4y2Ba 0gydF4y2Ba −gydF4y2Ba 问gydF4y2Ba ggydF4y2Ba ρgydF4y2Ba ggydF4y2Ba 0gydF4y2Ba =gydF4y2Ba ∂gydF4y2Ba ∂gydF4y2Ba tgydF4y2Ba VgydF4y2Ba ggydF4y2Ba tgydF4y2Ba ρgydF4y2Ba ggydF4y2Ba RgydF4y2Ba ,gydF4y2Ba 在哪里gydF4y2Ba ρgydF4y2Ba ggydF4y2Ba 0gydF4y2Ba 和gydF4y2Ba ρgydF4y2Ba ggydF4y2Ba RgydF4y2Ba 气体的密度在表面条件和储层条件,分别公斤/米gydF4y2Ba3gydF4y2Ba;和gydF4y2Ba VgydF4y2Ba ggydF4y2Ba tgydF4y2Ba 气体的体积在断裂系统在任何时候,mgydF4y2Ba3gydF4y2Ba。gydF4y2Ba

根据气体体积之间的关系gydF4y2Ba VgydF4y2Ba ggydF4y2Ba ,裂缝体积gydF4y2Ba VgydF4y2Ba fgydF4y2Ba 和水体积gydF4y2Ba VgydF4y2Ba wgydF4y2Ba 断裂系统,获得天然气压缩系数的定义gydF4y2Ba (5)gydF4y2Ba −gydF4y2Ba 问gydF4y2Ba ggydF4y2Ba BgydF4y2Ba ggydF4y2Ba =gydF4y2Ba VgydF4y2Ba ggydF4y2Ba cgydF4y2Ba ggydF4y2Ba ∂gydF4y2Ba PgydF4y2Ba fgydF4y2Ba ∂gydF4y2Ba tgydF4y2Ba −gydF4y2Ba ∂gydF4y2Ba VgydF4y2Ba wgydF4y2Ba ∂gydF4y2Ba tgydF4y2Ba +gydF4y2Ba ∂gydF4y2Ba VgydF4y2Ba fgydF4y2Ba ∂gydF4y2Ba tgydF4y2Ba ,gydF4y2Ba 在哪里gydF4y2Ba VgydF4y2Ba fgydF4y2Ba ,gydF4y2Ba VgydF4y2Ba ggydF4y2Ba ,gydF4y2Ba VgydF4y2Ba wgydF4y2Ba 裂缝体积、气体体积和裂缝系统的水量,分别mgydF4y2Ba3gydF4y2Ba。gydF4y2Ba

用天然气生产gydF4y2Ba GgydF4y2Ba pgydF4y2Ba 和水的生产gydF4y2Ba WgydF4y2Ba pgydF4y2Ba ,然后用链式法则,骨折派生的总量:gydF4y2Ba (6)gydF4y2Ba −gydF4y2Ba 1gydF4y2Ba VgydF4y2Ba figydF4y2Ba 问gydF4y2Ba ggydF4y2Ba BgydF4y2Ba ggydF4y2Ba +gydF4y2Ba 问gydF4y2Ba wgydF4y2Ba BgydF4y2Ba wgydF4y2Ba =gydF4y2Ba GgydF4y2Ba figydF4y2Ba −gydF4y2Ba GgydF4y2Ba pgydF4y2Ba BgydF4y2Ba ggydF4y2Ba cgydF4y2Ba ggydF4y2Ba GgydF4y2Ba figydF4y2Ba BgydF4y2Ba 胃肠道gydF4y2Ba /gydF4y2Ba 年代gydF4y2Ba 胃肠道gydF4y2Ba ∂gydF4y2Ba PgydF4y2Ba fgydF4y2Ba ∂gydF4y2Ba tgydF4y2Ba +gydF4y2Ba WgydF4y2Ba figydF4y2Ba −gydF4y2Ba WgydF4y2Ba pgydF4y2Ba BgydF4y2Ba wgydF4y2Ba cgydF4y2Ba wgydF4y2Ba WgydF4y2Ba figydF4y2Ba BgydF4y2Ba wigydF4y2Ba /gydF4y2Ba 年代gydF4y2Ba wigydF4y2Ba ∂gydF4y2Ba PgydF4y2Ba fgydF4y2Ba ∂gydF4y2Ba tgydF4y2Ba +gydF4y2Ba 1gydF4y2Ba VgydF4y2Ba figydF4y2Ba ∂gydF4y2Ba VgydF4y2Ba fgydF4y2Ba ∂gydF4y2Ba PgydF4y2Ba fgydF4y2Ba ∂gydF4y2Ba PgydF4y2Ba fgydF4y2Ba ∂gydF4y2Ba tgydF4y2Ba ,gydF4y2Ba 在哪里gydF4y2Ba BgydF4y2Ba wgydF4y2Ba 是水在初始条件,地层体积系数mgydF4y2Ba3gydF4y2Ba/ mgydF4y2Ba3gydF4y2Ba。gydF4y2Ba

早期的天然气生产出路),它是假定gydF4y2Ba BgydF4y2Ba wgydF4y2Ba ≈gydF4y2Ba BgydF4y2Ba wigydF4y2Ba ,最后的物质平衡方程,可以通过定义一个压缩性术语和总流体流量通过简化上述方程gydF4y2Ba (7)gydF4y2Ba ∂gydF4y2Ba PgydF4y2Ba fgydF4y2Ba ∂gydF4y2Ba tgydF4y2Ba =gydF4y2Ba −gydF4y2Ba 问gydF4y2Ba tgydF4y2Ba CgydF4y2Ba ~gydF4y2Ba tgydF4y2Ba VgydF4y2Ba figydF4y2Ba ,gydF4y2Ba 在哪里gydF4y2Ba 问gydF4y2Ba tgydF4y2Ba =gydF4y2Ba 问gydF4y2Ba ggydF4y2Ba BgydF4y2Ba ggydF4y2Ba +gydF4y2Ba 问gydF4y2Ba wgydF4y2Ba BgydF4y2Ba 总生产速度,mgydF4y2Ba3gydF4y2Ba/ d。gydF4y2Ba

鉴于单相稳态流可以描述使用连续性方程和达西定律,单相气体扩散方程在断裂系统是由张先生和冬季(gydF4y2Ba 51gydF4y2Ba]:gydF4y2Ba (8)gydF4y2Ba ∇gydF4y2Ba PgydF4y2Ba fgydF4y2Ba μgydF4y2Ba ggydF4y2Ba ZgydF4y2Ba ∇gydF4y2Ba PgydF4y2Ba fgydF4y2Ba =gydF4y2Ba ϕgydF4y2Ba fgydF4y2Ba KgydF4y2Ba fgydF4y2Ba PgydF4y2Ba fgydF4y2Ba ZgydF4y2Ba CgydF4y2Ba ~gydF4y2Ba tgydF4y2Ba ∂gydF4y2Ba PgydF4y2Ba fgydF4y2Ba ∂gydF4y2Ba tgydF4y2Ba ,gydF4y2Ba 在哪里gydF4y2Ba ZgydF4y2Ba 气体压缩系数,无因次;gydF4y2Ba φgydF4y2Ba fgydF4y2Ba 断裂系统的孔隙度,无量纲;和gydF4y2Ba ∇gydF4y2Ba 梯度算子。gydF4y2Ba

定义了拟压力和拟时间函数(gydF4y2Ba 46gydF4y2Ba,gydF4y2Ba 52gydF4y2Ba]gydF4y2Ba (9)gydF4y2Ba ψgydF4y2Ba PgydF4y2Ba fgydF4y2Ba =gydF4y2Ba ∫gydF4y2Ba 0gydF4y2Ba PgydF4y2Ba fgydF4y2Ba 2gydF4y2Ba PgydF4y2Ba fgydF4y2Ba μgydF4y2Ba ggydF4y2Ba ZgydF4y2Ba ∂gydF4y2Ba PgydF4y2Ba fgydF4y2Ba ,gydF4y2Ba (10)gydF4y2Ba tgydF4y2Ba 一个gydF4y2Ba =gydF4y2Ba ∫gydF4y2Ba 0gydF4y2Ba tgydF4y2Ba kgydF4y2Ba ggydF4y2Ba tgydF4y2Ba μgydF4y2Ba ggydF4y2Ba CgydF4y2Ba ~gydF4y2Ba tgydF4y2Ba ∂gydF4y2Ba tgydF4y2Ba ,gydF4y2Ba 在哪里gydF4y2Ba kgydF4y2Ba ggydF4y2Ba tgydF4y2Ba 气体的相对渗透率是在给定的时间,gydF4y2Ba μgydF4y2Ba米gydF4y2Ba2gydF4y2Ba;和gydF4y2Ba CgydF4y2Ba ~gydF4y2Ba tgydF4y2Ba 总有效压缩,爸爸gydF4y2Ba1gydF4y2Ba。gydF4y2Ba

然后,单相气体流动的控制方程在断裂系统gydF4y2Ba (11)gydF4y2Ba ∂gydF4y2Ba 2gydF4y2Ba ψgydF4y2Ba PgydF4y2Ba fgydF4y2Ba ∂gydF4y2Ba ygydF4y2Ba 2gydF4y2Ba =gydF4y2Ba ϕgydF4y2Ba fgydF4y2Ba KgydF4y2Ba fgydF4y2Ba ∂gydF4y2Ba ψgydF4y2Ba ∂gydF4y2Ba tgydF4y2Ba 一个gydF4y2Ba 。gydF4y2Ba

3.2。耦合与动态相对渗透率应力敏感性函数gydF4y2Ba

在方程(gydF4y2Ba 10gydF4y2Ba),gydF4y2Ba kgydF4y2Ba ggydF4y2Ba tgydF4y2Ba 随着时间的推移单相气体的变量函数。通过引入动态相对渗透率(DRP)函数gydF4y2Ba kgydF4y2Ba rggydF4y2Ba tgydF4y2Ba 单相气体的扩散方程转化为两阶段系统的扩散方程(gydF4y2Ba 53gydF4y2Ba]。gydF4y2Ba

本文的方法是类似于Ezulike和Dehghanpour的研究(gydF4y2Ba 53gydF4y2Ba),但由于地质条件形成属于四川南部不同于北美,动态气相的相对渗透率(DRP)函数获得通过使用合适的字段数据的方法。在先前的研究中,我们假定断裂系统是用天然气和水饱和。累积气体和水生产和时间之间的关系可以通过现场数据报告。因此,气水饱和的断裂系统可以与时间建立了气水产量的函数。鉴于这个区域异常超压的特点,这一地区的应力敏感性不容忽视(gydF4y2Ba 36gydF4y2Ba,gydF4y2Ba 54gydF4y2Ba,gydF4y2Ba 55gydF4y2Ba]。因此,物模型被认为是在动态的相对渗透率函数。具体的过程如图gydF4y2Ba 5gydF4y2Ba。gydF4y2Ba

过程估计动态相对渗透率(DRP)曲线考虑应力敏感性。gydF4y2Ba

在图gydF4y2Ba 6gydF4y2Ba的相对渗透率岩心压裂后通过实验测试研究领域。压力随时间变化的曲线在排液期间图所示gydF4y2Ba 7gydF4y2Ba。应力-应变试验曲线来描述执行无因次渗透率与有效应力之间通过稳定围压和减少内部压力。这是更符合实际的地层返排压力的变化。早期的页岩气井排液阶段后开井明显短;因此,可以认为,地层压力变化不大,但内部压力开始降低流体流出后的核心。gydF4y2Ba

气-水相对渗透率曲线的核心3和4。gydF4y2Ba

套管压力、有效应力改变随着时间的推移两井页岩气地层属于四川南部。3 (a)和(b) 4。gydF4y2Ba

取岩心的应力敏感性实验属于刺激面积3为例。矩阵和骨折最有可能包含在这些核心。我们需要尽可能多的实验测试数据获取物刺激地区压裂后尽可能准确。然而,受制于样本不足或水平井不同距离的样本,这些数据可能不能代表整个地区事实上,但仍然有很大的参考和研究价值。gydF4y2Ba

具体的实验步骤如下:(1)初始围压与原始地层压力是38 MPa,内部压力是23 MPa;(2)围压增加缓慢58 MPa的地层压力在开井排液,和内部压力增加到43 MPa同时保持有效应力常数;,(3)内部压力减少压力点增加有效应力不同,和样品的气体渗透率测量后每一个压力点是稳定的。gydF4y2Ba

无因次渗透率定义为gydF4y2Ba (12)gydF4y2Ba KgydF4y2Ba KgydF4y2Ba 0gydF4y2Ba =gydF4y2Ba 一个gydF4y2Ba egydF4y2Ba −gydF4y2Ba bgydF4y2Ba pgydF4y2Ba 我gydF4y2Ba −gydF4y2Ba pgydF4y2Ba fgydF4y2Ba ,gydF4y2Ba 在哪里gydF4y2Ba KgydF4y2Ba 气体的相对渗透率是在给定的时间,gydF4y2Ba μgydF4y2Ba米gydF4y2Ba2gydF4y2Ba;gydF4y2Ba KgydF4y2Ba 0gydF4y2Ba 是初始气相对渗透率,gydF4y2Ba μgydF4y2Ba米gydF4y2Ba2gydF4y2Ba;和gydF4y2Ba 一个gydF4y2Ba 是gydF4y2Ba bgydF4y2Ba 应力敏感性系数的实验中,无量纲。gydF4y2Ba

实验结果如图gydF4y2Ba 8gydF4y2Ba。根据图的计算方法gydF4y2Ba 5gydF4y2Ba,结果在图gydF4y2Ba 9gydF4y2Ba终于获得了。然后,根据拟合曲线,气相相对渗透率随时间变化的函数gydF4y2Ba (13)gydF4y2Ba kgydF4y2Ba rggydF4y2Ba tgydF4y2Ba =gydF4y2Ba αgydF4y2Ba tgydF4y2Ba −gydF4y2Ba βgydF4y2Ba 。gydF4y2Ba

应力-应变试验结果的3和4。gydF4y2Ba

动态相对渗透率函数随着时间的推移,考虑应力敏感性两井页岩气地层属于四川南部。(一)DRP 3的拟合曲线。4 (b)组成的拟合曲线。gydF4y2Ba

由于上述方程得到基于返排和核心实验相对渗透率数据,gydF4y2Ba αgydF4y2Ba 和gydF4y2Ba βgydF4y2Ba 被称为“返排系数和返排数据拦截,分别,这是由数据拟合曲线。gydF4y2Ba

 3.3。两相返排模型早期天然气生产出路)gydF4y2Ba

用方程(gydF4y2Ba 13gydF4y2Ba)方程(gydF4y2Ba 10gydF4y2Ba),气体单相流的控制方程在断裂系统变成两相流模型。这意味着当气水两相流发生在断裂系统,新实拟时间函数转化为gydF4y2Ba (14)gydF4y2Ba tgydF4y2Ba 一个gydF4y2Ba =gydF4y2Ba ∫gydF4y2Ba 0gydF4y2Ba tgydF4y2Ba kgydF4y2Ba rggydF4y2Ba tgydF4y2Ba μgydF4y2Ba ggydF4y2Ba CgydF4y2Ba ~gydF4y2Ba tgydF4y2Ba ∂gydF4y2Ba tgydF4y2Ba 。gydF4y2Ba

为了建立拟压力函数和拟时间函数之间的关系,我们做以下转变:gydF4y2Ba (15)gydF4y2Ba ∂gydF4y2Ba ψgydF4y2Ba ∂gydF4y2Ba tgydF4y2Ba 一个gydF4y2Ba =gydF4y2Ba ∂gydF4y2Ba PgydF4y2Ba fgydF4y2Ba ∂gydF4y2Ba tgydF4y2Ba ×gydF4y2Ba ∂gydF4y2Ba ψgydF4y2Ba ∂gydF4y2Ba PgydF4y2Ba fgydF4y2Ba ×gydF4y2Ba ∂gydF4y2Ba tgydF4y2Ba ∂gydF4y2Ba tgydF4y2Ba 一个gydF4y2Ba 。gydF4y2Ba

和gydF4y2Ba ∂gydF4y2Ba PgydF4y2Ba fgydF4y2Ba /gydF4y2Ba ∂gydF4y2Ba tgydF4y2Ba 推导出最终的物质平衡方程(方程(gydF4y2Ba 7gydF4y2Ba));以类似的方式,我们可以获得gydF4y2Ba ∂gydF4y2Ba ψgydF4y2Ba /gydF4y2Ba ∂gydF4y2Ba PgydF4y2Ba fgydF4y2Ba 和gydF4y2Ba ∂gydF4y2Ba tgydF4y2Ba /gydF4y2Ba ∂gydF4y2Ba tgydF4y2Ba 一个gydF4y2Ba 通过定义的拟压力(方程(gydF4y2Ba 9gydF4y2Ba))和新拟时间方程(方程(gydF4y2Ba 14gydF4y2Ba分别)。最后,方程(gydF4y2Ba 15gydF4y2Ba)成为gydF4y2Ba (16)gydF4y2Ba ∂gydF4y2Ba ψgydF4y2Ba ∂gydF4y2Ba tgydF4y2Ba 一个gydF4y2Ba =gydF4y2Ba −gydF4y2Ba 2gydF4y2Ba VgydF4y2Ba figydF4y2Ba 问gydF4y2Ba tgydF4y2Ba kgydF4y2Ba rggydF4y2Ba tgydF4y2Ba PgydF4y2Ba fgydF4y2Ba ZgydF4y2Ba 。gydF4y2Ba

等效气率定义为(gydF4y2Ba 48gydF4y2Ba]gydF4y2Ba (17)gydF4y2Ba 问gydF4y2Ba ggydF4y2Ba ∗gydF4y2Ba =gydF4y2Ba 1gydF4y2Ba kgydF4y2Ba rggydF4y2Ba tgydF4y2Ba 问gydF4y2Ba ggydF4y2Ba BgydF4y2Ba 胃肠道gydF4y2Ba +gydF4y2Ba 问gydF4y2Ba wgydF4y2Ba BgydF4y2Ba wgydF4y2Ba 。gydF4y2Ba

替换gydF4y2Ba PgydF4y2Ba fgydF4y2Ba /gydF4y2Ba ZgydF4y2Ba 在方程(gydF4y2Ba 16gydF4y2Ba使用真正的气体定律),然后把它与等效气体速度(方程(gydF4y2Ba 17gydF4y2Ba)):gydF4y2Ba (18)gydF4y2Ba ∂gydF4y2Ba ψgydF4y2Ba ∂gydF4y2Ba tgydF4y2Ba 一个gydF4y2Ba =gydF4y2Ba −gydF4y2Ba 2gydF4y2Ba VgydF4y2Ba figydF4y2Ba PgydF4y2Ba 我gydF4y2Ba ZgydF4y2Ba 我gydF4y2Ba 问gydF4y2Ba ggydF4y2Ba ∗gydF4y2Ba 。gydF4y2Ba

用gydF4y2Ba ∂gydF4y2Ba ψgydF4y2Ba /gydF4y2Ba ∂gydF4y2Ba tgydF4y2Ba 一个gydF4y2Ba (方程(gydF4y2Ba 18gydF4y2Ba在方程())gydF4y2Ba 11gydF4y2Ba),我们可以建立两阶段扩散方程:gydF4y2Ba (19)gydF4y2Ba ∂gydF4y2Ba 2gydF4y2Ba ψgydF4y2Ba PgydF4y2Ba fgydF4y2Ba ∂gydF4y2Ba ygydF4y2Ba 2gydF4y2Ba =gydF4y2Ba ϕgydF4y2Ba fgydF4y2Ba KgydF4y2Ba fgydF4y2Ba −gydF4y2Ba 2gydF4y2Ba VgydF4y2Ba figydF4y2Ba PgydF4y2Ba 我gydF4y2Ba ZgydF4y2Ba 我gydF4y2Ba 问gydF4y2Ba ggydF4y2Ba ∗gydF4y2Ba 。gydF4y2Ba

然后,下面的内外边界条件可以解决:gydF4y2Ba ygydF4y2Ba =gydF4y2Ba YgydF4y2Ba egydF4y2Ba ,gydF4y2Ba ∂gydF4y2Ba ψgydF4y2Ba PgydF4y2Ba fgydF4y2Ba /gydF4y2Ba ∂gydF4y2Ba ygydF4y2Ba =gydF4y2Ba 0gydF4y2Ba ,gydF4y2Ba ygydF4y2Ba =gydF4y2Ba 0gydF4y2Ba ,gydF4y2Ba ψgydF4y2Ba PgydF4y2Ba fgydF4y2Ba =gydF4y2Ba ψgydF4y2Ba PgydF4y2Ba wfgydF4y2Ba 。gydF4y2Ba

定义断裂存储系数gydF4y2Ba (20)gydF4y2Ba CgydF4y2Ba 圣gydF4y2Ba =gydF4y2Ba VgydF4y2Ba figydF4y2Ba ZgydF4y2Ba 我gydF4y2Ba 2gydF4y2Ba PgydF4y2Ba 我gydF4y2Ba 。gydF4y2Ba

然后,获得以下关系:gydF4y2Ba (21)gydF4y2Ba ψgydF4y2Ba PgydF4y2Ba ¯gydF4y2Ba fgydF4y2Ba −gydF4y2Ba ψgydF4y2Ba PgydF4y2Ba wfgydF4y2Ba 问gydF4y2Ba ggydF4y2Ba ∗gydF4y2Ba =gydF4y2Ba ϕgydF4y2Ba fgydF4y2Ba KgydF4y2Ba fgydF4y2Ba 1gydF4y2Ba CgydF4y2Ba 圣gydF4y2Ba YgydF4y2Ba egydF4y2Ba 2gydF4y2Ba 3gydF4y2Ba 。gydF4y2Ba

代入方程(gydF4y2Ba 20.gydF4y2Ba)方程(gydF4y2Ba 18gydF4y2Ba)和集成获得gydF4y2Ba (22)gydF4y2Ba tgydF4y2Ba 一个gydF4y2Ba =gydF4y2Ba CgydF4y2Ba 圣gydF4y2Ba ψgydF4y2Ba PgydF4y2Ba 我gydF4y2Ba −gydF4y2Ba ψgydF4y2Ba PgydF4y2Ba ¯gydF4y2Ba fgydF4y2Ba 问gydF4y2Ba ggydF4y2Ba ∗gydF4y2Ba 。gydF4y2Ba

结合方程(gydF4y2Ba 22gydF4y2Ba)和方程(gydF4y2Ba 15gydF4y2Ba),我们可以修改最后的两相流模型早期天然气生产阶段如下:gydF4y2Ba (23)gydF4y2Ba ψgydF4y2Ba PgydF4y2Ba 我gydF4y2Ba −gydF4y2Ba ψgydF4y2Ba PgydF4y2Ba wfgydF4y2Ba 问gydF4y2Ba ggydF4y2Ba ∗gydF4y2Ba =gydF4y2Ba 1gydF4y2Ba CgydF4y2Ba 圣gydF4y2Ba tgydF4y2Ba 一个gydF4y2Ba +gydF4y2Ba ϕgydF4y2Ba fgydF4y2Ba KgydF4y2Ba fgydF4y2Ba 1gydF4y2Ba CgydF4y2Ba 圣gydF4y2Ba YgydF4y2Ba egydF4y2Ba 2gydF4y2Ba 3gydF4y2Ba 。gydF4y2Ba

理论上,阴谋的rate-normalized拟压力(RNP)与拟时间应该产生直线关系。拟时间函数作为一个独立变量和方程左边作为因变量,斜率和截距可获得根据拟合曲线,和等效裂缝孔隙度之间的关系/半身的有效裂缝系统和有效裂缝系统的渗透率可以描述如下:gydF4y2Ba (24)gydF4y2Ba CgydF4y2Ba 圣gydF4y2Ba =gydF4y2Ba 1gydF4y2Ba 坡gydF4y2Ba =gydF4y2Ba VgydF4y2Ba figydF4y2Ba ZgydF4y2Ba 我gydF4y2Ba 2gydF4y2Ba PgydF4y2Ba 我gydF4y2Ba ,gydF4y2Ba (25)gydF4y2Ba YgydF4y2Ba DgydF4y2Ba =gydF4y2Ba ygydF4y2Ba 量gydF4y2Ba 轴gydF4y2Ba 拦截gydF4y2Ba 坡gydF4y2Ba =gydF4y2Ba ϕgydF4y2Ba fgydF4y2Ba KgydF4y2Ba fgydF4y2Ba YgydF4y2Ba egydF4y2Ba 2gydF4y2Ba 3gydF4y2Ba 。gydF4y2Ba

3.4。分析过程gydF4y2Ba

我们提出以下分析过程:gydF4y2Ba

获取和处理水和气体返排数据探索一个v型气水比的趋势(见图gydF4y2Ba 2gydF4y2Ba和gydF4y2Ba 3gydF4y2Ba与方程(gydF4y2Ba 2gydF4y2Ba))gydF4y2Ba

返排早期(出路是杰出的返排晚期(LGP)根据吉尼斯世界纪录的v型趋势曲线(见图gydF4y2Ba 3gydF4y2Ba)gydF4y2Ba

进行简单的出路裂缝网络系统模型(见图gydF4y2Ba 4gydF4y2Ba)gydF4y2Ba

计算有效压缩方程(gydF4y2Ba 3gydF4y2Ba)推导出两相物质平衡方程(MBE)(方程(gydF4y2Ba 7gydF4y2Ba))的断裂系统gydF4y2Ba

定义了拟压力和拟时间函数(方程(gydF4y2Ba 9gydF4y2Ba)和(gydF4y2Ba 10gydF4y2Ba)推导出气体扩散方程(方程(gydF4y2Ba 11gydF4y2Ba))的断裂系统gydF4y2Ba

气-水相对渗透率曲线(见图gydF4y2Ba 5gydF4y2Ba)和无因次渗透率与有效应力的曲线(见图gydF4y2Ba 8gydF4y2Ba通过实验)gydF4y2Ba

计算气体组成耦合应力敏感性(见图gydF4y2Ba 9gydF4y2Ba由图)之后的步骤gydF4y2Ba 6gydF4y2Ba将气体扩散方程(方程(gydF4y2Ba 11gydF4y2Ba)到两相扩散方程(方程(gydF4y2Ba 19gydF4y2Ba))gydF4y2Ba

修改最后的两相流模式(方程(gydF4y2Ba 23gydF4y2Ba通过结合方程())gydF4y2Ba 22gydF4y2Ba)和方程(gydF4y2Ba 15gydF4y2Ba)gydF4y2Ba

情节rate-normalized压力变化与拟时间(见图gydF4y2Ba 10gydF4y2Ba)gydF4y2Ba

计算SRV由方程(gydF4y2Ba 24gydF4y2Ba由方程()和有效裂缝系统渗透率gydF4y2Ba 25gydF4y2Ba)gydF4y2Ba

验证分析模型对微震的数据(见图gydF4y2Ba 11gydF4y2Ba)gydF4y2Ba

分析出路数据属于四川南部的两个井:rate-normalized压力与拟时间变化。3 (a)和(b) 4。gydF4y2Ba

(一)效果图3和4的微震的探测技术。(b)对比分析模型和微震的数据。gydF4y2Ba

返排的数据我们需要获得包括生产速度和压力和累积生产数据资料。然后,我们得到了一个v型趋势气水返排率曲线通过处理数据。因此,建立了一个分析模型的出路。最后,计算结果与微震的数据验证的数学结果。gydF4y2Ba

 4所示。结果和讨论gydF4y2Ba 4.1。刺激储集层体积gydF4y2Ba

我们上面提出的分析模型适用于分析的返排数据3和4。然而,有几个问题需要解决和讨论:gydF4y2Ba

很难获得一个适当的初始气体饱和度断裂系统从实际现场数据。与传统的数值模拟,这个参数是未知的实际字段数据。压裂操作完成后,有效的断裂系统连接井筒充满了压裂液约。因此,开好后返排操作在生产之前,恢复的体积压裂液在地上再度充满了气体。因此,在这篇文章中,我们需要一个合理的价值是经济复苏压裂液注入总额的百分比作为初始裂缝系统气体饱和度计算。根据徐的论文(gydF4y2Ba 48gydF4y2Ba),有一个小显著区别初始气体饱和度和复苏的实际值百分比。然而,这还不够严格。然后,其他可能的值代入计算我们得到新电脑。结果表明,该值对SRV估计的影响没有明显的误差范围内,和微震的不到1%的差异数据gydF4y2Ba

由于,压裂液可能泄漏到现有的活动自然骨折和矩阵在注射;返排液比例减少,但断裂系统的原始气体饱和度增加。因此,当注入液量的很大一部分不有助于创建骨折卷,这个值的初始气体饱和度计算实际上是低的。未来工作的一个方向是考虑使用结果迭代优化初始气体饱和度值。gydF4y2Ba

裂缝闭合的价值在总压缩系数不能准确表达。幸运的是,发现这个值没有明显的计算过程中对结果的影响。因为这个值是指的逆断裂刚度在处理骨折时,它被认为是给一个合理的价值在随后的工作从岩石断裂力学的角度gydF4y2Ba

最终的模型计算结果如图gydF4y2Ba 10gydF4y2Ba。自微震的检测技术被用来估计最后两井液压骨折刺激的影响。我们比较了计算结果与微震的数据,最后发现计算结果之间的差异的分析模型和微震的数据是次要的,如图gydF4y2Ba 11gydF4y2Ba。这个原因是,大量的骨折不连接有效时计入最终结果与断裂系统使用微震的检测技术。这些“死”骨折斜微震的结果。gydF4y2Ba

 4.2。有效的断裂系统刺激储层的渗透率gydF4y2Ba

如表所示gydF4y2Ba 1gydF4y2Ba后,获得rate-normalized压力变化的斜率与拟时间(见图gydF4y2Ba 10gydF4y2Ba),如果是通过方程(gydF4y2Ba 17gydF4y2Ba)。液压骨折的平均长度的基础上提供的微震数据的钻井平台3和4是位于核的平均孔隙度和压裂后通过实验,渗透率的有效刺激地区断裂系统计算通过使用方程(gydF4y2Ba 18gydF4y2Ba)。结果如表所示gydF4y2Ba 1gydF4y2Ba,绝对偏差指的是区别SRV计算的分析模型和微震的估计的数据,和相对偏差是指绝对偏差的比值SRV微震的估计数据。孔隙度的数值误差的有效渗透率的断裂系统负责错误。gydF4y2Ba

不同的结果的分析模型和微震的数据。gydF4y2Ba

参数名称gydF4y2Ba 3gydF4y2Ba 好4gydF4y2Ba
电脑计算的分析模型gydF4y2Ba 3.4811gydF4y2Ba ×gydF4y2Ba 10gydF4y2Ba 7gydF4y2Ba 米gydF4y2Ba 3gydF4y2Ba 3.8478gydF4y2Ba ×gydF4y2Ba 10gydF4y2Ba 7gydF4y2Ba 米gydF4y2Ba 3gydF4y2Ba
SRV微震的估计的数据gydF4y2Ba 3.5966gydF4y2Ba ×gydF4y2Ba 10gydF4y2Ba 7gydF4y2Ba 米gydF4y2Ba 3gydF4y2Ba 4.1305gydF4y2Ba ×gydF4y2Ba 10gydF4y2Ba 7gydF4y2Ba 米gydF4y2Ba 3gydF4y2Ba
绝对偏差gydF4y2Ba 1.1150gydF4y2Ba ×gydF4y2Ba 10gydF4y2Ba 6gydF4y2Ba 米gydF4y2Ba 3gydF4y2Ba 2.8270gydF4y2Ba ×gydF4y2Ba 10gydF4y2Ba 6gydF4y2Ba 米gydF4y2Ba 3gydF4y2Ba
相对偏差gydF4y2Ba 3.2%gydF4y2Ba 6.8%gydF4y2Ba
有效的断裂系统渗透率gydF4y2Ba 0.301医学博士gydF4y2Ba 0.244医学博士gydF4y2Ba
平均半身的液压骨折gydF4y2Ba 320米gydF4y2Ba 285米gydF4y2Ba
有效的断裂体系孔隙度gydF4y2Ba 15.5%gydF4y2Ba 13.8%gydF4y2Ba

如果SRV估计通过微震的数据作为比较的参考,具体计算过程如下:gydF4y2Ba (26)gydF4y2Ba dgydF4y2Ba 一个gydF4y2Ba =gydF4y2Ba VgydF4y2Ba 一个gydF4y2Ba −gydF4y2Ba VgydF4y2Ba 米gydF4y2Ba ,gydF4y2Ba (27)gydF4y2Ba dgydF4y2Ba RgydF4y2Ba =gydF4y2Ba VgydF4y2Ba 一个gydF4y2Ba −gydF4y2Ba VgydF4y2Ba 米gydF4y2Ba VgydF4y2Ba 米gydF4y2Ba ,gydF4y2Ba 在哪里gydF4y2Ba dgydF4y2Ba 一个gydF4y2Ba 和gydF4y2Ba dgydF4y2Ba RgydF4y2Ba 之间是绝对偏差和相对偏差计算的SRV分析模型和微震的估计的数据,分别mgydF4y2Ba3gydF4y2Ba;和gydF4y2Ba VgydF4y2Ba 一个gydF4y2Ba 和gydF4y2Ba VgydF4y2Ba 米gydF4y2Ba SRV计算的分析模型和微震的SRV估计的数据,mgydF4y2Ba3gydF4y2Ba。gydF4y2Ba

 4.3。不考虑应力敏感性分析模型gydF4y2Ba

如果不考虑应力敏感效应的过程中组成计算,组成曲线的偏差将发生,如图gydF4y2Ba 12gydF4y2Ba。和曲线有一个很大的错误被替换成步骤所示部分gydF4y2Ba 3.4gydF4y2Ba结果在图gydF4y2Ba 13gydF4y2Ba。因此,RNP曲线的斜率和截距都改变了。我们获得了新的曲线的斜率和截距和计算一个新的SRV根据方程(gydF4y2Ba 17gydF4y2Ba)和(gydF4y2Ba 18gydF4y2Ba)。然后,计算结果和相应的偏差,如图gydF4y2Ba 14gydF4y2Ba被获得。gydF4y2Ba

动态相对渗透率函数随着时间的推移,没有应力敏感性两个井属于四川南部。3 (a)组成。(b)组成的4。gydF4y2Ba

分析不考虑应力敏感性的出路数据属于四川南部的两个井页岩气:rate-normalized压力与拟时间变化。3 (a)和(b) 4。gydF4y2Ba

比较分析模型有或没有考虑应力敏感性和微震的数据。gydF4y2Ba

可以看出,相对偏差计算不考虑计算的应力敏感性效应大于考虑应力敏感性的影响。此外,由于真正的SRV高估了微震的数据,模型的计算结果不考虑应力敏感性的影响高于微震的数据,这也进一步表明,计算结果不考虑应力敏感性效应在高压区域将会出现一个很大的错误,这是不容忽视的。我们也计算了有效裂缝系统的渗透性,和3和4的有效渗透率是0.283 0.432医学和医学博士。可以看出,计算结果比考虑应力敏感性43.50%和15.75%的效果。gydF4y2Ba

 5。结论gydF4y2Ba

在这篇文章中,一个全面的模型来预测SRV在页岩气储层压裂后使用返排数据。包括应力敏感性。结果SRV与微震的数据进行比较。可以得出以下结论:gydF4y2Ba

返排的数据表明,气水比v型,即。,the early descending stage and the late ascending stage in this field, which can be used to estimate the volume of the effective fracture system

应力敏感性是一个关键因素影响的有效渗透率的断裂系统以及深水页岩气储层。一旦它被忽视,估计会高估了gydF4y2Ba

考虑应力敏感性预测深水救生艇在这个典型的块。结果表明,刺激了储集层体积的相对偏差计算模型和微震的数据小于10%,表明该方法可以提供合理的预测gydF4y2Ba

数据可用性gydF4y2Ba

测试数据用于支持本研究的结果包括在本文中。读者可以获得数据支持研究结果从测试数据表。gydF4y2Ba

 的利益冲突gydF4y2Ba

作者宣称没有利益冲突有关的出版。gydF4y2Ba

 确认gydF4y2Ba

这项研究是由中国国家科技重大项目(2017 zx05009 - 005),基础研究基金为中央大学(2652018209),和中国国家自然科学基金(51804282)。gydF4y2Ba

 王gydF4y2Ba G。gydF4y2Ba 贾gydF4y2Ba 一个。gydF4y2Ba 魏gydF4y2Ba Y。gydF4y2Ba 肖gydF4y2Ba C。gydF4y2Ba 瞬态压力分析multifractured水平井使用多重线性流模型在页岩气藏gydF4y2Ba GeofluidsgydF4y2Ba 2020年gydF4y2Ba 2020年gydF4y2Ba 20.gydF4y2Ba 8348205gydF4y2Ba 10.1155 / 2020/8348205gydF4y2Ba 周gydF4y2Ba W。gydF4y2Ba 巴纳吉gydF4y2Ba R。gydF4y2Ba 坡gydF4y2Ba b D。gydF4y2Ba SpathgydF4y2Ba J。gydF4y2Ba ThambynayagamgydF4y2Ba M。gydF4y2Ba Semianalytical生产复杂的水力裂缝网络的仿真gydF4y2Ba SPE杂志gydF4y2Ba 2014年gydF4y2Ba 19gydF4y2Ba 1gydF4y2Ba 6gydF4y2Ba 18gydF4y2Ba NeuzilgydF4y2Ba c, E。gydF4y2Ba 透水粘土和页岩如何?gydF4y2Ba 水资源研究gydF4y2Ba 1994年gydF4y2Ba 30.gydF4y2Ba 2gydF4y2Ba 145年gydF4y2Ba 150年gydF4y2Ba 10.1029/93 wr02930gydF4y2Ba 2 - s2.0 - 0028192646gydF4y2Ba 琼斯gydF4y2Ba r S。gydF4y2Ba Jr。gydF4y2Ba 的,gydF4y2Ba B。gydF4y2Ba 因特网gydF4y2Ba J。gydF4y2Ba 从早期的返排数据估算储层压力gydF4y2Ba 非常规资源技术会议gydF4y2Ba 2014年8月gydF4y2Ba 科罗拉多州丹佛市gydF4y2Ba 25gydF4y2Ba 27gydF4y2Ba ZanganehgydF4y2Ba B。gydF4y2Ba SoroushgydF4y2Ba M。gydF4y2Ba Williams-KovacsgydF4y2Ba j . D。gydF4y2Ba 克拉克森gydF4y2Ba c·R。gydF4y2Ba 参数影响负载恢复和石油油井水力压裂后突破时间很紧张gydF4y2Ba SPE / CSUR非常规资源会议gydF4y2Ba 2015年10月gydF4y2Ba 卡尔加里,阿尔伯塔省gydF4y2Ba 20.gydF4y2Ba 22gydF4y2Ba MohammadnejadgydF4y2Ba 一个。gydF4y2Ba 雪莱gydF4y2Ba r F。gydF4y2Ba 雷曼兄弟gydF4y2Ba l . V。gydF4y2Ba 沙阿gydF4y2Ba K。gydF4y2Ba GusaingydF4y2Ba D。gydF4y2Ba 康威gydF4y2Ba m . T。gydF4y2Ba 脆性页岩裂缝网络模型的发展gydF4y2Ba SPE水力压裂技术会议gydF4y2Ba 2013年2月gydF4y2Ba 林地,美国德克萨斯州gydF4y2Ba 4gydF4y2Ba 6gydF4y2Ba 徐gydF4y2Ba J。gydF4y2Ba 郭gydF4y2Ba C。gydF4y2Ba 腾gydF4y2Ba W。gydF4y2Ba 魏gydF4y2Ba M。gydF4y2Ba 江gydF4y2Ba R。gydF4y2Ba 致密油/气的生产性能分析水库考虑刺激水库卷使用椭圆流gydF4y2Ba 天然气的科学与工程》杂志上gydF4y2Ba 2015年gydF4y2Ba 26gydF4y2Ba 827年gydF4y2Ba 839年gydF4y2Ba 10.1016 / j.jngse.2015.06.057gydF4y2Ba 2 - s2.0 - 84937773191gydF4y2Ba 吴gydF4y2Ba M。gydF4y2Ba 丁gydF4y2Ba M。gydF4y2Ba 姚gydF4y2Ba J。gydF4y2Ba 徐gydF4y2Ba 年代。gydF4y2Ba 李gydF4y2Ba lgydF4y2Ba 李gydF4y2Ba X。gydF4y2Ba 多裂缝水平井的压力瞬态分析复合页岩气储层边界元法gydF4y2Ba 石油科学与工程》杂志上gydF4y2Ba 2018年gydF4y2Ba 162年gydF4y2Ba 84年gydF4y2Ba 101年gydF4y2Ba 10.1016 / j.petrol.2017.12.024gydF4y2Ba 2 - s2.0 - 85038258617gydF4y2Ba 程ydF4y2Ba Z。gydF4y2Ba 廖gydF4y2Ba X。gydF4y2Ba 赵gydF4y2Ba X。gydF4y2Ba 余gydF4y2Ba W。gydF4y2Ba SepehrnoorigydF4y2Ba K。gydF4y2Ba 工作流基于semianalytical模型来估计刺激的属性储层致密油水井的体积gydF4y2Ba 石油科学与工程》杂志上gydF4y2Ba 2019年gydF4y2Ba 178年gydF4y2Ba 892年gydF4y2Ba 903年gydF4y2Ba 10.1016 / j.petrol.2019.04.002gydF4y2Ba 2 - s2.0 - 85064178461gydF4y2Ba 同类gydF4y2Ba D。gydF4y2Ba 安德森gydF4y2Ba d . M。gydF4y2Ba StottsgydF4y2Ba g·w·J。gydF4y2Ba MattargydF4y2Ba lgydF4y2Ba BlasingamegydF4y2Ba t。gydF4y2Ba 生产数据分析、缺陷诊断gydF4y2Ba SPE Resrervoir评价&工程gydF4y2Ba 2013年gydF4y2Ba 13gydF4y2Ba 3gydF4y2Ba 538年gydF4y2Ba 552年gydF4y2Ba 同类gydF4y2Ba D。gydF4y2Ba 柯里gydF4y2Ba s M。gydF4y2Ba 西蒙斯gydF4y2Ba D。gydF4y2Ba 冲gydF4y2Ba j . A。gydF4y2Ba 布鲁萨德gydF4y2Ba n . J。gydF4y2Ba BlasingamegydF4y2Ba t。gydF4y2Ba 全面工作流程早期的分析和解释返排致密气和页岩储层的井数据系统gydF4y2Ba SPE年度技术会议和展览gydF4y2Ba 2010年9月gydF4y2Ba 意大利的佛罗伦萨gydF4y2Ba 19gydF4y2Ba 22gydF4y2Ba AbbasigydF4y2Ba m·A。gydF4y2Ba DehghanpourgydF4y2Ba H。gydF4y2Ba 霍克斯gydF4y2Ba r . V。gydF4y2Ba 返排分析断裂特征gydF4y2Ba SPE加拿大非常规资源会议gydF4y2Ba 2012年11月gydF4y2Ba 卡尔加里,阿尔伯塔省gydF4y2Ba AbbasigydF4y2Ba m·A。gydF4y2Ba EzulikegydF4y2Ba d . O。gydF4y2Ba DehghanpourgydF4y2Ba H。gydF4y2Ba 霍克斯gydF4y2Ba r . V。gydF4y2Ba 返排率和压力瞬态行为的比较研究在multifractured水平井在致密气和石油储层完成gydF4y2Ba 天然气的科学与工程》杂志上gydF4y2Ba 2014年gydF4y2Ba 17gydF4y2Ba 2gydF4y2Ba 82年gydF4y2Ba 93年gydF4y2Ba 10.1016 / j.jngse.2013.12.007gydF4y2Ba 2 - s2.0 - 84893428914gydF4y2Ba 张gydF4y2Ba Y。gydF4y2Ba Ehlig-EconomidesgydF4y2Ba C。gydF4y2Ba 占剩下的注入在页岩气井压裂液gydF4y2Ba SPE /中部凹陷非常规资源技术会议gydF4y2Ba 2014年8月gydF4y2Ba 美国科罗拉多州丹佛市gydF4y2Ba 25gydF4y2Ba 27gydF4y2Ba GhanbarigydF4y2Ba E。gydF4y2Ba 轴承gydF4y2Ba D。gydF4y2Ba AbbasigydF4y2Ba m·A。gydF4y2Ba DehghanpourgydF4y2Ba H。gydF4y2Ba 返排体积和化学分析评价负荷恢复及其对早期生产gydF4y2Ba SPE加拿大非常规资源会议gydF4y2Ba 2013年11月gydF4y2Ba 卡尔加里,加拿大阿尔伯塔省gydF4y2Ba 5gydF4y2Ba 7gydF4y2Ba AbbasigydF4y2Ba m·A。gydF4y2Ba 返排率和压力瞬态行为的比较研究在酸洗裂缝水平井,[博士。论文)gydF4y2Ba 2013年gydF4y2Ba 加拿大亚伯达省的埃德蒙顿,gydF4y2Ba 阿尔伯塔大学gydF4y2Ba 克拉克森gydF4y2Ba c·R。gydF4y2Ba 建模的两相返排multifractured水平井在页岩完成gydF4y2Ba SPE加拿大非常规资源会议gydF4y2Ba 2012年11月gydF4y2Ba 卡尔加里,加拿大阿尔伯塔省gydF4y2Ba Williams-KovacsgydF4y2Ba j . D。gydF4y2Ba 克拉克森gydF4y2Ba c·R。gydF4y2Ba 建模的两相返排酸洗裂缝水平紧气井刺激与氮增能压裂液gydF4y2Ba SPE加拿大非常规资源会议gydF4y2Ba 2013年11月gydF4y2Ba 卡尔加里,阿尔伯塔省gydF4y2Ba 道gydF4y2Ba Z。gydF4y2Ba 香坊gydF4y2Ba lgydF4y2Ba 总裁gydF4y2Ba Y。gydF4y2Ba 京gydF4y2Ba lgydF4y2Ba YonghuigydF4y2Ba W。gydF4y2Ba 越南盾gydF4y2Ba F。gydF4y2Ba 关井时间对返排率和生产率的影响的页岩气井gydF4y2Ba 天然气工业gydF4y2Ba 2017年gydF4y2Ba 37gydF4y2Ba 8gydF4y2Ba 48gydF4y2Ba 60gydF4y2Ba AdefidipegydF4y2Ba o . A。gydF4y2Ba 徐gydF4y2Ba Y。gydF4y2Ba DehghanpourgydF4y2Ba H。gydF4y2Ba 病毒传染gydF4y2Ba c·J。gydF4y2Ba 直接从页岩气井天然气生产:返排一个两阶段模型gydF4y2Ba SPE非常规资源联盟锦标赛gydF4y2Ba 2014年4月gydF4y2Ba 林地、德州gydF4y2Ba AdefidipegydF4y2Ba o . A。gydF4y2Ba 徐gydF4y2Ba Y。gydF4y2Ba DehghanpourgydF4y2Ba H。gydF4y2Ba 病毒传染gydF4y2Ba c·J。gydF4y2Ba 从早期生产数据估算有效裂缝体积:物质平衡方法gydF4y2Ba SPE / CSUR Conference-Canada非常规资源gydF4y2Ba 2014年10月gydF4y2Ba 卡尔加里,阿尔伯塔省gydF4y2Ba 徐gydF4y2Ba Y。gydF4y2Ba AdefidipegydF4y2Ba o . A。gydF4y2Ba DehghanpourgydF4y2Ba H。gydF4y2Ba 容量分析两相返排fata的断裂特性gydF4y2Ba SPE西部地区会议gydF4y2Ba 2015年4月gydF4y2Ba 加登格罗夫,加州gydF4y2Ba 徐gydF4y2Ba Y。gydF4y2Ba 傅gydF4y2Ba Y。gydF4y2Ba EzulikegydF4y2Ba d . O。gydF4y2Ba DehghanpourgydF4y2Ba H。gydF4y2Ba 病毒传染gydF4y2Ba c·J。gydF4y2Ba 使用开放的坦克模型建模两相返排数据gydF4y2Ba SPE / CSUR非常规资源会议gydF4y2Ba 2015年10月gydF4y2Ba 卡尔加里,阿尔伯塔省gydF4y2Ba EzulikegydF4y2Ba d . O。gydF4y2Ba DehghanpourgydF4y2Ba H。gydF4y2Ba 霍克斯gydF4y2Ba r . V。gydF4y2Ba 理解返排作为瞬态两相驱替过程:线性双孔隙度模型的一个扩展gydF4y2Ba 心肺复苏术非常规资源Conference-CanadagydF4y2Ba 2013年11月gydF4y2Ba 卡尔加里,阿尔伯塔省gydF4y2Ba EzulikegydF4y2Ba d . O。gydF4y2Ba DehghanpourgydF4y2Ba H。gydF4y2Ba 工作流analysis-creating返排数据价值的混乱gydF4y2Ba 2014年非常规资源技术会议gydF4y2Ba 2014年8月gydF4y2Ba 科罗拉多州丹佛市gydF4y2Ba EzulikegydF4y2Ba d . O。gydF4y2Ba DehghanpourgydF4y2Ba H。gydF4y2Ba 病毒传染gydF4y2Ba c·J。gydF4y2Ba flowback-guided生产数据分析方法在致密储层gydF4y2Ba 2014 SPE / CSUR Conference-Canada非常规资源gydF4y2Ba 2014年10月gydF4y2Ba 卡尔加里,阿尔伯塔省gydF4y2Ba 克拉克森gydF4y2Ba c·R。gydF4y2Ba Williams-KovacsgydF4y2Ba j . D。gydF4y2Ba 新方法建模多阶段油井排液酸洗裂缝水平紧来确定液压骨折属性gydF4y2Ba SPE年度技术会议和展览gydF4y2Ba 2013年10月gydF4y2Ba 新奥尔良,路易斯安那州gydF4y2Ba 克拉克森gydF4y2Ba c·R。gydF4y2Ba QanbarigydF4y2Ba F。gydF4y2Ba Williams-KovacsgydF4y2Ba j . D。gydF4y2Ba 创新使用rate-transient分析方法获得低渗透性储层表现出多相流的水力裂缝属性gydF4y2Ba 前缘gydF4y2Ba 2014年gydF4y2Ba 33gydF4y2Ba 10gydF4y2Ba 1108年gydF4y2Ba 1122年gydF4y2Ba 10.1190 / tle33101108.1gydF4y2Ba 2 - s2.0 - 84910149934gydF4y2Ba Williams-KovacsgydF4y2Ba j . D。gydF4y2Ba 克拉克森gydF4y2Ba c·R。gydF4y2Ba 随机建模的多阶段的排液酸洗裂缝水平致密油水井gydF4y2Ba SPE非常规资源Conference-CanadagydF4y2Ba 2013年11月gydF4y2Ba 卡尔加里,阿尔伯塔省gydF4y2Ba Williams-KovacsgydF4y2Ba j . D。gydF4y2Ba 克拉克森gydF4y2Ba c·R。gydF4y2Ba 随机建模的两相返排酸洗裂缝水平井水力裂缝属性估计和预测生产gydF4y2Ba SPE非常规资源会议美国gydF4y2Ba 2013年4月gydF4y2Ba 林地、德州gydF4y2Ba Williams-KovacsgydF4y2Ba j . D。gydF4y2Ba 克拉克森gydF4y2Ba c·R。gydF4y2Ba 修改方法建模阶段从酸洗裂缝水平页岩气井排液gydF4y2Ba 非常规资源技术会议gydF4y2Ba 2015年7月gydF4y2Ba 得克萨斯州圣安东尼奥市gydF4y2Ba Williams-KovacsgydF4y2Ba j . D。gydF4y2Ba 克拉克森gydF4y2Ba c·R。gydF4y2Ba 修改方法建模阶段从酸洗裂缝水平页岩气井排液gydF4y2Ba 天然气的科学与工程》杂志上gydF4y2Ba 2016年gydF4y2Ba 30.gydF4y2Ba 127年gydF4y2Ba 147年gydF4y2Ba 10.1016 / j.jngse.2016.02.003gydF4y2Ba 2 - s2.0 - 84958171408gydF4y2Ba 傅gydF4y2Ba Y。gydF4y2Ba EzulikegydF4y2Ba d . O。gydF4y2Ba DehghanpourgydF4y2Ba H。gydF4y2Ba 从早期的单相返排数据估算有效裂缝孔隙体积和裂缝设计参数联系起来gydF4y2Ba SPE / CSUR非常规资源会议gydF4y2Ba 2015年10月gydF4y2Ba 卡尔加里,阿尔伯塔省gydF4y2Ba AlkouhgydF4y2Ba 答:B。gydF4y2Ba McKettagydF4y2Ba 年代。gydF4y2Ba WattenbargergydF4y2Ba r。gydF4y2Ba 估计effective-fracture卷使用water-flowback和页岩气水井生产数据gydF4y2Ba 加拿大的石油技术杂志》上gydF4y2Ba 2014年gydF4y2Ba 53gydF4y2Ba 5gydF4y2Ba 293年gydF4y2Ba 303年gydF4y2Ba 王gydF4y2Ba X。gydF4y2Ba 朱gydF4y2Ba Y。gydF4y2Ba 傅gydF4y2Ba C。gydF4y2Ba 实验调查stress-dependent渗透率的Longmaxi页岩形成gydF4y2Ba 石油科学与工程》杂志上gydF4y2Ba 2019年gydF4y2Ba 175年gydF4y2Ba 932年gydF4y2Ba 947年gydF4y2Ba 10.1016 / j.petrol.2019.01.037gydF4y2Ba 2 - s2.0 - 85059911954gydF4y2Ba 杨ydF4y2Ba X。gydF4y2Ba 太阳gydF4y2Ba J。gydF4y2Ba 刘gydF4y2Ba D。gydF4y2Ba 数值模拟的页岩气多尺度渗流mechanism-coupled应力敏感性gydF4y2Ba 化学杂志gydF4y2Ba 2019年gydF4y2Ba 2019年gydF4y2Ba 13gydF4y2Ba 7387234gydF4y2Ba 10.1155 / 2019/7387234gydF4y2Ba 2 - s2.0 - 85064347902gydF4y2Ba 刘gydF4y2Ba W。gydF4y2Ba 刘gydF4y2Ba J。gydF4y2Ba 蔡gydF4y2Ba M。gydF4y2Ba 罗gydF4y2Ba C。gydF4y2Ba 史gydF4y2Ba X。gydF4y2Ba 张gydF4y2Ba J。gydF4y2Ba 页岩的孔隙演化特征Longmaxi形成,四川盆地gydF4y2Ba 石油研究gydF4y2Ba 2017年gydF4y2Ba 2gydF4y2Ba 4gydF4y2Ba 291年gydF4y2Ba 300年gydF4y2Ba 10.1016 / j.ptlrs.2017.03.003gydF4y2Ba 程ydF4y2Ba 年代。gydF4y2Ba 朱gydF4y2Ba Y。gydF4y2Ba 程ydF4y2Ba 年代。gydF4y2Ba 汉gydF4y2Ba Y。gydF4y2Ba 傅gydF4y2Ba C。gydF4y2Ba 生烃和页岩气积累Longmaxi形成,四川盆地南部中国gydF4y2Ba 海洋和石油地质学gydF4y2Ba 2017年gydF4y2Ba 86年gydF4y2Ba 248年gydF4y2Ba 258年gydF4y2Ba 王gydF4y2Ba C。gydF4y2Ba 张gydF4y2Ba B。gydF4y2Ba 陆gydF4y2Ba Y。gydF4y2Ba 蜀gydF4y2Ba Z。gydF4y2Ba 陆gydF4y2Ba Y。gydF4y2Ba 保gydF4y2Ba H。gydF4y2Ba 孟gydF4y2Ba Z。gydF4y2Ba 程ydF4y2Ba lgydF4y2Ba 页岩岩相分布特征及其控制因素的聚Formation-Member 1 Jiaoshiba Longmaxi形成的区域gydF4y2Ba 石油研究gydF4y2Ba 2018年gydF4y2Ba 3gydF4y2Ba 4gydF4y2Ba 306年gydF4y2Ba 319年gydF4y2Ba 10.1016 / j.ptlrs.2018.11.005gydF4y2Ba 越南盾gydF4y2Ba T。gydF4y2Ba 他gydF4y2Ba 年代。gydF4y2Ba 程ydF4y2Ba M。gydF4y2Ba 侯gydF4y2Ba Y。gydF4y2Ba 郭gydF4y2Ba X。gydF4y2Ba 魏gydF4y2Ba C。gydF4y2Ba 汉gydF4y2Ba Y。gydF4y2Ba 杨gydF4y2Ba R。gydF4y2Ba 石英类型和起源于古生代Wufeng-Longmaxi形成,四川盆地东部,中国:对页岩储层的孔隙度保存gydF4y2Ba 海洋和石油地质学gydF4y2Ba 2019年gydF4y2Ba 106年gydF4y2Ba 62年gydF4y2Ba 73年gydF4y2Ba 10.1016 / j.marpetgeo.2019.05.002gydF4y2Ba 2 - s2.0 - 85065229407gydF4y2Ba 程gydF4y2Ba Y。gydF4y2Ba 水动力学的影响在骨折采用水力压裂井在页岩气储层的性能gydF4y2Ba 加拿大的石油技术杂志》上gydF4y2Ba 2013年gydF4y2Ba 51gydF4y2Ba 2gydF4y2Ba 143年gydF4y2Ba 151年gydF4y2Ba 张gydF4y2Ba J。gydF4y2Ba KamenovgydF4y2Ba 一个。gydF4y2Ba 山gydF4y2Ba 答:D。gydF4y2Ba 朱gydF4y2Ba D。gydF4y2Ba 实验室测量页岩水力裂缝导率gydF4y2Ba SPE生产&操作gydF4y2Ba 2014年gydF4y2Ba 29日gydF4y2Ba 3gydF4y2Ba 216年gydF4y2Ba 227年gydF4y2Ba 10.2118 / 163839 - pagydF4y2Ba 2 - s2.0 - 84907509520gydF4y2Ba 风扇gydF4y2Ba lgydF4y2Ba 汤普森gydF4y2Ba j·W。gydF4y2Ba 罗宾逊gydF4y2Ba j . R。gydF4y2Ba 了解天然气生产机制和有效性的刺激在斯威尔页岩储层模拟gydF4y2Ba SPE非常规资源&国际石油会议gydF4y2Ba 2010年10月gydF4y2Ba 卡尔加里,阿尔伯塔省gydF4y2Ba KuchukgydF4y2Ba f·J。gydF4y2Ba 科夫gydF4y2Ba D。gydF4y2Ba 菲茨帕特里克gydF4y2Ba T。gydF4y2Ba 瞬态和不断地(双孔隙度)和离散的递减曲线分析天然裂缝性储层gydF4y2Ba SPE年度技术会议和展览gydF4y2Ba 2014年10月gydF4y2Ba 荷兰阿姆斯特丹gydF4y2Ba 拉赫曼gydF4y2Ba n·m·A。gydF4y2Ba 角gydF4y2Ba lgydF4y2Ba ZaoralgydF4y2Ba K。gydF4y2Ba 新方法计算伪时间真实气体流动物质平衡方程gydF4y2Ba 加拿大的石油技术杂志》上gydF4y2Ba 2006年gydF4y2Ba 45gydF4y2Ba 10gydF4y2Ba 36gydF4y2Ba 44gydF4y2Ba MoghadamgydF4y2Ba 年代。gydF4y2Ba JejegydF4y2Ba O。gydF4y2Ba MattargydF4y2Ba lgydF4y2Ba 先进的气体物质平衡简化格式gydF4y2Ba 加拿大的石油技术杂志》上gydF4y2Ba 2013年gydF4y2Ba 50gydF4y2Ba 1gydF4y2Ba 90年gydF4y2Ba 98年gydF4y2Ba 辛格gydF4y2Ba 诉K。gydF4y2Ba 物质平衡方程的概述(MBE)页岩气和非常规储层gydF4y2Ba SPE中东油气显示和会议gydF4y2Ba 2013年3月gydF4y2Ba 麦纳麦巴林gydF4y2Ba 徐gydF4y2Ba Y。gydF4y2Ba AdefidipegydF4y2Ba o . A。gydF4y2Ba DehghanpourgydF4y2Ba H。gydF4y2Ba 两相返排的流动物质平衡方程分析gydF4y2Ba 石油科学与工程杂志》上gydF4y2Ba 2016年gydF4y2Ba 142年gydF4y2Ba 170年gydF4y2Ba 185年gydF4y2Ba 10.1016 / j.petrol.2016.01.018gydF4y2Ba 2 - s2.0 - 84959234818gydF4y2Ba 马丁gydF4y2Ba j . C。gydF4y2Ba 简化的气驱油藏的流体方程和多相压力建立分析的理论基础gydF4y2Ba 交易的艾米gydF4y2Ba 2013年gydF4y2Ba 216年gydF4y2Ba 1gydF4y2Ba 321年gydF4y2Ba 323年gydF4y2Ba 克雷格gydF4y2Ba d . P。gydF4y2Ba 分析建模fracture-injection /下降sequenceand长垣候选人诊断测试的发展,[博士。论文)gydF4y2Ba 2006年gydF4y2Ba 德克萨斯州大学站gydF4y2Ba 德州农工大学gydF4y2Ba 张gydF4y2Ba D。gydF4y2Ba 冬天gydF4y2Ba C。gydF4y2Ba 弯矩方程方法单相流体在异构的水库gydF4y2Ba SPE杂志gydF4y2Ba 1999年gydF4y2Ba 4gydF4y2Ba 2gydF4y2Ba 118年gydF4y2Ba 127年gydF4y2Ba TabatabaiegydF4y2Ba s . H。gydF4y2Ba MattargydF4y2Ba lgydF4y2Ba Pooladi-DarvishgydF4y2Ba M。gydF4y2Ba 拟时间calculationin低渗透气藏gydF4y2Ba SPE加拿大非常规资源会议gydF4y2Ba 2013年11月gydF4y2Ba 卡尔加里,加拿大阿尔伯塔省gydF4y2Ba EzulikegydF4y2Ba d . O。gydF4y2Ba DehghanpourgydF4y2Ba H。gydF4y2Ba 返排造型作为瞬态两相损耗的过程gydF4y2Ba 天然气的科学与工程》杂志上gydF4y2Ba 2014年gydF4y2Ba 19gydF4y2Ba 258年gydF4y2Ba 278年gydF4y2Ba 10.1016 / j.jngse.2014.05.004gydF4y2Ba 2 - s2.0 - 84902462533gydF4y2Ba 赵gydF4y2Ba Y。gydF4y2Ba 王gydF4y2Ba Z。gydF4y2Ba 秦gydF4y2Ba X。gydF4y2Ba 李gydF4y2Ba J。gydF4y2Ba 杨gydF4y2Ba H。gydF4y2Ba Stress-dependent渗透率煤的裂隙网络:用晶格玻尔兹曼方法数值研究gydF4y2Ba 石油科学与工程》杂志上gydF4y2Ba 2019年gydF4y2Ba 173年gydF4y2Ba 1053年gydF4y2Ba 1064年gydF4y2Ba 10.1016 / j.petrol.2018.10.075gydF4y2Ba 2 - s2.0 - 85055863488gydF4y2Ba 俄文gydF4y2Ba Z。gydF4y2Ba 一个gydF4y2Ba K。gydF4y2Ba 胡gydF4y2Ba J。gydF4y2Ba 硫沉淀的影响在发展中含硫气藏的压敏效应gydF4y2Ba Geo-Energy研究进展gydF4y2Ba 2019年gydF4y2Ba 3gydF4y2Ba 3gydF4y2Ba 268年gydF4y2Ba 276年gydF4y2Ba 10.26804 / ager.2019.03.05gydF4y2Ba