-
适度出砂开采是目前国内外解决稠油生产过程中出砂问题的主要技术,稠油油井出砂的实时监测对适度出砂开采技术的应用具有非常重要的意义[1-4]。通过对油井生产过程中出砂状况进行持续有效的监测,获得实时出砂信息,可为适度出砂开采提供依据,最终指导油井生产,保障稠油高效开发。目前国内外主要采用振动监测方法[5-8]进行出砂监测,其主要由高频压电传感器、信号预处理系统、数据采集系统及计算机软件系统等组成。将高频压电传感器安装在出砂流动管道外壁上,生产过程中,油流携带砂粒以一定的速度撞击管壁,产生振动,高灵敏度传感器感应并接收这些振动信号,通过滤波、放大等预处理,最终传送到计算机信号分析系统,得到出砂量信息,达到实时监测目的,从而为评价出砂状况提供依据。
为了获得准确的出砂监测结果,需要将传感器安装在最有利于接收砂粒撞击管壁振动信号的位置,这需要模拟研究油水砂三相流通过90°弯管(油气井生产主要弯管形式)时的流动状态,确定砂粒经过弯管时撞击力度最大的地方,依此确定传感器的安装位置,以获取最强烈、最准确的出砂振动信号。利用Fluent软件对含砂油流进行了模拟,建立了油水砂三相流通过90°弯管时的流动模型,研究了砂粒撞击管壁的能量最大位置,同时建立了室内出砂监测模拟装置,进行了不同安装位置的室内对比实验,以验证理论模拟结果。
-
为更好地了解非牛顿流体油-水-固三相流流经弯管的流动特性,采用Fluent软件进行CFD计算,得到一定含水率下,一定流体黏度、含砂量、流速、砂粒粒径条件下弯管处多相流体流动的特性及固相砂粒运动轨迹。
Fluent软件Mixture模型可以模拟各相不同速度的多相流,也可以模拟各相以相同速度流动的多相流和有强耦合的各相同性多相流[9]。因此,选用Mixture模型研究弯管混合多相流,研究过程中,选用k-ε双方程湍流模型,模拟油水两相流在管道中的流动[10-12]。
-
利用有限元软件建立室内实验条件下的弯管实体模型,得到网格图如图1所示。模型进口边设为速度进口(velocity-inlet),模型出口边设为自由出口(out-flow),其余边设为管壁(wall)。弯管模型的局部放大图如图2所示。
-
固相砂粒在流场中的质量及动量载荷率非常低,连续相(油水相)基本不受离散相(砂粒)的影响,可以忽略砂粒对油水两相流的影响,因此选择非耦合计算;同时,考虑到砂粒与连续相湍流的相互作用及砂粒之间各自运动轨道的干扰,采用随机轨道模型计算。利用Fluent求解离散相问题主要包括以下步骤:(1)设定离散相轨道的迭代步数;(2)创建不同粒径砂粒的属性;(3)选择固相进入的平面及速度、粒径、质量流量等参数;(4)从计算好的油水两相流流场中查看颗粒轨迹。
为进一步确定油水砂多相流中固相砂粒冲击管壁产生最大动能的大致位置,即传感器最适合安装位置,采用Fluent流体力学分析软件对多相流中砂粒在弯管附近的流动状态进行分析[13-14]。
在多相流弯管模拟分析中,认为油-水两相流为连续相,油相均匀分布在水相中;砂粒之间的相互作用采用基于Lagrangian方法的小球模型;忽略黏性损耗,同时考虑重力的影响,得到连续性方程和动量方程
$$ \frac{\partial }{{\partial t}}({\alpha _k}{\rho _k}) + \nabla \cdot ({\alpha _k}{\rho _k}{u_k}) = 0 $$ (1) $$\begin{split} &\frac{\partial }{{\partial t}}({\alpha _k}{\rho _k}{u_k}) + \nabla \cdot ({\alpha _k}{\rho _k}{u_k}{u_k}) = \\& - {\alpha _k}\nabla p + \nabla \cdot ({\alpha _k}{\tau _k}) + {\alpha _k}{\rho _k}g - \frac{1}{{\Delta v}}\sum\nolimits_{i = 1}^n {{F_{{\rm{drag}},i}}} \end{split} $$ (2) $${\alpha _{\rm{o}}} + {\alpha _{\rm{w}}} + {\alpha _{\rm{s}}} = 1$$ (3) 式中,ρ为密度,kg/m³;α为各相体积分数,%;u为流体速度,m/s;p为压力,Pa;g为重力加速度,9.8 m/s²;τ为流体剪切压力,Pa;角标k代表不同的相,o为油相,w为水相,s为砂粒相;Fdrag为不同颗粒流动阻力,N;Δv为体积计算单元;n为网格内颗粒数量。
颗粒之间的运动方程由牛顿第二运动定律决定
$${m_{\rm{p}}}\frac{{{\rm{d}}{u_{\rm{p}}}}}{{{\rm{d}}t}} = {m_{\rm{p}}}g + {F_{{\rm{drag}}}}$$ (4) 式中,mp为颗粒质量,g;up为颗粒速度,m/s。
Fluent离散相模型对多相流中的砂粒流动特性的计算分析流程如图3所示。
-
采用润滑油模拟,油品黏度为330 mPa · s,液流含水20%,综合黏度200 mPa · s;固相砂粒密度为2 650 kg/m3,砂粒质量分数为0.1 %,粒径大小80~325目(表1),弯管外径30 mm,内径25 mm,弯管曲率半径约为2.5倍内径,流体入口流速为2.2 m/s,出口为自由出口。温度条件26 ℃,环境压力1个大气压,油样密度960 kg/m3。为了保证湍流流体在管路中充分发展,模型的流体入口和出口距离弯曲处均为32倍管内径。
表 1 砂粒目数与直径的对应关系
Table 1. Corresponding relationship between mesh and diameter of sand grain
目数 砂粒直径/μm 目数 砂粒直径/μm 80 180 200 75 100 150 325 45 150 106 -
分别在80目、100目、150目、200目、325目砂粒粒径条件下,模拟弯管内携砂流体流动,含砂量均为0.1%(质量含量),得到弯管处砂粒的运动轨迹如图4所示。
图 4 弯管处的不同粒径砂粒在油水砂三相流中的运动速度分布
Figure 4. Distribution of the movement velocity of the sand grains of different sizes in the oil-water-sand three-phase flow in the elbow
针对黏度较高的油水砂多相流体,进一步采用CFD计算方法对弯管处不同尺寸砂粒的速度和轨迹分布进行分析,结果如图4所示。图4中实线为不同尺寸砂粒在弯管内的速度和轨迹分布结果,对于油水砂三相流,不同尺寸砂粒(80~325目)均在弯管下游约2倍管径处速度达到最大值,此时动能最大,因此弯管下游2倍管径处也为监测砂粒撞击信号的传感器最优安装位置。图4中左上角方框内为某时刻,不同尺寸砂粒(80~325目)在弯管内的速度和状态分布图,结果表明无论是80目的大尺寸砂粒,还是325目的小尺寸砂粒均能冲击到管壁上。
-
在室内条件下,对安装在不同位置的传感器的出砂监测效果进行了对比实验研究,传感器的安装位置分别为弯管处、弯管下游约2倍管径处、弯管下游约3倍管径处。实验条件,采用综合黏度200 mPa · s、含水约20%的油流(主体油流为润滑油),砂样为100目的石英砂,含砂量0.1%,管流速度为2.2 m/s,与Fluent软件模拟条件一致,进行了出砂监测效果对比实验。多相含砂流体出砂监测模拟试验平台见图5,出砂监测传感器安装位置见图6。
图 5 多相含砂流体出砂监测模拟试验平台示意图
Figure 5. Schematic simulation experiment platform for monitoring the sand production of multi-phase sand containing fluid
图 6 室内条件下传感器不同安装位置示意图
Figure 6. Schematic installation position of the sensor under indoor conditions
出砂监测实验中,砂粒撞击管壁产生的振动信号实际上是确定性的受振系统在随机力的激励作用下产生振动的过程,因此,出砂振动信号可视为一种随机振动信号。采用振动信号的时频分析理论对实验数据进行了分析处理,流程见图7。图中,T为采样时间间隔,s;N为采样点数;k为傅里叶变化的第k个频谱;X(k)、Y(n)、Z(n)为傅里叶变换函数;x(n)为信号样本函数,−∞<n<+∞;W=exp(−2πj/N),j为虚单位,j2=1;P为功率谱函数。
对安装在不同位置传感器采集的出砂振动信号进行时频分析处理,得到出砂监测信号的时频特征值(表2)及功率谱图。分析可得:在其他实验条件均相同的情况下,安装在弯管下游2倍管径处得到的出砂监测信号功率谱幅值最大,信号能量最大,说明在此处砂粒撞击管壁振动最大,验证了多相管流条件下的理论模拟结果。因此,出砂监测传感器的最优安装位置为弯管下游约2倍管径处。
表 2 不同安装位置采集的出砂监测信号特征值对比
Table 2. Comparison between the characteristic values of the sand production monitoring signals acquired at different installation positions
参数 弯管处 弯管下游2倍处 弯管下游3倍处 最大时域幅值 0.292 3 0.633 7 0.425 9 最小时域幅值 −0.301 8 −0.641 8 −0.425 8 方差 2.420 0×10−4 8.894 3×10−4 4.468 1×10−4 平均值 −5.902 8×10−21 −9.635 0×10−21 −7.623 3×10−21 均方根值 0.015 6 0.029 8 0.021 1 功率谱值 2.42×10−4 8.90×10−4 4.47×10−4 -
(1)利用Fluent软件建立了含砂液流经过90°弯管流动的油水砂三相流模型,模拟了砂粒流经弯管的运动轨迹及速度分布,研究表明,砂粒流经弯管后,在下游方向约2倍管径的位置处撞击管壁的速度最大。
(2)建立了多相含砂流体出砂监测实验平台,进行了不同传感器安装位置的出砂监测实验,结合信号分析结果,表明安装在弯管下游约2倍管径处的监测效果最好,验证了模拟分析结果。因此,传感器安装在此处,可以获得最强的出砂振动信号,从而得到准确的出砂监测数据,指导适度出砂生产。
Numerical simulation and experimental verification on the installation position of sand production monitoring sensor in heavy oil wells
-
摘要: 非侵入式的稠油出砂监测技术是通过安装在管壁的传感器接收砂粒撞击信号实现出砂监测,为获得最优的砂粒撞击振动信号,保障出砂数据的可靠性和准确性,需要找到传感器在90°弯管处的最优安装位置。在室内条件下,建立了油水砂三相流的弯道流动模型,通过Fluent软件模拟了含砂油流在90°弯管中的运动轨迹,确定了砂粒在弯道流动的撞击管壁速度最大位置在弯管下游约管道直径2倍处,同时,通过建立的多相含砂流体出砂监测实验平台,对安装在不同位置的传感器的出砂监测信号进行了对比分析,验证了多相流条件下的理论模拟结果。研究表明,在弯管下游约2倍管径处安装传感器可以得到最可靠的出砂信号,该模拟研究结果为现场应用过程中传感器的安装提供了重要理论依据,保障了出砂监测数据的可靠性。Abstract: The non-intrusive sand production monitoring technology for the heavy oil production is to monitor the sand production by using the sensor installed on the pipe wall to receive the sand impinging signals. In order to get the optimal vibration signals of sand impinging and ensure the reliability and accuracy of sand production data, it is necessary to determine the optimal installation position of the sensor in the 90° elbow. In this paper, a flow-in-bend model for oil-water-sand three-phase flow was established under the indoor conditions. The movement trajectory of sand containing oil flow in the 90° elbow was simulated using Fluent software. It is determined that sand’s highest impact velocity on the pipe wall while flowing in the bend occurs at twice the pipeline diameter. What’s more, the sand production monitoring signals of the sensors installed at different positions were comparatively analyzed by means of the experimental platform for monitoring the sand production of multi-phase sand containing fluid. And based on this, the theoretical simulation results under the condition of multi-phase flow were verified. It is indicated that the most reliable sand production signal can be obtained by installing the sensor at twice the pipeline diameter in the downstream of the elbow. The simulation research results provide the important theoretical basis for the installation of the sensor in the process of field application, so as to ensure the reliability of sand production monitoring data.
-
表 1 砂粒目数与直径的对应关系
Table 1. Corresponding relationship between mesh and diameter of sand grain
目数 砂粒直径/μm 目数 砂粒直径/μm 80 180 200 75 100 150 325 45 150 106 表 2 不同安装位置采集的出砂监测信号特征值对比
Table 2. Comparison between the characteristic values of the sand production monitoring signals acquired at different installation positions
参数 弯管处 弯管下游2倍处 弯管下游3倍处 最大时域幅值 0.292 3 0.633 7 0.425 9 最小时域幅值 −0.301 8 −0.641 8 −0.425 8 方差 2.420 0×10−4 8.894 3×10−4 4.468 1×10−4 平均值 −5.902 8×10−21 −9.635 0×10−21 −7.623 3×10−21 均方根值 0.015 6 0.029 8 0.021 1 功率谱值 2.42×10−4 8.90×10−4 4.47×10−4 -
[1] SHAIEK S, ANRES S, VALDENAIRE S T. Sand management in subsea produced water separation unitreview of technologies and tests[R]. OMC 2015493, 2015. [2] HAMZA N, BAGHDADI F, YAHIA Z, et al. A new concept of sand management philosophy; revolutionary & cost effective solution[R]. SPE 186390, 2017. [3] GHARAIBAH E, ZHANG Y. Sand management and erosion prediction models for oil and gas applications-experimental and numerical validation studies[R]. SPE 181737, 2016. [4] 车争安, 孟召兰, 王尔钧, 等. 适度出砂技术在海上疏松低渗油田适应性实验评价[J]. 钻采工艺, 2017, 40(2):45-47. doi: 10.3969/J.ISSN.1006-768X.2017.02.14 CHE Zheng’an, MENG Zhaolan, WANG Erjun, et al. An experimental evaluation on the adaptability of sand management in loose and low permeability offshore oilfield[J]. Drilling and Production Technology, 2017, 40(2): 45-47. doi: 10.3969/J.ISSN.1006-768X.2017.02.14 [5] 武广瑷, 刘刚, 曹砚锋, 等. 稠油油田出砂地面实时监测技术[J]. 石油钻采工艺, 2016, 38(4):519-525. WU Guang’ai, LIU Gang, CAO Yanfeng, et al. Ground real-time monitoring on sand production of heavy oil reservoir[J]. Oil Drilling and Production Technology, 2016, 38(4): 519-525. [6] WANG Kai, LIU Zhiguo, LIU Gang, et al. Vibration sensor approaches for sand detection in oil-water-sand multiphase flow[J]. Powder Technology, 2015, 276: 183-192. doi: 10.1016/j.powtec.2015.02.041 [7] WANG Kai, LIU Zhiguo, LIU Gang, et al. Vibration sensor approaches for the monitoring of sand production in Bohai Bay[J]. Shock and Vibration, 2015: 1-6. [8] 刘刚, 贾宗文, 刘志国, 等. 基于高频振动信号分析的稠油出砂量监测方法[J]. 中国石油大学学报(自然科学版), 2011, 35(3):84-88. LIU Gang, JIA Zongwen, LIU Zhiguo, et al. Sand content monitoring method of heave oil wells based on high-frequency vibration signal analysis[J]. Journal of China University of Petroleum (Edition of Natural Science), 2011, 35(3): 84-88. [9] 董建. 油气混输泵压水室设计与流场模拟分析[D]. 兰州: 兰州理工大学 2010: 21-33. DONG Jian. Water-compressing chamber design and flow simulation analysis for oil-gas mixture pumps[D]. Lanzhou: Lanzhou University of Technology, 2010: 21-33. [10] 车得福, 李会雄. 多相流及其应用[M]. 西安: 西安交通大学出版社, 2007: 1-7. CHE Defu, LI Huixiong. Multiphase flow and its application[M]. Xi’an: Xi’an Jiaotong University, 2007: 1-7. [11] 于勇. 两相流动气体湍流变动模型和稠密两相湍流模型的研究[D]. 北京: 清华大学, 2004 YU Yong. Studies on gas turbulence modification in gas-particle flows and a dense two-phase turbulence model[D]. Beijing: Tsinghua University, 2004 [12] 曾卓雄. 稠密两相流动湍流模型及其应用[M]. 北京: 机械工业出版社, 2012: 95-133. ZENG Zhuoxiong. Turbulent model for deense two-phaseflow and its application[M]. Beijing: China Machine Press, 2012: 95-133 [13] 杜强, 李洋, 曾祥国. 数值模拟油气管线弯管处固液两相流场特性及冲刷腐蚀预测[J]. 腐蚀与防护, 2017, 38(10):15-19. DU Qiang, LI Yang, ZENG Xiangguo. Prediction of flow field characteristics of solid-liquid phases and erosion corrosion in elbows of oil-gas pipelines using numerical simulation[J]. Corrosion and Protection, 2017, 38(10): 15-19. [14] 杨从新, 董富弟, 程效锐, 等. 含沙水流中粒径对离心泵叶轮磨损特性影响的数值分析[J]. 兰州理工大学学报, 2014, 40(4):45-50. doi: 10.11662/j.issn.1673-5196.2014.04.011 YANG Congxin, DONG Fudi, CHENG Xiaorui, et al. Numerical analysis of effect of sand particle size of silt-laden water flow on wear characteristics of centrifugal pump impeller[J]. Journal of Lanzhou University of Technology, 2014, 40(4): 45-50. doi: 10.11662/j.issn.1673-5196.2014.04.011 -