葉 恒 李光越 劉 釗 張博越 劉家樂

1.中國石油天然氣股份有限公司北京油氣調(diào)控中心 2.中國地質(zhì)大學(xué)（北京） 3.中石油管道有限責(zé)任公司

0 引言

在管道投產(chǎn)實(shí)踐中，合理的注氮量是保證生產(chǎn)安全的關(guān)鍵。崔茂林等[1]根據(jù)自身多年參與管道投產(chǎn)經(jīng)驗(yàn)，提出中、大口徑天然氣管道投產(chǎn)中注氮量的經(jīng)驗(yàn)公式，為管道管容的6.00%～12.00%。付春麗[2]、付先惠等[3]各自研究了管道投產(chǎn)過程中氮?dú)饣鞖舛蔚倪\(yùn)移規(guī)律，發(fā)現(xiàn)混氣段長(zhǎng)度與管道長(zhǎng)度、管徑大小成正比。蒲麗珠等[4]統(tǒng)計(jì)了多條長(zhǎng)輸天然氣管道投產(chǎn)工程數(shù)據(jù)，一般現(xiàn)場(chǎng)注氮量為管道管容的7.00%～20.00%。陳傳勝等[5]根據(jù)川氣東送天然氣管道的投產(chǎn)經(jīng)驗(yàn)，推導(dǎo)出適合大口徑、高壓力管線投產(chǎn)的注氮量，同時(shí)提出在天然氣置換前管道內(nèi)氮?dú)夥獯婀苋荼葢?yīng)大于12.00%。在已進(jìn)行的各條管道投產(chǎn)實(shí)踐中，往往為了保證安全而注入過多氮?dú)?，?dǎo)致了大量的氮?dú)饫速M(fèi)，造成了不必要的經(jīng)濟(jì)損失[6]。

中俄東線天然氣管道是我國首條采用外徑1 422 mm大管徑、12 MPa壓力、X80高鋼級(jí)的長(zhǎng)輸天然氣管道，也是單管輸量最大的跨境長(zhǎng)輸天然氣管道。管道起自俄羅斯東西伯利亞，由布拉戈維申斯克進(jìn)入我國黑龍江省黑河，終點(diǎn)至上海市。其中，我國境內(nèi)段新建管道3 371 km，利用已建管道1 740 km，輸氣量將逐步提升至每年380h108m3。中俄東線北段已于2019年12月2日成功投產(chǎn)通氣，其中包含715 km的干線管道（外徑1 422 mm）以及 109 km長(zhǎng)嶺—長(zhǎng)春支線管道（外徑1 016 mm）[7-8]。為了保證安全投產(chǎn)同時(shí)避免氮?dú)獾睦速M(fèi)，開展對(duì)輸氣管道投產(chǎn)過程中氣體運(yùn)移規(guī)律以及注氮量的優(yōu)化研究，具有重要的意義。

為此，采用計(jì)算流體力學(xué)（CFD）數(shù)值模擬的方法，用中俄東線現(xiàn)場(chǎng)投產(chǎn)實(shí)際數(shù)據(jù)驗(yàn)證了該方法模型的可靠性，在此基礎(chǔ)上研究了外徑1 422 mm大口徑管道投產(chǎn)過程中各氣體的運(yùn)移規(guī)律及其影響因素，提出了不同條件下的理論最優(yōu)注氮管容比。

1 數(shù)值模型建立與邊界條件

1.1 數(shù)值模型的建立

采用CFD方法，模擬天然氣置換（以下簡(jiǎn)稱置換）過程，其二維幾何模型如圖1所示。管道長(zhǎng)L2=1 000 m，其中氮?dú)夥獯骈L(zhǎng)度為L(zhǎng)1，外徑為1 422 mm，壁厚為28 mm。天然氣從左側(cè)進(jìn)入，推著預(yù)先封存的氮?dú)飧綦x段向右側(cè)出口流動(dòng)。當(dāng)右側(cè)出口檢測(cè)到純天然氣時(shí)（甲烷體積分?jǐn)?shù)大于80%），認(rèn)為天然氣置換結(jié)束。由于該模擬尺寸較大，為保證計(jì)算收斂準(zhǔn)確的同時(shí)控制計(jì)算時(shí)間，整體采用四邊形結(jié)構(gòu)化網(wǎng)格，在縱向上采用漸變性網(wǎng)格，越靠近管壁處網(wǎng)格越密，管道中心區(qū)域相對(duì)稀疏，共劃分14h20 000個(gè)網(wǎng)格單元。由于本文研究的是外徑1 422 mm超大管徑，氮?dú)夂吞烊粴獾南鄬?duì)分子質(zhì)量相差較大，在建模加入了重力的影響因素。此外，中俄東線北段管道處于東北地區(qū)，沿途平坦，地形起伏低于200 m。因此本文考慮的重力為水平敷設(shè)管道垂直方向所受到的引力。

