文摘

阐明储层流动单元的动态演化特征在注水开发过程中,上层的新第三纪块M孤岛油田馆陶组地层是作为一个案例研究。基于测井资料、水注入剖面测试数据,subwell数据,和生产性能数据,其中,流动带指数(FZI静态)提出了流动单元的静态评价参数。累积注水(W之间的关系T)和FZI变化( )原来的样貌。因此, 由于注水结合流动单元的静态参数评价(FZI静态)流动单元的动态参数评价(FZI动态)。流的全面评估现实单位在不同时期的特点是公式 研究表明,流动单元的划分标准,FZI具有良好的相关性与吸水剖面的测试结果,如吸水强度和相对注水体积。使用FZI流动单元的静态参数评价,四种类型的流动单元划分如下:第一类流动单元, ;II型流动单元, ;类型III流动单元, ;IV型流动单元, 不同的流动单元的储层孔隙度和渗透率特征是高度相关的。此外,相对渗透率曲线证实,不同的流动单元有不同的渗流能力。虽然比较特征储层参数在不同时期,储层的物性参数更有利于流体流动与注水开发。增加相同的注水速度I和II型流动单元大于III和IV型流动单元。此外,当I型流动单元在大面积连续分布,高耗水量乐队成立,吸收大部分的水注水井注水。因此,注水效率很低。不同的流动单元的变化主要是由注入控制生产井网和WT。结合层间流动单元的相对变化特征,这些变化可分为增加类型变化和减少类型变化。最后,根据不同的流动单元的分布特征,含油饱和度和注水效果,策略,挖掘剩余油的潜力从三个方面(平面、层间和内层)制定。

1。介绍

储层流动单元研究是油藏描述的关键在石油和天然气领域的发展。储层流动单元的概念赫恩等人提出的类似油藏储层,在纵向和横向方向是连续的,和有一个内部对流体流动的影响1]。自那时以来,学者们对储层流动单元进行了广泛的研究。他们指导以后发展和残余油分布预测油气田的流动单元的分配机制。因此,他们已经逐渐形成一个共同的理解流单元的两个元素:(1)渗流特征存在于单个流单元相似,而不同的流动单元有不同的渗流特征;(2)有明显的渗流屏障接口和渗流差异流动单元之间的接口。

目前,最常用的流动单元划分方法包括流动带指数法(FZI) [2- - - - - -6),井间流量指数法(IFCI) [7,8),修改地层洛伦兹图方法(SMLP) (9),模糊聚类方法(10- - - - - -21),灰色系统理论方法(22,23),而非齐次综合指数法(24- - - - - -27]。这些方法大多是基于核心的渗透率和储层的形成。静态参数(如孔隙度、渗透率、泥内容,喉道半径、粒度中值、和水淹层解释结果)从测井和地震数据是用来划分流动单元。受到静态数据的瞬时性,获得的流量单位只反映了流体在储层渗流能力的特点在钻井数据采集。然而,在长期注水开发的砂岩储层,孔隙结构、渗透、泥沙含量,和其他特征的储层变化会引起渗流油藏连通性差(28- - - - - -35];因此,流动单元的类型也发生了改变。这个流单元的动态特性很难被评估使用上面的参数得到的静态数据。油藏描述的过程贯穿勘探和开发;从发现到最后放弃,它是反复进行的阶段和卷36]。流单位应该反映储层中的流体渗流特征的变化在不同时期,具有重要指导作用的油田开发调整措施,预测剩余油分布,和部署加密井。在当前的研究成果,许多学者的流动单元类型分为多个区域,以更好地了解水库之间的区别。然而,内部结构特性的变化引起的储层注水开发时往往忽略划分流动单元。然而,仍然有小的变化研究流动单元的类型。因此,储层流动单元的动态特性在不同时期的油田开发应该研究。

本研究认为上层部分馆陶组河流相储层的形成区M的孤岛油田为例,利用巨大的静态和动态数据在油田开发过程中积累起来的。之间的关系 并建立了注水体积。研究方法结合动态和静态储层流动单元的变化探讨了促进从储层流动单元的瞬态特征转变为现实的储层流动单元的特征。此外,储层流动单元的动态演化规律在不同时期在油田开发过程中澄清是指导油田开发调整措施和剩余油分布的预测。

2。地质背景

