基于生产动态数据的井间储层连通性识别方法

2019-02-21王任一

石油与天然气地质 2019年2期

王任一

(浙江海洋大学石化与能源工程学院,浙江舟山,316022)

目前注水开发仍是油田开采的最主要方式,而井间储层连通性认识,则是油田注水开发工程大厦的基石。井间储层连通性识别方法主要有3类： ① 数值模拟方法,所考虑的因素比较全面,所得结果比较直观,但由于考虑的因素较多,对所需资料的依赖性较大,尤其是对地质资料的要求较高,其预测结果可靠性,很大程度上取决于基于地质工程师小层对比及测井工程师储层参数解释结果所建立的三维地质模型可靠程度,并且周期长成本高； ② 油藏工程方法,主要有油水井生产动态资料相关分析法、试井法和示踪剂等方法。这些方法如果要实现定量化解释,一般都需建立一个理想化的复杂理论数学模型(与地下实际情况相差较大),还需要提供一个依赖于油藏静态认识成果的地质油藏参数(例如,小层对比成果和油藏参数解释结果等),一般均需要较长研究时间和较高实施成本费用； ③ 地球化学方法,主要有矿化度分析法和原油气相色谱指纹技术等。矿化度分析法,影响因素较多,不确定性较大。而气相色谱指纹技术也存在着一些自身难以克服的不足[1]。笔者提出一种利用油水井生产动态数据资料,通过多变量相空间重构方法,重构注采系统动力学特征,求取重构注采系统吸引子的关联维,可直接对井间储层连通性进行定量识别的方法。

1 方法原理

具有非线性渗流特征的低渗透油藏注水开发系统,是一个复杂的非线性多相流时变系统,油水井生产动态特征受多因素综合影响[2-3]。在现场生产上,一般简单通过注水井的注水量,与采油井的产液量(产水量或产油量等)之间的相关情况,可以大致判断油水井井间储层连通状况,也就是依据油水井之间哪个方向是增注增油,哪个方向是增注增水,哪个方向是既增油又增水,哪个方向是增注不增油等状况,为油田后续调整方案制定提供参考[4]。可见,在油水井注采系统的生产动态数据中,隐含着油水井井间储层连通性信息。

可以将油水井注采系统,简单视为由注水井、采油井和油水井井间储层,这三者组成的一个非线性多相流系统。多相流系统是一个非线性能量耗散动力学系统,它表现为一系列混沌运动特征[5]。根据Takens嵌入定理,系统中任一分量变化都是由与之相互作用着的其它分量所决定,单一变量时间序列隐含整个系统的运动规律[6]。可在不清楚影响这一注采系统各因素间相互作用机理的情况下,利用油水井生产动态数据,通过相空间重构方法,重构这一注采系统动力学特征,并与原注采系统动力学特征及吸引子拓扑结构完全相同,并可通过求取重构注采系统吸引子的关联维,来认识这一注采系统动力学特征[7-8]。因此,利用油水井生产动态数据,采用多变量相空间重构方法,就有可能从隐含着井间储层连通性的多变量生产动态数据中,提取出井间储层连通性信息。

1.1 注采系统的多变量相空间重构

依据Takens嵌入定理,对于无干扰噪音的足够长时间序列,只要时间延迟(τ)与嵌入维数(m)选择合适,单变量时间序列就足以重建注采系统的动力学特征。但在实际应用中,情况却不完全是这样,不能肯定基于特定观测变量一定可以正确地重构相空间。例如,对于Lorenz方程[9],基于z轴输出数据进行采样形成的时间序列,不能重构Lorenz系统动力学特征。因为系统x轴与y轴之间的对称关系无法反映到z轴数据列中。可见,对于任何给定的一维时间序列就可以重建动力学系统的看法并不总是正确[10]。由采油厂录取的油水井生产动态数据往往受多因素影响,数据波动很大,且生产动态数据长度有限。例如,每次修井或检泵周期的不确定性,或不同时间段注采压差等生产参数确定的人为经验性,或为了阶段性生产任务完成,人为性地加大措施强度等。如果仅仅通过油水井单变量生产动态数据(采油量或注水量),重构油水井注采系统动力学特征,有可能会产生一定误差。但是在多变量生产动态数据中,则包含了更多的反映注采系统动力学特征信息,采用多变量生产动态数据进行相空间重构,可以提高重构注采系统动力学特征的精度及可靠性。