圖1 天然氣置換過程的模擬二維幾何模型示意圖

天然氣置換過程涉及流動(dòng)、擴(kuò)散傳質(zhì)等多個(gè)過程，其中流動(dòng)包括單一和混合氣體的流動(dòng)，擴(kuò)散傳質(zhì)過程為氣體間的擴(kuò)散和流動(dòng)。每一個(gè)物理過程都可以用一個(gè)或一組控制方程來表示。氣體流動(dòng)過程可由動(dòng)量方程和連續(xù)性方程來表示，擴(kuò)散傳質(zhì)過程可由組分輸運(yùn)方程來描述。

1）連續(xù)方程：

式中v表示氣體的速度矢量，m/s；ρ表示密度，kg/m3；t表示時(shí)間，s；div表示散度算法。

2）動(dòng)量方程：

在坐標(biāo)系中i方向上的動(dòng)量守恒方程為（j方向同理）：

式中i、j為代表坐標(biāo)系中方向的單位向量；ui、uj表示流體在i、j方向的速度分量，m/s；p表示流體壓力，Pa；xi、xj表示i、j方向上的坐標(biāo)，m ；τij表示黏性應(yīng)力，Pa；μ表示動(dòng)力黏度，Pa·s；Fi表示i方向的單位質(zhì)量力，m/s2；δij表示克羅內(nèi)克符號(hào)，當(dāng)i=j時(shí)，δij=1，否則，δij＝0；λ表示流體的第二黏度系數(shù)，一般取－2/3。

3）組分方程：

式中φi表示i組分的體積分?jǐn)?shù)；Di表示i組分的擴(kuò)散系數(shù)，m2/s；Ri表示單位時(shí)間、體積下產(chǎn)生i組分的質(zhì)量，kg/(m3·s)；grad表示梯度算法。

1.2 邊界條件

模型初始條件：將氮?dú)庖砸欢ū壤獯嬖诠芏蔚那岸?，封存壓力?.02 MPa。天然氣置換過程中，天然氣從左側(cè)進(jìn)入管道置換氮?dú)?。將管線的出口設(shè)置成壓力出口與大氣連通，將天然氣進(jìn)入管線的進(jìn)口設(shè)置成速度入口。選擇基于壓力的分離式求解器，壓力速度耦合選擇SIMPLE算法求解，動(dòng)量方程、能量方程均采用二階迎風(fēng)離散格式，針對(duì)氣體間存在相互擴(kuò)散融合現(xiàn)象，采用Fluent軟件中組分輸運(yùn)模型進(jìn)行求解[9-11]，湍流模型選用k-ε模型，松弛因子介于0.3～0.5，步長(zhǎng)為10－3s，最大迭代步數(shù)為100。

2 數(shù)值模擬研究

2.1 模型的可靠性驗(yàn)證

考慮到利用CFD方法模擬實(shí)際管道長(zhǎng)度（715 km）的龐大計(jì)算量問題，對(duì)長(zhǎng)度為1 000 m、外徑為1 422 mm的管道進(jìn)行了模擬，同時(shí)為了驗(yàn)證1 000 m的數(shù)值模型的可靠性，利用中俄東線現(xiàn)場(chǎng)實(shí)際數(shù)據(jù)進(jìn)行驗(yàn)證。

理論最優(yōu)注氮管容比（本文的管容比均為物理管容比）定義為天然氣置換過程中，置換至管道末端時(shí)純氮?dú)舛蝿偤盟p為零的極限情況。數(shù)值上等于初始注氮管容比減去置換至末端時(shí)殘余純氮?dú)夤苋荼取?/p>

