通球過程皮碗球運動機理研究綜述

王學成，張雷雷，龔治海，張方曉，王國振，陳俊丞，趙磊*

通球過程皮碗球運動機理研究綜述

王學成1，張雷雷1，龔治海1，張方曉1，王國振1，陳俊丞2，趙磊2*

（1. 東莞市盛源石油化工有限公司，廣東 東莞 523000； 2. 遼寧石油化工大學，遼寧 撫順 113001）

通球過程利用氮氣推動皮碗球，將管道內的汽油或柴油輸送至儲罐，是連接碼頭與廠區完成油品運輸的有效手段。皮碗球是通球過程中的重要設備，在管道清理操作中也可作為清管器使用。 皮碗球的速度決定著通球過程的效率和安全性，受到后方氣體的推力，管壁的摩擦阻力和前方液體的阻力等因素的共同影響，皮碗球在管道中的運動是一個復雜的運動過程。本文對近10年通球過程皮碗球運動機理研究進行綜述，理論研究部分給出了氣體速度，皮碗球速度和阻力的計算公式，實驗研究部分介紹了旁通率對通球效率的影響，模擬研究部分介紹了通球過程中管道內的流場情況。研究結果為相關企業設計效率高，經濟性好的通球清管過程提供參考。同時發現，利用數值模擬方法解決工業生產中的復雜問題具有廣闊前景。

通球過程；皮碗球；油氣輸送；旁路清管；清管器

運輸汽油與柴油的船舶需要在碼頭或港口完成卸油操作，卸下的油品通過輸油管道輸送至廠區儲罐中進行儲存。在油品被注入管道之后，向管道中放置皮碗球設備，皮碗球放置完成后，利用輸氣軟管向管道內輸入氮氣，在氮氣的推動下，皮碗球向前移動，進而推動油品沿管道向儲罐方向移動，這一過程通常被稱作通球過程。除此之外，皮碗球設備也可作為清管器使用，當管道內存在殘存液體時，可以利用皮碗球將其推出管道；在管道除垢時，也可利用皮碗球推動清洗液，將垢層溶解在清洗液中一起推出管道。皮碗球在油品輸送與管道清管過程中均具有重要作用。受到后方氣體的推力，管壁的摩擦阻力和前方液體的阻力等因素的共同影響，皮碗球在管道中的運動是一個復雜的運動過程。在通球過程中，通入氣體的速度和皮碗球的速度決定著通球過程的效率和安全性，速度過低，油品的輸送時間變長；速度過高，管道內的壓力以及流體的沖擊力會對管道與儲罐造成破壞。皮碗球與管壁之間的摩擦阻力又會使皮碗球發生磨損。為了解決這些可能出現的問題，有必要對通球過程和清管過程中，推動氣體速度與皮碗球速度之間的關系，皮碗球受到的阻力以及影響通球和清管效率的因素進行分析，為相關企業設計效率高，經濟性好的通球清管過程。通過調研近十年通球過程的文獻發現，研究方法主要有三種，分別是理論研究，實驗研究與模擬研究，成果多集中于近五年，且目前并沒有關于這些研究的系統性綜述?；诖?，本文對上述三種研究方法的相關成果進行介紹，包括但不限于氣體速度，皮碗球速度和阻力的計算公式，旁通率對通球效率的影響，通球過程中管道內的流場情況等。同時發現，利用數值模擬方法解決工業生產中的復雜問題具有廣闊前景。

1 理論研究

M. Lesani等[1]推導和求解了皮碗球通過液體管道時的二維和三維動力學方程，采用動量方程和能量方程研究了流場對皮碗球運行軌跡的影響。M. Mirshamsi和M. Rafeeyan[2]利用動量方程和狀態方程計算氣體沿管道流動時的密度、速度和壓力等參數，導出了皮碗球運動的微分方程。同時，他們[3]在二維燃氣管道長清管器的動態分析中將清管器分成若干個小部分，依照牛頓第二定律對所有小部分分別列出方程，然后將所有方程相加，最終得到清管器總的運動動力學方程。在龍格庫塔法的基礎上，提出了求解運動微分方程的算法。