区M的孤岛油田位于孤岛低隆起旸萧条(数字1(一)1 (b))。上部的上第三系馆陶组地层是一组重要的工作区域(图中含油层系1 (c))。它主要由蜿蜒的河流相沉积和发展工微相,如河床滞留沉积、边滩、废弃河道、天然堤。此外,岩性主要为层间的中细砂岩、粉砂岩与灰色绿色和紫色红色泥岩。储层具有良好的物理性能,25.1% - -35.6%的孔隙度,渗透率 ,和储层非均质性强。上层成员馆陶组地层划分为五个砂组(从下到上Ng5-Ng1),其中Ng3和Ng4砂组是主要的含油层系。此外,NG3砂组进一步分为5个小层(Ng35-Ng31从下到上),和Ng4砂组进一步分为四个小层(Ng44-Ng41从下到上)。Ng35子层和Ng44子层是主要的含油储层,厚度约为7 ~ 13米。层间稳定发展,1 - 5米的厚度。Ng35子层和Ng44磷层是独立开发单位(37,38)(图1 (c))。

在6.4公里有127口井2研究的区域。自1976年以来,工作区域实现了两层的结合生产一个反向nine-point-well模式。在1983年,它最初是强化。然而,自1990年以来,它已经被强化,相反的行和列好模式。上部和下部的注塑生产流线角度系统模式大大不同。上层系统采用30°北偏西行井网,和较低的系统采用10°北偏东行井网。经过40多年的发展,已经进入了超高含水开发的后期阶段。2017年,研究人员进行了层系列井网交流,与整体优化改变了40°,关于大型注水分布差异和不均匀在该地区原油生产。交换后,某些发展成就。然而,问题依然存在,如不清楚储层特征和不确定极端水消费领域。 Therefore, it is important to study the change of reservoir flow unit under long-term water injection development in the study area. After more than 40 years of development, the well pattern density in the study area is large (the average well spacing is less than 200 m), and the logging data is complete. Moreover, a wealth of continuous time dynamic data has been accumulated, such as injection pressure, instantaneous water injection, cumulative water injection, and cumulative liquid production that have been accumulated. There are 18 wells with water injection profiles of test data in different periods and eight groups of subwells in different periods (well spacing is about 10–48 m). These data can provide data conditions for studying the dynamic change characteristics of flow units during water injection development.

3所示。方法

阐明储层变化在不同时期的特点,本研究提出以FZI (FZI静态)流动单元的静态评价参数基于测井数据,水注入剖面测试数据,subwell数据和生产性能数据。之间的关系T和FZI的变化( )在不同时期通过subwells安装钻。的 注水诱发的成交量结合流动单元的静态参数(FZI静态)评价的动态参数(FZI动态的)流动单元评估,两相结合实现综合评价,全面评估流动单元在不同时期的现实。

3.1。动态和静态评价流动单元

流动单元、储层流体流动能力的标准,改变随着时间的推移,注水,和其他因素。因此,在这项研究中,流动单元的参数分为静态和动态参数综合评价流动单元在不同时期。

3.1.1。流动单元的静态参数评估

FZI的方法确定了不同岩相的孔隙几何形状复杂变化通过核心数据。然后实现流量的分类单位根据储层孔隙度和渗透率的特点。这种方法具有广泛的实用性。与灰色系统理论法和模糊聚类法,必须全面的注水解释结果,这种方法最优反映了储层的瞬时特征。通过前面的实验分析和研究39- - - - - -43),注水开发引起的储层特征的变化可以直观地反映在储层孔隙度和渗透率的变化。该参数的选择提供了一个基础研究后续的动态参数流单元。流动单元的静态参数在本研究中采用了流动带指数(FZI)提出的“王et al。19]。公式(1)显示了计算公式。

公式:(我)FZI(静态)。流动带指数计算了静态参数,无量纲(2)K。渗透率(×103μ2)(3) 有效孔隙度(%)

由于随机误差的存在指数划分流动单元的流区,FZI相同的流动单元分布在其真正的平均值是直线段FZI的累积概率图。当有多个异构流量单位,FZI正态分布的叠加数正态分布。因此,它被表示为概率图上的多个直线段(22]。FZI概率累积曲线为研究区,流动和四个部门单位提供(图2):我流的单元类型, ;II型流动单元, ;类型III流动单元, ;IV型流动单元, I型和II流单元以细砂岩为主,第三类型流动单元是由细粉砂岩,和IV型流动单元包含细沙以粉砂岩为主。不同类型的流动单元有不同的岩性类型。

