基于NSGA-II算法的TEG脫水工藝能耗分析及參數優化*

王潤平 楊岳鵬 曹建峰

渤海鉆探第四鉆井工程分公司

天然氣脫水工藝是油氣田集輸過程中的核心環節，通過脫水處理，可以降低管輸過程中水合物的生成和堵塞問題，減緩天然氣中酸性物質對管材及設備的腐蝕速率，降低壓縮機的動力能耗，但脫水工藝會消耗大量的能源[1-2]。目前，以三甘醇（TEG）為代表的溶劑吸收法是油氣田系統應用最為廣泛的脫水工藝。因此，在保證脫水工藝適應和滿足現場需求的前提下，對TEG 工藝進行能耗分析和參數優化顯得至關重要。

迄今為止，已有許多學者對TEG 工藝的參數優化進行了研究，陳賡良對TEG 的再生工藝進行了完善，對比了惰性汽提、局部冷凝、減壓蒸餾和Drizo 工藝等對TEG 貧液質量分數和露點降的影響[3]；楊東磊等為減少閃蒸氣和再生氣對環境的影響，從廢氣回收、TEG 富液高壓汽提、TEG 富液低壓汽提等方面進行了工藝優化，在滿足露點的條件下，降低了再沸器能耗[4]；李天斌等分析了TEG 循環量、再沸器溫度和汽提氣流量對露點降的影響，得到了最佳操作點參數[5]。目前研究主要集中在工藝流程改進和系統能耗分析方面[6-7]，但脫水工藝涉及的影響因素較多，以上研究只選取了部分因素進行優化，未涉及主控因素，且干氣露點和能耗的變化趨勢通常相悖，如何實現多目標函數的優化也是函待解決的問題。基于此，采用ASPEN Hysys 軟件模擬天然氣脫水工藝流程，根據實際運行參數對脫水系統影響的強弱篩選優化參數，利用響應曲面法設計實驗方案并建立多目標回歸函數，采用第二代自適應非支配遺傳算法（NSGA-II）對函數進行求解，以期實現脫水深度和綜合能耗的有效平衡。

1 流程模擬與參數分析

1.1 TEG工藝流程模擬

以某氣田集氣站為例，其TEG 脫水工藝流程如圖1 所示。濕天然氣經過濾分離器分離出液體和固體雜質，氣田水進入污水處理系統；之后原料氣進入吸收塔底部，與頂部進入的TEG 貧液充分接觸，進行氣液傳質，脫除氣體中的水分；脫水后干氣進入干氣/貧液換熱器升溫，再經產品分離器分離后外輸。TEG富液從吸收塔底部流出，經過濾器過濾后進入再生塔頂部盤管加熱，然后進入閃蒸罐降壓，分離出的閃蒸氣去原料氣罐，閃蒸后的TEG富液從罐底流出，進入貧/富液換熱器，與再生塔底部流出的TEG 貧液換熱升溫，而后進入再生塔頂部。在再生塔中，TEG富液經精餾、提餾，在塔底重沸器的作用下將富甘醇中的水分蒸出。從重沸器流出的貧TEG 經貧液汽提柱汽提后溢流至TEG緩沖罐，再經換熱器降溫后進入循環泵，加壓至吸收塔頂部完成循環。

圖1 TEG脫水工藝流程Fig.1 TEG dehydration process flow

該TEG 工藝中加入了汽提塔對貧液進行再次提純，提高了流程中TEG 的質量分數，降低了重沸器出口溫度，有利于脫水深度的增加，因此在我國應用較為廣泛，屬于典型脫水工藝流程，具有一定的代表性。

由于ASPEN Hysys 軟件在油氣處理、天然氣凝液回收及液化方面均有較高的準確性，因此以Hysys軟件為平臺建立脫水工藝流程，采用PR方程計算氣液兩相熱力學參數，流程見圖2。

圖2 Hysys工藝流程Fig.2 Hysys process flow

現場隨機選擇一組運行數據輸入已建立的Hysys 模型中進行穩態求解，結果見表1。對比可知，出口干氣露點滿足設計要求（露點≤-10 ℃），與實際值的相對誤差為0.27%，貧TEG 進吸收塔溫度的相對誤差為2.77%，重沸器溫度的相對誤差為1.56%。造成誤差的原因一是該Hysys 模型為理想狀態的吸附-再生模型，二是現場測試參數存在一定的人工誤差。但總體看相對誤差較小，說明所建模型的合理性和可靠性較強，可以用于后續研究。