作为注采系统输入端的注水井生产动态数据变量有注入压力和注水量等，作为注采系统输出端的采油井生产动态数据变量有生产压差、产油量、产水量和产液量等。目前多变量时间序列相空间重构方法,多为经典G-P算法推广[11-12],本文采用文献[13]的重构方法。假设油水井多变量生产动态数据为xij,i=1,2,…,r,表示有r个变量生产动态数据；j=1,2,…,n,n为生产数据长度(时间)。

(1)

对于多变量生产动态数据(xij),取时间延迟τ=max(τi),i=1,2,…,r。τ为满足Takens嵌入定理共同时间延迟量。假定选择适当的嵌入维数为m,可重建r×m阶矩阵为

(2)

上式(2)中,i=1,2,…,M;M=n-(m-1)τ,为相空间中相点数。当r=1时,公式(2)就变成单变量时间序列相空间重构公式。由于任意两个r×m阶矩阵之差,仍然是一个r×m阶矩阵,则用多变量时间序列联合计算关联积分为:

(3)

其中

(4)

(5)

依据公式(5),在实际计算关联维D2时,由于注采系统动力学特征的复杂性,一般使用局部斜率方法,也即在某些尺度区间上,局部曲线段相对平滑,对应斜率也相对稳定的,这一局部尺度范围区间也称为无标度区。在m≥2D+1(D为系统维数)时,使得重构相空间中的“轨线”与原系统在微分同胚意义下是等价的[14],由此重构油水井注采系统与原系统具有相同动力学特征及吸引子拓扑结构。因此,通过求取反映注采系统吸引子特征不变量的关联维,就可提取出油水井井间储层连通性复杂程度的信息。

1.2 延迟时间与嵌入维数确定

相空间重构的思想是将一维时间序列变换到多维相空间,使原系统几何不变量保持不变[15]。相空间重构关键是确定嵌入维数(m)和时间延迟量(τ)。如果时间延迟量太小,则xi和xi+τ很靠近,相图运动轨迹压缩在对角线附近信息冗余,反之如果时间延迟量太大,相图吸引子则发生折叠产生几何畸变[16]。对于时间延迟量(τ)确定方法很多,本文采用互信息法确定。

互信息可以度量两个变量间相互依赖程度或两个变量间共有信息。不仅可以反映变量间线性相关性,也可以反映变量间的非线性相关性。对于确定的两个时间序列变量s=xi和q=xi+τ,二者互信息表达式为:

I(s,q)=H(s)+H(q)-H(s,q)

(6)

式中:I(s,q)为s和q这二个变量间互信息,H(s)和H(q)分别为时间序列s和q的信息熵,H(s,q)为二者联合信息熵。当s与q完全相关时互信息值最大,当s与q完全不相关时互信息值为零,本文互信息计算采用等间距格子算法[17]。通常取第一个互信息极小值点所在位置作为最佳延迟时间(τ)[18]。图1为利用研究区采油井WR531井76个月产油量数据(以月为计量单位,以下相同),求得在不同延迟时间τ下,I(s,q)互信息值变化关系曲线。从图中可看出,相空间重构最佳延迟时间为3个月(τ=3)。

图1 WR油田WR531井月产油量生产数据延迟时间与互信息关系Fig.1 Correlation between delay time and mutual information of monthly oil production data of Well WR531

嵌入维数(m)是相空间重构的另一个关键参数,本文采用比较成熟的CAO方法来确定嵌入维数(m)[19-20]。对于时间数据列xi,i=1,2,…,n,在m维相空间中重构向量为xi(m)=xi,xi+τ,…,xi+(m-1)τ,它们都存在一个距离最近向量xN(i),N=1,2,…,n-mτ。令

(7)

