天然氣輸氣站控排閘閥失效分析與泄漏路徑仿真

胡恒進 賈海平, 石 磊 王大為 郝 瑾, 尹志福

（1.陜西燃氣集團有限公司，陜西 西安 710016；2.陜西燃氣儲運及綜合利用工程研究中心，陜西 西安 710016；3.陜西燃氣集團富平能源科技有限公司，陜西 渭南 711700；4.陜西省表面工程與再制造重點實驗室，陜西 西安 710065）

0 引言

甲烷是僅次于二氧化碳的第二大溫室氣體，是溫室效應超過二氧化碳數倍的強效溫室氣體，天然氣中甲烷控排是減緩全球變暖的關鍵策略之一，閥門泄漏是造成天然氣輸氣站甲烷泄漏的關鍵因素。閥門是輸氣站場中使用最為廣泛的裝置部件之一，在裝置配管設計中正確選用、使用和維護閥門，對保證輸氣站各裝置正常運行十分重要。在長期生產運行中，因密閉性不好或出現故障，將會對安全平穩輸供氣帶來安全隱患，還會影響氣質達不到管輸要求。

針對輸油/輸氣站閥門排放控制失效問題，研究人員[1-3]研究從閥的結構原理、操作方式和維護保養等方面對排污閥故障原因進行分析以降低排污閥的故障率，并提出了預防閥門保養和維護等有關措施。劉健等[4]和馬玉寶等[5]通過分析輸氣站場閥門失效因素，建立了以閥門失效為事件的失效故障樹，采用多級模糊綜合評價法，開展了失效因素分析和可靠性評價。付安慶等[6]分析了西部某氣田天然氣管路球閥發生開裂泄漏事故的原因，認為螺紋根部存在應力集中和鑄造閥體存在缺陷，導致閥體承壓狀態下沿晶脆性斷裂。

本文對某天然氣輸氣站內甲烷控排失效閘閥進行理化性能檢測與分析，同通過有限元方法對閥套閥芯的腐蝕與磨損導致的密封問題進行了建模與計算，模擬結果對延長控排閘閥的安全使用壽命具有重要現實意義。

1 報廢閘閥失效分析

1.1 基本概況分析

本研究對象是陜西境內天然氣輸氣站內甲烷控排失效閘閥，采用的是中口徑閥門，公稱通徑DN50mm，采用閥體法蘭連接管道，公稱壓力為PN10MPa。閘閥是指啟閉體（閥板）由閥桿帶動閥座密封面作升降運動的閥門，可接通或截斷流體的通道。本項目涉及的閘閥閥體材料是碳素鋼（美國牌號WCB，國內牌號ZG250-485）鑄件。天然氣輸氣站內甲烷控排失效閘閥，內部的閥套閥芯存在腐蝕磨損，導致控排的密封不嚴。

1.2 閘閥檢測與失效分析

1.2.1 成分分析

閥體材質為WCB碳素鋼，一般應用于腐蝕極輕的水、油和氣的介質環境；閥芯閥座采用2Cr13馬氏體不銹鋼，具有高硬度、良好的耐磨性和強度，其化學成分測試結果如表1所示。閥體、閥芯閥座分別對照GB/T 12229-2005和GB/T 1220-2007，均符合標準要求；

1.2.2 顯微組織

通過閥座、閥芯的金相制樣和觀察，其顯微組織為較粗大馬氏體和鐵素體，如圖1所示。

1.2.3 顯微硬度

閥芯閥座材料2Cr13鋼的試驗加載負荷為200gf，保荷時間為15s，其維氏硬度度平均值為193；

1.2.4 形貌表征

（1）宏觀特征

針對永樂輸氣站閘閥進行了切割，觀察其內部結構（閥芯、閥套、閥座等）、腐蝕、磨損狀況，如圖2和圖3所示。

圖2 輸氣站控排閘閥的手輪、閥芯、閥座、閥體的宏觀特征

圖3 輸氣站控排閘閥的閥芯、閥套、閥座的宏觀特征

由圖2觀察，進氣口和排氣口的端面由封面墊圈與輸氣站管線或設備連接，存在一定程度腐蝕。內部進出口通道形成了較嚴重的沖刷腐蝕特征，尤其是彎道區域，但壁厚較大且天然氣為干氣介質，基本不會有服役安全問題。

