<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article article-type="brief-report" dtd-version="2.3" xml:lang="EN" xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Energy Res.</journal-id>
<journal-title>Frontiers in Energy Research</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Energy Res.</abbrev-journal-title>
<issn pub-type="epub">2296-598X</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="publisher-id">1135330</article-id>
<article-id pub-id-type="doi">10.3389/fenrg.2023.1135330</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Energy Research</subject>
<subj-group>
<subject>Brief Research Report</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>GRU-AGCN model for the content prediction of gases in power transformer oil</article-title>
<alt-title alt-title-type="left-running-head">Luo et al.</alt-title>
<alt-title alt-title-type="right-running-head">
<ext-link ext-link-type="uri" xlink:href="https://doi.org/10.3389/fenrg.2023.1135330">10.3389/fenrg.2023.1135330</ext-link>
</alt-title>
</title-group>
<contrib-group>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Luo</surname>
<given-names>Diansheng</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="corresp" rid="c001">&#x2a;</xref>
<uri xlink:href="https://loop.frontiersin.org/people/1741382/overview"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Chen</surname>
<given-names>Wengang</given-names>
</name>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Fang</surname>
<given-names>Jie</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/1730967/overview"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Liu</surname>
<given-names>Jianguo</given-names>
</name>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Yang</surname>
<given-names>Jinbiao</given-names>
</name>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Zhang</surname>
<given-names>Ke</given-names>
</name>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
</contrib>
</contrib-group>
<aff id="aff1">
<sup>1</sup>
<institution>College of Electrical and Information Engineering</institution>, <institution>Hunan University</institution>, <addr-line>Changsha</addr-line>, <country>China</country>
</aff>
<aff id="aff2">
<sup>2</sup>
<institution>State Grid Jincheng Power Supply Company</institution>, <addr-line>Jincheng</addr-line>, <country>China</country>
</aff>
<author-notes>
<fn fn-type="edited-by">
<p>
<bold>Edited by:</bold> <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/1259738/overview">Jian Zhao</ext-link>, Shanghai University of Electric Power, China</p>
</fn>
<fn fn-type="edited-by">
<p>
<bold>Reviewed by:</bold> <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/2051658/overview">Shu Zhang</ext-link>, Sichuan University, China</p>
<p>
<ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/1484038/overview">Tianguang Lu</ext-link>, Shandong University, China</p>
</fn>
<corresp id="c001">&#x2a;Correspondence: Diansheng Luo, <email>lhx20070322@sina.com</email>
</corresp>
<fn fn-type="other">
<p>This article was submitted to Process and Energy Systems Engineering, a section of the journal Frontiers in Energy Research</p>
</fn>
</author-notes>
<pub-date pub-type="epub">
<day>05</day>
<month>04</month>
<year>2023</year>
</pub-date>
<pub-date pub-type="collection">
<year>2023</year>
</pub-date>
<volume>11</volume>
<elocation-id>1135330</elocation-id>
<history>
<date date-type="received">
<day>31</day>
<month>12</month>
<year>2022</year>
</date>
<date date-type="accepted">
<day>27</day>
<month>03</month>
<year>2023</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#xa9; 2023 Luo, Chen, Fang, Liu, Yang and Zhang.</copyright-statement>
<copyright-year>2023</copyright-year>
<copyright-holder>Luo, Chen, Fang, Liu, Yang and Zhang</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/">
<p>This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.</p>
</license>
</permissions>
<abstract>
<p>Aiming at high accuracy of dissolved gas content prediction in transformer oil, a novel method based on gated recurrent unit and adaptive graph convolution network (GRU-AGCN) is proposed. For gated recurrent unit (GRU) can selectively choose the feature of time series, it is used to extract time series information of the gas content. Correlation among gases are extracted to improve the accuracy. The original adjacency matrix of the model is constructed according to the grey relational analysis (GRA), and the dynamic relation information between gases is extracted by adaptive graph convolution network (AGCN). The experimental result shows that the GRU-AGCN model can efficiently extract the temporal features and perceive the dynamic relationship of gases. The predictions error of the proposed method is lower than that of RNN, LSTM network and GRU network. The proposed method provides a reliable and accurate result for the prediction of dissolved gas content in transformer oil.</p>
</abstract>
<kwd-group>
<kwd>power transformer</kwd>
<kwd>dissolved gas</kwd>
<kwd>time series prediction</kwd>
<kwd>gated recurrent unit</kwd>
<kwd>graph convolutional network</kwd>
</kwd-group>
</article-meta>
</front>
<body>
<sec id="s1">
<title>1 Introduction</title>
<p>Transformer is an important equipment in power system, which is related to the safety and stability of power system (<xref ref-type="bibr" rid="B13">Wu et al., 2021</xref>; <xref ref-type="bibr" rid="B12">Su et al., 2022</xref>). During the operation of the power transformer, several gases are generated in oil by the influence of life cycle, load, oil temperature and other factors. When the transformer overheats or discharges, there will be abnormal changes in gases (<xref ref-type="bibr" rid="B5">G&#xf3;mez et al., 2014</xref>). Insulation will be broken down as the transformer fault develops, affecting the stability of the power system (<xref ref-type="bibr" rid="B6">Huang et al., 2021</xref>). Therefore, the state data of the power transformer can be obtained and its trend can be predicted by installing monitoring equipment on the transformer, which can provide a reliable judgment basis for the state maintenance of the transformer (<xref ref-type="bibr" rid="B14">Xie et al., 2020</xref>).</p>
<p>In recent years, a lot of research work has been done on the prediction of dissolved gas content in transformer oil. <xref ref-type="bibr" rid="B2">Chen et al. (2021)</xref> constructed an improved grey model for the time series prediction of gas in oil. K-nearest Neighbor was used to predict the furan level in transformers in (<xref ref-type="bibr" rid="B11">Shaban et al., 2016</xref>). In addition, many machine learning methods such as artificial neural networks (<xref ref-type="bibr" rid="B4">Ghunem et al., 2012</xref>) and support vector machines (<xref ref-type="bibr" rid="B1">Atherfold and van Zyl, 2020</xref>) are used for gas content prediction. The above methods process the data of each moment synchronously, so the feature significance of different moments cannot be considered, which means that their prediction accuracy is limited.</p>
<p>Recurrent neural network (RNN) extracts current and historical information by loop structure, which is suitable for processing time series. However, RNN is weak in learning information at too far time, which is called as &#x201c;Long-term dependency problem.&#x201d; The Long short-term Memory (LSTM) developed based on RNN introduces gating mechanism to alleviate the problem of RNN. Therefore, some prediction methods for gas content based on LSTM have been proposed (<xref ref-type="bibr" rid="B8">Lin et al., 2018b</xref>; <xref ref-type="bibr" rid="B9">Ma et al., 2021</xref>). But the complexity of LSTM leads to long training time. Gated Recurrent Unit (GRU) network was proposed as another variant of RNN, simplifies the gating mechanism, which can process long series efficiently and be trained fast (<xref ref-type="bibr" rid="B10">Minh et al., 2018</xref>).</p>
<p>There are complex dynamic relationships in gases in transformer oil (<xref ref-type="bibr" rid="B3">Fessler et al., 1989</xref>). When the hot spot temperature is low or the fault is slight, the insulating material produces hydrogen and alkanes through decomposition. When the fault is severe, the energy will be released by overheating or discharge, alkenes and alkynes will be produced by further decomposition. Therefore, the production of dissolved gases in oil has a certain correlation, and varies according to the fault type and development. Above methods are based on the prediction by single gas variable, which ignores the inherent correlations between gases, and the results are obtained with low accuracy and reliability. In the study of (<xref ref-type="bibr" rid="B9">Ma et al., 2021</xref>), temperature and seven gas content series were put into LSTM to predict gas content. <xref ref-type="bibr" rid="B7">Lin et al. (2018a)</xref> took the gas content and the ratios of various gas as the input features of neural network for the prediction. These methods consider the influence of various factors in the prediction, but the feature extraction methods of correlation factors are simple, and the promotion in prediction performance is limited.</p>
<p>For above problems, a content prediction method of gas in transformer oil based on gated recurrent unit and adaptive graph convolution network (GRU-AGCN) is proposed. GRU is used to extract time series information of gas content. To describe the correlation among gases, the original adjacency matrix of the model is constructed based on grey correlation analysis (GRA), and the adjacency matrix is real-time updated according to the input to ensure that the correlation obtained is accurate and reliable. The dynamic relation feature among gases is extracted by adaptive graph convolution network (AGCN). The accurate prediction of gas content can is achieved by integrating the historical information of gas itself and other relevant gas information.</p>
</sec>
<sec id="s2">
<title>2 Gated recurrent unit network</title>
<p>GRU network is one of the variants of RNN. In GRU, gating mechanism is introduced to RNN to selectively transmit temporal features.</p>
<p>The structure of GRU network is shown in <xref ref-type="fig" rid="F1">Figure 1A</xref>. Data [<bold>
<italic>x</italic>
</bold>
<sub>1</sub>, <bold>
<italic>x</italic>
</bold>
<sub>2</sub>, &#x2026;, <bold>
<italic>x</italic>
</bold>
<sub>
<italic>T</italic>
</sub>] with length T are input into corresponding units chronologically. The structure of GRU is shown in <xref ref-type="fig" rid="F1">Figure 1B</xref>. Take GRU of time <italic>t</italic> as an example (1 &#x3c; <italic>t</italic> &#x3c; <italic>T</italic>), the input of the unit includes current original data <bold>
<italic>x</italic>
</bold>
<sub>
<italic>t</italic>
</sub> and hidden feature <bold>
<italic>h</italic>
</bold>
<sub>
<italic>t</italic>&#x2212;1</sub> output by the last unit. And the output of the unit is the hidden feature at time <italic>t</italic> denoted as <bold>
<italic>h</italic>
</bold>
<sub>
<italic>t</italic>
</sub>. The workflow of GRU is as follows, Firstly, input is passed to reset gate and update gate to obtain reset state <bold>
<italic>r</italic>
</bold>
<sub>
<italic>t</italic>
</sub> and update state <bold>
<italic>u</italic>
</bold>
<sub>
<italic>t</italic>
</sub>. Secondly, the candidate hidden state <inline-formula id="inf1">
<mml:math id="m1">
<mml:mrow>
<mml:msub>
<mml:mover accent="true">
<mml:mi mathvariant="bold-italic">h</mml:mi>
<mml:mo>&#x5e;</mml:mo>
</mml:mover>
<mml:mi>t</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> at time <italic>t</italic> is calculated by reset state <bold>
<italic>r</italic>
</bold>
<sub>
<italic>t</italic>
</sub>. Finally, the hidden feature <bold>
<italic>h</italic>
</bold>
<sub>
<italic>t</italic>
</sub> at time <italic>t</italic> is updated by combining the updated state <bold>
<italic>u</italic>
</bold>
<sub>
<italic>t</italic>
</sub> with the hidden feature <bold>
<italic>h</italic>
</bold>
<sub>
<italic>t</italic>&#x2212;1</sub> at time (<italic>t</italic>&#x2212;1) and candidate hidden state <inline-formula id="inf2">
<mml:math id="m2">
<mml:mrow>
<mml:msub>
<mml:mover accent="true">
<mml:mi mathvariant="bold-italic">h</mml:mi>
<mml:mo>&#x5e;</mml:mo>
</mml:mover>
<mml:mi>t</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> at time <italic>t</italic>. The working principle is written as follows: <disp-formula id="e1">
<mml:math id="m3">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">r</mml:mi>
<mml:mi>t</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">s</mml:mi>
<mml:mi mathvariant="normal">i</mml:mi>
<mml:mi mathvariant="normal">g</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">x</mml:mi>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:msub>
<mml:mi mathvariant="bold-italic">W</mml:mi>
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi mathvariant="bold-italic">h</mml:mi>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>1</mml:mn>
<mml:msub>
<mml:mi mathvariant="bold-italic">W</mml:mi>
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi mathvariant="bold-italic">B</mml:mi>
<mml:mi>r</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(1)</label>
</disp-formula>
<disp-formula id="e2">
<mml:math id="m4">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">u</mml:mi>
<mml:mi>t</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">s</mml:mi>
<mml:mi mathvariant="normal">i</mml:mi>
<mml:mi mathvariant="normal">g</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">x</mml:mi>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:msub>
<mml:mi mathvariant="bold-italic">W</mml:mi>
<mml:mrow>
<mml:mi>u</mml:mi>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi mathvariant="bold-italic">h</mml:mi>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>1</mml:mn>
<mml:msub>
<mml:mi mathvariant="bold-italic">W</mml:mi>
<mml:mrow>
<mml:mi>u</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi mathvariant="bold-italic">B</mml:mi>
<mml:mi>u</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(2)</label>
</disp-formula>
<disp-formula id="e3">
<mml:math id="m5">
<mml:mrow>
<mml:msub>
<mml:mover accent="true">
<mml:mi mathvariant="bold-italic">h</mml:mi>
<mml:mo>&#x5e;</mml:mo>
</mml:mover>
<mml:mi>t</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mi mathvariant="italic">tan</mml:mi>
<mml:mo>&#x2061;</mml:mo>
<mml:mi mathvariant="normal">h</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">x</mml:mi>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:msub>
<mml:mi mathvariant="bold-italic">W</mml:mi>
<mml:mrow>
<mml:mi>h</mml:mi>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi mathvariant="bold-italic">r</mml:mi>
<mml:mi>t</mml:mi>
</mml:msub>
<mml:mo>&#x2297;</mml:mo>
<mml:msub>
<mml:mi mathvariant="bold-italic">h</mml:mi>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>1</mml:mn>
<mml:msub>
<mml:mi mathvariant="bold-italic">W</mml:mi>
<mml:mrow>
<mml:mi>h</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi mathvariant="bold-italic">B</mml:mi>
<mml:mi>h</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(3)</label>
</disp-formula>
<disp-formula id="e4">
<mml:math id="m6">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">h</mml:mi>
<mml:mi>t</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi mathvariant="bold-italic">u</mml:mi>
<mml:mi>t</mml:mi>
</mml:msub>
<mml:mo>&#x2297;</mml:mo>
<mml:msub>
<mml:mi mathvariant="bold-italic">h</mml:mi>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mn mathvariant="bold">1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi mathvariant="bold-italic">u</mml:mi>
<mml:mi>t</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x2297;</mml:mo>
<mml:msub>
<mml:mover accent="true">
<mml:mi mathvariant="bold-italic">h</mml:mi>
<mml:mo>&#x5e;</mml:mo>
</mml:mover>
<mml:mi>t</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
<label>(4)</label>
</disp-formula>
</p>
<fig id="F1" position="float">
<label>FIGURE 1</label>
<caption>
<p>Structure of GRU network, <bold>(A)</bold> is GRU network overall structure, <bold>(B)</bold> is GRU internal structure.</p>
</caption>
<graphic xlink:href="fenrg-11-1135330-g001.tif"/>
</fig>
<p>In above equations, <italic>f</italic>
<sub>sig</sub>(.) and <italic>f</italic>
<sub>tanh</sub>(.) are sigmoid function and hyperbolic tangent function, <bold>
<italic>W</italic>
</bold>
<sub>
<italic>r</italic>1</sub> and <bold>
<italic>W</italic>
</bold>
<sub>
<italic>r</italic>2</sub> are the weight of the reset gate, <bold>
<italic>B</italic>
</bold>
<sub>
<italic>r</italic>
</sub> is the bias of the reset gate, <bold>
<italic>W</italic>
</bold>
<sub>
<italic>u</italic>1</sub> and <bold>
<italic>W</italic>
</bold>
<sub>
<italic>u</italic>2</sub> are the weight of the update gate and <bold>
<italic>B</italic>
</bold>
<sub>
<italic>u</italic>
</sub> is the bias of the update gate. <bold>
<italic>W</italic>
</bold>
<sub>
<italic>h</italic>1</sub> and <bold>
<italic>W</italic>
</bold>
<sub>
<italic>h</italic>2</sub> are the weight and <bold>
<italic>B</italic>
</bold>
<sub>
<italic>h</italic>
</sub> is the bias in the candidate hidden state calculation. Hadamard Product operation is denoted as.</p>
<p>By above principle, historical information and current information can be filtered and processed flexibly in GRU network.</p>
</sec>
<sec id="s3">
<title>3 Adaptive graph convolutional network</title>
<sec id="s3-1">
<title>3.1 Grey relational analysis of gases in transformer oil</title>
<p>GRA is an analysis method based on grey system theory, which judges the correlation degree among factors according to the similarity of the shape of the factor curves. Therefore, the historical data of gas 1&#x223c; <italic>n</italic> with length <italic>M</italic> are analyzed by GRA.</p>
<p>Historical data for gas <italic>i</italic> is used as a reference sequence, which is denoted as <inline-formula id="inf3">
<mml:math id="m7">
<mml:mrow>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">X</mml:mi>
<mml:mi>M</mml:mi>
<mml:mi>i</mml:mi>
</mml:msubsup>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mfenced open="[" close="]" separators="|">
<mml:mrow>
<mml:msubsup>
<mml:mi>x</mml:mi>
<mml:mn>1</mml:mn>
<mml:mi>i</mml:mi>
</mml:msubsup>
<mml:mo>,</mml:mo>
<mml:msubsup>
<mml:mi>x</mml:mi>
<mml:mn>2</mml:mn>
<mml:mi>i</mml:mi>
</mml:msubsup>
<mml:mo>,</mml:mo>
<mml:mo>&#x22ef;</mml:mo>
<mml:mo>,</mml:mo>
<mml:msubsup>
<mml:mi>x</mml:mi>
<mml:mi>M</mml:mi>
<mml:mi>i</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>, and take the historical data of gas <italic>j</italic> as a comparison sequence, which is denoted as <inline-formula id="inf4">
<mml:math id="m8">
<mml:mrow>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">X</mml:mi>
<mml:mi>M</mml:mi>
<mml:mi>j</mml:mi>
</mml:msubsup>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mfenced open="[" close="]" separators="|">
<mml:mrow>
<mml:msubsup>
<mml:mi>x</mml:mi>
<mml:mn>1</mml:mn>
<mml:mi>j</mml:mi>
</mml:msubsup>
<mml:mo>,</mml:mo>
<mml:msubsup>
<mml:mi>x</mml:mi>
<mml:mn>2</mml:mn>
<mml:mi>j</mml:mi>
</mml:msubsup>
<mml:mo>,</mml:mo>
<mml:mo>&#x22ef;</mml:mo>
<mml:mo>,</mml:mo>
<mml:msubsup>
<mml:mi>x</mml:mi>
<mml:mi>M</mml:mi>
<mml:mi>j</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>, where <italic>i</italic>, <italic>j</italic>&#x2208;{1, 2, &#x2026;, <italic>n</italic>}, and <italic>i</italic>&#x2260;<italic>j</italic>. Grey correlation coefficient <italic>coeff</italic>
<sub>
<italic>ij</italic>
</sub> and grey relational degree <italic>grey</italic>
<sub>
<italic>ij</italic>
</sub> are calculated as follows,<disp-formula id="e5">
<mml:math id="m9">
<mml:mrow>
<mml:mfenced open="{" close="" separators="|">
<mml:mrow>
<mml:mtable columnalign="center">
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>c</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>f</mml:mi>
<mml:mi>f</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>m</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>min</mml:mi>
<mml:mrow>
<mml:mfenced open="|" close="|" separators="|">
<mml:mrow>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">X</mml:mi>
<mml:mi>M</mml:mi>
<mml:mi>i</mml:mi>
</mml:msubsup>
<mml:mo>&#x2212;</mml:mo>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">X</mml:mi>
<mml:mi>M</mml:mi>
<mml:mi>j</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>&#x3b1;</mml:mi>
<mml:mo>&#x2219;</mml:mo>
<mml:mi>max</mml:mi>
<mml:mrow>
<mml:mfenced open="|" close="|" separators="|">
<mml:mrow>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">X</mml:mi>
<mml:mi>M</mml:mi>
<mml:mi>i</mml:mi>
</mml:msubsup>
<mml:mo>&#x2212;</mml:mo>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">X</mml:mi>
<mml:mi>M</mml:mi>
<mml:mi>j</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:mrow>
<mml:mfenced open="|" close="|" separators="|">
<mml:mrow>
<mml:msubsup>
<mml:mi>x</mml:mi>
<mml:mi>m</mml:mi>
<mml:mi>i</mml:mi>
</mml:msubsup>
<mml:mo>&#x2212;</mml:mo>
<mml:msubsup>
<mml:mi>x</mml:mi>
<mml:mi>m</mml:mi>
<mml:mi>j</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>&#x3b1;</mml:mi>
<mml:mo>&#x2219;</mml:mo>
<mml:mi mathvariant="normal">m</mml:mi>
<mml:mi>a</mml:mi>
<mml:mi>x</mml:mi>
<mml:mrow>
<mml:mfenced open="|" close="|" separators="|">
<mml:mrow>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">X</mml:mi>
<mml:mi>M</mml:mi>
<mml:mi>i</mml:mi>
</mml:msubsup>
<mml:mo>&#x2212;</mml:mo>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">X</mml:mi>
<mml:mi>M</mml:mi>
<mml:mi>j</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>g</mml:mi>
<mml:mi>r</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>y</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mi>M</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:msubsup>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi>m</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mi>M</mml:mi>
</mml:msubsup>
<mml:msub>
<mml:mrow>
<mml:mi>c</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>f</mml:mi>
<mml:mi>f</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>m</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:math>
<label>(5)</label>
</disp-formula>
<italic>coeff</italic>
<sub>
<italic>ij</italic>
</sub>(<italic>t</italic>) is the correlation coefficient of reference sequence and comparison sequence at time <italic>t</italic>, the minimum and maximum values of the absolute difference between the reference sequence and other comparison sequences are denoted as <inline-formula id="inf5">
<mml:math id="m10">
<mml:mrow>
<mml:mi>min</mml:mi>
<mml:mrow>
<mml:mfenced open="|" close="|" separators="|">
<mml:mrow>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">X</mml:mi>
<mml:mi>M</mml:mi>
<mml:mi>i</mml:mi>
</mml:msubsup>
<mml:mo>&#x2212;</mml:mo>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">X</mml:mi>
<mml:mi>M</mml:mi>
<mml:mi>j</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>; <inline-formula id="inf6">
<mml:math id="m11">
<mml:mrow>
<mml:mi>max</mml:mi>
<mml:mrow>
<mml:mfenced open="|" close="|" separators="|">
<mml:mrow>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">X</mml:mi>
<mml:mi>M</mml:mi>
<mml:mi>i</mml:mi>
</mml:msubsup>
<mml:mo>&#x2212;</mml:mo>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">X</mml:mi>
<mml:mi>M</mml:mi>
<mml:mi>j</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>, <italic>&#x3b1;</italic>&#x2208;(0, 1) is resolution factor, which determines the distribution of <italic>grey</italic>
<sub>
<italic>ij</italic>
</sub>. To avoid the distribution of GRA being too compact or too sparse, it is appropriate to set <italic>&#x3b1;</italic> at 0.5.</p>
<p>If the grey relational degree <italic>grey</italic>
<sub>
<italic>ij</italic>
</sub>&#x2208;(0, 1) obtained is close to 1, it indicates that gas <italic>j</italic> and gas <italic>i</italic> have a strong correlation. Otherwise, the correlation between gas <italic>j</italic> and gas <italic>i</italic> is weak. GRA preliminarily determined the correlation information among gases in oil.</p>
</sec>
<sec id="s3-2">
<title>3.2 Graph convolutional network</title>
<p>Graph convolutional network is designed to realize information fusion and feature extraction of graph data with complex association by means of adjacency matrix and graph kernel. Therefore, in this paper, the gas in oil is taken as the graph node, the result of GRA constitute edges, and the weighted gas adjacency matrix <bold>
<italic>A</italic>
</bold>
<sub>0</sub> is constructed based on the grey relational degree to extract the gas correlation features.</p>
<p>Traditional adjacency matrix can only describe the connection between nodes, that is, if there is an edge between node <italic>i</italic> and node <italic>j</italic>, the element <italic>a</italic>
<sub>
<italic>ij</italic>
</sub> in adjacency matrix is 1, otherwise, it is 0, which cannot express the strength of connections. Based on the results of GRA in the last section, the element <italic>a</italic>
<sub>
<italic>ij</italic>
</sub> of gas adjacency matrix is redefined as follows,<disp-formula id="e6">
<mml:math id="m12">
<mml:mrow>
<mml:msub>
<mml:mi>a</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mfenced open="{" close="" separators="|">
<mml:mrow>
<mml:mtable columnalign="center">
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>g</mml:mi>
<mml:mi>r</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>y</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:mi>i</mml:mi>
<mml:mo>&#x2260;</mml:mo>
<mml:mi>j</mml:mi>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mn>0</mml:mn>
<mml:mo>,</mml:mo>
<mml:mi>i</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mi>j</mml:mi>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(6)</label>
</disp-formula>
</p>
<p>The gas original adjacency matrix <bold>
<italic>A</italic>
</bold>
<sub>0</sub> (dimension <italic>n</italic> &#xd7; <italic>n</italic>) is obtained by above equation.</p>
<p>To make the gas node introduce its own features in the graph convolution operation, self-loop is added to the gas node in the graph, and the graph convolution network is as follows,<disp-formula id="e7">
<mml:math id="m13">
<mml:mrow>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">G</mml:mi>
<mml:mi mathvariant="normal">C</mml:mi>
<mml:mi mathvariant="normal">N</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi mathvariant="bold-italic">H</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">s</mml:mi>
<mml:mi mathvariant="normal">i</mml:mi>
<mml:mi mathvariant="normal">g</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msup>
<mml:mi mathvariant="bold-italic">D</mml:mi>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mn mathvariant="bold">1</mml:mn>
<mml:mo>/</mml:mo>
<mml:mn mathvariant="bold">2</mml:mn>
</mml:mrow>
</mml:msup>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi mathvariant="bold-italic">I</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi mathvariant="bold-italic">A</mml:mi>
<mml:mn>0</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:msup>
<mml:mi mathvariant="bold-italic">D</mml:mi>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mn mathvariant="bold">1</mml:mn>
<mml:mo>/</mml:mo>
<mml:mn mathvariant="bold">2</mml:mn>
</mml:mrow>
</mml:msup>
<mml:mi mathvariant="bold-italic">H</mml:mi>
<mml:msub>
<mml:mi mathvariant="bold-italic">W</mml:mi>
<mml:mi mathvariant="normal">G</mml:mi>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi mathvariant="bold-italic">B</mml:mi>
<mml:mi>G</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(7)</label>
</disp-formula>
</p>
<p>Where <bold>
<italic>H</italic>
</bold> (dimension <italic>n</italic> &#xd7; <italic>C</italic>
<sub>hidden</sub>) is the input feature of gas node 1&#x223c;<italic>n</italic>, denoted as <bold>
<italic>H</italic>
</bold> &#x3d; [<bold>
<italic>h</italic>
</bold>
<sup>1</sup>, <bold>
<italic>h</italic>
</bold>
<sup>2</sup>, &#x2026;, <bold>
<italic>h</italic>
</bold>
<sup>
<italic>n</italic>
</sup>], <bold>
<italic>W</italic>
</bold>
<sub>G</sub> (dimension <italic>C</italic>
<sub>hidden</sub> &#xd7; <italic>C</italic>
<sub>out</sub>) is the graph convolution kernel, <bold>
<italic>B</italic>
</bold>
<sub>G</sub> (dimension <italic>n</italic> &#xd7; <italic>C</italic>
<sub>out</sub>) is the graph convolution bias, and <bold>
<italic>D</italic>
</bold> (dimension <italic>n</italic> &#xd7; <italic>n</italic>) is degree matrix, which is the diagonal matrix composed of the number of edges of gas node 1&#x223c; <italic>n</italic> in the graph with self-loop. The self-loop refers to the edge where node links to itself, which avoid the omission of node itself caused by only extracting neighbor node information in feature update. In gas prediction, the self-loop enables the model to consider the original features of the gas when extracting related gas features.</p>
</sec>
<sec id="s3-3">
<title>3.3 Reconstruction of adjacency matrix and adaptive graph convolution networks</title>
<p>To capture real-time dynamic gas state information, the original graph convolution algorithm is improved. In the improved model, the original static adjacency matrix <bold>
<italic>A</italic>
</bold>
<sub>0</sub> is modified to a variable that changes adaptively according to the input feature <bold>
<italic>H</italic>
</bold>. A learnable distance measure is introduced into the adjacency matrix. The distance <italic>dist</italic>
<sub>
<italic>ij</italic>
</sub> between input feature <bold>
<italic>h</italic>
</bold>
<sup>
<italic>i</italic>
</sup> (dimension 1&#xd7;<italic>C</italic>
<sub>hidden</sub>) of gas node <italic>i</italic> and input feature <bold>
<italic>h</italic>
</bold>
<sup>
<italic>j</italic>
</sup> (dimension 1 &#xd7; <italic>C</italic>
<sub>hidden</sub>) of gas node <italic>j</italic> is proposed as follows, <disp-formula id="e8">
<mml:math id="m14">
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>s</mml:mi>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msqrt>
<mml:mrow>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msup>
<mml:mi mathvariant="bold-italic">h</mml:mi>
<mml:mi>i</mml:mi>
</mml:msup>
<mml:msub>
<mml:mi mathvariant="bold-italic">W</mml:mi>
<mml:mi mathvariant="normal">d</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:msup>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msup>
<mml:mi mathvariant="bold-italic">h</mml:mi>
<mml:mi>j</mml:mi>
</mml:msup>
<mml:msub>
<mml:mi mathvariant="bold-italic">W</mml:mi>
<mml:mi mathvariant="normal">d</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mi mathvariant="normal">T</mml:mi>
</mml:msup>
</mml:mrow>
</mml:msqrt>
</mml:mrow>
</mml:math>
<label>(8)</label>
</disp-formula>
</p>
<p>Where <bold>
<italic>W</italic>
</bold>
<sub>d</sub> (dimension <italic>C</italic>
<sub>hidden</sub> &#xd7; <italic>C</italic>
<sub>hidden</sub>) is a learnable weight matrix.</p>
<p>Gaussian function could improve data differentiation, map distance <italic>dist</italic>
<sub>
<italic>ij</italic>
</sub> non-linear to [0,1], and enhance the expression ability of gas correlation. Here, Gaussian function is used to process distance measurement, and reconstructed adjacency matrix elements are obtained. Specific formula is as follows, <disp-formula id="e9">
<mml:math id="m15">
<mml:mrow>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">G</mml:mi>
<mml:mi mathvariant="normal">a</mml:mi>
<mml:mi mathvariant="normal">u</mml:mi>
<mml:mi mathvariant="normal">s</mml:mi>
<mml:mi mathvariant="normal">s</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>s</mml:mi>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mi>exp</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mfrac>
<mml:msup>
<mml:msub>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>s</mml:mi>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mrow>
<mml:mn mathvariant="bold">2</mml:mn>
<mml:msup>
<mml:mi>&#x3c3;</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(9)</label>
</disp-formula>
</p>
<p>In the formula, exp(.) is an exponential function based on the constant e, <italic>&#x3c3;</italic> is a bandwidth, which adjusts the influence range of the Gaussian.</p>
<p>The reconstructed adjacency matrix <bold>
<italic>A</italic>
</bold>
<sub>rec</sub> (dimension <italic>n</italic> &#xd7; <italic>n</italic>) is obtained by above steps, which is expressed as:<disp-formula id="e10">
<mml:math id="m16">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">A</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">r</mml:mi>
<mml:mi mathvariant="normal">e</mml:mi>
<mml:mi mathvariant="normal">c</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">G</mml:mi>
<mml:mi mathvariant="normal">a</mml:mi>
<mml:mi mathvariant="normal">u</mml:mi>
<mml:mi mathvariant="normal">s</mml:mi>
<mml:mi mathvariant="normal">s</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msqrt>
<mml:mrow>
<mml:mi mathvariant="bold-italic">H</mml:mi>
<mml:msub>
<mml:mi mathvariant="bold-italic">W</mml:mi>
<mml:mi mathvariant="normal">d</mml:mi>
</mml:msub>
<mml:msup>
<mml:mrow>
<mml:msup>
<mml:msub>
<mml:mi mathvariant="bold-italic">W</mml:mi>
<mml:mi mathvariant="normal">d</mml:mi>
</mml:msub>
<mml:mi mathvariant="normal">T</mml:mi>
</mml:msup>
<mml:mi mathvariant="bold-italic">H</mml:mi>
</mml:mrow>
<mml:mi mathvariant="normal">T</mml:mi>
</mml:msup>
</mml:mrow>
</mml:msqrt>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(10)</label>
</disp-formula>
</p>
<p>The reconstructed adjacency matrix <bold>
<italic>A</italic>
</bold>
<sub>rec</sub> is introduced into the graph convolution formula, which is used as a modification to the original adjacency matrix <bold>
<italic>A</italic>
</bold>
<sub>0</sub>, and adaptive graph convolution network formula is shown as follows:<disp-formula id="e11">
<mml:math id="m17">
<mml:mrow>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">A</mml:mi>
<mml:mi mathvariant="normal">G</mml:mi>
<mml:mi mathvariant="normal">C</mml:mi>
<mml:mi mathvariant="normal">N</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi mathvariant="bold-italic">H</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">s</mml:mi>
<mml:mi mathvariant="normal">i</mml:mi>
<mml:mi mathvariant="normal">g</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msup>
<mml:mi mathvariant="bold-italic">D</mml:mi>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mrow>
<mml:mn mathvariant="bold">1</mml:mn>
<mml:mo>/</mml:mo>
<mml:mn mathvariant="bold">2</mml:mn>
</mml:mrow>
</mml:mrow>
</mml:msup>
<mml:msub>
<mml:mi mathvariant="bold-italic">A</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">a</mml:mi>
<mml:mi mathvariant="normal">d</mml:mi>
<mml:mi mathvariant="normal">a</mml:mi>
<mml:mi mathvariant="normal">p</mml:mi>
<mml:mi mathvariant="normal">t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msup>
<mml:mi mathvariant="bold-italic">D</mml:mi>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mrow>
<mml:mn mathvariant="bold">1</mml:mn>
<mml:mo>/</mml:mo>
<mml:mn mathvariant="bold">2</mml:mn>
</mml:mrow>
</mml:mrow>
</mml:msup>
<mml:mi mathvariant="bold-italic">H</mml:mi>
<mml:msub>
<mml:mi mathvariant="bold-italic">W</mml:mi>
<mml:mi mathvariant="normal">G</mml:mi>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi mathvariant="bold-italic">B</mml:mi>
<mml:mi>G</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(11)</label>
</disp-formula>
</p>
<p>Where <bold>
<italic>A</italic>
</bold>
<sub>adapt</sub> &#x3d; <bold>
<italic>I</italic>
</bold> &#x2b; <bold>
<italic>A</italic>
</bold>
<sub>0</sub> &#x2b; <bold>
<italic>A</italic>
</bold>
<sub>rec</sub>. It means that the model can learn the current optimal gas relationship from the gas relationship information obtained by GRA. For the learning range is small, the model converges quickly through the training.</p>
</sec>
</sec>
<sec id="s4">
<title>4 Realization of GRU-AGCN model for gas content prediction</title>
<sec id="s4-1">
<title>4.1 GRU-AGCN model design</title>
<p>Traditional methods focus on the single gas sequence but concerns little on the correlation information among gases, and their ability to perceive transformer operation status is limited. In addition, some prediction models with multiple gases also lack effective ways to analyze and extract correlation information of input features. The above problems lead to low stability and reliability of gas content prediction in transformer oil.</p>
<p>Therefore, a GRU-AGCN model for gas content prediction is designed in this paper. The model structure is shown in <xref ref-type="fig" rid="F2">Figure 2</xref>. The model takes the content series of <italic>T</italic> time points of gas 1&#x223c;<italic>n</italic> as input, the temporal features of each gas are extracted through several parallel GRU networks. And temporal features of gases from GRU networks are input into AGCN to extract the gas dynamic correlation features. AGCN firstly calculates the distance measure between gases according to the time series features, and updates the adaptive adjacency matrix based on the original adjacency matrix obtained by GRA. Then, the updated graph adjacency matrix and temporal features are passed into the graph convolution for feature fusion and feature transformation. Finally, the predicted values of gas contents at time (<italic>T</italic> &#x2b; 1) are output by linear layer and activation function.</p>
<fig id="F2" position="float">
<label>FIGURE 2</label>
<caption>
<p>Structure of GRU-AGCN model for gas content prediction.</p>
</caption>
<graphic xlink:href="fenrg-11-1135330-g002.tif"/>
</fig>
</sec>
<sec id="s4-2">
<title>4.2 Prediction process based on GRU-AGCN model</title>
<p>Content prediction process of gases in transformer oil based on GRU-AGCN model are as follows,<list list-type="simple">
<list-item>
<p>1) Obtain historical data of gas content in transformer oil.</p>
</list-item>
<list-item>
<p>2) Remove the outliers in the gas content data and fill in the gaps by interpolation.</p>
</list-item>
<list-item>
<p>3) Normalize the content data of each gas.</p>
</list-item>
<list-item>
<p>4) The normalized data are rolled and intercepted into multiple groups of samples based on the time length <italic>T</italic>, and sample set is separated into training set and test set.</p>
</list-item>
<list-item>
<p>5) Based on the training set, the relational degrees of gases are obtained by GRA, and original gas topology and adjacency matrix <bold>
<italic>A</italic>
</bold>
<sub>0</sub> are constructed.</p>
</list-item>
<list-item>
<p>6) The smoothed mean absolute error is set to be loss function. And Adam optimization algorithm is used to make the model learn training samples. After several iterations of training and parameters adjusting, the prediction requirement could be satisfied by the model.</p>
</list-item>
<list-item>
<p>7) The test set samples are input to model obtained in Step (6) for prediction test, and the predicted results are analyzed.</p>
</list-item>
</list>
</p>
</sec>
</sec>
<sec id="s5">
<title>5 Analysis of experimental examples</title>
<p>To verify the superiority of the proposed model, a GRU-AGCN model is built by Python language based on Spyder platform and PyTorch library.</p>
<p>Mean absolute error (MAE) and mean absolute percentage error (MAPE) are used to evaluate the prediction performance of the model. Suppose that gas content series in time 1&#x223c;<italic>t</italic> is used for prediction, denoted as <bold>
<italic>X</italic>
</bold> &#x3d; [<italic>x</italic>
<sub>1</sub>, <italic>x</italic>
<sub>2</sub>, &#x2026;, <italic>x</italic>
<sub>
<italic>t</italic>
</sub>], the corresponding predicted value series is denoted as <bold>
<italic>Y</italic>
</bold> &#x3d; [<italic>y</italic>
<sub>1</sub>, <italic>y</italic>
<sub>2</sub>, &#x2026;, <italic>y</italic>
<sub>
<italic>t</italic>
</sub>]. The MAE and MAPE are calculated as follows:<disp-formula id="e12">
<mml:math id="m18">
<mml:mrow>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">M</mml:mi>
<mml:mi mathvariant="normal">A</mml:mi>
<mml:mi mathvariant="normal">E</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi mathvariant="bold-italic">X</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi mathvariant="bold-italic">Y</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mrow>
<mml:mstyle displaystyle="true">
<mml:munderover>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mi>t</mml:mi>
</mml:munderover>
</mml:mstyle>
<mml:mrow>
<mml:mfenced open="|" close="|" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>x</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>y</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(12)</label>
</disp-formula>
<disp-formula id="e13">
<mml:math id="m19">
<mml:mrow>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">M</mml:mi>
<mml:mi mathvariant="normal">A</mml:mi>
<mml:mi mathvariant="normal">P</mml:mi>
<mml:mi mathvariant="normal">E</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi mathvariant="bold-italic">X</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi mathvariant="bold-italic">Y</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mrow>
<mml:mstyle displaystyle="true">
<mml:munderover>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mi>t</mml:mi>
</mml:munderover>
</mml:mstyle>
<mml:mfrac>
<mml:mrow>
<mml:mfenced open="|" close="|" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>x</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>y</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:msub>
<mml:mi>y</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mfrac>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(13)</label>
</disp-formula>
</p>
<p>A converter transformer in Henan Province is taken as the research object. It was put into operation in 2014 and has been operating normally ever since. Dissolved gases monitored include hydrogen (H<sub>2</sub>), methane (CH<sub>4</sub>), ethane (C<sub>2</sub>H<sub>6</sub>), ethylene (C<sub>2</sub>H<sub>4</sub>), acetylene (C<sub>2</sub>H<sub>2</sub>), carbon monoxide (CO) and carbon dioxide (CO<sub>2</sub>), and their content is within the normal range stipulated by IEC. Eight hundred and twenty sets of data of the transformer from May to October 2019 were selected for the experiment. The time period is 4&#xa0;h. The content of C<sub>2</sub>H<sub>2</sub> is always 0 in the monitoring interval, which is not used as model input. The input of the model is set as the content series of six gases in 20 consecutive time points, and the output is set as the predicted value of gases at the next time point. The original gas content data are processed according to the above steps and 800 groups of samples were obtained. The size of training set and test set are divided by 3:1. The first 600 groups are set as the training set, and the last 200 groups are the test set.</p>
<p>The grey relational degree table of six gas content is obtained by GRA. As shown in <xref ref-type="table" rid="T1">Table 1</xref>, the correlations rang is distributed in [0.634, 1.000], which indicates that there are clearly exists correlations among six gases, but their intensity are different. Regardless of the self-correlation, CO and CO<sub>2</sub> are the most correlated, while the correlation between H<sub>2</sub> and CO<sub>2</sub> is the weakest.</p>
<table-wrap id="T1" position="float">
<label>TABLE 1</label>
<caption>
<p>Grey relational degree of gases in transformer oil (<italic>&#x3b1;</italic> &#x3d; 0.5).</p>
</caption>
<table>
<thead valign="top">
<tr>
<th rowspan="2" align="center">References</th>
<th colspan="6" align="center">Comparison</th>
</tr>
<tr>
<th align="center">H<sub>2</sub>
</th>
<th align="center">CH<sub>4</sub>
</th>
<th align="center">C<sub>2</sub>H<sub>6</sub>
</th>
<th align="center">C<sub>2</sub>H<sub>4</sub>
</th>
<th align="center">CO</th>
<th align="center">CO<sub>2</sub>
</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="center">H<sub>2</sub>
</td>
<td align="center">1.000</td>
<td align="center">0.675</td>
<td align="center">0.713</td>
<td align="center">0.694</td>
<td align="center">0.646</td>
<td align="center">0.634</td>
</tr>
<tr>
<td align="center">CH<sub>4</sub>
</td>
<td align="center">0.675</td>
<td align="center">1.000</td>
<td align="center">0.763</td>
<td align="center">0.727</td>
<td align="center">0.815</td>
<td align="center">0.807</td>
</tr>
<tr>
<td align="center">C<sub>2</sub>H<sub>6</sub>
</td>
<td align="center">0.713</td>
<td align="center">0.763</td>
<td align="center">1.000</td>
<td align="center">0.767</td>
<td align="center">0.749</td>
<td align="center">0.733</td>
</tr>
<tr>
<td align="center">C<sub>2</sub>H<sub>4</sub>
</td>
<td align="center">0.694</td>
<td align="center">0.727</td>
<td align="center">0.767</td>
<td align="center">1.000</td>
<td align="center">0.665</td>
<td align="center">0.661</td>
</tr>
<tr>
<td align="center">CO</td>
<td align="center">0.646</td>
<td align="center">0.815</td>
<td align="center">0.749</td>
<td align="center">0.665</td>
<td align="center">1.000</td>
<td align="center">0.854</td>
</tr>
<tr>
<td align="center">CO<sub>2</sub>
</td>
<td align="center">0.634</td>
<td align="center">0.807</td>
<td align="center">0.733</td>
<td align="center">0.661</td>
<td align="center">0.854</td>
<td align="center">1.000</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>Based on <xref ref-type="table" rid="T1">Table 1</xref>, the correlation degrees are introduced into AGCN of the model as original adjacency matrix. And parameter <italic>&#x3c3;</italic>
<sup>2</sup> in the AGCN is set to 0.5, because it achieves the best result in model training.</p>
<p>Then, GRU-AGCN model is trained for 200 times, and the initial learning rate is set to 0.0005. The sample loss curve in <xref ref-type="fig" rid="F3">Figure 3</xref> is obtained. As shown in <xref ref-type="fig" rid="F3">Figure 3</xref>, the model parameters converge to the global optimum quickly and smoothly during training.</p>
<fig id="F3" position="float">
<label>FIGURE 3</label>
<caption>
<p>Loss curve in training.</p>
</caption>
<graphic xlink:href="fenrg-11-1135330-g003.tif"/>
</fig>
<p>To verify the accuracy and reliability of the proposed model in gas content prediction task, the GRU-AGCN model and other traditional models are used for single-step prediction experiment. The models used for comparison include RNN and its variants GRU Net and LSTM Net. In addition, original GCN is introduced into GRU Net and LSTM Net to form new models, denoted as GRU-GCN and LSTM-GCN respectively. Specific parameter settings of models are shown in <xref ref-type="table" rid="T2">Table 2</xref>.</p>
<table-wrap id="T2" position="float">
<label>TABLE 2</label>
<caption>
<p>Parameters of prediction models.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="center">Model type</th>
<th align="center">Input dimension</th>
<th align="center">Output dimension</th>
<th align="center">Number of layers</th>
<th align="center">Hidden dimension</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="center">GRU-AGCN</td>
<td align="center">6 &#xd7; 20</td>
<td align="center">6 &#xd7; 1</td>
<td align="center">GRU:1 AGCN:1</td>
<td align="center">[20, 10]</td>
</tr>
<tr>
<td align="center">GRU-GCN</td>
<td align="center">6 &#xd7; 20</td>
<td align="center">6 &#xd7; 1</td>
<td align="center">GRU:1 GCN:1</td>
<td align="center">[20, 10]</td>
</tr>
<tr>
<td align="center">LSTM-GCN</td>
<td align="center">6 &#xd7; 20</td>
<td align="center">6 &#xd7; 1</td>
<td align="center">GRU:1 GCN:1</td>
<td align="center">[20, 10]</td>
</tr>
<tr>
<td align="center">GRU Net</td>
<td align="center">20</td>
<td align="center">1</td>
<td align="center">1</td>
<td align="center">20</td>
</tr>
<tr>
<td align="center">LSTM Net</td>
<td align="center">20</td>
<td align="center">1</td>
<td align="center">1</td>
<td align="center">20</td>
</tr>
<tr>
<td align="center">RNN</td>
<td align="center">20</td>
<td align="center">1</td>
<td align="center">2</td>
<td align="center">[20, 10]</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>Except for the model proposed in this paper and models introduced GCN, all other models are designed for univariate prediction. The input time scale, prediction range and training conditions of all models remained consistent. Hidden features of models are output through a linear layer and ReLU activation function. All models were trained 200 times.</p>
<p>Limited by space, this section mainly analyzes the result of CO<sub>2</sub> content prediction experiment. <xref ref-type="table" rid="T3">Table 3</xref> records the training time and prediction effect of CO<sub>2</sub> content prediction by above models.</p>
<table-wrap id="T3" position="float">
<label>TABLE 3</label>
<caption>
<p>Effects of models on CO<sub>2</sub> content prediction.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="center">Model type</th>
<th align="center">MAE (%)</th>
<th align="center">MAPE (%)</th>
<th align="center">Training time (s)</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="center">GRU-AGCN</td>
<td align="center">1.95</td>
<td align="center">3.18</td>
<td align="center">289.9</td>
</tr>
<tr>
<td align="center">GRU-GCN</td>
<td align="center">3.17</td>
<td align="center">4.50</td>
<td align="center">275.8</td>
</tr>
<tr>
<td align="center">LSTM-GCN</td>
<td align="center">3.32</td>
<td align="center">4.59</td>
<td align="center">1188.8</td>
</tr>
<tr>
<td align="center">GRU Net</td>
<td align="center">3.31</td>
<td align="center">4.67</td>
<td align="center">269.3</td>
</tr>
<tr>
<td align="center">LSTM Net</td>
<td align="center">3.56</td>
<td align="center">4.99</td>
<td align="center">1149.6</td>
</tr>
<tr>
<td align="center">RNN</td>
<td align="center">4.43</td>
<td align="center">6.04</td>
<td align="center">221.3</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>Firstly, RNN and its variants GRU Net and LSTM Net are compared. Obviously, RNN has the fastest training, but its prediction accuracy is far lower than other networks. Both GRU Net and LSTM Net introduce gating mechanism, which selectively extract features at different moments. Compared with RNN, MAE percentage of the two models decreased by 1.12 and 0.87 respectively. However, the training time of LSTM Net is much longer than that of GRU Net. It is proved that GRU Net can extract temporal features of gas content efficiently and is easier to train than LSTM Net.</p>
<p>Then, GRU Net and LSTM Net are compared with GRU-GCN and LSTM-GCN respectively. The prediction accuracy of model introduced GCN is significantly higher than the original model. The reason lies in that the information of the related gases has been used for the prediction by GRU in GRU-GCN model and LSTM-GCN model.</p>
<p>Finally, GRU-AGCN model is compared with GRU Net and GRU-GCN model. The prediction result by above models are shown in <xref ref-type="fig" rid="F4">Figure 4</xref>. Based on the analysis of <xref ref-type="fig" rid="F4">Figure 4</xref> and <xref ref-type="table" rid="T3">Table 3</xref>, it can be found that GRU Net model only accurately predict the stable development of CO<sub>2</sub>, as shown in 1st&#x223c;125th point of <xref ref-type="fig" rid="F4">Figure 4</xref>, which shows a poor prediction effect on series with large fluctuation. The GRU- GCN utilized the correlations between the gases. But the model cannot obtain satisfactory prediction effect in long time span, because the real-time changes in gas relationships cannot be represented by static adjacency matrices. The proposed GRU-AGCN model has the smallest MAE and MAPE in prediction, which is only 20.6&#xa0;s slower than GRU Net during the training. The proposed method obtains the best effect with an acceptable time cost. The reason is that GRU part of the proposed model can efficiently extract the temporal features of gas series, while the AGCN part adaptively adjust the graph adjacency matrix and dynamically acquire the gas relationship from learning samples, which makes the prediction result more accurate and reasonable than other models. The proposed model keeps high accuracy and reliability for series in stable situation or large fluctuations.</p>
<fig id="F4" position="float">
<label>FIGURE 4</label>
<caption>
<p>The effect of three models on CO<sub>2</sub> content prediction.</p>
</caption>
<graphic xlink:href="fenrg-11-1135330-g004.tif"/>
</fig>
</sec>
<sec sec-type="conclusion" id="s6">
<title>6 Conclusion</title>
<p>The change of gas content in transformer oil reflects the health status of transformer, which needs to consider both the temporal information from gas itself and the dynamic correlation information from related gases. Base on above two aspects, a prediction method based on GRU-AGCN model is proposed. From the experiment comparison and analysis, it can be concluded that:<list list-type="simple">
<list-item>
<p>1) GRU Net can extract key information of important moments from gas series, and its prediction is more accurate than RNN. Compared with LSTM Net, GRU Net is simpler, training faster, and extraction efficiency for the temporal features is higher.</p>
</list-item>
<list-item>
<p>2) Compared with the single gas variable prediction model, the multi-gas variable combined prediction model with GCN consider the correlation between gases, which significantly improves the prediction accuracy.</p>
</list-item>
<list-item>
<p>3) The AGCN proposed in this paper extracts the gas relations dynamically, which makes the prediction result more accurate and reasonable. The proposed model can maintain high accuracy for different time period.</p>
</list-item>
</list>
</p>
</sec>
</body>
<back>
<sec sec-type="data-availability" id="s7">
<title>Data availability statement</title>
<p>The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.</p>
</sec>
<sec id="s8">
<title>Author contributions</title>
<p>DL and WC contributed to the development and design of the model. JF implemented the programming of the proposed model algorithm. JL collected experimental data. JY and KZ designed the experiment to prove the superiority of the proposed method. All authors wrote parts of the manuscript. All authors have read and agreed to the submission of the manuscript.</p>
</sec>
<sec id="s9">
<title>Funding</title>
<p>The research is supported by the Science and Technology Project of State Grid Shanxi Electric Power Co., LTD., China (No. SGSXJC00YCJS2100679).</p>
</sec>
<sec sec-type="COI-statement" id="s10">
<title>Conflict of interest</title>
<p>Authors WC, JL, JY, and KZ were employed by State Grid Jincheng Power Supply Company.</p>
<p>The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.</p>
</sec>
<sec sec-type="disclaimer" id="s11">
<title>Publisher&#x2019;s note</title>
<p>All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.</p>
</sec>
<ref-list>
<title>References</title>
<ref id="B1">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Atherfold</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>van Zyl</surname>
<given-names>T. L.</given-names>
</name>
</person-group> (<year>2020</year>). &#x201c;<article-title>A method for dissolved gas forecasting in power transformers using LS-SVM</article-title>,&#x201d; in <source>2020 IEEE 23rd international conference on information fusion (FUSION)</source> (<publisher-loc>Rustenburg, South Africa</publisher-loc>, <fpage>1</fpage>&#x2013;<lpage>8</lpage>. <comment>6-9 July 2020</comment>. <pub-id pub-id-type="doi">10.23919/FUSION45008.2020.9190216</pub-id>
</citation>
</ref>
<ref id="B2">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Chen</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Hu</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Dong</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Zheng</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>Wei</surname>
<given-names>J.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Transformer operating state monitoring system based on wireless sensor networks</article-title>. <source>IEEE Sensors J.</source> <volume>21</volume>, <fpage>25098</fpage>&#x2013;<lpage>25105</lpage>. <pub-id pub-id-type="doi">10.1109/JSEN.2021.3050763</pub-id>
</citation>
</ref>
<ref id="B3">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Fessler</surname>
<given-names>W. A.</given-names>
</name>
<name>
<surname>Rouse</surname>
<given-names>T. O.</given-names>
</name>
<name>
<surname>McNutt</surname>
<given-names>W. J.</given-names>
</name>
<name>
<surname>Compton</surname>
<given-names>O. R.</given-names>
</name>
</person-group> (<year>1989</year>). <article-title>A refined mathematical model for prediction of bubble evolution in transformers</article-title>. <source>IEEE Trans. Power Deliv.</source> <volume>4</volume>, <fpage>391</fpage>&#x2013;<lpage>404</lpage>. <pub-id pub-id-type="doi">10.1109/61.19229</pub-id>
</citation>
</ref>
<ref id="B4">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ghunem</surname>
<given-names>R. A.</given-names>
</name>
<name>
<surname>Assaleh</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>El-hag</surname>
<given-names>A. H.</given-names>
</name>
</person-group> (<year>2012</year>). <article-title>Artificial neural networks with stepwise regression for predicting transformer oil furan content</article-title>. <source>IEEE Trans. Dielectr. Electr. Insulation</source> <volume>19</volume>, <fpage>414</fpage>&#x2013;<lpage>420</lpage>. <pub-id pub-id-type="doi">10.1109/TDEI.2012.6180233</pub-id>
</citation>
</ref>
<ref id="B5">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>G&#xf3;mez</surname>
<given-names>N. A.</given-names>
</name>
<name>
<surname>Wilhelm</surname>
<given-names>H. M.</given-names>
</name>
<name>
<surname>Santos</surname>
<given-names>C. C.</given-names>
</name>
<name>
<surname>Stocco</surname>
<given-names>G. B.</given-names>
</name>
</person-group> (<year>2014</year>). <article-title>Dissolved gas analysis (DGA) of natural ester insulating fluids with different chemical compositions</article-title>. <source>IEEE Trans. Dielectr. Electr. Insulation</source> <volume>21</volume>, <fpage>1071</fpage>&#x2013;<lpage>1078</lpage>. <pub-id pub-id-type="doi">10.1109/TDEI.2014.6832250</pub-id>
</citation>
</ref>
<ref id="B6">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Huang</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Feng</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>X.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Multi-field collaborative of oil-immersed transformer for distributed energy resources temperature rise considering the influence of heat transfer oil</article-title>. <source>Front. Energy Res.</source> <volume>9</volume>. <pub-id pub-id-type="doi">10.3389/fenrg.2021.701873</pub-id>
</citation>
</ref>
<ref id="B7">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lin</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Sheng</surname>
<given-names>G. H.</given-names>
</name>
<name>
<surname>Yan</surname>
<given-names>Y. J.</given-names>
</name>
<name>
<surname>Dai</surname>
<given-names>J. J.</given-names>
</name>
<name>
<surname>Jiang</surname>
<given-names>X. C.</given-names>
</name>
</person-group> (<year>2018a</year>). <article-title>Prediction of dissolved gas concentrations in transformer oil based on the KPCA-FFOA-GRNN model</article-title>. <source>ENERGIES</source> <volume>11</volume>, <fpage>225</fpage>. <pub-id pub-id-type="doi">10.3390/en11010225</pub-id>
</citation>
</ref>
<ref id="B8">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lin</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Su</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Yan</surname>
<given-names>Y. J.</given-names>
</name>
<name>
<surname>Sheng</surname>
<given-names>G. H.</given-names>
</name>
<name>
<surname>Xie</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Jiang</surname>
<given-names>X. C.</given-names>
</name>
</person-group> (<year>2018b</year>). <article-title>Prediction method for power transformer running state based on LSTM_DBN network</article-title>. <source>ENERGIES</source> <volume>11</volume>, <fpage>1880</fpage>. <pub-id pub-id-type="doi">10.3390/en11071880</pub-id>
</citation>
</ref>
<ref id="B9">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ma</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Hu</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Shang</surname>
<given-names>Y.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>A new method for transformer fault prediction based on multifeature enhancement and refined long short-term memory</article-title>. <source>IEEE Trans. Instrum. Meas.</source> <volume>70</volume>, <fpage>1</fpage>&#x2013;<lpage>11</lpage>. <pub-id pub-id-type="doi">10.1109/TIM.2021.3098383</pub-id>
</citation>
</ref>
<ref id="B10">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Minh</surname>
<given-names>D. L.</given-names>
</name>
<name>
<surname>Sadeghi-Niaraki</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Huy</surname>
<given-names>H. D.</given-names>
</name>
<name>
<surname>Min</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Moon</surname>
<given-names>H.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>Deep learning approach for short-term stock trends prediction based on two-stream gated recurrent unit network</article-title>. <source>IEEE Access</source> <volume>6</volume>, <fpage>55392</fpage>&#x2013;<lpage>55404</lpage>. <pub-id pub-id-type="doi">10.1109/ACCESS.2018.2868970</pub-id>
</citation>
</ref>
<ref id="B11">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Shaban</surname>
<given-names>K. B.</given-names>
</name>
<name>
<surname>El-Hag</surname>
<given-names>A. H.</given-names>
</name>
<name>
<surname>Benhmed</surname>
<given-names>K.</given-names>
</name>
</person-group> (<year>2016</year>). <article-title>Prediction of transformer furan levels</article-title>. <source>IEEE Trans. Power Deliv.</source> <volume>31</volume>, <fpage>1778</fpage>&#x2013;<lpage>1779</lpage>. <pub-id pub-id-type="doi">10.1109/TPWRD.2016.2521320</pub-id>
</citation>
</ref>
<ref id="B12">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Su</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Huang</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Qin</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Zhao</surname>
<given-names>W.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Transformer vibration detection based on YOLOv4 and optical flow in background of high proportion of renewable energy access</article-title>. <source>Front. Energy Res.</source> <volume>10</volume>. <pub-id pub-id-type="doi">10.3389/fenrg.2022.764903</pub-id>
</citation>
</ref>
<ref id="B13">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wu</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Zhou</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Lin</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Chen</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Huang</surname>
<given-names>Y.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Improved genetic algorithm and XGBoost classifier for power transformer fault diagnosis</article-title>. <source>Front. Energy Res.</source> <volume>9</volume>. <pub-id pub-id-type="doi">10.3389/fenrg.2021.745744</pub-id>
</citation>
</ref>
<ref id="B14">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Xie</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Zhao</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Hong</surname>
<given-names>T.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Transformer monitoring and protection in dynamic power systems &#x2013; a review</article-title>. <source>Front. Energy Res.</source> <volume>8</volume>. <pub-id pub-id-type="doi">10.3389/fenrg.2020.00150</pub-id>
</citation>
</ref>
</ref-list>
</back>
</article>