LNG到港交接計量的影響因素及誤差核定

仲從響

（中國石油天然氣股份有限公司天然氣銷售江蘇分公司，江蘇南京 211102）

能源物資是現代社會的根本驅動力，液化天然氣（LNG）是一種清潔能源而被廣泛關注。我國“多煤、缺油、少氣”的資源結構決定了天然氣需要依賴大規模的進口。2020年，我國天然氣進口量超過1億t，其中液化天然氣進口量占比超過7成[1]。在液化天然氣貿易及運輸過程中LGN接收站扮演著重要的角色，其不僅需要對LGN進行接收，更需要對其進行計量核算。一般規則下，LNG的國際貿易遵循到港交接的方式，這更進一步提高了LNG接收站計量誤差的要求。基于此，本文從交接計量方式出發，探究交接中的誤差計算方式，為精準計量提供必要保障。

1 液化天然氣的計量方式與計算規則

1.1 液化天然氣的計量

基于能源資源的不均衡分布，液化天然氣的國際貿易較為頻繁。在長期的國際貿易實踐中，液化天然氣貿易建立了一整套基于接收站的到港貿易流程與模式，并具有相對可靠且得到廣泛認可的計量體系。目前在交易過程中，雙反基本遵循國標組織ISO6578關于冷凍碳氫化合物液體靜態計量程序的要求進行實操[2]。在具體的計量過程中大體分為三個基本過程。

①對卸載體積進行計量。LNG的體積隨著溫度壓強的變化而變化，在具體的卸載體積核定過程中采用的是差量法對其進行計算，即用卸載前體積與卸載后的剩余體積差值來表征卸載體積，在實操層面上是采用隨船安裝的交接計量系統來進行。值得注意的是該過程應該在貨運方、收貨方、審核方以及進出口管理部門等多方的共同見證下進行。②對返氣熱量進行計算。在LNG卸載過程中采用壓力自留的方式進行卸載，進而需要對艙內的壓力進行平衡，此部分平衡均通過LNG蒸發氣來進行。常見的方式是直接計算返倉天然氣的熱量，計算中需要對氣態樣本進行采樣以求得不同組分所占比例，再根據體積比例按照不同物質的熱值求得返倉氣態的熱量[3]。③對卸載熱量進行計算。卸載熱量的計算是根據卸載體積求出卸載液化天然氣在燃燒情況下能夠產生的總熱量，再減去返倉氣體的熱量進行求解。同時，在卸載過程中為了維持LNG倉溫度與壓力，船體需要啟動發動機，發動機的能源消耗也需要計算到熱量扣除量中。

1.2 液化天然氣計量的計算規則

液化天然氣的計算分為三個方面，進而在具體的計算規則上也包含了三個基本步驟。

（1）卸載體積的計算

卸載體積應該根據溫度、壓力以及船體的傾斜角度來進行綜合計算。在實操過程中LNG運輸船屬于專用運輸船的范疇，隨船安裝交接計量系統，該系統通過內置方程算法的方式對溫度、壓力以及船體傾斜角度進行了公式參數設定。進而在實際的交接過程中，監督單位需要對上述指標的系統參數設定進行監督，而不需要進行具體計算。計算過程可以簡化為卸貨前的LNG體積與卸貨后的體積差值。

其中：VLNG表示卸載的實際LNG體積；V1表示卸載前的LNG體積；V2表示卸載后的LNG體積。

（2）返氣熱量的計算

不同物質的熱值不同，在具體熱量計算中需要按照檢測樣品的不同氣態組分所占的質量比例，按照不同物質的熱值對返氣熱量進行求解。為了計算的精準性，一般按照摩爾分數來進行計算，在特定溫度下，返倉氣體的體積按照公式（2）計算[4]。

其中：DLNG代表溫度T下的氣體體積，kg/m3；M為組分的摩爾質量，kg/kmol；V為組分的摩爾分數，%；k1、k2為液體體積縮減的校正系數；Xn為氮氣的摩爾分數，%；Xc為甲烷的摩爾分數，%。值得注意的是溫度應該取液化天然氣卸載過程中的平均溫度。

液化天然氣的組分相對復雜，除了主要成分甲烷之外，還包括乙烷、丙烷、異戊烷等，需要按照上述公式對可檢出的組分進行分別求解，最后采用加權的方式計算出最終的返氣熱量，見式（3）。

其中：G表示發熱量；xi、Mi、Gi分別代表了i組分的摩爾分數、摩爾質量與發熱量。單位分別為%、kg/kmol與 kJ/kg。

（3）交接天然氣熱量的計算

交接天然氣的熱量為交接總熱量扣除返倉天然氣熱量及交接過程中發動機的消耗熱量，即可以表示為如下公式。

其中：Q表示交接天然氣熱量，MJ；VLNG、DLNG、G氣分別表示液化天然氣的交接體積、摩爾質量以及高位熱量；E返倉、E發動機分別代表返倉氣體與發動機的能量損耗。