定义a(i,m)的平均值为：

(8)

为了确定a(i,m)平均值E(m)随m变化情况,定义

(9)

如果时间序列具有非线性特征时,则m大于某个m0值时,E1(m)就趋于不变,相空间重构的最小嵌入维数就是m0。图2也是利用WR531井月产油量数据,按CAO方法计算求得的E1(m)值随嵌入维数(m)变化曲线图(τ=3)。从图2中可看出,当嵌入维数(m)大致大于15时,E1(m)基本趋于不变。也说明WR531井产油量数据具有很强非线性特征。

实际上根据经典G-P关联维算法[11],也可求得最小嵌入维数。对于非线性时间序列而言,一般随着嵌入维数m增大,lnε-lnC(ε)曲线簇变化差异变小。如图3所示,当m>15时,lnε-lnC(ε)曲线簇变化趋于重叠不变。稳定直线段区间(无标度区)的斜率,即为关联维数D2也趋于稳定值。这时所对应的嵌入维数(m),就是最小嵌入维数(m0)。图3也表明,利用WR531井月产油量数据,进行相空间重构合适嵌入维数为15(m=15,τ=3)。

图2 WR油田WR531井月产油量数据E1(m)值随嵌入维数(m)变化曲线Fig.2 Plot showing the variation of monthly oil production data of Well WR531 E1(m)value with the change of embedding dimensions(m)

图3 WR油田WR531井月产油量数据的lnε-ln C(ε)曲线图(τ=3个月)Fig.3 Plot showing the lnε-lnC(ε) of monthly oil production data of Well WR531(τ=3 months)

1.3 井间储层连通程度的确定

如果反映非线性动力学系统某一侧面的二变量时间序列为xa和xb,则可通过3种方式进行相空间重构,重建这一非线性动力学系统。与之相对应,也可分别求出3种不同方式所重构这一非线性动力学系统吸引子的关联维。也就是根据公式(5),可分别求出联合利用二变量时间序列(xa和xb)所重构动力学系统吸引子的关联维D2(xa,xb)、利用单变量时间序列(xa)所重构动力学系统吸引子的关联维D2(xa),以及利用单变量时间序列(xb)所重构动力学系统吸引子的关联维D2(xb)。文献[21]已经证明,如果D2(xa,xb)≈D2(xa)≈D2(xb),则变量xa和xb时间序列来源于同一非线性动力学系统。如果D2(xa,xb)=D2(xa)+D2(xb),则变量xa和xb时间序列来源于不同的非线性动力学系统,变量xa和xb属于非线性无关,彼此不能提供任何信息。

对于由水井、油井以及二者间未知连通性的储层,这三者组成注采系统而言,假定有二变量生产动态数据序列,例如注水井注水量(xa)和采油井产水量(xb)。则可通过3种方式进行相空间重构,重建这一注采系统动力学特征。与之相对应,也可分别求出这3种不同方式重构注采系统吸引子的关联维。如果利用注水井注水量单变量生产动态数据(xa),所重构注采系统吸引子的关联维为D2(xa)；利用采油井产水量单变量生产动态数据(xb),所重构注采系统吸引子的关联维为D2(xb)；联合利用注水井注水量和采油井产水量二变量生产动态数据(xa,xb),所重构注采系统吸引子的关联维为D2(xa,xb),则油水井井间储层连通性可通过如下方法确定：

1)D2(xa,xb)≈D2(xa)≈D2(xb)。说明注水井注水量(xa)和采油井产水量(xb)来自同一注采系统,二者具有很强的非线性相关性。3种不同方式所重构注采系统,具有相同动力学特征及吸引子拓扑结构。也说明油水井井间储层完全连通,具有同一个注采系统的动力学响应特征。

2)D2(xa,xb)≈D2(xa)+D2(xb)。说明注水井注水量(xa)和采油井产水量(xb)来自完全不同的注采系统,二者非线性无关,彼此间不能提供任何信息。3种不同方式所重构注采系统,具有完全不同的动力学特征及吸引子拓扑结构,没有任何相似性,油水井井间储层是完全不连通的,具有不同注采系统的动力学响应特征。