根据孔隙度渗透率曲线不同类型的流动单元,不同的流动单元的核心孔隙度和渗透率特征有很高的相关系数(图3)。因此,这种方法可以有效地描述复杂的高孔隙度、高渗透率储层特征。

结合储层流动单元的类型被流动带指数取心井被分成流动单元类型。自底向上的流动单元的质量逐渐变化后通过岩心分析孔隙度和渗透率测试的特点(图4)。流动单元的水平越高,驱油效率越高。因此,流动单元的类型确定基于FZI指数可以有效地反映储层特征的程度。

结合分析的相对渗透率曲线的特点划分流动单元类型、相对渗透率端点和曲线形状的不同类型的流量单位是不同的(图5)。此外,流动单元的划分可以反映储层渗流特征的差异。从我流单元IV型流动单元类型,束缚水饱和度和残余油增加,和两阶段的时间间隔coinfiltration区降低(表1)。

3.1.2。流动单元的动态参数评估

流动单元的动态评估的困难选择动态参数。通过比较分析和测试数据超过十取心井的研究区域内的三个时期,可以看出随着注水的发展,泥质含量,和碳酸盐含量减少,稍有增加,孔隙度和渗透率大大变化(表2)。

然而,许多先前的实验结果表明,在注水开发过程中,水库水量的增加提高了储层的孔隙度和渗透率(39- - - - - -43]。因此,WT水库诱发的水库在一定程度上的变化。计算WT每个流单元和建立之间的关系T数量和流动单元的特征是评估的关键流的动态参数研究单位。

使用FZI流动单元的静态参数评估,以达到注水分割基于流动单元分类、储层吸水能力之间的相关性和FZI决心通过水注入剖面测试数据。此外,FZI引入油藏工程方法实现油藏的注水分裂。通过 类似的储层在不同时期对subwells钻,W之间的关系T 原来的样貌。此外, 注水引起的静态参数结合流动单元评价(FZI静态)流动单元的动态参数评估(FZI动态;公式(2))。

自从WT数据是时间连续动态数据,从理论上讲,(FZI动态)可以计算在不同时期的表现 研究的动态变化特征流动单元在注水开发过程中。

(1)注水分割方法。动态数据,每个流单元的水流的水注入井和石油生产井测量是很困难的。因此,它不能直接用于计算的计算累积注水。获得每个流的累积注水数据单位,应该研究注水分裂方法计算每个流单元的内部水流和获得累积注水。

储层流动单元特征的影响不能被忽略在注水分裂。此外,相对注水和吸水强度动态监测数据反映了不同储层的吸水能力和储层流动单元的不同部分。根据81年的相关分析结果吸水剖面测试数据点FZI和水的吸收强度和相对18井注水,吸水强度之间的相关系数单位厚度和FZI为0.7608(图6(一))。此外,之间的相关系数相对单位厚度和注水FZI为0.845(图6 (b))。这表明FZI的流动单元划分的基础,反映了储层的吸水能力在某种程度上。

为了阐明水流水注水井之间的关系和石油生产井,根据渗流力学和油藏工程理论,渗流阻力的歧视模型的约束下注射生产关系是通过达西公式和水电相似原理建立了。

根据渗流阻力的大小,每喷射水流生产关系等参数进行分析,参照井间流动单元类型、有效厚度、流动带指数和井距。全面概括为三种模式:飞机系列,平面平行,垂直多层并行(图7)。

步骤1计算层的渗流阻力,我的水注入井 受影响的。以注水井为中心,计算每个小层的渗流阻力系数方向油井周围的水注入井:

步骤2我计算的总平面渗流阻力层的水注入井(水电相似原理):

步骤3计算注水的比例系数 - - - - - -水注入井的层 - - - - - -th油井的影响:

步骤4计算注水的总电阻:

步骤5比例计算分解水注入我的水注入井层系数:

步骤6计算从注水井注水体积 - - - - - -th层:

步骤7注水的 - - - - - -水注入井的层 - - - - - -油井th有效: 的公式:(我) 两个点之间的电阻(m3/不)(2) 两个点之间的间距(米)(3) 流动单元的厚度(m)(iv) 流动带指数(v) 总平面层的渗流阻力 水的注入井(vi) 比例系数的注水 - - - - - -水注入井的层 - - - - - -th有效油井(七) 总阻力的水注入井(八) 比例系数的分解水注入层 水的注入井(第九) 水从水注入井层体积 (x) 注入水的数量 - - - - - -有效的油井 - - - - - -th层