其中返倉氣體的總消耗按照公式3進行計算；發動機的總消耗按照交接過程中消耗的氣體總量與氣體熱值的乘積來計算，消耗的氣體總量可以根據發動機燃燒爐的進氣儀表直接測量。

2 液化天然氣的誤差核定

通過對液化天然氣到港交接的計算過程分析，不難發現其計算過程較為復雜，且存在多個流程與步驟，進而計算后的結果必然存在一定的誤差。按照不同的來源，誤差可以分為兩種，即天然誤差與測量誤差。天然誤差主要是指摩爾質量、常數等參數在引入公式的過程中精確的有效數字位數帶來的誤差，此種誤差是無法避免的。測量誤差主要是指對參數進行測定過程中可能產生的誤差。如對天然氣倉內液面的測量、發動機燃燒室用氣流量的測量、對氣體組分的測量等。這些過程中由于測量精度的問題其誤差必然存在，但是可以通過多次測量以及提高設備與人工精度的方式來盡量降低。

無論是何種誤差均無法完全消除，進而需要對可能誤差進行經驗性求解，從而達到修正的效果，進一步使得計算值與實際交接值更為貼近。值得注意的是，液化天然氣交接環節較為復雜，計量方式經過了多步驟計算，而原則上每次測量與計算均會產生誤差，便會形成誤差的累積。在多冪累積下誤差會得到放大，這也是針對交接誤差的計算十分必要的原因所在。

在具體的計算過程中，誤差采用聯動累積的方式予以計算，此種計算方式符合液化天然氣交接的計算規律，也與客觀事實相符，具體如公式（5）。

其中：Y為間接測量值誤差；Xi為直接測量值誤差。直接測量值誤差可以采用經驗帶入的方式進行，也可以通過多次測量后進行誤差計算的方式予以確定。在實際操作中固定設備的直接測量誤差變化較小，一般多采用經驗值代入的方式來進行。

多元函數中的增量可以通過全微分進行表示，進而公式（5）中的增量可以表示為不同直接測量過程中全微分加權值，公式（5）可以變形為公式（6）。

在公式（6）中，可以對誤差較小變量的微分量進行替代，從而達到縮減運算量的根本目的，在經過了縮減后，利用誤差傳遞公式可以對不同誤差過程中的比例進行求解，從而獲得變量與最終計算結果（熱量）之間的相關關系。在物理層面上的意義表現為不同直接測量過程對交接熱量的影響比例大小。

從這一角度，誤差的累積最終計算公式如下[5]：

經過整理與參數帶入結果見式（8）：

3 液化天然氣的誤差實證與影響因素

為了進一步探討與認證本文對誤差計算方式的正確性，并對具體的誤差計算結果應用方式進行探究。選擇2021年7月到港的某艘液化天然氣到港船只的實際交接過程為具體的數據來源，對各步驟的計算方式進行實證計算。可直接測量參數如表1所示。

表1 直接測量指標結果統計表

通過直接測量的方式能夠獲得包括溫度、卸貨體積等直接數據，將所獲得的數據代入到上述公式體系中可以對實際的交接結算量進行計算。本實證計算則更關注不同部位的誤差擬合計算的最終結果。通過現場判斷發現誤差主要來源于直接測量數值，具體包括液態體積誤差、返倉溫度誤差、返倉壓力誤差、發動機損耗誤差、甲烷含量誤差等，具體的誤差數值及產生原因如表2所示。

表2 誤差來源及結果統計表

將誤差參數代入公式（8）中，對總誤差進行計算，并分別計算各個誤差在總體誤差中的比例。同時，熱量誤差為最終的誤差結果，在計算過程中以熱量誤差為100%，對其他各誤差影響因素在誤差中的占比進行計算，其結果如表3所示。

表3 誤差所占比例結果統計表

從表2中看出，案例船只的總交易發熱量為35.57億MJ。發熱量的總誤差為57.45萬MJ，占交易量的比例為0.016 2%。

其中各個部分的誤差占比存在較大的差異，按照誤差占比進行排序，其中液態體積測量過程中的誤差占比最高，能夠解釋總誤差中56.07%的變量，甲烷含量測量誤差占比為31.51%。上述兩項指標累積解釋的誤差變量超過85%，在實際操作過程中，為進一步簡便計算流程與結算規則，可以在上述兩者誤差之和的基礎上浮10%的方式來表征總體誤差。其他誤差占比均較小，如返倉氣體溫度誤差占比3.69%，發動機損耗誤差占比1.96%等。

4 總結

液化天然氣到港交接受到溫度、壓力等影響，且計算過程較為復雜。在本文的研究中，最終的熱量誤差在0.015%以上，這對于大宗交易已推高了成本及交易損耗。在此種情況下，本文在總結了計量方式及計算規則的基礎上，對誤差的來源及影響因素進行了系統分析。并通過實例計算的方式探討不同誤差對總誤差的占比，發現液態體積測量過程以及甲烷含量測量過程中產生的誤差最大，分別達到86.07%和31.51%，二者可以共同解釋誤差變量大于85%適宜于作為誤差計算的標志性指標。