井筒失效数据采集分析方法及其装置、存储介质与流程

1.本发明涉及石油开发技术，尤其涉及一种井筒失效数据采集分析方法及其装置、存储介质。


背景技术：

2.气井的油管作为地层高温高压气体产出通道，所处环境和条件极为复杂，在应力、高温、高压和流体的腐蚀作用下，会发生完整性失效，如丝扣泄漏、腐蚀穿孔和应力开裂、乃至断裂等。油管完整性失效后，油管内的高温高压气体进入油管和套管之间的空间，并上浮聚集在井口处，进而改变了原有的井筒内部流体分布，引起井口高压，即引起持续环空压力。当井口的压力释放后，在压差的作用下气体会沿着气窜通道重新补充进环形空间，增加了生产成本和管理风险。
3.因此，有必要分析油管完整性失效后井筒的压力变化规律和流体分布特征。


技术实现要素：

4.本发明提供一种井筒失效数据采集分析方法及其装置、存储介质，用以解决分析油管完整性失效后井筒的压力变化规律和流体分布特征的问题。
5.一方面，本发明提供一种井筒失效数据采集分析方法，包括：
6.获取井筒中的油管和套管所形成的环形空间的流体温度变化平均值；
7.获取所述环形空间的环空压力；
8.检测所述环形空间的参数数据组，所述参数数据组包括所述环形空间的液柱初始体积、所述环形空间的液体等压膨胀系数、所述环形空间的横截面积、所述环形空间的液体等温压缩系数、所述环形空间的长度；
9.根据所述参数数据组、所述流体温度变化平均值和所述环空压力，确定所述环形空间的液柱长度和气柱长度；
10.获取油管泄漏点的内外侧压力比；
11.若所述内外侧压力比未达到预设压力比，则返回执行步骤所述获取所述井筒中油管和套管之间的流体温度变化平均值；
12.若所述内外侧压力比达到所述预设压力比，则根据获取到的多组所述环形空间的液柱长度和气柱长度，分析井筒压力变化规律和流体分布特征。
13.其中一项实施例中，所述根据所述参数数据组、所述流体温度变化平均值和所述环空压力，确定所述环形空间的液柱长度和气柱长度，包括：
14.根据所述液柱初始体积、所述流体温度变化平均值、所述液体等压膨胀系数、所述环空压力、所述液体等温压缩系数和所述横截面积确定所述液柱长度；
15.根据所述环形空间的长度和所述液柱长度的差值确定所述气柱长度。
16.其中一项实施例中，所述根据所述液柱初始体积、所述流体温度变化平均值、所述液体等压膨胀系数、所述环空压力、所述液体等温压缩系数和所述横截面积确定所述液柱
长度，包括：
17.根据所述液柱初始体积、所述流体温度变化平均值、所述液体等压膨胀系数、所述环空压力和所述液体等温压缩系数确定第一参数；
18.根据所述第一参数和所述横截面积的比值确定所述液柱长度。
19.其中一项实施例中，所述根据所述第一参数和所述横截面积的比值确定所述液柱长度包括：
20.根据公式确定所述液柱长度，其中，hy代表所述环形空间的液柱长度；v
al
代表所述环形空间的液柱初始体积；δt
al
代表所述流体温度变化平均值；α
p
代表所述液体等压膨胀系数；p
an
代表所述环空压力；k
t
代表所述液体等温压缩系数；aa代表所述环形空间的横截面积。
21.其中一项实施例中，所述获取井筒中的油管和套管所形成的环形空间的流体温度变化平均值，包括：
22.自井底至井口将所述油管等间距划分为i个分段，并用井筒分段编号值标记每个分段；
23.根据公式计算所述流体温度变化平均值；
24.其中，δt
al
代表所述流体温度变化平均值；代表第i个分段对应的环形空间的流体温度；代表第i个分段对应的油管的初始地层温度；n代表所述环形空间的流体终止深度对应的分段井筒编号，m代表所述环形空间的流体起始深度对应的井筒分段编号。
25.其中一项实施例中，满足公式：
[0026][0027]
其中，代表所述第i个分段对应的环形空间的流体温度；代表第i个分段对应的油管内的气体温度；代表第i个分段对应的油管的油管中心到井筒外侧的径向热阻；代表所述第i个分段对应的油管的初始地层温度；td代表无因次地层温度；λe代表地层导热系数；代表第i个分段对应的油管的油管中心至油管外侧的径向热阻；代表第i个分段对应的油管的油管中心至套管内侧的径向热阻；π代表圆周率。
[0028]
其中一项实施例中，所述获取所述环形空间的环空压力包括：
[0029]
根据公式v
ag
v
al
δv
al
＝v
an
计算所述环形空间的环空压力，其中，v
an
代表所述环形空间的总体积；其中，δv
al
代表所述环形空间的液柱体积变化量；v
al
代表所述环形空间的液柱初始体积；v
ag
代表所述环形空间的气柱体积；
[0030]
其中，v
ag
满足公式：
[0031]
n代表所述环形空间的流体终止深度对应的分段井筒编号，m代表所述环形空间的流体起始深度对应的井筒分段编号；v
gx
代表第x个时间步长时气体自所述油管泄漏入所述
环形空间的总体积；v
gx-1
代表第x-1个时间步长时气体自所述油管泄漏入所述环形空间的总体积；ps代表标况压力；p
al
代表气体泄漏点处的外侧压力；ts代表标况温度；t
ai
代表所述环形空间的流体温度；
[0032]vgx
满足公式：
[0033]
其中，q
lgx-1
代表第x-1个时间步长时油管内的气体在完整性失效点的泄漏速率；v
gx
代表第x个时间步长时气体自所述油管泄漏入所述环形空间的总体积；v
gx-1
代表第x-1个时间步长时气体自所述油管泄漏入所述环形空间的总体积；δt代表时间步长；
[0034]qlg
满足公式：
[0035][0036]
cre满足公式：
[0037][0038]
p
al
满足公式：
[0039][0040]
其中，
[0041]qlg
代表所述油管内的气体在完整性失效点的泄漏速率；co代表流量系数；p
fl
代表气体泄漏点处的内侧压力；a
l
代表油管上泄漏点的当量面积；ρ
ga
代表所述环形空间的气体密度；kg代表气体的绝热指数；mg代表气体摩尔质量；zg代表油管内气体的压缩因子；r代表气体常数；t
fl
代表气体泄漏点处的温度；p
al
代表气体泄漏点处的外侧压力；cre代表临界压力比；p
an
代表所述环形空间的环空压力；g代表重力加速度；ρ
l
代表所述环形空间的液体密度；h
l
自气体泄漏点至井口的距离，即泄漏点所处的深度；hg代表所述环形空间的气柱长度。
[0042]
其中一项实施例中，t
fl
满足公式：
[0043][0044]
其中，t
fl
代表气体泄漏点处的温度；代表第i个分段对应的油管内的气体温度；i代表所述油管的失效位置对应的分段编号；
[0045]
zg满足公式：
[0046][0047]
其中，代表第i个分段对应的油管内气体的压缩因子；zg代表气体压缩因子；i代表所述油管的失效位置对应的分段编号；
[0048]
p
fl
满足公式：
[0049][0050]
其中，p
fl
代表气体泄漏点处的内侧压力；代表第i个分段对应的油管内气体压力；i代表所述油管的失效位置对应的分段编号；
[0051]
满足公式：
[0052][0053]
其中，代表第i个分段对应的油管内气体压力；代表第i-1个分段对应的油管内气体压力；δz代表每个所述分段对应的油管的长度；代表第i个分段对应的油管内气体压力；g代表重力加速度；θ代表井斜角；fi代表第i个分段的油管内的气体摩擦系数；代表第i个分段对应的油管内气体的流速；代表第i个分段对应的油管的内径；代表第i个分段对应的油管内气体的流速变化值；
[0054]
满足公式：
[0055][0056]
其中，代表第i个分段对应的油管内气体的流速变化值；代表第i个分段对应的油管内气体的流速；代表第i-1个分段对应的油管内气体的流速；
[0057]fi
满足公式：
[0058][0059]
其中，fi代表第i个分段的油管内的气体摩擦系数；rai代表第i个分段对应的油管的粗糙度；代表第i个分段对应的油管的内径；μi代表第i个分段对应的油管内的气体粘度；代表第i个分段对应的油管内气体压力；代表第i个分段对应的油管内气体的流速变化值；
[0060]
满足公式：
[0061][0062]
其中，第i个分段对应的油管内气体压力；代表第i个分段对应的油管内气体的压缩因子；mg代表气体摩尔质量；r代表气体常数；代表第i个分段对应的油管内的气体温度；
[0063]
满足公式：
[0064][0065]
其中，wf代表油管内的气体质量流量；代表第i个分段对应的油管内气体的密度；代表第i个分段对应的油管内横截面的面积；
[0066]hl
满足公式：
[0067]hl
＝h
b-iδz；
[0068]hl
代表自气体泄漏点至井口的距离，即泄漏点所处的深度；hb代表自井底至井口的距离，即井深；δz代表每个所述分段对应的油管的长度。
[0069]
其中一项实施例中，μi满足公式：
[0070][0071]
其中，μi代表第i个分段对应的油管内的气体粘度；μ0代表气体测量粘度；t0代表测量粘度所对应的温度；代表第i个分段对应的油管内的气体温度；b代表常数；
[0072]
其中：
[0073][0074]
其中，代表第i个分段对应的油管内气体的压缩因子；a1、a2、a3、a4、a5、a6、a7、a8、a9、a
10
均代表常数；pr代表拟对比压力；tr代表拟对比温度；γg代表气体相对密度；代表第i-1个分段对应的油管内气体压力；代表第i个分段对应的油管的拟对比温度。
[0075]
另一方面，本发明提供一种井筒失效数据采集分析装置，包括：
[0076]
获取模块，用于获取井筒中的油管和套管所形成的环形空间的流体温度变化平均值；
[0077]
所述获取模块还用于获取所述环形空间的环空压力；
[0078]
检测模块，用于检测所述环形空间的参数数据组，所述参数数据组包括所述环形空间的液柱初始体积、所述环形空间的液体等压膨胀系数、所述环形空间的横截面积、所述环形空间的液体等温压缩系数、所述环形空间的长度；
[0079]
计算模块，用于根据所述参数数据组、所述流体温度变化平均值和所述环空压力，确定所述环形空间的液柱长度和气柱长度；
[0080]
所述获取模块还用于获取油管泄漏点的内外侧压力比；
[0081]
若所述内外侧压力比未达到预设压力比，则返回步骤所述获取所述井筒中油管和套管之间的流体温度变化平均值，直到所述内外侧压力达到预设压力比；
[0082]
所述获取模块若所述内外侧压力比未达到预设压力比，则返回执行步骤所述获取所述井筒中油管和套管之间的流体温度变化平均值。
[0083]
另一方面，本发明还提供一种计算机设备，包括存储器，处理器和收发器，所述存储器用于存储指令，所述收发器用于和其他设备通信，所述处理器用于执行所述存储器中存储的指令，以使所述计算机设备执行如上实施例所述的方法。
[0084]
另一方面，本发明还提供一种计算机可读存储介质，所述计算机可读存储介质中存储有计算机执行指令，当所述指令被执行时，使得计算机执行如上实施例所述的方法。
附图说明
[0085]
此处的附图被并入说明书中并构成本说明书的一部分，示出了符合本公开的实施例，并与说明书一起用于解释本公开的原理。
[0086]
图1为本发明的一个实施例提供的井筒失效数据采集方法的应用示意图。
[0087]
图2为本发明的一个实施例提供的井筒失效数据采集方法的流程示意图。
[0088]
图3为本发明的又一个实施例提供的井筒失效数据采集方法的流程示意图。
[0089]
图4为本发明的一个实施例提供的井筒失效数据采集装置的示意图。
[0090]
图5为本发明的一个实施例提供的计算机设备的示意图。
[0091]
通过上述附图，已示出本公开明确的实施例，后文中将有更详细的描述。这些附图和文字描述并不是为了通过任何方式限制本公开构思的范围，而是通过参考特定实施例为本领域技术人员说明本公开的概念。
具体实施方式
[0092]
这里将详细地对示例性实施例进行说明，其示例表示在附图中。下面的描述涉及附图时，除非另有表示，不同附图中的相同数字表示相同或相似的要素。以下示例性实施例中所描述的实施方式并不代表与本公开相一致的所有实施方式。相反，它们仅是与如所附权利要求书中所详述的、本公开的一些方面相一致的装置和方法的例子。
[0093]
本发明提供的井筒失效数据采集分析方法，应用于气井，图1为本发明提供的井筒失效数据采集分析方法的应用示意图，图中包括油管、套管、由油管和套管形成的环形空间、套管外层的水泥环等。
[0094]
气井的油管作为地层高温高压气体产出通道，所处环境和条件极为复杂，在应力、高温、高压和流体的腐蚀作用下，会发生完整性失效，如丝扣泄漏、腐蚀穿孔和应力开裂、乃至断裂等。油管完整性失效后，油管内的高温高压气体进入油管和套管之间的空间，并上浮聚集在井口处，进而改变了原有的井筒内部流体分布，引起井口高压，即引起持续环空压力。当井口的压力释放后，在压差的作用下气体会沿着气窜通道重新补充进环形空间，增加了生产成本和管理风险。因此，有必要分析油管完整性失效后井筒的压力变化规律和流体分布特征。
[0095]
本发明提供的井筒失效数据采集分析方法及其装置、存储介质旨在解决现有技术的如上问题。
[0096]
下面以具体地实施例对本发明的技术方案以及本技术的技术方案如何解决上述技术问题进行详细说明。下面这几个具体的实施例可以相互结合，对于相同或相似的概念
或过程可能在某些实施例中不再赘述。下面将结合附图，对本发明的实施例进行描述。
[0097]
请参见图2，本发明提供一种井筒失效数据采集分析方法，包括：
[0098]
s100，获取井筒中的油管和套管所形成的环形空间的流体温度变化平均值。
[0099]
该流体温度变化平均值，指的是流体温度的变化值的平均值。在一个实施例中，用δt
al
代表该流体温度变化平均值，单位为℃。
[0100]
s200，获取所述环形空间的环空压力。
[0101]
在一个实施例中，用p
an
代表所述环形空间的环空压力，单位为pa。
[0102]
s300，检测所述环形空间的参数数据组，所述参数数据组包括所述环形空间的液柱初始体积、所述环形空间的液体等压膨胀系数、所述环形空间的横截面积、所述环形空间的液体等温压缩系数、所述环形空间的长度。
[0103]
在一个实施例中，用v
al
代表所述环形空间的液柱初始体积，单位为m3；用α
p
代表所述液体等压膨胀系数，单位为℃-1
；用aa代表所述环形空间的横截面积，单位为m2；用k
t
代表所述液体等温压缩系数，单位为mpa-1
；用ha代表该环形空间的长度，单位为m。
[0104]
s400，根据所述参数数据组、所述流体温度变化平均值和所述环空压力，确定所述环形空间的液柱长度和气柱长度。
[0105]
在一个实施例中，用hy代表所述环形空间的液柱长度，单位为m；用hg代表所述环形空间的气柱长度，单位为m。
[0106]
s500，获取油管泄漏点的内外侧压力比。
[0107]
s600，若所述内外侧压力比未达到预设压力比，则返回执行步骤所述获取所述井筒中油管和套管之间的流体温度变化平均值。
[0108]
s700，若所述内外侧压力比达到所述预设压力比，则根据获取到的多组所述环形空间的液柱长度和气柱长度，分析井筒压力变化规律和流体分布特征。
[0109]
在一个实施例中，该预设压力比未99％-100％，达到该预设压力比证明油管的泄漏点已不再泄露。在油管泄漏点的内外侧压力比还没有达到该预设压力比的时候，需要不断得采集该环形空间的液柱长度和气柱长度，最后通过采集到的多组该环形空间的液柱长度和气柱长度来分析井筒压力变化规律和流体分布特征。
[0110]
本发明提供一种井筒失效数据采集分析方法，在油管泄漏点的内外侧压力比达到预设压力比之前，不断采集油管和套管之间，即环形空间的液柱长度数据和气柱长度数据，通过多组液柱长度数据和多组气柱长度数据分析井筒压力变化规律和流体分布特征。本发明提供的井筒失效数据采集分析方法可以对井筒失效后的压力变化规律和流体分布特征进行分析，为气井的设计、生产、管理和维修作业提供指导。
[0111]
请参见图3，在本发明的一个实施例中，s400包括：
[0112]
s410，根据所述液柱初始体积、所述流体温度变化值、所述液体等压膨胀系数、所述环形压力、所述液体等温压缩系数和所述横截面积确定所述液柱长度。
[0113]
其中，s410包括：
[0114]
s411，根据所述液柱初始体积、所述流体温度变化平均值、所述液体等压膨胀系数、所述环空压力和所述液体等温压缩系数确定第一参数。
[0115]
s412，根据所述第一参数和所述横截面积的比值确定所述液柱长度。
[0116]
其中，s412包括：根据公式确定所述液柱长度。其中，hy代表所述环形空间的液柱长度，单位m；v
al
代表所述环形空间的液柱初始体积，为已知量，单位为m3；δt
al
代表所述流体温度变化平均值，单位为℃；α
p
代表所述液体等压膨胀系数，为已知量，单位为℃-1
；p
an
代表所述环空压力，单位为pa；k
t
代表所述液体等温压缩系数，为已知量，单位为mpa-1
；aa代表所述环形空间的横截面积，为已知量，单位为m2。该第一参数指的是v
al
[1-(δt
al
α
p-10-6
p
ankt
)]。
[0117]
s420，根据所述环形空间的长度和所述液柱长度的差值确定所述气柱长度。
[0118]
s420包括：根据公式hg＝h
a-hy确定该气柱长度。其中，hg代表该气柱长度，单位为m；ha代表该环形空间的长度，为已知量，单位为m；hy代表所述环形空间的液柱长度，单位m。
[0119]
在本发明的一个实施例中，s100包括：
[0120]
s110，自井底至井口将所述油管等间距划分为i个分段，并用井筒分段编号值标记每个分段。
[0121]
每个分段的长度为δz，δz能被井筒长度整除，在一个实施例中，δz的取值范围为1米至10米。该井筒分段编号指的是自井底至井口，按照1，2，3，4等进行编号，例如，包含井底的第一个分段编号为1。在本实施例中，井筒被等间距划分后一共分为i个分段，本实施例中用第i个分段描述是第几个分段。
[0122]
s120，根据公式计算所述流体温度变化平均值；
[0123]
s130，其中，δt
al
代表所述流体温度变化平均值，单位为℃；代表第i个分段对应的环形空间的流体温度，单位为℃；代表第i个分段对应的油管的初始地层温度，单位为℃；n代表所述环形空间的流体终止深度对应的分段井筒编号，m代表所述环形空间的流体起始深度对应的井筒分段编号。
[0124]
其中，满足公式：
[0125][0126]
其中，代表所述第i个分段对应的环形空间的流体温度，单位为℃；代表第i个分段对应的油管内的气体温度，单位为℃；代表第i个分段对应的油管的油管中心到井筒外侧的径向热阻，为已知量，单位为m
·
℃/w；代表所述第i个分段对应的油管的初始地层温度，为已知量，单位为℃；td代表无因次地层温度，为已知量；λe代表地层导热系数，为已知量；代表第i个分段对应的油管的油管中心至油管外侧的径向热阻，为已知量，单位为m
·
℃/w；代表第i个分段对应的油管的油管中心至套管内侧的径向热阻，为已知量，单位为m
·
℃/w；π代表圆周率。
[0127]
在本发明的一个实施例中，s200包括：
[0128]
根据公式v
ag
v
al
δv
al
＝v
an
计算所述环形空间的环空压力，其中，v
an
代表所述环形空间的总体积，v
an
为定值，单位为m3；其中，δv
al
代表所述环形空间的液柱体积变化量，单
位为m3；v
al
代表所述环形空间的液柱初始体积，单位为m3；v
ag
代表所述环形空间的气柱体积，单位为m3；
[0129]
其中，v
ag
满足公式：
[0130]
n代表所述环形空间的流体终止深度对应的分段井筒编号，m代表所述环形空间的流体起始深度对应的井筒分段编号；v
gx
代表标准状况下第x个时间步长时气体自所述油管泄漏入所述环形空间的总体积，单位为m3；v
gx-1
代表标准状况下第x-1个时间步长时气体自所述油管泄漏入所述环形空间的总体积，单位为m3；在一个实施例中，取时间步长δt，叠加获取气体泄漏进入该环形空间的体积，由此计算气体自漏入该环形空间的总体积，第x个时间步长，指的是第x个时间步长δt。比如x为3，则v
gx
是第3个时间步长δt时气体自该油管泄漏入该环形空间的总体积。
[0131]
ps代表标况压力，为已知量；p
al
代表气体泄漏点处的外侧压力，为已知量，单位为pa；ts代表标况温度，为已知量，单位为℃；t
ai
代表所述环形空间的流体温度，为已知量，单位为℃。
[0132]vgx
满足公式：
[0133]
其中，q
lgx-1
代表标准状况下，第x-1个时间步长时油管内的气体在完整性失效点的泄漏速率，单位为m3/s；v
gx
代表标准状况下第x个时间步长时气体自所述油管泄漏入所述环形空间的总体积，单位为m3；v
gx-1
代表标准状况下第x-1个时间步长时气体自所述油管泄漏入所述环形空间的总体积，单位为m3；δt代表时间步长，单位为s；
[0134]qlg
满足公式：
[0135][0136]
cre满足公式：
[0137][0138]
p
al
满足公式：
[0139][0140]
其中，
[0141]qlg
代表标准状况下所述油管内的气体在完整性失效点的泄漏速率，单位为m3/s；co代表流量系数，为已知量；p
fl
代表气体泄漏点处的内侧压力，单位为pa；a
l
代表油管上泄漏点的当量面积，单位为m2；ρ
ga
代表所述环形空间的气体密度，单位为kg/m3；kg代表气体的
绝热指数，为已知量；mg代表气体摩尔质量，为已知量，单位为kg/mol；zg代表油管内气体的压缩因子，为已知量；r代表气体常数，为已知量；t
fl
代表气体泄漏点处的温度，为已知量；p
al
代表气体泄漏点处的外侧压力，为已知量，单位为pa；cre代表临界压力比；p
an
代表所述环形空间的环空压力，单位为pa；g代表重力加速度；ρ
l
代表所述环形空间的液体密度，单位为kg/m3；h
l
自气体泄漏点至井口的距离，即泄漏点所处的深度，单位为m；hg代表所述环形空间的气柱长度，单位为m。
[0142]
可以看到，v
ag
、v
al
、δv
al
都是关于该环形空间的环空压力p
an
的变量，而v
an
是一个定值，因此根据v
an
便可以推算出p
an
的值。
[0143]
在本发明的一个实施例中，t
fl
满足公式：
[0144][0145]
其中，t
fl
代表气体泄漏点处的温度，单位为℃；代表第i个分段对应的油管内的气体温度，单位为℃；i代表所述油管的失效位置对应的分段编号。
[0146]
zg满足公式：
[0147][0148]
其中，代表第i个分段对应的油管内气体的压缩因子；zg代表气体压缩因子；i代表所述油管的失效位置对应的分段编号；
[0149]
pf
l
满足公式：
[0150][0151]
其中，p
fl
代表气体泄漏点处的内侧压力，单位为pa；代表第i个分段对应的油管内气体压力，单位为pa；i代表所述油管的失效位置对应的分段编号；
[0152]
满足公式：
[0153][0154]
其中，代表第i个分段对应的油管内气体压力，为已知量，单位为pa；代表第i-1个分段对应的油管内气体压力，为已知量，单位为pa；δz代表每个所述分段对应的油管的长度，为已知量，单位为m；代表第i个分段对应的油管内气体的密度，单位为kg/m3；g代表重力加速度；θ代表井斜角；fi代表第i个分段的油管内的气体摩擦系数；代表第i个分段对应的油管内气体的流速，单位为m/s；代表第i个分段对应的油管的内径，单位为m；代表第i个分段对应的油管内气体的流速变化值，单位为m/s；
[0155]
满足公式：
[0156][0157]
其中，代表第i个分段对应的油管内气体的流速变化值，单位为m/s；代表第i
个分段对应的油管内气体的流速，单位为m/s；代表第i-1个分段对应的油管内气体的流速，单位为m/s；
[0158]fi
满足公式：
[0159][0160]
其中，fi代表第i个分段的油管内的气体摩擦系数，无因次；rai代表第i个分段对应的油管的粗糙度，为已知量，单位为m；代表第i个分段对应的油管的内径，为已知量，单位为m；μi代表第i个分段对应的油管内的气体粘度，单位为pa
·
s；代表第i个分段对应的油管内气体的密度，单位为kg/m3；代表第i个分段对应的油管内气体的流速变化值，单位为m/s；
[0161]
其中，μi满足公式：
[0162][0163]
其中，μi代表第i个分段对应的油管内的气体粘度，单位为pa
·
s；μ0代表气体测量粘度，单位为pa
·
s；t0代表测量粘度所对应的温度，单位为℃；代表第i个分段对应的油管内的气体温度，单位为℃；b代表常数，取110.4；
[0164]
其中：
[0165][0166]
其中，代表第i个分段对应的油管内气体的压缩因子，为已知量，无因次；a1、a2、a3、a4、a5、a6、a7、a8、a9、a
10
均代表常数，为已知量，无因次；pr代表拟对比压力，无因次；tr代表拟对比温度，无因次；γg代表气体相对密度，为已知量，无因次；代表第i-1个分段对应的油管内气体压力，单位为pa；代表第i个分段对应的油管的拟对比温度，无因次。
[0167]
满足公式：
[0168][0169]
其中，第i个分段对应的油管内气体压力，单位为pa；代表第i个分段对应的油管内气体的压缩因子，无因次；mg代表气体摩尔质量，单位为kg/mol；r代表气体常数，无因次；代表第i个分段对应的油管内的气体温度，单位为℃；
[0170]
满足公式：
[0171][0172]
其中，wf代表油管内的气体质量流量，单位为kg/s；代表第i个分段对应的油管内气体的密度，单位为kg/m3；代表第i个分段对应的油管内横截面的面积，单位为m2；
[0173]hl
满足公式：
[0174]hl
＝h
b-iδz；
[0175]hl
代表自气体泄漏点至井口的距离，即泄漏点所处的深度，单位为m；hb代表自井底至井口的距离，即井深，单位为m；δz代表每个所述分段对应的油管的长度，单位为m。
[0176]
请参见图4，本发明还提供一种井筒失效数据采集分析装置10，包括：
[0177]
获取模块11，用于获取井筒中的油管和套管所形成的环形空间的流体温度变化平均值。
[0178]
所述获取模块11还用于：
[0179]
自井底至井口将所述油管等间距划分为i个分段，并用井筒分段编号值标记每个分段；
[0180]
根据公式计算所述流体温度变化平均值；
[0181]
其中，δt
al
代表所述流体温度变化平均值；代表第i个分段对应的环形空间的流体温度；代表第i个分段对应的油管的初始地层温度；n代表所述环形空间的流体终止深度对应的分段井筒编号，m代表所述环形空间的流体起始深度对应的井筒分段编号。
[0182]
满足公式：
[0183][0184]
其中，代表所述第i个分段对应的环形空间的流体温度；代表第i个分段对应的油管内的气体温度；代表第i个分段对应的油管的油管中心到井筒外侧的径向热阻；代表所述第i个分段对应的油管的初始地层温度；td代表无因次地层温度；λe代表地层导热系数；代表第i个分段对应的油管的油管中心至油管外侧的径向热阻；代表第i个分段对应的油管的油管中心至套管内侧的径向热阻；π代表圆周率。
[0185]
所述获取模块11还用于获取所述环形空间的环空压力。
[0186]
所述获取模块还用于根据公式v
ag
v
al
δv
al
＝v
an
计算所述环形空间的环空压力，其中，v
an
代表所述环形空间的总体积；其中，δv
al
代表所述环形空间的液柱体积变化量；v
al
代表所述环形空间的液柱初始体积；v
ag
代表所述环形空间的气柱体积；
[0187]
其中，v
ag
满足公式：
[0188]
n代表所述环形空间的流体终止深度对应的分段井筒编号，m代表所述环形空间的
流体起始深度对应的井筒分段编号；v
gx
代表第x个时间步长时气体自所述油管泄漏入所述环形空间的总体积；v
gx-1
代表第x-1个时间步长时气体自所述油管泄漏入所述环形空间的总体积；ps代表标况压力；p
al
代表气体泄漏点处的外侧压力；ts代表标况温度；t
ai
代表所述环形空间的流体温度；
[0189]vgx
满足公式：
[0190]
其中，q
lgx-1
代表第x-1个时间步长时油管内的气体在完整性失效点的泄漏速率；v
gx
代表第x个时间步长时气体自所述油管泄漏入所述环形空间的总体积；v
gx-1
代表第x-1个时间步长时气体自所述油管泄漏入所述环形空间的总体积；δt代表时间步长；
[0191]qlg
满足公式：
[0192][0193]
cre满足公式：
[0194][0195]
p
al
满足公式：
[0196][0197]
其中，
[0198]qlg
代表所述油管内的气体在完整性失效点的泄漏速率；co代表流量系数；p
fl
代表气体泄漏点处的内侧压力；a
l
代表油管上泄漏点的当量面积；ρ
ga
代表所述环形空间的气体密度；kg代表气体的绝热指数；mg代表气体摩尔质量；zg代表油管内气体的压缩因子；r代表气体常数；t
fl
代表气体泄漏点处的温度；p
al
代表气体泄漏点处的外侧压力；cre代表临界压力比；p
an
代表所述环形空间的环空压力；g代表重力加速度；ρ
l
代表所述环形空间的液体密度；h
l
自气体泄漏点至井口的距离，即泄漏点所处的深度；hg代表所述环形空间的气柱长度。
[0199]
t
fl
满足公式：
[0200][0201]
其中，t
fl
代表气体泄漏点处的温度；代表第i个分段对应的油管内的气体温度；i代表所述油管的失效位置对应的分段编号；
[0202]
zg满足公式：
[0203][0204]
其中，代表第i个分段对应的油管内气体的压缩因子；zg代表气体压缩因子；i代
表所述油管的失效位置对应的分段编号；
[0205]
p
fl
满足公式：
[0206][0207]
其中，pf
l
代表气体泄漏点处的内侧压力；代表第i个分段对应的油管内气体压力；i代表所述油管的失效位置对应的分段编号；
[0208]
满足公式：
[0209][0210]
其中，代表第i个分段对应的油管内气体压力；代表第i-1个分段对应的油管内气体压力；δz代表每个所述分段对应的油管的长度；代表第i个分段对应的油管内气体压力；g代表重力加速度；θ代表井斜角；fi代表第i个分段的油管内的气体摩擦系数；代表第i个分段对应的油管内气体的流速；代表第i个分段对应的油管的内径；代表第i个分段对应的油管内气体的流速变化值；
[0211]
满足公式：
[0212][0213]
其中，代表第i个分段对应的油管内气体的流速变化值；代表第i个分段对应的油管内气体的流速；代表第i-1个分段对应的油管内气体的流速；
[0214]fi
满足公式：
[0215][0216]
其中，fi代表第i个分段的油管内的气体摩擦系数；rai代表第i个分段对应的油管的粗糙度；代表第i个分段对应的油管的内径；μi代表第i个分段对应的油管内的气体粘度；代表第i个分段对应的油管内气体压力；代表第i个分段对应的油管内气体的流速变化值；
[0217]
满足公式：
[0218][0219]
其中，第i个分段对应的油管内气体压力；代表第i个分段对应的油管内气体的压缩因子；mg代表气体摩尔质量；r代表气体常数；代表第i个分段对应的油管内的气体温度；
[0220]
满足公式：
[0221][0222]
其中，wf代表油管内的气体质量流量；代表第i个分段对应的油管内气体的密度；代表第i个分段对应的油管内横截面的面积；
[0223]hl
满足公式：
[0224]hl
＝h
b-iδz；
[0225]hl
代表自气体泄漏点至井口的距离，即泄漏点所处的深度；hb代表自井底至井口的距离，即井深；δz代表每个所述分段对应的油管的长度。
[0226]
μi满足公式：
[0227][0228]
其中，μi代表第i个分段对应的油管内的气体粘度；μ0代表气体测量粘度；t0代表测量粘度所对应的温度；代表第i个分段对应的油管内的气体温度；b代表常数；
[0229]
其中：
[0230][0231]
其中，代表第i个分段对应的油管内气体的压缩因子；a1、a2、a3、a4、a5、a6、a7、a8、a9、a
10
均代表常数；pr代表拟对比压力；tr代表拟对比温度；γg代表气体相对密度；代表第i-1个分段对应的油管内气体压力；代表第i个分段对应的油管的拟对比温度。
[0232]
检测模块12，用于检测所述环形空间的参数数据组，所述参数数据组包括所述环形空间的液柱初始体积、所述环形空间的液体等压膨胀系数、所述环形空间的横截面积、所述环形空间的液体等温压缩系数、所述环形空间的长度；
[0233]
计算模块13，用于根据所述参数数据组、所述流体温度变化平均值和所述环空压力，确定所述环形空间的液柱长度和气柱长度；
[0234]
所述计算模块13还用于根据所述液柱初始体积、所述流体温度变化平均值、所述液体等压膨胀系数、所述环空压力、所述液体等温压缩系数和所述横截面积确定所述液柱长度；根据所述环形空间的长度和所述液柱长度的差值确定所述气柱长度。所述计算模块13还用于根据所述液柱初始体积、所述流体温度变化平均值、所述液体等压膨胀系数、所述环空压力和所述液体等温压缩系数确定第一参数；根据所述第一参数和所述横截面积的比值确定所述液柱长度。
[0235]
所述计算模块13还用于根据公式确定所述液柱长度，其中，hy代表所述环形空间的液柱长度；v
al
代表所述环形空间的液柱初始体积；δt
al
代表所述流体温度变化平均值；α
p
代表所述液体等压膨胀系数；p
an
代表所述环空压力；k
t
代表所述液体等温压缩系数；aa代表所述环形空间的横截面积。
[0236]
所述获取模块11还用于获取油管泄漏点的内外侧压力比；
[0237]
所述获取模块11还用关于若所述内外侧压力比未达到预设压力比，则返回步骤所述获取所述井筒中油管和套管之间的流体温度变化平均值；
[0238]
所述获取模块11还用于若所述内外侧压力比达到预设压力比，则返回执行步骤所述获取所述井筒中油管和套管之间的流体温度变化平均值；
[0239]
所述计算模块13还用于若所述内外侧压力比达到所述预设压力比，则根据获取到的多组所述环形空间的液柱长度和气柱长度，分析井筒压力变化规律和流体分布特征。
[0240]
请参见图5，本发明还提供一种计算机设备20，包括存储器21，处理器22和收发器23，该存储器21用于存储指令，该收发器23用于和其他设备通信，该处理器22用于执行该存储器21中存储的指令，以使该计算机设备执行如上任一项所述的井筒失效数据采集分析方法。
[0241]
本发明还提供一种计算机可读存储介质，该计算机可读存储介质中存储有计算机执行指令，当该指令被执行时，使得计算机执行指令被处理器执行时用于实现如上任一项实施例提供的该井筒失效数据采集分析方法。本发明还提供一种另一种计算机可读存储介质，该计算机可读存储介质中存储有计算机执行指令，当该指令被执行时，使得计算机执行如上任一项实施例提供的接口报文转换方法。
[0242]
需要说明的是，上述计算机可读存储介质可以是只读存储器(read only memory，rom)、可编程只读存储器(programmable read-only memory，prom)、可擦除可编程只读存储器(erasable programmable read-only memory，eprom)、电可擦除可编程只读存储器(electrically erasable programmable read-only memory，eeprom)、磁性随机存取存储器(ferromagnetic random access memory，fram)、快闪存储器(flash memory)、磁表面存储器、光盘、或只读光盘(compact disc read-only memory，cd-rom)等存储器；也可以是包括上述存储器之一或任意组合的各种电子设备，如移动电话、计算机、平板设备、个人数字助理等。
[0243]
需要说明的是，在本文中，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者装置不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、物品或者装置所固有的要素。在没有更多限制的情况下，由语句“包括一个
……”
限定的要素，并不排除在包括该要素的过程、方法、物品或者装置中还存在另外的相同要素。
[0244]
上述本发明实施例序号仅仅为了描述，不代表实施例的优劣。
[0245]
通过以上的实施方式的描述，本领域的技术人员可以清楚地了解到上述实施例方法可借助软件加必需的通用硬件平台的方式来实现，当然也可以通过硬件，但很多情况下前者是更佳的实施方式。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质
(如rom/ram、磁碟、光盘)中，包括若干指令用以使得一台终端设备(可以是手机，计算机，服务器，空调器，或者网络设备等)执行本发明各个实施例所描述的方法。
[0246]
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
[0247]
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
[0248]
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
[0249]
以上仅为本发明的优选实施例，并非因此限制本发明的专利范围，凡是利用本发明说明书及附图内容所作的等效结构或等效流程变换，或直接或间接运用在其他相关的技术领域，均同理包括在本发明的专利保护范围内。