上述研究給出了皮碗球的總運動方程，為后續氣體與皮碗球之間關系的研究打下了基礎。在此之后，鄢華士和余恒東[4]對清管器的運行狀態進行了分析：清管器在運行過程中，速度呈現加速、減速的交替變化。清管器運行的最大阻力出現在初始階段，即清管器從初始速度為零加速至勻速這一過程，此時沿程摩阻最大，清管器所需的氣體驅動力最大，對應的氣體流量最大。因此，以清管器初始運行階段作為研究對象，計算出驅動清管器所需的氣體流量，即可作氣源管道設計的依據。清管器所需的驅動力由（1）-（4）聯立計算可得，驅動清管器所需的氣體質量流量則由式（5）計算：

式中：F—清管器受到的氣體驅動力，N；

F—清管器受到的液體阻力（包括管道摩阻及液位高差引起的阻力），N；

F—清管器與管壁的摩擦力，N；

—清管器的質量，kg；

—清管器的重量，N；

Δ—清管器加速時間，s；

—清管器勻速運行速度，m·s-1。

式中：λ—管道的摩阻系數；

ρ—液體化工品的密度，kg·m-3；

L—通球管線的長度（包括管線長度及管件當量摩阻長度），m；

—重力加速度，m·s-2；

H—收球端與發球端之間的液位差壓，Pa。

式中：—清管器皮碗與管內壁的摩擦系數；

—管道內徑，m；

—清管器皮碗數量；

1—清管器中心軸直徑，m；

—清管器的過盈率；

—清管器材料的彈性模量；

—清管器材料的泊桑系數。

式中：P—清管器初始受到的氣體最大驅動壓力，Pa。

式中：Q —氣體質量流量，kg·s-1；

ρ0—氣體標準狀態下的密度，kg·m-3;

P0—氣體標準壓力，取 98 000 Pa;

0—氣體標準溫度，K;

T—氣體實際溫度，K。

王鋒等研究發現[5]皮碗清管器的磨損量與滑動距離、法向載荷、皮碗過盈量、材料的屈服極限等有關。通球管道內的輸送介質不同，所使用的公式也存在差異。Chen等[6]通過對天然氣管道通球實驗所得數據的整理分析，提出了表征皮碗球速度與后方推進氣體速度之間關系的公式。Zhang等[7]研究了皮碗球在丘陵天然氣管道中運行過程的動力學模型。根據響應面法（RSM）的模擬結果，提出了一個預測水平輸氣管道皮碗球穩態速度的公式

在明確了皮碗球與推動氣體之間的關系后，就可以針對管道的實際情況進行皮碗球的選型。白港生等[8]推導出了清管器通過凹陷變形直管段、三通、厚壁及閥門、彎頭時需要滿足的長度、骨架鋼體直徑、皮碗密封間距等基本設計參數，并給出具體的指導公式和經驗值。何畏等[9]通過分析清管器運行過程中碟形皮碗相較于直板皮碗的優勢，優選碟形皮碗作為清管器的皮碗類型；對碟形皮碗進行力學分析并采用有限元軟件 Abaqus 分析得出，過盈量和皮碗唇部厚度與接觸長度呈正相關，皮碗唇部長度與接觸長度呈負相關；接觸應力隨著過盈量和唇部厚度增加而逐漸增大，隨著唇部長度增加而逐漸減小。

2 實驗研究

Chen等[10]采用水平，豎直和傾斜管道相結合的實驗模型，以皮碗球的旁通率為自變量進行多組實驗，對通球過程中管內壓力波動、皮碗球速度以及管道末端排出液體流量進行監測。實驗結果表明，當管內液體體積由2L提高到8L時，旁路管道壓力峰值增加14.4%，無旁路管道壓力峰值增加34%，使用帶旁路的皮碗球裝置可以顯著降低通球過程中管內的壓力波動，同時對管內液體加載量的變化具有更好的適應性。皮碗球兩端的壓差與旁路分數無關，與皮碗球運動過程中的沿程阻力有關。皮碗球在管道中的平均速度隨旁路分數的的增大線性減小。在此之后，Chen等[11]在只包含水平管道的實驗裝置上繼續使用不同旁通率的皮碗球進行實驗，旁通率從0%提高到2%，液體流出速率降低了75%。