3) 对于介于上述两者之间的情况,可以通过定义归一化指标来确定

(10)

其中,M(xa,xb)=D2(xa)+D2(xb)-D2(xa,xb)。R称为井间储层连通性系数。当油水井间储层为完全连通时,R=1；当油水井间储层为完全不连通时,R=0；当连通性处于二者之间情况,R=1～0。这样就可以利用公式(10),对油水井间储层通性进行定量解释与评价。

2 应用实例

如图4所示,WR油田为一简单断背斜边水砂岩油藏,储层为一套扇三角洲前缘沉积体,含油层段为三间房组(J2s)。主力油层单一(J2s2(3)),平均孔隙度为16%,渗透率为33.4×10-3μm2。油田采用注水开发方式,共有油水井32口,其中油井20口,注水井12口,综合含水达到48.7%。在油田西南部,靠近南部断层高部位黑色虚线标注的长椭圆形区域,为实例研究区。将利用上述方法,对研究区内2口采油井(WR52和WR531井)与2口注水井(WR51和WR53)之间储层连通情况进行分析。原综合油藏数值模拟结果和生产动态资料,对这4口井井间储层连通性认识结论为：① 注水井WR53井和采油井WR531井之间不连通；② 注水井WR53井和采油井WR52井之间连通；③ 注水井WR51和采油WR52井之间不连通。下面利用上述新方法,对这3组油水井井间储层连通性进行分析。

图4 WR油田J2s2(3)油层注采井网构造图Fig.4 Structural map showing the injection-production pattern of J2s2(3) oil layer in WR oilfield

1) WR53井和WR531井井间储层连通性

从注水井WR53井注水量和采油井WR531井产水量随时间变化曲线可看出(图5a)： ① 在前31个月,WR53井注水量呈持续快速单边下降的趋势,而采油井WR531的产水量较为稳定,波动不大,但整体处于较低水平。WR53井注水量和采油井WR531的产水量间没有相关性； ② 在第31～51个月区间段,WR53井注水量整体仍呈继续快速下降的趋势,而采油井WR531井产水量也随之下降,二者呈正相关性； ③ 在第51个月时,在对WR53井进行提压增注措施之后,WR53井注水量得到明显恢复,WR531井产水量也随之相应增加,似乎井间储层具有一定的连通性。

图5b为WR53注水量与WR531井产水量lnε-lnC(ε)曲线图(说明：为了图中显示简洁明了,在lnε-lnC(ε)曲线簇叠合图中,仅显示τ=3,m=15时曲线,以下相同)。从图5b可看出,利用WR53井注水量、WR531井产水量,以及联合利用WR53注水量与WR531井产水量二变量生产动态数据,通过这3种不同生产动态数据相空间重构方法,所重构这一注采系统动力学特征完全一致。相空间重构lnε-lnC(ε)曲线非常相似,三者几乎粘合在一起,所重构注采系统吸引子的关联维,即在无标度区间三者的斜率基本一致,在tan(45°)左右,即D2(xa)≈D2(xb)≈D2(xa,xb)≈1。因此,无论利用WR53井注水量、WR531井产水量单变量生产动态数据,还是联合利用WR53井注水量和WR531井产水量二变量生产动态数据,所重构注采系统吸引子均具有相同的拓扑结构,属于同一注采系统的动力学响应特征。根据公式(10),WR53井和WR531井井间储层连通性系数R为1,连通性很好,而原认识结论却是不连通的。

图5 WR油田WR53井和WR531井井间储层连通关系Fig.5 Analysis of the interwell reservoir connectivity between Well WR53 and Well WR531a.WR53井注水量与WR531井产水量随时间变化曲线；b.WR53井注水量和WR531井产水量的lnε-lnC(ε)曲线

2) WR53井和WR52井井间储层连通性

