【《注入井井筒物理模型分析概述》1500字】

注入井井筒物理模型分析概述目录TOC\o"1-3"\h\u11537注入井井筒物理模型分析 141421.1注入井井筒结构物理模型 1136081.2 注CO2井筒温度分布计算方法 2210471.3注CO2井筒温度分布计算方法 31.1注入井井筒结构物理模型依据目前大部分的单井注气进行研究可以发现，大部分的井筒结构都是由套管及水泥环组成（如图1所示）。因此提出了如下的前提条件和假设：假设1-管柱的密封良好；假设2-大部分的地层中物理参数恒定；假设3-水泥环外缘到地层之间的热量传递属于一维的非稳态；假设4-在实际数学模型运算中考虑环空流体以及水泥环的热阻。图2-1中，为注入CO2温度，℃；为井的实际深度，m；为地表实际检测温度，℃；为地层温度，℃；为套管的外半径，m；为套管的内半径，m；为水泥环外半径，m。图2-1CO2注入井井筒结构物理模型Figure2-1PhysicalmodelofwellborestructureofCO2injectionwell1.2 注CO2井筒温度分布计算方法根据能量守恒定理，结合传热学理论可以知道，目前的井筒CO2变化产生的热量与水泥环外壁向CO2传递的热量之间存在相互恒定的情况，因此本人把无因次时间的函数纳入其中，整合到实际的温度梯度计算方程内。 （2-1）其中： （2-2）式中，为注入CO2温度，℃；为地温梯度，℃/m；为总传热系数，；为地表温度，℃；为地层热导率，；为注入流体定压比热容，；为注入质量流量，；为无因次时间函数；系数表示温度梯度与注入CO2的物性和CO2的注入参数密不可分。而基于无因次的时间函数，具体的可以参考Hasan-Kabir方法： （2-3）其中： （2-4）式中，为CO2注入的实际检测时间，；为地层在实际运作过程中的热扩散系数。在注气井中，环空气体：比如说CO2、空气等的实际热阻不易检测，因此通常会使用方程进行考虑，同时由于CO2属于结构对称三原子，因此它不存在辐射情况，鉴于套管壁和油管的厚度不大，因此可以不计算辐射换热系数。在考虑方程以后基于复合多层圆管壁热阻串联原理，同时结合热力学理论，计算出总传热系数： （2-5）式中，为水泥环热导率，；为环空流体的实际流换热系数，。环空流体的对流换热系数与总传热系数挂钩，因为换热系数与环形空间的温降呈现正相关，与此同时，环空气体为二氧化碳，因此整体的热传导性不高，与地层之间的温降不大，所以可以采用平均求值的方法进行优化。1.3注CO2井筒温度分布计算方法以井筒其中一个微元段为研究对象，由一般流体连续性方程可得： （2-6）根据气体稳定流动能量方程，可以得到井筒压力梯度方程为： （2-7）式中，为流速，m/s；为氮体密度，kg/m3；为压强，Pa；为重力加速度，9.8m/s2；为注气井井筒的倾角，在实际计算中，需要充分考虑水平位移以及垂深两大因素，使用条件一般为直井或者稍微有一点斜角的定向井，°；为摩阻系数；为油管内径，m。上式右侧中的三项分别为重力势能，摩阻和动能损失，毫无疑问动能损失是因为CO2的密度变化，井深增加后，一系列的参数如温度、压力、黏度、密度、压缩因子等都会有所不同，因此在模型推导过程中必须要进行有效迭代。此外，借助压缩因子来修正CO2的密度方程后可以表示为： （2-8）式中，为压缩因子；为温度，℃。气体压缩因子需要使用查图法、查表法、和状态方程法等确定，而本文文献[33]提出的的CO2压缩因子计算方程式，根据此次试验CO2井筒的压力范围，合理推导出公式： （2-9）其中：

CO2的运动黏度是压力和温度的两个相关值，从牛顿的内摩擦定律出发，以“二氧化碳黏度表”的数据作为参考，回归出CO2运动和黏度的方程，同时把雷诺数结合其中，两种公式合并后可得： （2-10）式中，为雷诺数；为CO2体积流量，m3/h；为油管内径，mm；为温度，K；雷诺数是气体流动特性评估的关键数据，通过雷诺数可以判断管道内的流体状态和实际情况。此次试验的CO2雷诺数大于3000，流动状态为湍流流动状态，因此实际的摩阻系数，经过对比研究后决定使用尼古拉斯方程计算[33]： （2-11）式中，为达西摩阻系