一种气井井喷着火情况下的地层压力反演系统与反演方法与流程

技术特征：
1.一种气井井喷着火情况下的地层压力反演系统，其特征在于，包括地面井喷火焰高度测量系统和数据信息反演系统；地面井喷火焰高度测量系统包括：红外热成像摄像机、可见光摄像机、小型气象站、测距传感器；红外热成像摄像机和可见光摄像机用于拍摄气井井喷火焰现场图像，用来研究、分析和确定火焰的高度；小型气象站用来连续记录以下气象条件：风向、风速和环境温度；测距传感器用来计算布置红外热成像摄像机和可见光摄像机距井喷火焰距离；数据信息反演系统包括：数据采集管线、计算机、显示器；数据采集管线连接计算机和地面井喷火焰高度测量系统，向计算机传递测量数据；计算机根据获得的测量数据进行地层信息反演，得到地层压力；显示器连接计算机，用于直观输出地层压力反演结果。2.根据权利要求1所述的气井井喷着火情况下的地层压力反演系统，其特征在于，火焰高度确定由红外热成像摄像机、可见光摄像机拍摄到的图像，由计算机处理识别，以温度800k作为井喷火焰边界。3.根据权利要求1所述的气井井喷着火情况下的地层压力反演系统，其特征在于，所述的测距传感器为防爆激光测距传感器。4.一种利用权利要求1-3任意一项权利要求所述的气井井喷着火情况下的地层压力反演系统的反演方法，其特征在于，包括步骤如下：在确保人员生命安全和设备正常运行的情况下，采用测距传感器确定安全距离；将现场的红外热成像摄像机和可见光摄像机设置为记录模式，测得的火焰图像和测距传感器测得的距离通过数据采集管线导入到计算机，通过地层压力反演进行处理分析，反演得到地层压力。5.一种气井井喷着火情况下的地层压力反演系统的反演方法，其特征在于，具体步骤包括如下：(1)根据现场钻井施工资料，确定地层压力反演所需要的工程参数；(2)根据红外热成像摄像机和可见光摄像机测得的火焰图像计算火焰高度，由于环境风的存在，火焰高度计算公式为：式(1)中l0为井喷火焰长度，m；h为测量的火焰高度，m；θ为井喷火焰与水平线的夹角(0
°
～90
°
)；u
w
为当地风速，m/s；(3)由于失控井喷流射流速度通常处于超音速区域，属于气体射流，因此采用spalart—allmaras模型对失控井喷射流进行描述：式中t表示时间，s；ρ表示混合气体密度，kg/m3；x
i
，x
j
分别表示不同位置的张量符号；表示湍流黏度系数；u
i
表示时均匀速度分量；g
v
表示湍流黏度产生项；y
v
表示由于壁面阻挡与黏性阻尼引起的湍流黏度减少项；和c
b2
为常数；μ表示分子黏度，pa
·
s；(4)在超音速流动的射流过程中，由于真实气体具有可压缩性，故利用泊松式和pv＝
znrt得：pρ-γ
＝c
ꢀꢀꢀꢀ
(3)再结合欧拉方程推导出真实气体伯努利方程：式中p为井喷出口位置压力，pa；v为气体体积，m3；z为天然气压缩系数；n为气体的物质的量；mol；r为摩尔气体常数；j/(mol
·
k)；t表示绝对温度，k；γ为比热比；v为井喷气体速度，m/s；g为重力加速度m/s2；h为火焰高度，m；c为常数；(5)由真实气体伯努利方程建立火焰高度与流速的表达公式：式中p0为井底压力，pa；ρ0为井底气体密度；p为井喷出口位置压力，pa；f为阻力系数，m/s；l为套管长度，m；r为套管半径，m；表示雷诺平均后的平均速度，m/s；(6)根据现场施工资料确定气体组分，可由公式(6)计算井喷气体相对分子质量：式中m
g
表示天然气相对分子质量；x
i
表示天然气组分i的摩尔数；n表示组分数；m
i
为组分i的相对分子质量；(7)将气体理想状态方程(7)与公式(5)联立可求出气体井喷着火井口处的气体流速表达式(8)：达式(8)：式中ρ
g
表示气体密度kg/m3；t0表示相对温度，k；r为摩尔气体常数；j/(mol
·
k)；(8)在获得井口流速的基础上，计算井喷气体产量q公式(9)为：式中k表示井喷气体产量与流量的转换系数，该系数由气体质量流量转换系数表查得；(9)在对地层压力求解过程中，需要先对井底流动压力p
s
进行求解，得到井底流动压力后，方可求出地层压力；在对井底流动压力p
s
的求解方法如下：9-1、基本方程气体从井底沿井筒流动到井口过程中，假设流体为稳定流，取长度dl的管段为控制体进行研究，则可以根据能量方程写成方程式(10)：d
p
ρvgdv dw dl
w
＝0
ꢀꢀꢀꢀ
(10)式中d
p
为管长dl内对应的总压降；ρ为流动状态下的气体密度，kg/m3；g为重力加速度m/s2；v为气体速度，m/s；dv为管长dl内对应的速度变化；dw为外界对气体所做的功；dl
w
摩擦引起的压力损失；
对于垂直管气体流动：(1)从井底到井口没有功的变化，即dw＝0；(2)动能损失忽略不计，即vdv＝0；则公式(10)可以写成：式中d为管径，m，f为阻力系数可以用colebrook公式：式中e套管内壁绝对粗糙度，m；标准状态下，取标准状况下大气压p
sc
＝0.01325mpa，标准状况下温度t
sc
＝293k，q
sc
为标准状态下气体流量，m3/d，则在管内任意一点(p，t)下：/d，则在管内任意一点(p，t)下：式中，b
g
为气体非达西流动系数；a为达西流动系数；γ
g
为井喷气体摩尔分数；将公式(13)、(14)代入公式(12)，并除以重力加速度，整理后得到公式(15)：考虑到井内气体向上流动时，由于压力差的作用，气体的压力是逐渐减小的，可以对上式进行变形，并写成积分形式：式中p2计算段终点流动压力，mpa；p1计算段起流动压力，mpa；h2计算段终点的高度，m；h1计算段起点的高度，m；9-2、求解方法在公式(16)中，p、t为l的函数，z又是t和p的函数，因此，公式(16)的求解可用迭代法进行计算，整个井筒的长度h平均分为n段，每段长度为δh；在所在的每一段中z、t用相应段的平均值，即公式(16)可以写成如下形式：
将公式(17)，简化后得公式(18)：式中，p2为计算段终点流动压力，mpa；p1为计算段起流动压力，mpa；在在条件下得压缩系数；为流动管柱δh段内得气体平均温度，k；从已知的井口压力开始，逐段往下计算，直到井底即可求出井底流动压力p
s
；(10)在求解得到井底流动压力后，将井喷气体产量代入产气二项式方程即可求得地层压力；井喷后的地层压力可以根据预测的天然气的井喷量、井喷时的井底流动压力应用产气二项方程式(19)求得；p
f
＝(p
s2
aq bq2)
0.5
ꢀꢀꢀꢀ
(19)式中p
f
为地层压力，pa；p
s
为发生井喷时的井底流动压力，pa；a为达西流动系数；b为非达西流动系数。6.根据权利要求5所述的气井井喷着火情况下的地层压力反演系统的反演方法，其特征在于，步骤(1)所述工程参数是指根据现场施工资料确定所需计算参数，包括：井眼尺寸，套管压力、地面温度、阻力系数、气体的黏度、气体的比热比、产出气体的摩尔数、产出气体相对密度、产出气的相对分子质量、放喷管线内壁绝对粗糙度、井喷气体压缩系数、达西流动系数、非达西流动系数。7.根据权利要求5所述的气井井喷着火情况下的地层压力反演系统的反演方法，其特征在于，步骤(2)中，火焰高度确定由红外热成像摄像机图像800k为火焰边界和可见光摄像机确定火焰高度。

技术总结
本发明涉及一种气井井喷着火情况下的地层压力反演系统与反演方法，属于石油与天然气工程技术领域，基于计算的火焰长度建立井口参数计算模型，通过对求解计算模型得到井喷着火工况下井口的气体流量，最后根据计算得到的井口参数反演地层压力。借助本发明的气井井喷着火情况下的地层压力反演系统精准测量井喷着火工况下火焰高度，本发明气井井喷着火情况下的地层压力反演系统与方法可以精准确定井喷着火井口气体流速及流量，为井喷着火条件下地层压力反演及抢险作业提供指导。层压力反演及抢险作业提供指导。层压力反演及抢险作业提供指导。


技术研发人员：葛磊 王志远 王海斌 崔海林 刘晓
受保护的技术使用者：中石化胜利石油工程有限公司塔里木分公司
技术研发日：2022.03.07
技术公布日：2022/6/1