由圖3觀察，閥蓋下表面存在不同程度的腐蝕，閥芯密封圈良好，閥芯下端部表面有輕微的腐蝕磨損特征，閥套外表面有較嚴重的磨損和腐蝕現象，是由天然氣中含有的極微量的H2O、較高CO2和微量H2S等侵蝕性氣體以及氣體含有極少量的沙土沖擊引起。靠近閥套進出氣口附近區域腐蝕更為明顯，閥座局部發生了顯著的減薄特征，可看到有明顯的沖刷腐蝕痕跡；

（2）微觀特征

圖4為閥座失效樣品兩個部位的微觀形貌（×200倍）。閥座外側面和內表面進行微觀觀察，外表面存在磨損條紋狀，內表面腐蝕程度較重。

圖4 輸氣站控排閘閥的閥座試樣微觀形貌

圖5為閥芯失效樣品兩個部位的微觀形貌（×200倍）。閥芯的內表面存在磨損劃痕，局部腐蝕較嚴重；其外側面損傷劃痕處腐蝕較嚴重。其表面腐蝕產物為碳酸亞鐵和鐵的氧化物等主要物質。

圖5 輸氣站控排閘閥的閥芯試樣微觀形貌

2 閥門泄漏路徑仿真

由于閥體是鑄造成形，閥芯和閥座機構也比較復雜，直接用有限元軟件建立模型比較困難。本文先通過三維軟件建立功能部分模型，再導入到有限元仿真軟件進行計算。

2.1 閘閥泄漏路徑建模

2.1.1 三維結構建模

本輸氣站控排閘閥的實際結構數據如表2所示。針對其主要功能部分的閥體、閥套和閥芯建立三維結構模型如圖6所示，裝配后閥體內部沿中面剖開的剖視圖如圖7所示。

表2 閥芯和閥套模型參數

圖6 控排閘閥閥體、閥套、閥芯的三維建模

圖7 控排閘閥閥體剖面的三維建模

2.1.2 有限元模型建立

為了探究控排閘閥在流固耦合條件下的性能，假設天然氣流速10m/s條件下，基于abaqus軟件的歐拉-拉格朗日流固耦合算法（CEL），對實際工況下天然氣流經對閥套和閥芯腐蝕和磨損的影響規律，將三維閥體模型導入到有限元仿真軟件中，進行各部分材料參數設置與裝配。

（1）材料參數設置

為了保證計算效率，閥體部分采用剛體建模，連接部、法蘭和腔體等與質量點與相耦合，歐拉體形狀根據閥體決定。剛體部分不參與計算，可以不設置材料參數；閥芯和閥套設置為2Cr13不銹鋼材料，具體材料參數如表2所示；

（2）單元類型的選擇及網格劃分

由于計算結果對流體域網格大小較為敏感，為此對整個流體域分區進行網格劃分。閥體、閥芯和閥套采用四面體網格劃分，歐拉體部分用六面體網格劃分。在穩態條件下，隨著網格數量的增加，相對誤差逐漸減小；

（3）接觸類型和邊界條件

在仿真天然氣流經閥體時，氣體對閥體、閥套和閥芯必然會產生一定力學作用，天然氣和閥各部分的作用通過接觸面來實現，需要設置天然氣和閥體之間的接觸方式，同時天然氣自身也相互作用，可選擇通用接觸算法，為了防止氣體穿透選擇硬接觸方式。邊界條件需要固定閥體、閥套和閥芯，給定流體初始速度假設10m/s，定義重力場為9.8m/s2。

2.2 閘閥泄漏路徑仿真分析

在閥門完全打開狀態下，天然氣以高流速流入閥體，沖擊閥套和閥芯后再流出閥體，長時間沖蝕作用會造成閥芯和閥套的磨損和腐蝕，最終造成閥門失效，如圖8所示。

圖8 天然氣流經閘閥模擬過程

（1）閥套磨損與腐蝕分析

天然氣從閥門入口處沖入，首先沖擊到閥套，逐漸形成沖擊應力，隨著時間推移沖擊應力逐漸增大，在0.02s時達到最大；從0.03s開始，最大應力值并沒有出現在閥套壁上，而是出現在閥套底部與閥芯接觸部分，具體受力云圖如圖9所示。

圖9 不同時刻閥套受力云圖

