含水饱和度的确定方法、装置、电子及存储介质制造方法及图纸

本申请提供的一种含水饱和度的确定方法、装置、电子设备及存储介质，通过预先建立含水饱和度计算式，在获取目标井的目标生产时间后，可以输入到含水饱和度计算式中，快速计算出含水饱和度值，提高了计算效率。提高了计算效率。提高了计算效率。

【技术实现步骤摘要】
含水饱和度的确定方法、装置、电子及存储介质


[0001]本申请涉及页岩气藏开发领域，特别地涉及一种含水饱和度的确定方法、装置、电子设备及存储介质。

技术介绍

[0002]页岩气井通过压裂水平井技术获得工业油气流，在返排和生产初期压裂液与气体共同产出，处于气
‑
水同流的状态。含水饱和度是研究两相渗流的基础，亦是页岩气田数值模拟及动态分析不可缺少的重要资料。目前对于含水饱和度的获取主要通过在数值模拟技术获取，但是该方法费时费力，而通过生产数据计算含水饱和度的研究相对较少。如何采用油田现场的生产数据快速、准确地计算含水饱和度，对于开展数值模拟和动态分析工作、进一步调整开发方案和提高页岩气藏开采效率具有重要意义。
[0003]目前，页岩气藏含水饱和度的核心问题在于计算压裂过程中压裂液对储层含水饱和度的影响。主要结合渗吸理论和数值模拟方法计算含水饱和度分布，渗吸量计算方法包括Handy模型、Shembreet方法和K.Makhanov方法等。其中，Handy模型使用广泛，适用于气水两相流，渗吸量与渗吸时间的平方根成正比，没有考虑毛细管间的横向流动特征。Shembreet方法在Handy模型的基础上，考虑了毛细管压力等因素，推出了渗吸量与时间平方根的正比关系。数值模拟方法需结合渗吸实验测试结果和气—液两相流动的连续性方程，建立压裂过程中压裂液渗吸侵入页岩储层的数值模型进行计算，能反映页岩气井压裂后井周围含水饱和度的重新分布过程复杂，但是耗时耗力，难以提高页岩气开采效率。

技术实现思路