中俄東線北段現(xiàn)場(chǎng)投產(chǎn)實(shí)踐中，初始注氮管容比為8.00%，封存壓力為0.02 MPa。實(shí)際置換速度的平均值為7 m/s，現(xiàn)場(chǎng)投產(chǎn)實(shí)踐中，采用調(diào)整投產(chǎn)調(diào)節(jié)閥開度的方法控制置換速度，盡量保持速度恒定。當(dāng)檢測(cè)到氮?dú)怏w積分?jǐn)?shù)低于95.00%時(shí)認(rèn)為純氮?dú)舛蜗?，天然氣—氮?dú)饣鞖獾竭_(dá)，現(xiàn)場(chǎng)數(shù)據(jù)分析得到理論最優(yōu)注氮管容比為3.50%。

采用上述1 000 m的數(shù)值模型進(jìn)行模擬，初始注氮管容比設(shè)置為8.00%，天然氣置換速度設(shè)置為7 m/s。當(dāng)?shù)獨(dú)怏w積分?jǐn)?shù)低于95.00%時(shí)認(rèn)為純氮?dú)舛蜗В玫搅瞬煌脫Q距離時(shí)（純氮?dú)忸^前端的位置），管道內(nèi)純氮?dú)舛喂艽姹鹊淖兓厔?shì)（圖2）。如圖2所示，置換至管道末端時(shí)殘余純氮?dú)舛喂苋荼葹?.02%。根據(jù)定義，理論最優(yōu)注氮管容比為3.98%，相比于同工況下中俄東線現(xiàn)場(chǎng)數(shù)據(jù)3.50%，絕對(duì)誤差為0.48%，說明模型相對(duì)合理。

根據(jù)《天然氣管道運(yùn)行規(guī)范：SY/T 5922ü2012》[12]，為保證天然氣置換過程中的安全，置換過程中氮?dú)飧綦x段的氮?dú)怏w積分?jǐn)?shù)應(yīng)大于98.00%，大于該值則認(rèn)為是純氮?dú)舛?。若管道?nèi)不存在氮?dú)怏w積分?jǐn)?shù)大于98.00%的氮?dú)舛?，則認(rèn)為天然氣已經(jīng)突破。由于該標(biāo)準(zhǔn)比中俄東線北段現(xiàn)場(chǎng)實(shí)踐中氮?dú)怏w積分?jǐn)?shù)大于95.00%的標(biāo)準(zhǔn)更加嚴(yán)格，因此考慮安全因素，下面的數(shù)值模擬研究將按照98.00%的標(biāo)準(zhǔn)進(jìn)行。

圖2 外徑1 422 mm管道置換過程純氮?dú)舛喂苋荼茸兓厔?shì)模擬結(jié)果圖

2.2 重力對(duì)置換過程氮?dú)夥植家?guī)律的影響

為研究重力的影響，模擬在氮?dú)夥獯娑喂苋荼葹?.00%、置換速度為5 m/s條件下，4種常見的管徑（外徑711 mm、外徑1 016 mm、外徑1 219 mm、外徑1 422 mm）的置換過程，模擬結(jié)果如圖3所示，分別表示置換至管道模型前段（0～120 m）、中段（350～600 m）和末段（750～1 000 m）對(duì)應(yīng)的置換前期（20 s）、置換中期（100 s）、置換后期（180 s）時(shí)管道內(nèi)氮?dú)鉂舛确植紝?duì)比云圖。

如圖3所示，對(duì)于外徑711 mm、外徑1 016 mm管道，重力的影響較小，置換過程中，天然氣總體表現(xiàn)出沿管道中心線“錐形”突進(jìn)的現(xiàn)象。對(duì)于外徑1 219 mm、外徑1 422 mm管道，重力影響較大，天然氣置換過程中，天然氣沿著管段的頂部突進(jìn)，而在管段的底部發(fā)生“延后”現(xiàn)象。基于此特點(diǎn)，對(duì)于中俄東線外徑1 422 mm的管道，重力的影響不可忽略。因此，后續(xù)模擬均考慮重力的影響。

2.3 不同氮?dú)夥獯媪繉?duì)置換過程氮?dú)夥植家?guī)律的影響

模擬外徑1 422 mm管道不同氮?dú)夥獯媪康那闆r，采用單因素變量法，設(shè)計(jì)了初始注氮管容比（L1）分別為6.00%、7.00%、8.00%、9.00%、10.00%的5種工況進(jìn)行模型。參考《天然氣管道運(yùn)行規(guī)范：SY/T 5922ü2012》[12]，天然氣置換速度不宜超過5 m/s，將各工況天然氣置換速度均設(shè)置為5 m/s。