根据这个注水分裂方法,每个注水井的注水体积的方向周围油井和累积注水体积(WT每个流单元可以计算的)。注水体积分割不同的流动单元的每个方向可以通过结合FZI之间的相关性和吸水能力确认。

(2)计算流动带指数的变化 )。从上面的分析, 为确定动态参数是关键参数(FZI动力学)流的单位评估。subwells,井网相对较近,和原来的储层流动单元类型有类似的特征。假设两个井的钻井时间不同,后来的钻井数据可以用来评估流单元第一钻井和注水开发后变化特征(图8(一个))。此外,后来的钻井数据可用于计算 上节中有八双subwells Guanguan的研究领域。钻井时差范围从39到367个月,和WT差异范围从 通过W之间的相关分析T 在subwells,观察,WT不同类型的流动与单位有很好的相关性 (图8 (b))。其中,I型和II流单位有更好的储层特征,和注水体积越大,越大FZI的增加。相反,III和IV型流动单元较低增加FZI注水体积的增加。

在上面的分析中,W之间的拟合关系T四种类型的单位和流动 可获得双subwells之间。因此, 可以估计W的大小T在流动单元。因此,流动单元的动态参数评价(FZI动态)可以确定。

此外,通过上面的拟合公式,FZI流动单元在不同时期的预测。通过计算FZI postdrilling验证,符合率达到88.3%。因此,这种方法可以预测储层流动单元的变化引起的注水开发。

4所示。结果

4.1。流动单元评价

9显示的结合动态和静态评价储层流动单元可以使用FZI进行静态和FZI动态参数实现从流的瞬时特征表示单位转换到现实的特征表示。

9(一个)的计划是1983年Ng35层流动单元分布,用FZI静态作为评价参数,而图9 (b)的计划是2021年Ng35层流动单元分布,用FZI动态作为评价参数。图9表明,流动单元的平面分布评价FZI静态和FZI动态参数是不同的。相反,I型和II流单元的范围在图9 (b)明显更大。这表明,注水开发有积极的趋势改善流动单元的类型,它与上面的认可。

在不同的时间,有653注水资料从44井测井资料的研究领域。总体错误率低于8%,验证计算结果的动态流动单元使用相对单井注水在不同的时间。因此,相对变化量注水的注水剖面测试数据证实了动态静态相结合的合理性评价流动单元(图10)。随着水的不断发展注入,FZI玫瑰,和相应的注水相对上升。然而, 不同类型的流量单位是不同的, I和II型流动单元大大增加。

4.2。流动单元在注水开发的变化特征
4.2.1。准备周期性划分流动单元进化

流的变化特征单元累积注水是密切相关的。他们受到的改进模式和注水开发的过程。根据井网加密调整的特点和大规模的注水开发调整措施的研究区域,Ng35层的开发过程分为五个阶段(表3)研究的动态变化特征在注水开发过程中储层流动单元的生产特性曲线Ng35层(图11),五个生产周期明显不同特征,代表流动单元的变化特点进行了讨论。

4.2.2。流动单元的平面变化特征

FZI(静态)和FZI(动态的)被用来评估Ng35层的流动单元在五个时期,画的流动单元平面分布图5期(图12),研究了动态变化在注水开发过程中储层流动单元的特征。

(1)在不同时期流动单元的平面分布特征。阶段1:流动单元主要是II型的类型,类型我本地流量单位开发和IV型流单位大多是在南方(图12(a))。这一时期是油田生产的早期阶段,较低的水注入,石油产量高,平均含水率低于40%,初始储层流动单元的特点反映在一定程度的研究区域。(我)第二阶段。I型流动单元的比例增加(图12(b)),水积累层( );液体生产大大增加,含水率显著增加的形成。流动单元类型的变化在这个时期主要是流动单元分布特征的变化引起的井网加密调整(2)期3。WT是关于 ;含水率大于90%。受到大量的注水开发研究西北地区,当地II型流单位变成了I型流单元(图12(c))(3)期4。每月注水降低;WT是关于 开发时间长、WT是大。研究区进入高含水期。受长期注水,I型流动单元的分布格局逐渐符合水注入井的方向(图网络12(d))。一些类型III和IV型流动单元在东南注入水注入井的面积,导致大的水流,增加了FZI,增加流动单元的水平(iv)期5。在此期间,井网优化调整,精简的角度改变了大约40°,水注水井的分布位置大大改变,和注水浓度的变化导致I型流动单元分布形状的变化受井网流线的变化(图12(e))