从注水井WR53井注水量和采油井WR52井产水量随时间变化曲线可看出(图6a)：① 在前44个月,WR53井注水量呈持续快速单边下降的趋势,而WR52井产水量整体却呈显著上升趋势,二者呈明显负相关性；② 在第44～51个月区间段,WR53井注水量很低,注入能力大幅下降。而采油井WR52产水量也随之明显下降；③ 在注水第51个月时,在对WR53井进行提压增注措施之后,注水井WR53井注水量得到明显恢复,并与采油井WR52产水量间整体呈较好正相关性,显示井间储层具有较好的连通性。

图6b为WR53注水量和WR52井产水量的lnε-lnC(ε)曲线图,从图6b可看出,利用WR53注水井注水量、WR52采油井产水量,以及联合WR53注水量和WR52井产水量二变量生产动态数据,通过这3种不同生产动态数据相空间重构方法,所重构这一注采系统动力学特征完全一致,相空间重构lnε-lnC(ε)曲线也非常相似,三者近乎平行,所重构注采系统吸引子的关联维,即在无标度区间三者的斜率基本一致,在tan 42°左右,也即D2(xa)≈D2(xb)≈D2(xa,xb)≈0.9。

图6 WR油田WR53井和WR52井井间储层连通关系分析Fig.6 Analysis of the interwell reservoir connectivity between Well WR53 and Well WR52a.WR53井注水量与WR52井产水量随时间变化曲线；b.WR53井注水量和WR52井产水量的lnε-lnC(ε)关系曲线

因此,无论利用WR53井注水量、WR52油井产水量单变量生产动态数据,还是联合利用WR53井注水量和WR52井产水量二变量生产动态数据,所重构注采系统吸引子均具有相同的拓扑结构,属于同一注采系统动力学响应特征。根据公式(10),WR53井和WR52井井间储层连通性系数R为1,井间连通性很好,并与原认识结论一致。

另外,由于实例(2)中所重构注采系统吸引子的关联维为0.9,要小于前面实例(1)关联维1.0,也说明井间储层连通性的复杂程度可能要略小于前者。

3) WR51井和WR52井井间储层连通性

从注水井WR51井注水量和采油井WR52井产水量随时间变化曲线可看出(图7a),在前6个月,WR51井注水量处于高位运行,第7个月之后,WR51井注水量大幅下降,处于低位运行,并且整体一直保持相对稳定不变状态。而采油井WR52井产水量整体却呈明显逐渐上升趋势,二者间没有任何相关性。

图7 WR油田WR51井和WR52井井间储层连通关系Fig.7 Analysis of the interwell reservoir connectivity between Well WR51 and Well WR52a.WR51井注水量与WR52井产水量随时间变化曲线; b.WR51井注水量和WR52井产水量的lnε-lnC(ε)关系曲线

图7b为WR51注水量和WR52井产水量的lnε-lnC(ε)曲线图,从图7b可看出,利用WR51注水井注水量、WR52采油井产水量,以及联合WR51注水量与WR52井产水量二变量生产动态数据,通过这3种不同生产动态数据相空间重构方法,所重构这一注采系统动力学特征完全不同,相空间重构的lnε-lnC(ε)曲线间相差很大,三者各不相同。利用注水井WR51井注水量单变量生产动态数据,所重构注采系统吸引子的关联维,即无标度区间斜率D2(xa)为1.74;利用采油井WR52井产水量单变量生产动态数据,所重构注采系统征吸引子的关联维,即无标度区间斜率D2(xb)为0.27；联合利用WR51注水量与WR52井产水量二变量生产动态数据,所重构注采系统吸引子的关联维,即无标度区间斜率D2(xa,xb)为2.04。由于D2(xa)+D2(xb)≈D2(xa,xb)(1.74+0.27≈2.04),因此,无论利用WR51井注水量、WR52井产水量单变量生产动态数据,还是联合利用WR53井注水量和WR52油井产水量二变量生产动态数据,所重构注采系统吸引子均具有完全不同的拓扑结构,不属于同一注采系统的动力学响应特征。根据公式(10),WR51井和WR52井井间储层连通性系数R接近零(0.014),二者为非线性不相关,彼此间不能提供任何信息,井间储层完全不连通,与原认识结论一致。