此外，相關企業結合已有線路對海上油氣田通球過程進行了有效的實驗研究[12-15]。

3 模擬研究

左施施等[16]利用fluent對清管器模型內部流場進行仿真，獲得不同形狀的旁通閥在清管器工作時的壓力云圖和湍流動能圖。田宏軍等[17]利用fluent進行數值模擬。結果表明: 在旁通狀態下，高速流體通過旁路，對管道內壁產生強剪切作用，流過清管器后，在下游形成渦流區。壓力場分為上游高壓區、清管器中壓區、頭部低壓區和下游穩定區; 旁路上和軸線上最大流體速度隨間隙量的增加而減小，流場壓降及所受的驅動力隨間隙量的增加而急劇減小。Asghar Talbizadeh 和Mohammad Mehdi Keshtkar[18]對通球過程中皮碗球周圍的全湍流流動進行建模，對影響皮碗球兩端壓降的參數進行了研究，確定了這些基本參數在通球過程中的最優狀態。

Li等[19]提出了一種基于移動網格法的皮碗球模型，將皮碗球與流體的接觸面視為離散網格節點，采用精細的運動網格對皮碗球周圍的流場進行了研究，模擬著重考慮了氣體的可壓縮性和通球過程中皮碗球受到的摩擦力。他們[20]為了快速準確地描述皮碗球在通球過程中的運動狀態，首先建立了一個適用于瞬態與穩態條件下的分析模型，然后建立了一個包含管道進口氣體質量流量、皮碗球摩擦力和旁路參數等變量的模型，用以計算通球過程所需時間。在考慮氣體可壓縮性和摩擦剖面變化的基礎上，在模擬程序中添加摩擦段，研究了皮碗球在通球過程中可能出現的卡死和滑移現象。在設計的摩擦段內正摩擦偏差是造成突然減速然后突然停止的主要原因，而當皮碗球進入低摩擦區域時，氣體可壓縮性則成為皮碗球速度波動或卡死的主要原因。皮碗球與管壁之間的配合關系、皮碗球表面的的摩擦系數、管材特性(內徑、內壁粗糙度)和皮碗球速度對摩擦力均有影響。

Liu等[21]采用多相流瞬態模擬軟件對中國普光氣田P102集氣站通球清管過程進行了研究。左施施等人[22]以可控速清管器為例，建立了清管器的流體區域泄流流場模型，通過仿真得到了旁通閥不同開度下清管器的流場壓力分布，得到了旁通閥相對開度與壓差的關系曲線（圖1）。隨著旁通閥開度增大，清管器前后壓差先急劇下降再趨于平緩，對附近流場的湍動能影響逐漸變小，清管器在管道中的運行速度逐漸趨于平穩。為得到一個合理的開度區間來平穩地調節清管器運行速度，擬合了開度與壓差的方程式（6），根據擬合函數可對閥門不同開度時的壓差進行估算。

式中，α定義為旁通閥開度，由式（7）計算可得：

式中:S—旁通閥任意狀態下的泄流面積;

1—旁通閥最大泄流面積。

Liu等[23]采用OLGA和CFD模擬相結合的方法，分析了不同進口質量流量和旁通率對提升管底部壓力、皮碗球速度和液相分布形態等特征參數的影響。模擬結果表明:立管底部最大壓力隨著皮碗球旁通率的增加先減小后基本保持不變，最小壓力和通球時間增大，壓力變化范圍和通球速度減小; 當旁通率從0%增加到2%時，壓力和皮碗球速度波動下降最顯著，而進口質量流量對立管底部壓力的影響不大。過高的旁通率或不足的進口質量流量會降低通球操作的效率，并可能不足以提供皮碗球運動所需的壓差。