[0004]针对上述问题，本申请提供一种含水饱和度的确定方法、装置、电子设备及存储介质。
[0005]本申请提供了一种含水饱和度的确定方法，包括：
[0006]获取目标井的目标生产时间；
[0007]将所述目标生产时间输入至预先建立的含水饱和度计算式中，得到所述目标生产时间对应的含水饱和度值，其中，所述含水饱和度计算式包括生产时间与含水饱和度之间的关系。
[0008]在一些实施例中，所述方法还包括：
[0009]获取所述目标井的生产时间和所述生产时间对应的日产气量、日产水量、气相相对渗透率数据、水相相对渗透率数据、含水饱和度；
[0010]基于所述生产时间、所述日产气量、日产水量拟合水气比与生产时间的第一关系曲线；
[0011]基于生产时间、气相相对渗透率数据、水相相对渗透率数据确定相对渗透率比值与时间的第一关系；
[0012]基于所述第一关系和所述第一关系曲线确定水气比与相对渗透率比值之间的转
换关系；
[0013]基于所述气相相对渗透率数据、水相相对渗透率数据及含水饱和度拟合相对渗透率比值与含水饱和度之间的第二关系；
[0014]至少基于所述转换关系和所述第二关系确定含水饱和度与生产时间之间的关系，得到含水饱和度计算式。
[0015]在一些实施例中，基于所述生产时间、所述日产气量、日产水量拟合水气比与生产时间的第一关系曲线，包括：
[0016]基于所述所述日产气量和日产水量确定水气比；
[0017]从所述水气比中挑选水气比大于预设水气比、且生产时间连续的目标水气比，得到第一关系曲线。
[0018]在一些实施例中，所述第一关系曲线为：WGR＝aln(t)+b，其中，a和b为无因次量。
[0019]在一些实施例中，所述第二关系为：其中，k
rw
为水相相对渗透率；k
rg
为气相相对渗透率，c和d为无因次量，S
w
为含水饱和度。
[0020]在一些实施例中，所述含水饱和度计算式为：
[0021]其中，ρ
w
为水密度；ρ
g
为气密度；μ
w
为水粘度；μ
g
为气体粘度；B
w
为水的体积系数；B
g
为气的体积系数。
[0022]本申请实施例提供一种含水饱和度的确定装置，包括：
[0023]第二获取模块，用于获取所述目标井的生产时间和所述生产时间对应的日产气量、日产水量、气相相对渗透率数据、水相相对渗透率数据、含水饱和度；
[0024]第一确定模块，用于基于所述生产时间、所述日产气量、日产水量拟合水气比与生产时间的第一关系曲线；
[0025]第二确定模块，用于基于生产时间、气相相对渗透率数据、水相相对渗透率数据确定相对渗透率比值与时间的第一关系；
[0026]第三确定模块，用于基于所述第一关系和所述第一关系曲线确定水气比与相对渗透率比值之间的转换关系；
[0027]第四确定模块，用于基于所述气相相对渗透率数据、水相相对渗透率数据及含水饱和度拟合相对渗透率比值与含水饱和度之间的第二关系；
[0028]第五确定模块，用于至少基于所述转换关系和所述第二关系确定含水饱和度与生产时间之间的关系，得到含水饱和度计算式。
[0029]本申请实施例提供一种电子设备，包括存储器和处理器，所述存储器上存储有计算机程序，该计算机程序被所述处理器执行时，执行上述任意一项所述含水饱和度的确定方法。
[0030]本申请实施例提供一种存储介质，该存储介质存储的计算机程序，能够被一个或多个处理器执行，能够用来实现上述任一项所述含水饱和度的确定方法。
[0031]本申请提供的一种含水饱和度的确定方法、装置、电子设备及存储介质，通过预先建立含水饱和度计算式，在获取目标井的目标生产时间后，可以输入到含水饱和度计算式中，快速计算出含水饱和度值，提高了计算效率。
附图说明
[0032]在下文中将基于实施例并参考附图来对本申请进行更详细的描述。
[0033]图1为本申请实施例提供的一种含水饱和度的确定方法的实现流程示意图；
[0034]图2为本申请实施例提供的另一种含水饱和度的确定方法的实现流程示意图；
[0035]图3为本申请实施例提供的一种生产数据随时间的变化示意图；
[0036]图4为本申请实施例提供的一种相对渗透率比值与含水饱和度的关系曲线示意图；
[0037]图5为本申请实施例提供的电子设备的组成结构示意图。
[0038]在附图中，相同的部件使用相同的附图标记，附图并未按照实际的比例绘制。
具体实施方式
[0039]为了使本申请的目的、技术方案和优点更加清楚，下面将结合附图对本申请作进一步地详细描述，所描述的实施例不应视为对本申请的限制，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其它实施例，都属于本申请保护的范围。
[0040]在以下的描述中，涉及到“一些实施例”，其描述了所有可能实施例的子集，但是可以理解，“一些实施例”可以是所有可能实施例的相同子集或不同子集，并且可以在不冲突的情况下相互结合。
[0041]如果申请文件中出现“第一\第二\第三”的类似描述则增加以下的说明，在以下的描述中，所涉及本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种含水饱和度的确定方法，其特征在于，所述方法包括：获取目标井的目标生产时间；将所述目标生产时间输入至预先建立的含水饱和度计算式中，得到所述目标生产时间对应的含水饱和度值，其中，所述含水饱和度计算式包括生产时间与含水饱和度之间的关系。2.根据权利要求1所述的方法，其特征在于，所述方法还包括：获取所述目标井的生产时间和所述生产时间对应的日产气量、日产水量、气相相对渗透率数据、水相相对渗透率数据、含水饱和度；基于所述生产时间、所述日产气量、日产水量拟合水气比与生产时间的第一关系曲线；基于生产时间、气相相对渗透率数据、水相相对渗透率数据确定相对渗透率比值与时间的第一关系；基于所述第一关系和所述第一关系曲线确定水气比与相对渗透率比值之间的转换关系；基于所述气相相对渗透率数据、水相相对渗透率数据及含水饱和度拟合相对渗透率比值与含水饱和度之间的第二关系；至少基于所述转换关系和所述第二关系确定含水饱和度与生产时间之间的关系，得到含水饱和度计算式。3.根据权利要求2所述的方法，其特征在于，基于所述生产时间、所述日产气量、日产水量拟合水气比与生产时间的第一关系曲线，包括：基于所述所述日产气量和日产水量确定水气比；从所述水气比中挑选水气比大于预设水气比、且生产时间连续的目标水气比，得到第一关系曲线。4.根据权利要求3所述的方法，其特征在于，所述第一关系曲线为：WGR＝aln(t)+b，其中，a和b为无因次量。5.根据权利要求4所述的方法，其特征在于，所述第二关系为：其中，k
rw
为水相相对渗透率；k
rg
为气相相对渗透率，c和d为无因次量，S
w
为含水饱和度。6.根据权利要求5所述的方法，其特征在于，所述含水饱和度计算式为：其中，ρ
w
为水密度；ρ

【专利技术属性】
技术研发人员：卢婷，王鸣川，龙胜祥，张殿伟，彭泽阳，乔辉，杨帆，
申请(专利权)人：中国石油化工股份有限公司石油勘探开发研究院，
类型：发明
国别省市：