表1 實際運行值與Hysys模擬值對比Tab.1 Comparison between actual running values and Hysys simulated values

1.2 關鍵參數能耗分析

通過對工藝流程進行分析，發現脫水系統的主要耗能設備為貧液循環泵（Q-100）、再生塔塔底重沸器（Q-101）和汽提氣加熱器（Q-102）。結合前人研究的脫水工藝影響因素[8-9]，在表1 中運行參數的基礎上進行定量分析和篩選優化參數。從Q-100 的能耗看，貧液循環泵的能耗隨原料氣溫度、再生塔底壓力的增加而降低，隨TEG 循環量、汽提氣量和重沸器溫度的增加而增加，隨原料氣壓力的增加先增后減再增。原料氣溫度較高時會使水蒸氣在TEG 溶液中的溶解度降低，進而影響吸收效果和能耗；原料氣壓力較高時會引發吸收塔內氣液接觸不充分，進而增加塔板壓降，影響脫水效果，但壓力對貧液循環泵能耗的影響較小（圖3）。以無量綱化處理后的參數為自變量，以貧液循環泵能耗為因變量，得到敏感性分析結果（圖4）。斜率越大表示該參數的敏感性越強，故影響貧液循環泵能耗的關鍵參數為TEG 循環量，循環量每增加1 kmol/h，能耗增加約950 kJ/h。

圖3 運行參數對貧液循環泵的能耗影響Fig.3 Influence of operating parameters on energy consumption of lean liquid circulating pump

圖4 貧液循環泵能耗敏感性分析Fig.4 Sensitivity analysis of lean liquid circulating pump energy consumption

從Q-101 的能耗看，重沸器的能耗隨汽提氣量、再生塔底壓力的增加而降低，隨原料氣壓力、TEG循環量和重沸器溫度的增加而增加，隨原料氣溫度的增加先緩慢下降后快速下降最后保持穩定（圖5）。同理，以無量綱化處理后的參數為自變量，以重沸器能耗為因變量，得到敏感性分析結果（圖6）。影響重沸器能耗的關鍵參數為重沸器溫度，溫度每增加2 ℃，能耗增加約14 900 kJ/h。

圖5 運行參數對重沸器的能耗影響Fig.5 Influence of operating parameters on energy consumption of reboiler

圖6 重沸器能耗敏感性分析Fig.6 Sensitivity analysis of reboiler energy consumption

從Q-102的能耗看，汽提氣加熱器的能耗隨原料氣溫度、TEG 循環量和重沸器溫度的增加而降低，隨原料氣壓力、汽提氣量和再生塔底壓力的增加而增加（圖7）。同理，以無量綱化處理后的參數為自變量，以汽提氣加熱器能耗為因變量，得到敏感性分析結果（圖8）。影響汽提氣加熱器能耗的關鍵參數為汽提氣量，氣量每增加0.1 kmol/h，能耗增加約740 kJ/h。

圖7 運行參數對汽提氣加熱器的能耗影響Fig.7 Influence of operating parameters on energy consumption of stripping gas heater

圖8 汽提氣加熱器能耗敏感性分析Fig.8 Sensitivity analysis of stripping gas heater energy consumption

2 多目標函數

2.1 建立函數

通過上述分析可知，耗能設備中重沸器的能耗占比最大，且能耗與脫水效果呈正比，即脫水深度越深，能耗越高。為保證結果科學有效，將三個設備的能耗進行歸一處理，引入等量功的概念進行衡量[10]：

式中：Weq為等量功，kJ/h；η為渦輪機效率，取0.85；Qreb為重沸器能耗，kJ/h；Wpump為貧液循環泵能耗，kJ/h；Whot為汽提氣加熱器能耗，kJ/h。

以TEG 循環量x1、重沸器溫度x2和汽提氣量x3為自變量，建立以等量功Weq和出口干氣露點Tdew的多目標優化函數，借助Design Expert 軟件進行三因素分布分析（表2），采用BBD 法構建實驗設計方案，并利用Hysys軟件得到模擬結果（表3）。

表2 因素分布Tab.2 Distribution of factors