(2)控制的平面变化因素流单元。经过很长时间的注水开发,I型和II型流动单元的比例增加,流动单元的分布趋势是大致相同的生产井。影响注水,流动单元的类型生产水注水井位置大大改变,大部分变成了I型流动单元。井网交换后,流动单元的分布趋势与流线型的变化改变了生产井排方向。IV型流动单元的总体变化很小,说明低级流单元的渗流能力弱,储层特征的变化是小注入过程中生产的发展。同时,由于注水波和它的影响,中部地区的流动单元类型相对大大改变了。I型流单位显著增加与最初的发展阶段,形成了主要在飞机上高水消费区。这个地区大多数注水,吸收液体高生产、高含水量和其他地区的水驱效率低。不同类型的流动单元的变化主要是由注入控制生产井网和累积注水。

4.3。层间流动单元的相对变化特征

小层之间的stable-developed夹层确保每一层的独立开采。在注入和生产相结合,不同的流动单元层间吸水不均匀引起的。结合流动单元的变化特点,分析了层间干扰的变化,主要表现在两个特征。

4.3.1。增加类型改变

这些特征主要反映在合并后的注入和生产井。水库的水吸收强度指数FZI相关;即吸水能力指数流单元的高流区更大比指数流单元的低流区。与长期注入和生产相结合,水库之间的差异进一步扩大,增加吸水强度的差异,层间干扰增强。油井NG35和NG44开发一组稠油层(图13)于1988年钻。他们最初是用作生产井,然后转化为水注水井。在水里注入剖面测试在2000年进行,两层的相对注水比为1.77,和水吸收层之间的差别很小。在测试在2006年和2015年,相对注水比为4.26和13.2,分别。长期注水开发提高了储层非均质性和形成越来越层间差异。

4.3.2。振幅降低类型变化

这些特点主要是减少造成的层间差异单层发展。例如,穿孔的位置是在小的上层。较低的小层没有穿孔,诱导FZI的增加和减少FZI的区别和较低的小层较低的底层流单元的小上层的影响下注水开发。随后,体积相对注水层之间的差异在随后的注入和生产相结合的过程变得更小。油井钻探时在1991年,只有上层NG35穿孔在最初的阶段(图14),低于Ng44 FZI。2008年,WTNg35是 ,Ng44层不主注塑生产线,水流很小,和FZI的增加远远低于Ng35,随后,well-patched Ng44层。结合后注水剖面的分析测试结果,相对注水层之间的差异很小,由于减少层间流动带指数的差异。这个特性是amplitude-reducing层间差异引起的不完美的注塑生产关系。

一般来说,长期注水开发将增加FZI的水库。注水体积的差异会导致不同的变化FZI的水库,然后影响油藏的注水开发。

5。讨论

5.1。生产问题流单元的控制之下

这些特点主要是减少造成的层间差异单层发展。例如(44,45),穿孔的位置是在小的上层。较低的小层没有穿孔,诱导FZI的增加和减少FZI的区别和较低的小层较低的底层流单元的小上层的影响下注水开发。随后,体积相对注水层之间的差异在随后的注入和生产相结合的过程变得更小。油井钻探时在1991年,只有上层NG35穿孔在最初的阶段(图14),低于Ng44 FZI。2008年,WTNg35是 ,Ng44层不主注塑生产线,水流很小,和FZI的增加远远低于Ng35,随后,well-patched Ng44层。结合后注水剖面的分析测试结果,相对注水层之间的差异很小,由于减少层间流动带指数的差异。这个特性是amplitude-reducing层间差异引起的不完美的注塑生产关系。