為了進一步分析天然氣對閥套底部進氣口的沖擊和腐蝕，將閥套沿中間截面剖開，進氣口處應力云圖，如圖10所示。從0.03s開始，最大應力值出現在底部進氣口，隨時間推移，進氣口的應力最大值越來越大，應力分布的范圍有擴大趨勢，說明閥門在工作過程中，閥套進氣口受到天然氣沖擊和腐蝕越來越嚴重；

圖10 不同時刻閥套進氣口應力云圖

（2）閥芯磨損與腐蝕分析

閥芯與閥套進氣口是控制天然氣是否能通過閥門的主要功能部件，當閥芯進入閥套進氣孔時，將阻斷天然氣通過閥門，當閥芯遠離閥套進氣孔時，天然氣通過進氣孔流入閥套內部，再流出到下一個管道中，閥門導通。

根據閥門的結構，盡管閥芯遠離進氣孔，閥門打開，但是高速流動的天然氣通過進氣孔時第一時間沖擊到了閥芯錐形體，不同時刻天然氣對閥芯沖擊產生的應力，如圖11所示。在0～0.02s之間時，天然氣還沒通過進氣孔，此時閥芯錐形體上沒有應力產生，但閥芯的底座與閥套接觸，此時閥套正受到天然氣沖擊產生一定變形，從而影響閥芯底座上產生較小的應力。在0.03～0.05s之間時，天然氣通過閥套上的進氣孔進入閥套內部，并直接沖擊上閥芯上的錐形體，隨時間推移，沖擊產生的應力有沿錐形體逐漸擴大的趨勢。

圖11 不同時刻的閥芯受力云圖

為了進一步分析天然氣對閥芯錐形體的沖擊和腐蝕，將錐形體局部放大，如圖12所示。從0.03s時，最大應力值出現在錐形體底部，且應力分布不均勻，隨時間推移，進氣口的天然氣逐漸變得比較穩定，對錐形體的底部的沖擊逐漸均勻，應力最大值越來越大，應力分布的范圍有擴大趨勢。說明閥門在工作過程中，閥芯受到天然氣沖擊和腐蝕越來越嚴重，最終將造成閥門不能完全阻斷天然氣。

2.3 閘閥仿真和實際情況對比分析

仿真獲得閥座進氣孔、閥套內壁和閥芯錐形體的沖擊腐蝕情況與實際閥門中的相應組件進行對比，以及閥體進氣孔進行對比，分別如圖13和圖14所示。在進氣孔璧面、閥套內壁、閥芯錐形體上均有比較嚴重的沖擊和腐蝕，仿真結果與實際情況基本吻合，說明仿真獲得的結果比較復合生產實際，數據可靠。

圖13 閥座進氣孔空的實際與仿真比較

圖14 閥套底部內壁、閥芯錐形體的實際與仿真比較

為了延長閥門的使用壽命，可以通過兩個途徑來實現：①改變進氣孔和閥芯錐形體的結構，使天然氣不能直接沖擊到進氣孔內壁和閥芯錐形體外壁，減小沖擊和腐蝕；②在天然氣直接沖擊的進氣孔和閥芯錐形體上進行表面合金化或表面涂層技術處理，增加其抗腐蝕和抗沖擊損傷能力。

3 結語

（1）進氣口和排氣口的端面存在一定程度腐蝕，內部進出口通道形成了較嚴重的沖刷腐蝕特征，閥芯下端部表面有輕微的腐蝕磨損特征，閥套外表面有較嚴重的磨損和腐蝕現象，閥座局部發生了顯著的減薄特征，是由天然氣中含有的極微量的H2O、較高CO2和微量H2S等侵蝕性氣體以及氣體含含有極少量的沙土沖擊引起；

（2）天然氣對閥套底部進氣口的沖擊和腐蝕，從0.03s開始，最大應力值出現在閥套底部進氣口，隨時間推移其進氣口的應力最大值越來越大，應力分布范圍有擴大趨勢；

（3）在0～0.02s時，閥芯錐形體上沒有應力產生，但閥座與閥套接觸影響到閥座產生了較小應力；在0.03～0.05s時，天然氣進入閥套內部并直接沖擊閥芯錐形體，隨時間推移，沖擊產生應力有沿錐形體逐漸擴大趨勢。