表3 BBD實驗方案和模擬結果匯總Tab.3 Summary of BBD experimental scheme and simulation results

將實驗數據采用多元二次回歸方程擬合，得到等量功Weq和出口干氣露點Tdew的回歸函數：

對回歸模型進行方差分析（表4）可知。兩個參數的模型p值均小于0.01，說明回歸模擬的擬合精度較高；信噪比均大于4，說明模型的可靠性和顯著性較強；三個相關系數均較大，說明模型預測值與實際值的吻合性較好，誤差不顯著。

表4 回歸方程方差分析結果Tab.4 Results of regression equation variance analysis

2.2 參數優化

NSGA-Ⅱ算法是帶有精英保留策略的自適應非支配多目標遺傳算法，可保證找到的最優解不會被拋棄，提高了算法收斂性和潛在并行解的搜索性能，算法耗時較短[11]。當等量功Weq達到最優時，出口干氣露點Tdew往往無法達到最優，即兩個目標函數存在相互沖突，不可能有唯一確定解，因此得到折衷方案，此時的解稱為Pareto 集，將Pareto 集構成的曲面稱為Pareto 前沿。對公式（2）和（3）進行多目標優化求解，就是采用NSGA-Ⅱ算法對其中的Pareto集進行求解。

算法步驟如下：①在決策變量的范圍空間內，產生隨機數，形成參數初始值；②輸入目標函數，計算種群個體的適應度，對個體進行快速非支配排序，計算種群的擁擠程度，進行交叉和變異操作，得到下一代種群個體；③以總的進化次數為限制條件，判斷輸出結果是否滿足終止條件，如不滿足，返回步驟①，如滿足，輸出Pareto前沿。

設置種群為100，進化次數10000，交叉因子0.8，得到最佳的Pareto前沿（圖9）。圖中的紅點A為當前操作工況，由Pareto 前沿可以得到該工況下的優化操作點B 和C，以及對應的運行參數和目標函數值（表5）。在保持與當前水露點接近的條件下，Weq降低了8373.19 kJ/h，與優化前相比降低了4.18%；在保持與當前等量功接近的條件下，Tdew降低了1.92 ℃，與優化前相比降低了17.47%。在實際的生產中，可根據圖6 選擇不同的解集以滿足不同工況需求。

表5 優化前后的最優解Tab.5 Optimal solutions before and after optimization

圖9 Pareto前沿Fig.9 Pareto frontier

2.3 結果對比

為驗證NSGA-Ⅱ算法的科學性，利用Hysys 軟件自帶的優化器求解，根據上述條件建立目標函數、決策變量和約束條件，采用Box（黑箱）和Mix（混合）算法求解，兩種算法的適用范圍見表6，決策變量的步長取0.5，優化結果見表7。經過優化，TEG循環量和汽提氣量均有所減小，但重沸器溫度未得到優化，Weq有所降低，同時Tdew有所上升，但能耗優化結果不及NSGA-Ⅱ算法。綜上可見，NSGA-Ⅱ算法在能耗降低及參數優化上具有優越性，3 個運行參數均得到優化，說明NSGA-Ⅱ算法得到是全局最優解。

表6 優化算法適用范圍Tab.6 Application scope of optimization algorithm

表7 Hysys自帶優化器的優化結果Tab.7 Optimization results of Hysys built-in optimizer

3 結論

以Hysys 軟件為平臺，建立了天然氣地面脫水系統模型，對影響脫水系統的運行參數進行了篩選，引入等量功和干氣露點作為目標函數，通過BBD 實驗設計和NSGA-Ⅱ算法實現了工藝參數的多目標優化求解，得到如下結論：

（1）影響貧液循環泵能耗、重沸器能耗和汽提氣加熱器能耗的關鍵參數分別為TEG 循環量、重沸器溫度和汽提氣量。

（2）在保持與當前工況水露點接近的條件下，Weq降低了8373.19 kJ/h，降低了4.18%；在保持與當前工況等量功接近的條件下，Tdew降低了1.92 ℃，降低了17.47%。

（3）采用Hysys 軟件自帶的優化器求解時，TEG循環量和汽提氣量均有所減小，但重沸器溫度未得到優化；NSGA-Ⅱ算法在能耗降低及參數優化上具有優越性，可以得到全局最優解。