注射生产流线型的变化角度实现2017年研究区(37)可以减轻缺陷,如水库引起的低生产效率差异,在一定程度上提高原油采收率。然而,由于流动单元类型的变化,研究成果和注水剖面的测试结果,仍有许多注采井网流线的变化,和注射生产关系是位于高层流动单元。当地水吸收太严重了。根据生产数据的综合分析,含油饱和度的计划,和动态流动单元分布特征(图15),井网地区整体含油饱和度研究区北部的很低,,平面非均质性严重。几个水注水井I型流动单位,唯一的油井附近的油井是我流动单元类型。自井网交换,平均每日液体生产122.9米3远高于其他油井。I型流动单元的分布诱导形成的主要渗流通道之间,附近的水注水井。这诱导低注入和附近的油井的生产效率和困难在挖掘剩余油的潜力。

整个研究区高含水。然而,中间部分仍然有含油饱和度高,表明流动单元层的差异也很大,和异质性很强。类型III、IV和一些II型流单位仍有很高的储备。因此,挖掘剩余油的潜力的原因在高含水油藏井网调整,建立一个新的注射生产关系,并调整飞机的发展模式,一个intralayer,层间水以避免极端消费的主要渗流通道。

5.2。剩余油分布的控制下流动单元及其潜在的开发策略

根据流动单元的分布特征,含油饱和度,和注水的结果在研究区,利用策略的剩余油潜力飞机(46,47)、夹层和intralayer分布特征的基础上提出了动态流动单元(图16)。

5.2.1。潜在的平面剩余油开发策略

大多数平面剩余油的形成是由于水驱效率的不平衡引起的飞机不同的储层流动单元。目前井网状态下,中间两个水井和一个油井位于I型流动单元的分布范围内,形成一个有利的注塑生产渠道稳定力量。由于流动单元类型的差异,为其他油井水流很小,且位移参差不齐。因此,根据动态流动单元的分布特征,“高注入和低生产”的策略可以采用;即I型流动单元井作为主要的水注水井,和其他类型的流量单位是用作生产井。这个井网功能可以避免非常优势渗流通道的形成在一定程度上,提高注入和生产效率,挖掘剩余油的潜力。

5.2.2。潜在的层间剩余油开发策略

在生产相结合,层间流动单元之间存在较大的差异,导致很小的整体相对吸水和驱替效率低。这个特性常常诱导极富残余含油量和低注水学位地平线与小相对吸水。这个特性,“高堵塞和低注入”的方法可以用于生产,也就是说,阻塞的洞地平线的吸水,专注于利用地平线与相对较小的水吸收。以前的研究证实,注水的发展增加了FZI和水水库的吸收能力。此外,当FZIes的两层之间的差异很小,联合注入和生产可以提高驱油效率。

5.2.3。开发策略形成的剩余油潜力

层的剩余油主要是水库的节奏特点,导致高FZI的底部砂体,严重注水,注水水平低或没有注水砂体的顶部,和相对浓缩剩余的石油。关于这个特性,剩余油潜力可以挖掘使用集中生产的水平井钻井砂体顶部或压裂砂体的顶部增加储层的流体流动能力。

6。结论

(1)使用FZI流动单元的静态参数评价,四种类型的流动单元划分:I型流动单元, ;II型流动单元, ;类型III流动单元, ;IV型流动单元, 储层孔隙度和渗透率的特点,不同类型的流量单位是高度相关的(2)实现注水分割基于流动单元分类、储层吸水能力之间的相关性和FZI决心通过水注入剖面测试数据。随后,流动带指数引入油藏工程方法来划分储层注水。W之间的关系T 是通过双井钻安装在不同的时期。的 造成注水成交量结合流动单元的静态参数评价(FZI静态)流动单元的动态参数评价(FZI动态)(3)储层流动单元在不同时期的演变结果表明,注水开发可以提高流动带指数,但增加注水I和II型流单位远远大于III和IV型流动单元。因此,注水开发将进一步提高储层的非均质性。不同类型的流动单元的变化主要是由注采井网控制和WT。两种层间流动的相对变化特征单元进行了总结(4)根据流动单元的分布特征,含油饱和度,和注水的结果在研究区,四种剩余油潜力开发策略在一个平面上,一个intralayer,和一个层间提出了:在平面上(高注入和低生产)、层间(高堵塞和低注入),层间(水平井生产),压裂

数据可用性

使用的数据来支持本研究的结果包括在本文中。

的利益冲突

作者宣称没有利益冲突。

确认

这项工作是由黑龙江省自然科学基金资助(没有。YQ2019D002)和中国国家自然科学基金(41772149和41772149号)。作者感谢东北石油大学提供的指导和中石化胜利油田公司。