在長距離氣液管道通球過程中，皮碗球和管道內壁之間的摩擦會導致皮碗球皮碗的直徑小于管道的內徑，從而產生間隙，間隙的存在對液體的流動有重要影響。為了探究流體在管道中的流動特性，Chen等人[24]基于CFD方法進行了一系列數值模擬。結果表明：間隙越大，管道內殘留液體越多，輸送流體效率越低，此時需要適當增大后方推動氣體速度，進而增大皮碗球速度。

圖1 旁通閥相對開度與壓差的關系[21]

4 結束語

1）利用理論研究得到的公式，可以計算在管道輸送不同介質時，皮碗球的速度，后方推動氣體的流量，皮碗球在運動過程中受到的阻力等。通過分析可知，皮碗球速度的主要影響因素是管道中的氣體速度，其次是旁通面積和皮碗球受到的摩擦阻力，而提高皮碗球清管器的旁通率，是提升通球效率的有效手段。

2）通過搭建實驗模型，利用壓力監測設備和速度監測設備，可以較好的得到通球過程中管道內的壓力變化情況和皮碗球的速度，同時使用透明管道可以清晰看到管內輸送介質的流動情況。實驗結果表明，使用帶旁路的皮碗球裝置可以顯著降低通球過程中管內的壓力波動，同時對管內液體加載量的變化具有更好的適應性；皮碗球兩端的壓差與旁路分數無關，與皮碗球運動過程中的沿程阻力有關；皮碗球在管道中的平均速度隨旁路分數的的增大線性減小。

3）使用數值模擬軟件對通球過程進行仿真模擬可以得到通球過程中管內速度場和壓力場的分布情況，利用流固耦合功能可以分析皮碗球的變形與應力，使用瞬態模型還可以分析某一時刻的相關結果，仿真模擬所得結果與理論和實驗的結果吻合程度好，這些優勢使得模擬方法成為研究通球過程的最佳手段，相關成果在近五年的文獻中得到大量介紹。

［1］LESANI M, RAFEEYAN M, SOHANKAR A. Dynamic Analysis of Small Pig through Two and Three-Dimensional Liquid Pipeline [J]. Journal of Applied Fluid Mechanics, 2012, 5(2): 75-83.

［2］MIRSHAMSI M, RAFEEYAN M. Dynamic Analysis of Pig through Two and Three Dimensional Gas Pipeline [J]. Journal of Applied Fluid Mechanics, 2015, 8(1): 43-54.

［3］MIRSHAMSI M, RAFEEYAN M. Dynamic analysis and simulation of long pig in gas pipeline [J]. Journal of Natural Gas Science and Engineering, 2015, 23:294-303.

［4］鄢華士, 余恒東. 液體化工品管道通球動力氣源設計計算方法研究 [J]. 化學工程與裝備, 2015, (07): 50-52.

［5］王鋒, 聞亞星, 陳金忠. 油氣管道皮碗清管器材料磨損機理研究 [J]. 管道技術與設備, 2019, (01): 40-42.

［6］CHEN J H, HE L M, LUO X M, et al. Characterization of bypass pig velocity in gas pipeline: An experimental and analytical study [J]. Journal of Natural Gas Science and Engineering, 2020, 73.

［7］ZHANG L, ZHOU J W, HE H G. Modeling and Simulation of Pigging for a Gas Pipeline Using a Bypass Pig [J]. Mathematical Problems in Engineering, 2020, 2020.

［8］白港生, 李衛全, 張元, et al. 油氣管道清管器的通過性能設計 [J]. 油氣田地面工程, 2021, 40(04): 71-76.

［9］何畏, 孫磊, 黃飛楊, et al. 清管器聚氨酯皮碗接觸性能分析與結構優選 [J]. 工程塑料應用, 2021, 49(03): 118-123+129.

[10] CHEN J H, LUO X M, ZHANG H L, et al. Experimental study on movement characteristics of bypass pig [J]. Journal of Natural Gas Science and Engineering, 2018, 59:212-223.

[11] CHEN J H, HE L M, LUO X M, et al. Bypass pigging technology on amelioration of pigging-induced liquid flow: An experimental and modelling study [J]. Ocean Engineering, 2020, 198.