圖3 不同管徑管道置換過程氮?dú)鉂舛确植紝?duì)比云圖

圖4 不同注氮量條件下置換過程不同時(shí)間、不同位置的氮?dú)鉂舛确植紝?duì)比云圖

圖5 不同注氮量條件下置換過程純氮?dú)舛喂苋荼茸兓厔?shì)圖

2.3.1 不同工況下置換過程中氮?dú)夥植家?guī)律分析

對(duì)上述5種工況的模擬結(jié)果，選取天然氣置換至管道模型前段（0～120 m）、中段（350～600 m）和末段（750～1 000 m）對(duì)應(yīng)的置換前期（20 s）、置換中期（100 s）、置換后期（180 s）時(shí)管道內(nèi)氮?dú)鉂舛确植紝?duì)比云圖（圖4）。

如圖4所示，在置換前期（20 s），各工況均存在100.00%氮?dú)舛?；在置換中期（100 s），初始注氮管容比6.00%工況下氮?dú)庾畲篌w積分?jǐn)?shù)開始下降至99.60%，其余各工況均存在100.00%氮?dú)舛?；在置換后期（180 s），初始注氮管容比6.00%、7.00%工況下氮?dú)庾畲篌w積分?jǐn)?shù)分別降至91.60%和96.40%，低于98.00%的標(biāo)準(zhǔn)，說明天然氣已經(jīng)突破。其余各工況均存在100.00%氮?dú)舛巍?/p>

2.3.2 不同工況下置換過程中純氮?dú)舛巫兓厔?shì)研究

不同注氮量條件下隨著置換的進(jìn)行純氮?dú)舛喂苋荼鹊淖兓厔?shì)如圖5所示。

如圖5所示，各工況下曲線趨勢(shì)幾乎一致，隨著置換時(shí)間的增加，純氮?dú)舛喂苋荼入S之遞減，整體遞減速率先增大后減小。對(duì)于各工況，初始注氮管容比為6.00%、7.00%工況下純氮?dú)舛畏謩e在置換140 s、160 s時(shí)耗盡，而初始注氮管容比為8.00%、9.00%、10.00%的工況下在置換結(jié)束（188 s）后仍保有純氮?dú)舛?，分別為0.40%、2.10%、3.50%。因此，在5 m/s的置換速度下，對(duì)應(yīng)理論最優(yōu)注氮管容比為7.60%，約等于8.00%。低于該值的初始注氮管容比是不安全的，會(huì)導(dǎo)致天然氣與空氣的接觸，造成極大的安全隱患[13-15]。

2.4 不同置換速度對(duì)置換過程氮?dú)夥植家?guī)律的影響

天然氣在進(jìn)入管道后流速會(huì)發(fā)生變化。當(dāng)管道出口壓力恒定（大氣壓）時(shí)，管道入口某一恒定的置換速度對(duì)應(yīng)一種固定的流速動(dòng)態(tài)分布，因此，可用置換速度來代表某種流速動(dòng)態(tài)分布的情況。此外，若已知置換速度可求出對(duì)應(yīng)的流量，從而指導(dǎo)現(xiàn)場(chǎng)進(jìn)行調(diào)節(jié)閥的設(shè)置，工程上也具有實(shí)際意義。綜合考慮，根據(jù)單因素變量法，設(shè)置初始注氮管容比為8.00%，針對(duì)不同的天然氣置換速度（3 m/s、5 m/s、7 m/s、9 m/s）進(jìn)行模擬研究，其他條件不變。

2.4.1 不同工況下置換過程中氮?dú)夥植家?guī)律分析

圖6為不同工況置換100 s時(shí)氮?dú)鉂舛确植紝?duì)比云圖。如前所述，由于重力的影響，天然氣沿著管段的上部突進(jìn)，而管段底部發(fā)生“延后”現(xiàn)象。隨著置換速度的增加，這種不均勻分布的現(xiàn)象越來越弱，逐步呈現(xiàn)出“錐形”突進(jìn)的現(xiàn)象。該規(guī)律也是外徑1 422 mm管道投產(chǎn)置換時(shí)的重要特點(diǎn)。

圖6 置換100 s后氮?dú)鉂舛确植紝?duì)比云圖

2.4.2 不同工況下置換過程純氮?dú)舛巫兓厔?shì)研究

各工況置換速度不同導(dǎo)致各工況置換完成的時(shí)間不同，當(dāng)純氮?dú)舛吻熬壉恢脫Q至對(duì)應(yīng)距離時(shí)，不同工況下純氮?dú)舛喂艽姹鹊内厔?shì)變化規(guī)律如圖7所示。

由圖7可以看出，當(dāng)置換速度為3 m/s時(shí)，純氮?dú)舛卧谥脫Q至600 m時(shí)消失，天然氣與空氣接觸，應(yīng)避免該工況。置換速度越快，置換至相同距離時(shí)管道內(nèi)純氮?dú)舛喂苋荼仍酱?，單調(diào)遞增。為找出是否存在拐點(diǎn)，增加了天然氣置換速度為15 m/s和30 m/s的兩種工況，模擬結(jié)果表明，單調(diào)遞增的趨勢(shì)不變，但這兩種工況純氮?dú)舛喂苋荼茸兓厔?shì)已趨于重合。

由于更大的流速是不經(jīng)濟(jì)的，因此可以得到結(jié)論：投產(chǎn)時(shí)，對(duì)于外徑1 422 mm管道，在經(jīng)濟(jì)流速范圍內(nèi)，置換相同距離時(shí)，純氮?dú)舛喂苋荼入S著天然氣置換速度的增加而增加，且增加速度越來越慢，最后幾乎趨于一個(gè)極大值。為保證天然氣投產(chǎn)置換過程的安全，防止氮?dú)獾睦速M(fèi)，應(yīng)在合理范圍內(nèi)盡量增加天然氣置換速度；對(duì)于單一工況（單一置換速度），隨著置換的進(jìn)行，純氮?dú)舛喂苋荼炔粩嘞陆?，但下降速度越來越慢；置換速度為5 m/s、7 m/s、9 m/s、15 m/s、30 m/s的工況下，置換結(jié)束仍殘留有純氮?dú)舛?，?duì)應(yīng)的理論最優(yōu)注氮量為管容的7.60%、5.00%、4.50%、4.00%、4.00%。

圖7 不同工況下純氮?dú)舛喂苋荼茸兓厔?shì)圖

2.4.3 不同工況下置換過程混氣段變化趨勢(shì)研究

混氣段分為天然氣—氮?dú)饣鞖舛魏偷獨(dú)狻諝饣鞖舛?，分別發(fā)生在天然氣和氮?dú)?、氮?dú)夂涂諝獾慕唤缣帯?/p>

由置換前端往后推算，各個(gè)氣體段的判斷標(biāo)準(zhǔn)如下：當(dāng)氧氣體積分?jǐn)?shù)低于18.00%，則認(rèn)為氮?dú)狻諝饣鞖舛我训竭_(dá)；當(dāng)?shù)獨(dú)怏w積分?jǐn)?shù)大于98.00%，則認(rèn)為純氮?dú)舛我训竭_(dá)；當(dāng)?shù)獨(dú)夂康陀?8.00%時(shí)，認(rèn)為天然氣—氮?dú)饣鞖舛我训竭_(dá)；當(dāng)甲烷體積分?jǐn)?shù)大于80.00%則認(rèn)為純天然氣段已到達(dá)。

針對(duì)上述天然氣未突破的工況，對(duì)外徑1 422 mm管道天然氣置換過程中的天然氣—氮?dú)饣鞖舛我约暗獨(dú)狻諝饣鞖舛蔚淖兓厔?shì)進(jìn)行分析，結(jié)果如圖8、9所示。結(jié)果表明，當(dāng)置換相同距離時(shí)，天然氣—氮?dú)饣鞖舛喂苋荼入S著天然氣置換速度增加，均呈現(xiàn)出遞減的趨勢(shì)，且遞減速度越來越慢，最后幾乎趨于一個(gè)極小值；而氮?dú)狻諝饣鞖舛喂苋荼扰c天然氣置換速度的關(guān)系不明顯，各工況條件下接近；對(duì)于單一工況，隨著置換的進(jìn)行，天然氣—氮?dú)饣鞖舛我约暗獨(dú)狻諝饣鞖舛蔚墓苋荼炔粩嗌仙?，上升速度表現(xiàn)為先增加后減少的趨勢(shì)。