[12] 宮鵬群. 海上無人平臺油氣海管內環路通球問題及解決方法 [J]. 石化技術, 2021, 28(09): 42-43.

[13] 呂長亮. 某海上油田混輸柔性軟管首次在線通球分析研究 [J]. 石油和化工設備, 2021, 24(09): 20-22.

[14] 馬駿, 陳啟明, 鄭武. 某海上油田清管球收球筒工藝流程優化改造 [J]. 天津科技, 2022, 49(05): 18-21.

[15] 王玉偉, 郝敬雷. 海洋石油平臺立管陸地通球施工問題探討 [J]. 遼寧化工, 2018, 47(09): 933-935.

[16] 左施施, 張恕遠, 呂光遠, et al. 旁通閥形狀對清管器管內流場的影響研究 [J]. 管道技術與設備, 2019, (04): 59-62.

[17] 田宏軍, 魏云港, 曹宇光, et al. 管道清管器旁通流場三維數值模擬分析 [J]. 石油機械, 2019, 47(01): 120-129.

[18] TALBIZADEH A, KESHTKAR M M. Numerical and experimental study on a bypass pig motion in oil transmission pipeline: a case study [J]. Journal of Petroleum Exploration and Production Technology, 2020, 10(7): 3007-3023.

[19] LI X W, HE L M, LUO X M, et al. Numerical simulation and experimental study of bypass pigging slippage [J]. Ocean Engineering, 2021, 230.

[20] LI X W, HE L M, LUO X M, et al. Transient pigging dynamics in gas pipelines: Models, experiments, and simulations [J]. Ocean Engineering, 2021, 232.

[21] LIU W, HUI L X, SHANG J F, et al. Study on Pigging Process Analyses and Control Methods for Moisture Gathering Lines: Study in Puguang Gas Field, China [J]. Journal of Pipeline Systems Engineering and Practice, 2021, 12(3): 04021011.

[22] 左施施, 張恕遠, 呂光遠, 等. 旁通閥開度對清管器泄流流場影響研究 [J]. 管道技術與設備, 2021, (01): 33-36.

[23] LIU Y T, ZHU X X, SONG J P, et al. Research on bypass pigging in offshore riser system to mitigate severe slugging [J]. Ocean Engineering, 2022, 246:110606.

[24] CHEN S T, ZHANG Y B, WANG X L, et al. Numerical simulation and experiment of the gas-liquid two-phase flow in the pigging process based on bypass state [J]. Ocean Engineering, 2022, 252:111184.

Research on the Movement Mechanism of the Bowl Ball During the Passing Process

1,1,1,1,1,2,2*

（1. Dongguan Shengyuan Petrochemical Co., Ltd., Dongguan Guangdong 523000, China;2. Liaoning Petrochemical University, Fushun Liaoning 113001, China）

The ball passing process uses nitrogen to push the cup ball to deliver the gasoline or diesel in the pipeline to the storage tank.It is an effective means to connect the wharf and the factory to complete the transportation of oil products.The ball is an important equipment in the process of passing the ball, and it can also be used as a pig in the pipeline cleaning operation.The speed of the bowl ball determines the efficiency and safety of the passing process.Affected by the thrust of the gas behind, the frictional resistance of the pipe wall and the resistance of the liquid in the front, the movement of the ball in the pipe is a complex movement process. In this paper, the research on the movement mechanism of the ball during the passing process in the past 10 years was summarized. In the theoretical part, the calculation formulas of gas velocity, cup ball velocity and resistance were given; in the experimental part, the influence of bypass rate on ball passing efficiency was introduced; and in the simulation part, the flow field in the pipeline during ball passing was introduced. The research results can provide some reference for related enterprises to design high efficiency and good economy pigging process. At the same time, it is found that the numerical simulation method has broad prospects for solving complex problems in industrial production.

Ball pass process; Bowl ball; Oil and gas transportation; Bypass pigging; Pig

2022-06-27

王學成（1966-），男，新疆人，研究方向：成品油庫運營管理。

趙磊（1982-），男，講師，博士，儲運動態系統控制與仿真。

TE83

A

1004-0935（2023）01-0117-04