圖8 不同工況下天然氣ü氮?dú)饣鞖舛喂苋荼茸兓厔?shì)圖

圖9 不同工況下氮?dú)猕箍諝饣鞖舛喂苋荼茸兓厔?shì)圖

2.5 機(jī)理分析

上述現(xiàn)象的機(jī)理分析如下。

1）隨著置換速度的增加，一方面，分子擴(kuò)散與對(duì)流擴(kuò)散作用增強(qiáng)，導(dǎo)致天然氣—氮?dú)饣鞖舛喂苋荼仍龆?；另一方面，重力引發(fā)頂端突進(jìn)而造成天然氣非活塞式前進(jìn)的現(xiàn)象減弱，導(dǎo)致天然氣—氮?dú)饣鞖舛喂苋荼葴p少。由于重力影響是外徑1 422 mm管道的投產(chǎn)時(shí)影響氣體運(yùn)移規(guī)律的重要特征。重力的影響更大，總體表現(xiàn)出隨著置換速度的增加，天然氣—氮?dú)饣鞖舛喂苋荼入S之減小。

2）由于氮?dú)夂涂諝庵亓ο嘟ǚ肿淤|(zhì)量分別為28、29），兩者界面處幾乎沒有重力影響，且氮?dú)狻諝饣鞖舛慰拷髿鈮毫σ粋?cè)，不同置換速度對(duì)其影響較小。因此，不同工況下比例接近，變化不明顯。

3）純氮?dú)舛喂苋荼鹊乃p主要是由于兩個(gè)混氣段的增加。置換速度增加，天然氣—氮?dú)饣鞖舛喂苋荼葴p少，氮?dú)狻諝饣鞖舛喂苋荼葞缀醪蛔?。因此，表現(xiàn)出純氮?dú)舛喂苋荼入S置換速度增加而增加的趨勢(shì)。

3 現(xiàn)場(chǎng)投產(chǎn)數(shù)據(jù)分析研究

中俄東線北段于2019年底進(jìn)行了干線715 km的天然氣置換工作，通過投產(chǎn)調(diào)節(jié)閥保持7 m/s的置換速度，總置換時(shí)間為35 h。天然氣置換前進(jìn)行了氮?dú)夥獯?，注氮管容比?.00%，封存壓力為0.02 MPa。置換過程中，根據(jù)投產(chǎn)方案，在各閥室和站場(chǎng)進(jìn)行純氮?dú)舛我约盎鞖舛蔚臋z測(cè)，而純天然氣的檢測(cè)只在5個(gè)站場(chǎng)進(jìn)行。當(dāng)置換終點(diǎn)黑河壓氣站檢測(cè)到純天然氣（甲烷體積分?jǐn)?shù)高于80.00%）時(shí)，天然氣置換工作完畢。需要說明，本次投產(chǎn)實(shí)踐中，當(dāng)檢測(cè)到氮?dú)怏w積分?jǐn)?shù)低于95.00%時(shí)認(rèn)為純氮?dú)舛蜗?，天然氣—氮?dú)饣鞖舛蔚竭_(dá)。與《天然氣管道運(yùn)行規(guī)范：SY/T 5922ü2012》[12]要求的98.00%標(biāo)準(zhǔn)略有出入。

通過現(xiàn)場(chǎng)實(shí)測(cè)氣頭數(shù)據(jù)分析，得到結(jié)論如下。

1）天然氣—氮?dú)饣鞖舛巫兓?guī)律（圖10）：隨著置換時(shí)間的增加，天然氣—氮?dú)饣鞖舛喂苋荼戎鸩皆黾?，最終穩(wěn)定在0.60%左右。從僅有的5個(gè)點(diǎn)的數(shù)據(jù)分析，隨著天然氣置換的進(jìn)行，天然氣—氮?dú)饣鞖舛喂苋荼戎饾u增加，但增加速度變慢最終趨于定值的規(guī)律，與模擬趨勢(shì)一致。

圖10 天然氣ü氮?dú)饣鞖舛喂苋荼茸兓厔?shì)圖

2）純氮?dú)舛巫兓?guī)律（圖11）：開始置換后，由于氮?dú)舛吻岸伺c大氣接觸，所以壓力降低，氣體膨脹；同時(shí)，氮?dú)夥獯鏁r(shí)管道內(nèi)氮?dú)鉂舛却笥?9.00%，而現(xiàn)場(chǎng)氣頭檢測(cè)時(shí)認(rèn)為天然氣—氮?dú)饣鞖忸^到達(dá)的標(biāo)準(zhǔn)為氮?dú)怏w積分?jǐn)?shù)低于95.00%，相當(dāng)于把一部分天然氣—氮?dú)饣鞖舛握J(rèn)為是純氮?dú)舛危瑥亩鴮?dǎo)致純氮?dú)舛喂苋荼仍龃?。綜合導(dǎo)致開始置換后，檢測(cè)到的純氮?dú)舛喂苋荼却笥诘獨(dú)夥獯鏁r(shí)的管容比。因此，忽略初始注氮管容比8.00%，以置換2 h（在第一個(gè)檢測(cè)點(diǎn)檢測(cè)到天然氣/氮?dú)饣鞖忸^）對(duì)應(yīng)的純氮?dú)舛喂苋荼?2.00%作為初始值。隨著置換時(shí)間的增加，純氮?dú)舛蔚墓苋荼认陆?，前期下降較快，11 h降至9.70%，后期相對(duì)平緩，在接下來的18 h僅下降1.30%，最終放空時(shí)，仍殘存有8.50%的純氮?dú)?。與模擬趨勢(shì)一致。

圖11 純氮?dú)舛喂苋荼茸兓厔?shì)圖

3）氮?dú)狻諝饣鞖舛巫兓?guī)律（圖12）：現(xiàn)場(chǎng)測(cè)試氧氣體積分?jǐn)?shù)，低于18.00%認(rèn)為氮?dú)狻諝饣鞖舛我堰_(dá)到。測(cè)試時(shí)由于儀表精度等問題，存在較大誤差，反映出曲線不停振蕩，初始管容比為2.10%，5 h內(nèi)降至1.00%以內(nèi)，隨后在1.00%以內(nèi)振蕩，無較好的規(guī)律。但考慮到后期天然氣—氮?dú)饣鞖舛我约凹兊獨(dú)舛蜗鄬?duì)穩(wěn)定，可以判斷，氮?dú)狻諝饣鞖舛魏笃谝矐?yīng)較為穩(wěn)定。

圖12 氮?dú)猕箍諝饣鞖舛喂苋荼茸兓厔?shì)圖

4）理論最優(yōu)注氮管容比：根據(jù)定義，中俄東線投產(chǎn)實(shí)踐在7 m/s的條件下，對(duì)應(yīng)理論最優(yōu)注氮管容比為3.50%。

4 結(jié)論

1）中俄東線（外徑1 422 mm）投產(chǎn)置換過程中，大管徑對(duì)應(yīng)的重力因素不可忽略。在重力的影響下，天然氣置換時(shí)，天然氣沿著管段的頂部突進(jìn)，而在管段的底部發(fā)生“延后”，造成天然氣非活塞式前進(jìn)現(xiàn)象。作為對(duì)比，外徑711 mm、外徑1 016 mm等相對(duì)較小管徑置換時(shí)，整體則表出現(xiàn)天然氣沿管道中心線“錐形”突進(jìn)現(xiàn)象。

2）對(duì)于中俄東線（外徑1 422 mm），置換速度是影響氣體運(yùn)移規(guī)律的主要因素，在3～30 m/s的置換速度范圍內(nèi)，隨著置換速度的增加，置換至相同距離時(shí)，氮?dú)狻諝饣鞖舛喂苋荼葞缀醪蛔?，天然氣—氮?dú)饣鞖舛喂苋荼认陆?，總體表現(xiàn)出純氮?dú)舛喂苋荼仍黾拥内厔?shì)，且增加速率越來越慢。實(shí)際投產(chǎn)時(shí)為降低氮?dú)獾膿p失，應(yīng)適當(dāng)提高天然氣的置換速度。

3）對(duì)于中俄東線（外徑1 422 mm）， 氮?dú)夥獯鎵毫?.02 MPa條件下，置換速度分別為5 m/s、7 m/s、9 m/s、15 m/s、30 m/s對(duì)應(yīng)的理論最優(yōu)注氮管容比分別為7.60%、5.00%、4.50%、4.00%、4.00%。