<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article article-type="research-article" 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">775259</article-id>
<article-id pub-id-type="doi">10.3389/fenrg.2021.775259</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Energy Research</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Techno-Economic Evaluation of Biological and Fluidised-Bed Based Methanation Process Chains for Grid-Ready Biomethane Production</article-title>
<alt-title alt-title-type="left-running-head">Gantenbein et al.</alt-title>
<alt-title alt-title-type="right-running-head">Techno-Economic Evaluation Biological/Fluidised-Bed Methanation</alt-title>
</title-group>
<contrib-group>
<contrib contrib-type="author">
<name>
<surname>Gantenbein</surname>
<given-names>Andreas</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/1474266/overview"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Kr&#xf6;cher</surname>
<given-names>Oliver</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Biollaz</surname>
<given-names>Serge M. A.</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Schildhauer</surname>
<given-names>Tilman J.</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/999936/overview"/>
</contrib>
</contrib-group>
<aff id="aff1">
<sup>1</sup>
<institution>Paul Scherrer Institut (PSI)</institution>, <addr-line>Villigen</addr-line>, <country>Switzerland</country>
</aff>
<aff id="aff2">
<sup>2</sup>
<institution>Institute for Chemical Sciences and Engineering, &#x00C9;cole polytechnique f&#x00E9;d&#x00E9;rale de Lausanne (EPFL)</institution>, <addr-line>Lausanne</addr-line>, <country>Switzerland</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/1264063/overview">Ilkka Hannula</ext-link>, International Energy Agency, France</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/277673/overview">Sunil A. Patil</ext-link>, Indian Institute of Science Education and Research Mohali, India</p>
<p>
<ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/1367088/overview">Kristian Melin</ext-link>, LUT University, Finland</p>
</fn>
<corresp id="c001">&#x2a;Correspondence: Tilman J. Schildhauer, <email>tilman.schildhauer@psi.ch</email>
</corresp>
<fn fn-type="other">
<p>This article was submitted to Bioenergy and Biofuels, a section of the journal Frontiers in Energy Research</p>
</fn>
</author-notes>
<pub-date pub-type="epub">
<day>24</day>
<month>03</month>
<year>2022</year>
</pub-date>
<pub-date pub-type="collection">
<year>2021</year>
</pub-date>
<volume>9</volume>
<elocation-id>775259</elocation-id>
<history>
<date date-type="received">
<day>13</day>
<month>09</month>
<year>2021</year>
</date>
<date date-type="accepted">
<day>23</day>
<month>12</month>
<year>2021</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#xa9; 2022 Gantenbein, Kr&#xf6;cher, Biollaz and Schildhauer.</copyright-statement>
<copyright-year>2022</copyright-year>
<copyright-holder>Gantenbein, Kr&#xf6;cher, Biollaz and Schildhauer</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>Three different power-to-methane process chains with grid injection in two scales (1&#xa0;MW<sub>el</sub> and 6&#xa0;MW<sub>el</sub>) were analysed regarding their investment and operation cost. The process chains were based on biological or catalytic bubbling fluidised bed methanation in combination with proton exchange membrane or solid oxide electrolyser cells. A bottom-up techno-economic analysis showed a cost benefit of around 17&#x2013;19% lower biomethane production cost for the bubbling fluidised bed technology as less than a third of the reactor volumes is required for catalytic methanation. This cost benefit is only given in combination with PEM electrolysis, as the high-temperature electrolyser stacks currently result in high investment cost. Based on electricity cost of 5&#xa0;&#x20AC;-ct/kWh<sub>el</sub> and a plant size of 6&#xa0;MW<sub>el</sub>, biomethane production cost of 13.95&#xa0;&#x20ac;-ct./kWh for catalytic and 17.30&#xa0;&#x20ac;-ct/kWh for biological methanation could be obtained, both including PEM electrolysis. A significant efficiency increase by integrating the heat of catalytic methanation reaction with the high-temperature electrolysis can be achieved; however investment cost have to decrease below 1000&#xa0;&#x20ac;/kW<sub>el</sub> to obtain economically feasible production cost of biomethane. Under current economic and technological circumstances, CO<sub>2</sub> methanation using the bubbling fluidised bed technology is the most cost effective.</p>
</abstract>
<kwd-group>
<kwd>biogas</kwd>
<kwd>power-to-gas</kwd>
<kwd>methanation</kwd>
<kwd>cost</kwd>
<kwd>biological methanation</kwd>
<kwd>catalytic fluidised bed methanation</kwd>
</kwd-group>
<contract-sponsor id="cn001">Bundesamt f&#xfc;r Energie<named-content content-type="fundref-id">10.13039/501100005380</named-content>
</contract-sponsor>
</article-meta>
</front>
<body>
<sec id="s1">
<title>1 Introduction</title>
<p>The integration of biomass in the future energy system bears a large potential for the decarbonisation of society. Switzerland&#x2019;s gas providers agreed on increasing the share of renewable gas in the distribution grid to 30% of the residential heating consumption by 2030 (<xref ref-type="bibr" rid="B38">Schmid and Decurtins, 2018</xref>). The largest potential in unused domestic bio-resources is associated with biogas from digestion processes (<xref ref-type="bibr" rid="B43">Thees et al., 2017</xref>). In order to make use of this potential, biogas upgrading and utilisation processes on a relatively small scale are necessary. Methanation of such biogenic gases also enables a flexible way of integration and seasonal storage of renewable electricity in the future energy system. The Power-to-Gas technology can provide a large, yet unused potential for such chemical energy storage (<xref ref-type="bibr" rid="B20">G&#xf6;tz et al., 2016</xref>; <xref ref-type="bibr" rid="B5">B&#xf6;hm et al., 2020</xref>). The technology allows the coupling between different energy sectors (i.e. electricity, gas, and heat) and has also a large potential in reducing the carbon footprint of the existing gas grid (<xref ref-type="bibr" rid="B25">Kober et al., 2019</xref>). The underlying methanation reaction is currently being tested in a variety of lab-scale and demonstration-scale projects (<xref ref-type="bibr" rid="B44">Thema et al., 2019</xref>). The technologies investigated involve mostly fixed-bed and biological methanation concepts (<xref ref-type="bibr" rid="B44">Thema et al., 2019</xref>).</p>
<p>In previous work, Witte et al. compared the cost structure of catalytic fixed bed and catalytic bubbling fluidised bed (BFB) methanation. It could be shown, that also fluidised bed methanation and subsequent membrane upgrading can provide a cost-efficient way to produce biomethane (<xref ref-type="bibr" rid="B53">Witte et al., 2018b</xref>). The maturity of the technology could be demonstrated by stable operation for more than 1,000&#xa0;h with subsequent grid injection (<xref ref-type="bibr" rid="B51">Witte et al., 2019</xref>), but no demonstration project is in operation yet.</p>
<p>While biogas proved to be a cost-efficient carbon source for direct methanation of biogas, any other CO<sub>2</sub> source can be used (<xref ref-type="bibr" rid="B52">Witte et al., 2018a</xref>; <xref ref-type="bibr" rid="B8">Calbry-Muzyka and Schildhauer, 2020</xref>). Direct air capture and similar technologies to obtain CO<sub>2</sub> are associated with high cost. However, also waste streams of almost pure CO<sub>2</sub> are available and can be used for methanation. These include for example the off-gas streams of conventional biogas upgrading units. The Store&#x26;Go demonstration plants make use of different technologies, such as direct air capture, waste-CO<sub>2</sub> from a waste water treatment plants, as well as a CO<sub>2</sub> stream from a bioethanol plant (<xref ref-type="bibr" rid="B40">Store and Go Project, 2016b</xref>). The carbon source, as well as the source of electricity for the electrolysis highly influence the profitability of the overall process, as a higher selling price for biomethane can be achieved by renewable sources (<xref ref-type="bibr" rid="B13">Energie360&#xb0;, 2021</xref>).</p>
<p>Due to the high share in operation cost, the efficiency of the electrolyser technology has a large influence on the overall production cost of the biomethane (<xref ref-type="bibr" rid="B52">Witte et al., 2018a</xref>). Alkaline electrolysis (AEL) is the most commonly used technology (<xref ref-type="bibr" rid="B32">Posdziech et al., 2019</xref>). Nevertheless, proton exchange membrane (PEM) technology allows for higher current densities and higher operation pressures than conventional alkaline electrolysis. This technology is increasing its industrial maturity and commercial applications are now available on the market. The solid oxide electrolysis cell (SOE) technology is still under development, mainly because stack lifetime issues due to material degradation need to be solved (<xref ref-type="bibr" rid="B32">Posdziech et al., 2019</xref>). SOE technology allows an increase of electrolyser efficiency by integration of an external heat source, for example catalytic methanation (<xref ref-type="bibr" rid="B21">Gruber et al., 2018</xref>).</p>
<p>This work provides a bottom-up techno-economical evaluation of three different Power-to-Gas (PtG) process chains in two different scales. For each of these processes the size and performance of the main equipment is determined, and its capital and operation costs are evaluated. The study compares the cost structure of the bubble column-based biological methanation and catalytic bubbling fluidised bed methanation. Furthermore, the influence of electrolyser technology and cost is investigated.</p>
</sec>
<sec id="s2">
<title>2 Methods</title>
<sec id="s2-1">
<title>2.1 Process Design</title>
<p>The processes evaluated in this work are divided into two scales corresponding to the electrolyser power of two demonstration projects. The first scale corresponds to the 1&#xa0;MW<sub>el</sub> BioCat plant in Aved&#xf8;re (Denmark), based on biological methanation (<xref ref-type="bibr" rid="B12">Electrochaea.dk ApS, 2014</xref>). The second scale refers to the 6&#xa0;MW<sub>el</sub> Audi e-gas project in Werlte (Germany), which operates alkaline electrolysers with cooled, fixed-bed methanation reactors (<xref ref-type="bibr" rid="B20">G&#xf6;tz et al., 2016</xref>). Furthermore, two methanation technologies are compared, namely bubbling fluidised bed methanation and biological methanation in a stirred bubble column.</p>
<p>As electrolyser technologies, proton exchange membrane, as well as high-temperature electrolysis by solid oxide electrolysis cells are considered. Based on the available electrolysis power and subsequent H<sub>2</sub>-feed, a CO<sub>2</sub> feed stream is calculated according to the respective needs of the methanation technology. For the cases including the SOE technology, the hydrogen output of the electrolyser is kept the same as in the PEM-cases. Therefore, plant costs are comparable to the PEM-BFB case.</p>
<p>An overview on the considered cases is given in <xref ref-type="table" rid="T1">Table 1</xref>.</p>
<table-wrap id="T1" position="float">
<label>TABLE 1</label>
<caption>
<p>Power-to-Gas process chains for CO<sub>2</sub> methanation considered in this study. The processes abbreviated with a &#x201c;B&#x201d; refer to the biological methanation, those abbreviated with &#x201c;C&#x201d; to the catalytic case. The numbers indicate the 1&#xa0;MW and 6&#xa0;MW-scale respectively. The variations marked with &#x201c;bg&#x201d; use biogas as feed instead of pure CO<sub>2</sub>. The processes based on SOE instead of PEM electrolysis are indicated with &#x201c;S&#x201d;.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="left">Case</th>
<th align="center">Electrolysis power</th>
<th align="center">Electrolyser technology</th>
<th align="center">Methanation technology</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="left">B1</td>
<td align="center">1&#xa0;MW</td>
<td align="center">PEM</td>
<td align="center">Biological</td>
</tr>
<tr>
<td align="left">B6</td>
<td align="center">6&#xa0;MW</td>
<td align="center">PEM</td>
<td align="center">Biological</td>
</tr>
<tr>
<td align="left">C1/C1bg</td>
<td align="center">1&#xa0;MW</td>
<td align="center">PEM</td>
<td align="center">Catalytic BFB</td>
</tr>
<tr>
<td align="left">C6/C6bg</td>
<td align="center">6&#xa0;MW</td>
<td align="center">PEM</td>
<td align="center">Catalytic BFB</td>
</tr>
<tr>
<td align="left">C1S/C1Sbg</td>
<td align="center">1&#xa0;MW<sub>eq</sub>
</td>
<td align="center">SOE</td>
<td align="center">Catalytic BFB</td>
</tr>
<tr>
<td align="left">C6S/C6Sbg</td>
<td align="center">6&#xa0;MW<sub>eq</sub>
</td>
<td align="center">SOE</td>
<td align="center">Catalytic BFB</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>The combination of biological methanation with solid oxide electrolysis is not considered in this work since the temperature level of excess heat from biological methanation is too low for integration with high temperature electrolysis. Such a process would either require the addition of an external high temperature heat source or result in a decreased electrical efficiency of the SOE.</p>
<p>The two catalytic processes were furthermore evaluated for the case of direct methanation of biogas. In this case, the electrolyser power is fixed to the respective CO<sub>2</sub>-case and a biogas-feed containing the necessary CO<sub>2</sub> is provided to the process. The raw biogas is assumed to have a methane content of 60 vol.-%. The cases with biogas as feed are marked with &#x201c;bg.&#x201d;</p>
<sec id="s2-1-1">
<title>2.1.1 Cases B1 and B6: Biological Methanation With PEM Electrolysis (PEM-Bio)</title>
<p>The process in cases B1 and B6 consist of a biological methanation unit fed by a PEM electrolyser in the 1&#xa0;MW and 6&#xa0;MW-scale, as shown in <xref ref-type="fig" rid="F1">Figure 1</xref>. A purified CO<sub>2</sub> feed stream is compressed to the operation pressure of the reactor and mixed with hydrogen from the electrolyser. The CO<sub>2</sub> feed stream to the plant and the water feed stream to the electrolyser are chosen according to the electrolyser power and the H<sub>2</sub>/CO<sub>2</sub>-ratio required in the feed of the reactor. The water feed stream is compressed to the operation pressure of the PEM-electrolyser, which is 30&#xa0;bar. The valorisation of oxygen produced by the electrolyser is optional, but not part of this work.</p>
<fig id="F1" position="float">
<label>FIGURE 1</label>
<caption>
<p>Process chain for the biological methanation of CO<sub>2</sub>.</p>
</caption>
<graphic xlink:href="fenrg-09-775259-g001.tif"/>
</fig>
<p>The feed gas mixture is then pre-heated to the reactor inlet temperature, which is 60&#xb0;C in the case of biological methanation. As this temperature level is low and temperature differences of the feed gases relatively small compared to the required temperature, no heat exchanger equipment is considered in the upstream part of the plant. Furthermore, the feed gas is injected into a large volume of liquid, which allows for fast adjustment of the temperature.</p>
<p>The biological methanation reactor is a tall continuously stirred tank reactor (CSTR) unit with multiple stirrer stages. The design of the sizing of the unit is given by the data in the work of Inkeri et al. (<xref ref-type="bibr" rid="B23">Inkeri et al., 2018</xref>), summarised in <xref ref-type="table" rid="T2">Table 2</xref>. The rate-based model developed in this study was validated by comparison of model predictions with literature data. The biological methanation kinetics itself are based on the growth of the archaeon <italic>Methanobacter thermoautotrophicum</italic> obtained from the work of Schill et al. (<xref ref-type="bibr" rid="B35">Schill et al., 1996</xref>). The up-scaled reactor system operates at a temperature of 60&#x2013;65&#xb0;C.</p>
<table-wrap id="T2" position="float">
<label>TABLE 2</label>
<caption>
<p>Data used to estimate the performance of the biological methanation system according to Inkeri et al. (<xref ref-type="bibr" rid="B23">Inkeri et al., 2018</xref>). The methane production rate was corrected according to <xref ref-type="disp-formula" rid="e2">Eqs 2</xref>, <xref ref-type="disp-formula" rid="e3">3</xref>.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="left">Electrolyser power P<sub>ael</sub> [MW<sub>el</sub>]</th>
<th align="center">Reactor volume V<sub>re</sub> [m<sup>3</sup>]</th>
<th align="center">Reactor diameter d<sub>re</sub> [m]</th>
<th align="center">Reactor height h [m]</th>
<th align="center">Impeller diameter d<sub>im</sub> [m]</th>
<th align="center">Stirring power P<sub>stirrer</sub> [kW]</th>
<th align="center">Product gas flow (methane) V<sub>CH4,&#xa0;Prod.</sub>&#xa0;[m<sup>3</sup>/&#xa0;h]</th>
<th align="center">Heat&#xa0;production Q [MW]</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="left">1.0</td>
<td align="char" char=".">8</td>
<td align="char" char=".">1.0</td>
<td align="char" char=".">10</td>
<td align="char" char=".">0.5</td>
<td align="char" char=".">11</td>
<td align="char" char=".">54</td>
<td align="char" char=".">0.16</td>
</tr>
<tr>
<td align="left">6.0</td>
<td align="char" char=".">115</td>
<td align="char" char=".">2.4</td>
<td align="char" char=".">24.5</td>
<td align="char" char=".">1.2</td>
<td align="char" char=".">56</td>
<td align="char" char=".">324</td>
<td align="char" char=".">0.97</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>The reactor geometry is fixed for both scales to 1:2:20 for the stirrer diameter (<inline-formula id="inf1">
<mml:math id="m1">
<mml:mrow>
<mml:msub>
<mml:mi>d</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>m</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>), reactor diameter (<inline-formula id="inf2">
<mml:math id="m2">
<mml:mrow>
<mml:msub>
<mml:mi>d</mml:mi>
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mi>e</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>) and reactor height (<inline-formula id="inf3">
<mml:math id="m3">
<mml:mi>h</mml:mi>
</mml:math>
</inline-formula>). The stirrer furthermore consists of 10 impeller stages along the height. The study of Inkeri et al. furthermore provides data for the reactor volume (<inline-formula id="inf4">
<mml:math id="m4">
<mml:mrow>
<mml:msub>
<mml:mi>V</mml:mi>
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mi>e</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>), the power required for stirring (<inline-formula id="inf5">
<mml:math id="m5">
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mi>t</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>r</mml:mi>
<mml:mi>r</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>r</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>), the methane product gas flow rate (<inline-formula id="inf6">
<mml:math id="m6">
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mover accent="true">
<mml:mi>V</mml:mi>
<mml:mo>&#x2d9;</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mrow>
<mml:mi>C</mml:mi>
<mml:msub>
<mml:mi>H</mml:mi>
<mml:mn>4</mml:mn>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:mi>P</mml:mi>
<mml:mi>r</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>d</mml:mi>
<mml:mo>.</mml:mo>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>) and the heat production of the reactor system (<inline-formula id="inf7">
<mml:math id="m7">
<mml:mi>Q</mml:mi>
</mml:math>
</inline-formula>).</p>
<p>In biological methanation, the reaction is catalysed by one or multiple strains of microorganisms. In addition to the suitable temperature and pressure conditions, those organisms also have to be provided with nutrients in order to foster their growth. Therefore, these nutrients are provided as solution to the reactor. To avoid the accumulation of potentially toxic metabolites in the liquid phase, a part of the slurry has to be discarded as wastewater. This wastewater stream furthermore contains the water produced by the methanation reaction according to <xref ref-type="disp-formula" rid="e1">Eq. 1</xref>. To ensure a proper distribution of nutrients and microorganisms along the reactor height, a part of the slurry is pumped from the bottom of the reactor to the top. Inkeri et al. provides a liquid feed rate D<sub>L</sub> of 0.05&#xa0;h<sup>&#x2212;1</sup> and a recycle ratio of 0.75, which allows the calculation of the required nutrient feed.<disp-formula id="e1">
<mml:math id="m8">
<mml:mrow>
<mml:mn>4</mml:mn>
<mml:mo>&#xa0;</mml:mo>
<mml:msub>
<mml:mi>H</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>C</mml:mi>
<mml:msub>
<mml:mi>O</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#x2194;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mi>C</mml:mi>
<mml:msub>
<mml:mi>H</mml:mi>
<mml:mn>4</mml:mn>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>2</mml:mn>
<mml:mo>&#xa0;</mml:mo>
<mml:msub>
<mml:mi>H</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:mi>O</mml:mi>
<mml:mo>,</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mi>&#x394;</mml:mi>
<mml:msub>
<mml:mi>H</mml:mi>
<mml:mi>R</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>165</mml:mn>
<mml:mo>&#xa0;</mml:mo>
<mml:mi>k</mml:mi>
<mml:mi>J</mml:mi>
<mml:mo>&#xa0;</mml:mo>
<mml:mi>m</mml:mi>
<mml:mi>o</mml:mi>
<mml:msup>
<mml:mi>l</mml:mi>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msup>
</mml:mrow>
</mml:math>
<label>(1)</label>
</disp-formula>
</p>
<p>The wastewater flow also contains a small amount of dissolved gases. In order to avoid insufficient CO<sub>2</sub> supply at the top of the reactor, 1% excess CO<sub>2</sub> is fed to the system (<xref ref-type="bibr" rid="B23">Inkeri et al., 2018</xref>). For the calculations, a H<sub>2</sub>/CO<sub>2</sub> ratio of 3.96 is chosen. The methane production rates published by Inkeri et al. are slightly above the stoichiometric possible values. They were corrected for this work by <xref ref-type="disp-formula" rid="e2">Eqs 2</xref>, <xref ref-type="disp-formula" rid="e3">3</xref> and listed in <xref ref-type="table" rid="T2">Table 2</xref>. We assumed that the CO<sub>2</sub> is fully dissolved in the liquid phase and the raw biomethane leaving the reactor only contains CH<sub>4</sub> and unreacted H<sub>2</sub>. This assumption can be confirmed by measurements at the Solothurn pilot plant (<xref ref-type="bibr" rid="B22">Huettenrauch et al., 2020</xref>).<disp-formula id="e2">
<mml:math id="m9">
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mover accent="true">
<mml:mi>n</mml:mi>
<mml:mo>&#x2d9;</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>H</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>e</mml:mi>
<mml:mi>l</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>c</mml:mi>
<mml:mi>t</mml:mi>
<mml:mi>r</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>l</mml:mi>
<mml:mi>y</mml:mi>
<mml:mi>s</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>r</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x22c5;</mml:mo>
<mml:msub>
<mml:mi>&#x3b7;</mml:mi>
<mml:mrow>
<mml:mi>P</mml:mi>
<mml:mi>E</mml:mi>
<mml:mi>M</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>L</mml:mi>
<mml:mi>H</mml:mi>
<mml:mi>V</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mi>L</mml:mi>
<mml:mi>H</mml:mi>
<mml:msub>
<mml:mi>V</mml:mi>
<mml:mrow>
<mml:msub>
<mml:mi>H</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:mrow>
</mml:msub>
<mml:mo>&#x22c5;</mml:mo>
<mml:mi>M</mml:mi>
<mml:msub>
<mml:mi>W</mml:mi>
<mml:mrow>
<mml:msub>
<mml:mi>H</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
<label>(2)</label>
</disp-formula>
<disp-formula id="e3">
<mml:math id="m10">
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mover accent="true">
<mml:mi>n</mml:mi>
<mml:mo>&#x2d9;</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mrow>
<mml:mi>C</mml:mi>
<mml:msub>
<mml:mi>H</mml:mi>
<mml:mn>4</mml:mn>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:mi>b</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>m</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>t</mml:mi>
<mml:mi>h</mml:mi>
<mml:mi>a</mml:mi>
<mml:mi>n</mml:mi>
<mml:mi>e</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mover accent="true">
<mml:mi>n</mml:mi>
<mml:mo>&#x2d9;</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>H</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:mrow>
</mml:msub>
<mml:mfrac>
<mml:mn>1</mml:mn>
<mml:mrow>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>x</mml:mi>
<mml:mrow>
<mml:mi>C</mml:mi>
<mml:msub>
<mml:mi>H</mml:mi>
<mml:mn>4</mml:mn>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:mi>b</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>m</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>t</mml:mi>
<mml:mi>h</mml:mi>
<mml:mi>a</mml:mi>
<mml:mi>n</mml:mi>
<mml:mi>e</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>x</mml:mi>
<mml:mrow>
<mml:mi>C</mml:mi>
<mml:msub>
<mml:mi>H</mml:mi>
<mml:mn>4</mml:mn>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:mi>b</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>m</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>t</mml:mi>
<mml:mi>h</mml:mi>
<mml:mi>a</mml:mi>
<mml:mi>n</mml:mi>
<mml:mi>e</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>4</mml:mn>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
<label>(3)</label>
</disp-formula>
</p>
<p>The hydrogen production rate <inline-formula id="inf8">
<mml:math id="m11">
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mover accent="true">
<mml:mi>n</mml:mi>
<mml:mo>&#x2d9;</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>H</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is calculated from the nominal power of the electrolyser <inline-formula id="inf9">
<mml:math id="m12">
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>e</mml:mi>
<mml:mi>l</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>c</mml:mi>
<mml:mi>t</mml:mi>
<mml:mi>r</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>l</mml:mi>
<mml:mi>y</mml:mi>
<mml:mi>s</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>r</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> and its efficiency <inline-formula id="inf10">
<mml:math id="m13">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3b7;</mml:mi>
<mml:mrow>
<mml:mi>P</mml:mi>
<mml:mi>E</mml:mi>
<mml:mi>M</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>. The efficiency is given by 65%, based on the lower heating value (LHV), a value applying for alkaline as well as PEM technology (<xref ref-type="bibr" rid="B4">Bertuccioli et al., 2014</xref>; <xref ref-type="bibr" rid="B53">Witte et al., 2018b</xref>; <xref ref-type="bibr" rid="B23">Inkeri et al., 2018</xref>).</p>
<p>The product gas stream contains around 98% CH<sub>4</sub> and fulfils grid injection limitations. Nevertheless, it is saturated with water and contaminated with volatile metabolites from the biological process. Therefore, a temperature swing adsorption (TSA) process for drying and a sorption step for desulphurization is considered. The TSA consists of two vessels operated alternatingly, one in the main product gas stream in adsorption mode and the other in regeneration mode. The adsorption process occurs at ambient temperature, whereas for desorption the vessel is electrically heated to 200&#xb0;C and flushed with a part of the dried product gas stream. This water-saturated stream is then fed to a condenser and recompressed to be fed to the main stream again. This procedure allows to regenerate the sorbent without losing product gas or contaminating the product gas with other components, such as air or nitrogen. In the TSA, silica gel is used as a sorbent.</p>
<p>The desulphurisation step consists also of two vessels which are filled with a specialised sorbent for biogas polishing. The first vessel is in operation mode and fed with the main gas stream whereas the second one serves as back up in case of maintenance or sorbent replacement.</p>
<p>The now cleaned and dried biomethane can be expanded to the required pressure and injected to the gas distribution grid.</p>
</sec>
<sec id="s2-1-2">
<title>2.1.2 Cases C1 and C6: Bubbling Fluidised Bed Methanation With PEM Electrolysis (PEM-BFB)</title>
<p>The process consisting of the bubbling fluidised bed methanation is designed as shown in the process flow diagram in <xref ref-type="fig" rid="F2">Figure 2</xref>. It consists of the same upstream equipment as the process in PEM-bio case and is based on a previous study by Witte et al. (<xref ref-type="bibr" rid="B53">Witte et al., 2018b</xref>) on direct methanation of biogas. Again, cleaned CO<sub>2</sub> and water for electrolysis are fed to the plant and compressed to the operation pressures of the electrolyser and reactor. Hydrogen is mixed to the main CO<sub>2</sub> stream. This feed gas has to be pre-heated to the minimal feed temperature of the reactor, which is set to 280&#xb0;C. This temperature is the minimum useful operation temperature of the nickel-based catalyst. The catalyst has to be protected from contamination with sulphur compounds and subsequent deactivation (<xref ref-type="bibr" rid="B51">Witte et al., 2019</xref>). Since biogenic CO<sub>2</sub> still contains a small amount of such contaminants a gas polishing step is introduced in the CO<sub>2</sub> feed stream (<xref ref-type="bibr" rid="B7">Calbry-Muzyka et al., 2019</xref>).</p>
<fig id="F2" position="float">
<label>FIGURE 2</label>
<caption>
<p>Process chain for bubbling fluidised bed methanation of CO<sub>2</sub> using a PEM electrolyser as hydrogen source.</p>
</caption>
<graphic xlink:href="fenrg-09-775259-g002.tif"/>
</fig>
<p>In order to avoid coking and subsequent deactivation of the catalyst, an over-stoichiometric ratio of H<sub>2</sub>/CO<sub>2</sub> of 4.01 is fed to the reactor. The product gas of the reactor is led to a heat exchanger and a condenser, which allow the valorisation of the off gas heat and removal of product water. The heat removed in these units is used to pre-heat the feed streams of CO<sub>2</sub> and H<sub>2</sub>.</p>
<p>The heterogeneous methanation reaction is limited by the thermodynamic equilibrium, which results in hydrogen concentrations of more than 2 vol.-% in the dry product gas. In order to fulfil grid injection limitations, excess H<sub>2</sub> has to be removed to allow the use of the biomethane in CNG cars whose tanks are certified for 2% H<sub>2</sub> fraction at maximum. As previous work showed, membrane separation can be a cost efficient way to do so (<xref ref-type="bibr" rid="B52">Witte et al., 2018a</xref>). We implemented a counter-current hollow fibre membrane upgrading system to obtain full grid compliance. The resulting H<sub>2</sub>-rich recycle stream, which also contains CH<sub>4</sub> and CO<sub>2</sub>, is fed to the CO<sub>2</sub> feed stream before the main compressor. This recycle stream directly influences the gas composition in the reactor feed, which demands an adjustment of the CO<sub>2</sub> feed stream in such way that the required H<sub>2</sub>/CO<sub>2</sub> is fulfilled.</p>
</sec>
<sec id="s2-1-3">
<title>2.1.3 Cases C1S and C6S: Bubbling Fluidised Bed Methanation With High Temperature Electrolysis (SOE-BFB)</title>
<p>This case represents a modification of the previous system: the PEM electrolyser is replaced by a high-temperature solid oxide electrolyser system. This allows the comparison of the cost structure, when BFB methanation is combined with an SOE system. Instead of liquid water, steam is used for electrolysis, which allows the provision of heat by an external heat source. The stack itself can be operated at a thermoneutral voltage, where the heat flux from the stack is zero (<xref ref-type="bibr" rid="B19">Giglio et al., 2015</xref>). This allows a very high electrical efficiency of the stack, as a part of the energy for the electrolysis is provided by the thermal energy of the steam (<xref ref-type="bibr" rid="B6">Buttler and Spliethoff, 2018</xref>). In this study, a stack efficiency <inline-formula id="inf11">
<mml:math id="m14">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3b7;</mml:mi>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mi>t</mml:mi>
<mml:mi>a</mml:mi>
<mml:mi>c</mml:mi>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> of 89% (LHV) and a steam conversion efficiency of 90% were assumed (<xref ref-type="bibr" rid="B21">Gruber et al., 2018</xref>).</p>
<p>Such SOE systems are under experimental investigation in pilot scale, such as in the <xref ref-type="bibr" rid="B24">Helmeth project (2021)</xref>. In this project, a hot-box is used to keep the temperature of the stack constant and to allow heat integration between the feed streams and hot product gas streams. Furthermore, this box allows the pressurisation of the sensitive high temperature components and enables an operation pressure up to 15&#xa0;bar <xref ref-type="bibr" rid="B24">Helmeth project (2021)</xref>. For this work, a similar design is assumed, where heat exchangers allow the pre-heating of the feed water and superheating of the steam by the product gas streams (H<sub>2</sub> &#x2b; H<sub>2</sub>O and O<sub>2</sub>). The heat required for evaporation of the feed water is provided by the off-heat from the methanation reactor.</p>
<p>Instead of scaling this process by the electrical input, the electrolyser was scaled by the same hydrogen output as in the previous cases. This allows a direct comparison of the plant cost structure and indicates the benefits of including a SOE system in the methanation process.</p>
<p>All other process units are designed the same way as in case PEM-BFB.</p>
</sec>
</sec>
<sec id="s2-2">
<title>2.2 Main Models</title>
<sec id="s2-2-1">
<title>2.2.1 Bubbling Fluidized Bed Reactor</title>
<p>This work considers an internally cooled, bubbling fluidised bed reactor with a nickel-based methanation catalyst. The model of this unit is based on a pseudo-homogeneous two-phase approach according to Kopyscinski et al. (<xref ref-type="bibr" rid="B26">Kopyscinski et al., 2011</xref>), which applies hydrodynamic correlations of freely bubbling fluidised beds. This model was already modified and used in previous studies by Witte et al. (<xref ref-type="bibr" rid="B53">Witte et al., 2018b</xref>). A more detailed description can be found in the <xref ref-type="sec" rid="s9">Supplementary Material</xref>.</p>
<p>The reaction kinetics implemented in the reactor model are expressed as Langmuir-Hinshelwood-type rate equations. The kinetics of CO<sub>2</sub> methanation are expressed as a water-gas-shift (WGS) reaction coupled with CO methanation as shown in <xref ref-type="disp-formula" rid="e4">Eqs 4</xref>, <xref ref-type="disp-formula" rid="e5">5</xref>:<disp-formula id="e4">
<mml:math id="m15">
<mml:mrow>
<mml:msub>
<mml:mi>r</mml:mi>
<mml:mrow>
<mml:mi>W</mml:mi>
<mml:mi>G</mml:mi>
<mml:mi>S</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>k</mml:mi>
<mml:mrow>
<mml:mi>W</mml:mi>
<mml:mi>G</mml:mi>
<mml:mi>S</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mi>&#x3b1;</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mrow>
<mml:mi>C</mml:mi>
<mml:mi>O</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mrow>
<mml:mi>H</mml:mi>
<mml:mn>2</mml:mn>
<mml:mi>O</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mrow>
<mml:mi>C</mml:mi>
<mml:mi>O</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mrow>
<mml:mi>H</mml:mi>
<mml:mn>2</mml:mn>
<mml:mi>O</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mrow>
<mml:mi>e</mml:mi>
<mml:mi>q</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>W</mml:mi>
<mml:mi>G</mml:mi>
<mml:mi>S</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:msubsup>
<mml:mi>p</mml:mi>
<mml:mrow>
<mml:mi>H</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mn>0.5</mml:mn>
</mml:mrow>
</mml:msubsup>
<mml:msup>
<mml:mrow>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mi>C</mml:mi>
</mml:msub>
<mml:msubsup>
<mml:mi>p</mml:mi>
<mml:mrow>
<mml:mi>C</mml:mi>
<mml:mi>O</mml:mi>
</mml:mrow>
<mml:mi>d</mml:mi>
</mml:msubsup>
<mml:msubsup>
<mml:mi>p</mml:mi>
<mml:mrow>
<mml:mi>H</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
<mml:mi>e</mml:mi>
</mml:msubsup>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mrow>
<mml:mi>O</mml:mi>
<mml:mi>H</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mrow>
<mml:mi>H</mml:mi>
<mml:mn>2</mml:mn>
<mml:mi>O</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msubsup>
<mml:mi>p</mml:mi>
<mml:mrow>
<mml:mi>H</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>0.5</mml:mn>
</mml:mrow>
</mml:msubsup>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
<label>(4)</label>
</disp-formula>
<disp-formula id="e5">
<mml:math id="m16">
<mml:mrow>
<mml:msub>
<mml:mi>r</mml:mi>
<mml:mrow>
<mml:mi>M</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>t</mml:mi>
<mml:mi>h</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>k</mml:mi>
<mml:mrow>
<mml:mi>M</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>t</mml:mi>
<mml:mi>h</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mi>C</mml:mi>
</mml:msub>
<mml:msubsup>
<mml:mi>p</mml:mi>
<mml:mrow>
<mml:mi>C</mml:mi>
<mml:mi>O</mml:mi>
</mml:mrow>
<mml:mi>a</mml:mi>
</mml:msubsup>
<mml:msubsup>
<mml:mi>p</mml:mi>
<mml:mrow>
<mml:mi>H</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
<mml:mi>b</mml:mi>
</mml:msubsup>
<mml:msup>
<mml:mrow>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mrow>
<mml:mi>C</mml:mi>
<mml:mi>H</mml:mi>
<mml:mn>4</mml:mn>
</mml:mrow>
</mml:msub>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mrow>
<mml:mi>H</mml:mi>
<mml:mn>2</mml:mn>
<mml:mi>O</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mrow>
<mml:mi>e</mml:mi>
<mml:mi>q</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>M</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>t</mml:mi>
<mml:mi>h</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mrow>
<mml:mi>C</mml:mi>
<mml:mi>O</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msubsup>
<mml:mi>p</mml:mi>
<mml:mrow>
<mml:mi>H</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
<mml:mn>3</mml:mn>
</mml:msubsup>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mi>c</mml:mi>
</mml:msup>
</mml:mrow>
<mml:mrow>
<mml:msup>
<mml:mrow>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mi>C</mml:mi>
</mml:msub>
<mml:msubsup>
<mml:mi>p</mml:mi>
<mml:mrow>
<mml:mi>C</mml:mi>
<mml:mi>O</mml:mi>
</mml:mrow>
<mml:mi>d</mml:mi>
</mml:msubsup>
<mml:msubsup>
<mml:mi>p</mml:mi>
<mml:mrow>
<mml:mi>H</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
<mml:mi>e</mml:mi>
</mml:msubsup>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mrow>
<mml:mi>O</mml:mi>
<mml:mi>H</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mrow>
<mml:mi>H</mml:mi>
<mml:mn>2</mml:mn>
<mml:mi>O</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msubsup>
<mml:mi>p</mml:mi>
<mml:mrow>
<mml:mi>H</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>0.5</mml:mn>
</mml:mrow>
</mml:msubsup>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
<label>(5)</label>
</disp-formula>
</p>
<p>The temperature dependence of rate and adsorption coefficients are based on the Arrhenius and Van&#x2019;t Hoff approach according to <xref ref-type="disp-formula" rid="e6">Eqs 6</xref>, <xref ref-type="disp-formula" rid="e7">7</xref>:<disp-formula id="e6">
<mml:math id="m17">
<mml:mrow>
<mml:msub>
<mml:mi>k</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi>k</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>T</mml:mi>
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>f</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:msub>
<mml:mi>e</mml:mi>
<mml:mi>x</mml:mi>
<mml:mi>p</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mrow>
<mml:mi>A</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mi>R</mml:mi>
<mml:msub>
<mml:mi>T</mml:mi>
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>f</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>T</mml:mi>
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>f</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mi>T</mml:mi>
</mml:mfrac>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>,</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mi>i</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mi>W</mml:mi>
<mml:mi>G</mml:mi>
<mml:mi>S</mml:mi>
<mml:mo>,</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mi>M</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>t</mml:mi>
<mml:mi>h</mml:mi>
</mml:mrow>
</mml:math>
<label>(6)</label>
</disp-formula>
<disp-formula id="e7">
<mml:math id="m18">
<mml:mrow>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mi>j</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mrow>
<mml:mi>j</mml:mi>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>T</mml:mi>
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>f</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:msub>
<mml:mi>e</mml:mi>
<mml:mi>x</mml:mi>
<mml:mi>p</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mi>&#x394;</mml:mi>
<mml:msub>
<mml:mi>H</mml:mi>
<mml:mi>j</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mi>R</mml:mi>
<mml:msub>
<mml:mi>T</mml:mi>
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>f</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>T</mml:mi>
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>f</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mi>T</mml:mi>
</mml:mfrac>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>,</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mi>j</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mi>&#x3b1;</mml:mi>
<mml:mo>,</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mi>C</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>O</mml:mi>
<mml:mi>H</mml:mi>
<mml:mo>&#x2026;</mml:mo>
</mml:mrow>
</mml:math>
<label>(7)</label>
</disp-formula>
</p>
<p>The values of kinetic parameters are listed in <xref ref-type="sec" rid="s9">Supplementary Table S1</xref>.</p>
<p>Hydrodynamics in the reactor is the second main influence on the performance. This influence is expressed through bubble size correlations, from which bubble rise velocities and the bubble hold-up <inline-formula id="inf28">
<mml:math id="m35">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3b5;</mml:mi>
<mml:mi>b</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> can be determined. These values allow the estimation of the total bubble surface area available for inter-phase mass transfer.</p>
<p>A variety of bubble size correlations exist in literature, but they are only valid for conditions with no heat exchanger internals present in the reactor. Therefore, the bubble growth correlation with the smallest overall bubble size <inline-formula id="inf29">
<mml:math id="m36">
<mml:mrow>
<mml:msub>
<mml:mi>d</mml:mi>
<mml:mi>B</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> by Werther et al. (<xref ref-type="bibr" rid="B49">Werther, 1976</xref>) was used in the model and it was further modified by a factor of 0.5 in order to account for the internals present in the reactor that significantly inhibit the bubble growth in the reactor (<xref ref-type="bibr" rid="B29">Maurer et al., 2016</xref>; <xref ref-type="bibr" rid="B36">Schillinger et al., 2017</xref>).<disp-formula id="e8">
<mml:math id="m37">
<mml:mrow>
<mml:msub>
<mml:mi>d</mml:mi>
<mml:mi>B</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>0.5</mml:mn>
<mml:mo>&#x22c5;</mml:mo>
<mml:mn>0.835</mml:mn>
<mml:msup>
<mml:mrow>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>0.272</mml:mn>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>u</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mrow>
<mml:mi>m</mml:mi>
<mml:mi>f</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>/</mml:mo>
<mml:mn>3</mml:mn>
</mml:mrow>
</mml:msup>
<mml:msup>
<mml:mrow>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>0.0684</mml:mn>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:mn>1.21</mml:mn>
</mml:mrow>
</mml:msup>
</mml:mrow>
</mml:math>
<label>(8)</label>
</disp-formula>
</p>
<p>As shown in <xref ref-type="disp-formula" rid="e8">Eq. 8</xref>, the bubble diameter <inline-formula id="inf30">
<mml:math id="m38">
<mml:mrow>
<mml:msub>
<mml:mi>d</mml:mi>
<mml:mi>B</mml:mi>
</mml:msub>
<mml:mo>&#xa0;</mml:mo>
</mml:mrow>
</mml:math>
</inline-formula> is dependent on the gas velocity <inline-formula id="inf31">
<mml:math id="m39">
<mml:mi>u</mml:mi>
</mml:math>
</inline-formula> in the reactor and the gas velocity required for minimal fluidisation of the bed <inline-formula id="inf32">
<mml:math id="m40">
<mml:mrow>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mrow>
<mml:mi>m</mml:mi>
<mml:mi>f</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>.</p>
<p>The reactor model is based on the following assumptions:<list list-type="simple">
<list-item>
<label>-</label>
<p>Steady-state conditions and ideal gas behaviour</p>
</list-item>
<list-item>
<label>-</label>
<p>No reaction in the bubble phase</p>
</list-item>
<list-item>
<label>-</label>
<p>Laminar boundary layers around the catalyst particles and influence of pore diffusion are neglected. Therefore, the gas concentration in and on the particles equals the dense phase concentrations.</p>
</list-item>
<list-item>
<label>-</label>
<p>Radial gas concentrations are neglected</p>
</list-item>
</list>
</p>
<p>Catalyst deactivation is not considered in the model. The economic analysis includes an annual replacement of the full catalyst load in order to compensate for minor deactivation effects. Field experiments showed that a stable operation of fluidised bed methanation in biogas is possible for more than 1000&#xa0;h (<xref ref-type="bibr" rid="B7">Calbry-Muzyka et al., 2019</xref>; <xref ref-type="bibr" rid="B51">Witte et al., 2019</xref>). In order to account for deactivation effects and the additional reactor height needed for particle settling, the reactor height and catalyst mass was multiplied by a design factor of 1.5 prior to cost calculations.</p>
</sec>
<sec id="s2-2-2">
<title>2.2.2 Biological Methanation Reactor</title>
<p>Data published by Inkeri et al. (<xref ref-type="bibr" rid="B23">Inkeri et al., 2018</xref>) was used to estimate the process performance in the work presented here. The study chose the scale of well-known demonstration projects for methanation, namely the Danish BioCat project (<xref ref-type="bibr" rid="B12">Electrochaea.dk ApS, 2014</xref>), the Audi e-gas project (<xref ref-type="bibr" rid="B56">Specht et al., 2016</xref>) in Germany, as well as a large electrolysis plant in Finland, which is not considered in this work. The BioCat project consists of a biological methanation reactor, which is fed by a 1&#xa0;MW<sub>el</sub> electrolysis unit. The Audi e-gas project is a pilot plant consisting of a fixed bed methanation reactor in the 3&#xa0;MW<sub>CH4</sub> scale.</p>
<p>The work of Inkeri et al. estimates the gas-liquid mass transfer by a lumped mass transfer coefficient, taking into account the reactor height, stirring energy and stirrer geometry. Furthermore, the influence of physical properties of gas and liquid phase components on the bubble size and subsequently the gas-liquid mass transfer are approximated. This approach, based on simplified fundamental equations describing physical phenomena allows a scale-up of the reactor model.</p>
<p>Inkeri et al. verified the scale-up of their model using estimated dimensions of the plant in Aved&#xf8;re (Denmark). In the meantime another pilot plant by the same technology was operated in Solothurn (Switzerland) as part of the Store&#x26;Go-project (<xref ref-type="bibr" rid="B40">Store and Go Project, 2016b</xref>). This plant was designed for a hydrogen input corresponding to a 700&#xa0;kW<sub>el</sub> electrolysis unit, which required a vessel volume of around 7&#xa0;m<sup>3</sup>. Furthermore, a second biological methanation plant of 2.5&#xa0;MW<sub>el</sub> is currently under construction in Switzerland. This plant uses the technology of microbEnergy (<xref ref-type="bibr" rid="B30">microbEnergy GmbH, 2019</xref>), which is also based on a stirred bubble column, but makes use of strains of microorganisms already present in the waste water treatment plant. Limited technical information is accessible, but based on the available illustrations, we estimate the reactor size to around 45&#xa0;m<sup>3</sup> (<xref ref-type="bibr" rid="B27">Limeco, 2020</xref>). These values indicate that the dimensions given by Inkeri et al. are in a valid range for such systems.</p>
<p>Biological systems, where substantial gas formation occurs in the liquid phase, tend to form large amounts of foam. This foam formation has to be minimised, as it can cause severe problems in the downstream part of the plant. Usually, anti-foaming agents are added to the slurry and technical equipment for de-foaming is used (<xref ref-type="bibr" rid="B41">Store and Go Project, 2016a</xref>). These measures are beyond the scope of this work, nevertheless, we assumed an additional 2&#xa0;m of reactor height for foam protection.</p>
<p>The data used is summarized in <xref ref-type="table" rid="T2">Table 2</xref>.</p>
</sec>
<sec id="s2-2-3">
<title>2.2.3 Hollow Fibre Gas Separation Membrane</title>
<p>For the gas separation unit required for the removal of excess hydrogen after BFB methanation, a rate-based model from literature was implemented. Makaruk et al. published an algorithm simulating a counter-current hollow fibre membrane (<xref ref-type="bibr" rid="B28">Makaruk and Harasek, 2009</xref>). This algorithm is based on an iterative finite difference Gau&#xdf;&#x2013;Seidel method. A more detailed description can be found in the <xref ref-type="sec" rid="s9">Supplementary Material</xref>.</p>
<p>The counter-current membrane model is based on the following assumptions:<list list-type="simple">
<list-item>
<label>-</label>
<p>Ideal gas behaviour</p>
</list-item>
<list-item>
<label>-</label>
<p>Axial dispersion is neglected</p>
</list-item>
<list-item>
<label>-</label>
<p>Pressure drop on bore and shell side is neglected</p>
</list-item>
<list-item>
<label>-</label>
<p>The influence of the gas composition on the permeability is neglected. Therefore, the permeability parameter is set constant.</p>
</list-item>
</list>
</p>
<p>In this work, a constant value for membrane permeance is introduced, based on the commercial polymer type Matrimid. Permeability data for the relevant species were taken from Zhang et al. (<xref ref-type="bibr" rid="B54">Zhang et al., 2008</xref>). Since there was no permeability determined for water, we assumed it to be equal to the permeability of hydrogen. This is justified by other commercial biogas upgrading membranes, which indicate a higher permeability for water than for hydrogen (<xref ref-type="bibr" rid="B14">Evonik Industries AG SEPURAN, 2021</xref>). In order to obtain the permeance parameter, the permeability has to be divided by the effective thickness of the separation layer. The actual separation layer is only a fraction of the wall thickness of the membrane fibre and has to be estimated for this case. These thicknesses range from 500 to even less than 100&#xa0;nm (<xref ref-type="bibr" rid="B2">Baker, 2002</xref>). In this work, we assumed the effective thickness to be 100&#xa0;nm.</p>
</sec>
<sec id="s2-2-4">
<title>2.2.4 Auxiliary Units</title>
<p>For all auxiliary process units, short-cut models are used, which are based on correlations (<xref ref-type="bibr" rid="B53">Witte et al., 2018b</xref>).</p>
<p>For the PEM electrolysers, a specific electricity consumption of 4.6&#xa0;kWh/m<sup>3</sup> (H<sub>2</sub>, NTP) is considered, which is an average value for PEM and alkaline technology (<xref ref-type="bibr" rid="B6">Buttler and Spliethoff, 2018</xref>). This value corresponds to an efficiency of 65% (LHV). The data used to calculate the stack performance of SOE is given in <xref ref-type="table" rid="T3">Table 3</xref>.</p>
<table-wrap id="T3" position="float">
<label>TABLE 3</label>
<caption>
<p>Input parameters for the bubbling fluidised bed-based processes for CO<sub>2</sub> methanation, including the parameters for the sub-systems reactor, membrane upgrading and electrolyser.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="left"/>
<th align="left"/>
<th align="left"/>
<th align="center">C1</th>
<th align="center">C1S</th>
<th align="center">C6</th>
<th align="center">C6S</th>
<th align="center">C1bg</th>
<th align="center">C1Sbg</th>
<th align="center">C6bg</th>
<th align="center">C6Sbg</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td colspan="11" align="left">Reactor</td>
</tr>
<tr>
<td align="left">&#x2003;Reactor temperature</td>
<td align="center">T<sub>r</sub>
</td>
<td align="center">&#xb0;C</td>
<td align="center">360</td>
<td align="char" char=".">360</td>
<td align="center">360</td>
<td align="char" char=".">360</td>
<td align="center">360</td>
<td align="char" char=".">360</td>
<td align="center">360</td>
<td align="char" char=".">360</td>
</tr>
<tr>
<td align="left">&#x2003;Temperature of the gas feed</td>
<td align="center">T<sub>r,feed</sub>
</td>
<td align="center">&#xb0;C</td>
<td align="center">280</td>
<td align="char" char=".">280</td>
<td align="center">280</td>
<td align="char" char=".">280</td>
<td align="center">280</td>
<td align="char" char=".">280</td>
<td align="center">280</td>
<td align="char" char=".">280</td>
</tr>
<tr>
<td align="left">&#x2003;System pressure</td>
<td align="center">p<sub>react</sub>
</td>
<td align="center">bara</td>
<td align="center">10</td>
<td align="char" char=".">10</td>
<td align="center">10</td>
<td align="char" char=".">10</td>
<td align="center">10</td>
<td align="char" char=".">10</td>
<td align="center">10</td>
<td align="char" char=".">10</td>
</tr>
<tr>
<td align="left">&#x2003;Diameter of heat exchanger tubes</td>
<td align="center">d<sub>HEXtubes</sub>
</td>
<td align="center">m</td>
<td align="center">0.02</td>
<td align="char" char=".">0.02</td>
<td align="center">0.02</td>
<td align="char" char=".">0.02</td>
<td align="center">0.02</td>
<td align="char" char=".">0.02</td>
<td align="center">0.02</td>
<td align="char" char=".">0.02</td>
</tr>
<tr>
<td align="left">&#x2003;Feed ratio H<sub>2</sub>/CO<sub>2</sub>
</td>
<td align="center">r<sub>H2/CO2</sub>
</td>
<td align="center">mol/mol</td>
<td align="center">4.01</td>
<td align="char" char=".">4.01</td>
<td align="center">4.01</td>
<td align="char" char=".">4.01</td>
<td align="center">4.01</td>
<td align="char" char=".">4.01</td>
<td align="center">4.01</td>
<td align="char" char=".">4.01</td>
</tr>
<tr>
<td align="left">&#x2003;Condenser temperature</td>
<td align="center">T<sub>cond</sub>
</td>
<td align="center">&#xb0;C</td>
<td align="center">30</td>
<td align="char" char=".">30</td>
<td align="center">30</td>
<td align="char" char=".">30</td>
<td align="center">30</td>
<td align="char" char=".">30</td>
<td align="center">30</td>
<td align="char" char=".">30</td>
</tr>
<tr>
<td align="left">Membrane module (counter-current)</td>
<td align="left"/>
<td align="left"/>
<td align="left"/>
<td align="left"/>
<td align="left"/>
<td align="left"/>
<td align="left"/>
<td align="left"/>
<td align="left"/>
<td align="left"/>
</tr>
<tr>
<td align="left">&#x2003;Number of modules (4 inch)</td>
<td align="left"/>
<td align="center">&#x2014;</td>
<td align="center">2</td>
<td align="char" char=".">2</td>
<td align="center">11</td>
<td align="char" char=".">11</td>
<td align="center">4</td>
<td align="char" char=".">4</td>
<td align="center">22</td>
<td align="char" char=".">22</td>
</tr>
<tr>
<td colspan="11" align="left">Electrolysis</td>
</tr>
<tr>
<td align="left">&#x2003;Temperature of the SOE</td>
<td align="center">T<sub>soec</sub>
</td>
<td align="center">&#xb0;C</td>
<td align="center">n/a</td>
<td align="char" char=".">800</td>
<td align="center">n/a</td>
<td align="char" char=".">800</td>
<td align="center">n/a</td>
<td align="char" char=".">800</td>
<td align="center">n/a</td>
<td align="char" char=".">800</td>
</tr>
<tr>
<td align="left">&#x2003;El. Efficiency stack (LHV)</td>
<td align="center">
<inline-formula id="inf44">
<mml:math id="m55">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3b7;</mml:mi>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mi>t</mml:mi>
<mml:mi>a</mml:mi>
<mml:mi>c</mml:mi>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>
</td>
<td align="center">&#x2014;</td>
<td align="center">0.65</td>
<td align="char" char=".">0.89</td>
<td align="center">0.65</td>
<td align="char" char=".">0.89</td>
<td align="center">0.65</td>
<td align="char" char=".">0.89</td>
<td align="center">0.65</td>
<td align="char" char=".">0.89</td>
</tr>
<tr>
<td align="left">&#x2003;Steam conversion</td>
<td align="center">
<inline-formula id="inf45">
<mml:math id="m56">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3b7;</mml:mi>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mi>c</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>
</td>
<td align="center">&#x2014;</td>
<td align="center">n/a</td>
<td align="char" char=".">0.9</td>
<td align="center">n/a</td>
<td align="char" char=".">0.9</td>
<td align="center">n/a</td>
<td align="char" char=".">0.9</td>
<td align="center">n/a</td>
<td align="char" char=".">0.9</td>
</tr>
<tr>
<td align="left">&#x2003;Electrolyser pressure</td>
<td align="center">p<sub>electrolysis</sub>
</td>
<td align="center">bara</td>
<td align="center">30</td>
<td align="char" char=".">15</td>
<td align="center">30</td>
<td align="char" char=".">15</td>
<td align="center">30</td>
<td align="char" char=".">15</td>
<td align="center">30</td>
<td align="char" char=".">15</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>The TSA unit is based on a zeolite sorbent with a maximum capacity of 0.15&#xa0;kg<sub>H2O</sub>/kg<sub>sorbent</sub>. A cycle time of 8&#xa0;h and a regeneration temperature of 200&#xb0;C is assumed.</p>
<p>For the desulphurisation after the biological methanation and in the feed of the BFB methanation, we assume an inlet concentration of 50&#xa0;ppm H<sub>2</sub>S. Vessel size calculations are based on a conservative maximum load of 7&#xa0;wt.-% H<sub>2</sub>S on the sorbent. Further data is based on the commercially available sorbent Sulfatrap R7 (<xref ref-type="bibr" rid="B7">Calbry-Muzyka et al., 2019</xref>).</p>
<p>Pumps and compressors are based on thermodynamic correlations, with an electrical efficiency of 0.8.</p>
<p>The heat exchanger equipment was sized based on heat transfer coefficients available as correlations for different systems (<xref ref-type="bibr" rid="B48">Ulrich and Vasudevan, 2004</xref>). The exchangers are assumed to be shell and tube systems in counter-current operation.</p>
<p>All volumetric flow rates reported in this work refer to standard conditions at 1&#xa0;atm and 0&#xb0;C according to DIN 1343 (<xref ref-type="bibr" rid="B11">DIN, 1990</xref>).</p>
</sec>
</sec>
<sec id="s2-3">
<title>2.3 Cost Calculations</title>
<p>The comparison of the cost structure of the three processes is performed by determining the operational cost (OPEX) and capital cost. Based on the technical performance and size of the units, cost functions for the units are applied, which were obtained from Ulrich and Vasudevan (<xref ref-type="bibr" rid="B48">Ulrich and Vasudevan, 2004</xref>) and already part of previous work (<xref ref-type="bibr" rid="B52">Witte et al., 2018a</xref>). These functions have the following structure:<disp-formula id="e9">
<mml:math id="m41">
<mml:mrow>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mrow>
<mml:mi>B</mml:mi>
<mml:mi>M</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>P</mml:mi>
</mml:msub>
<mml:mo>&#x22c5;</mml:mo>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mrow>
<mml:mi>B</mml:mi>
<mml:mi>M</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mtext>,&#xa0;with&#xa0;</mml:mtext>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mrow>
<mml:mi>B</mml:mi>
<mml:mi>M</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mi>M</mml:mi>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mi>p</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(9)</label>
</disp-formula>
</p>
<p>The bare module cost <inline-formula id="inf33">
<mml:math id="m42">
<mml:mrow>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mrow>
<mml:mi>B</mml:mi>
<mml:mi>M</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> consists of the purchased equipment cost <inline-formula id="inf34">
<mml:math id="m43">
<mml:mrow>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>P</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>, which is multiplied by a bare module factor <inline-formula id="inf35">
<mml:math id="m44">
<mml:mrow>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mrow>
<mml:mi>B</mml:mi>
<mml:mi>M</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>. This factor is a function of the material factor <inline-formula id="inf36">
<mml:math id="m45">
<mml:mrow>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mi>M</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> and the pressure factor <inline-formula id="inf37">
<mml:math id="m46">
<mml:mrow>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mi>p</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>, which are given in literature and consider the cost increase by the use of special materials and elevated pressures. The accuracy of such cost evaluations is within the range of &#xb1;20% (<xref ref-type="bibr" rid="B48">Ulrich and Vasudevan, 2004</xref>).</p>
<p>In order to obtain the investment cost <inline-formula id="inf38">
<mml:math id="m47">
<mml:mrow>
<mml:msubsup>
<mml:mi>C</mml:mi>
<mml:mi>k</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>n</mml:mi>
<mml:mi>v</mml:mi>
</mml:mrow>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> for a unit <inline-formula id="inf39">
<mml:math id="m48">
<mml:mi>k</mml:mi>
</mml:math>
</inline-formula>, the cost for piping, instrumentation and further surrounding equipment have to be considered. Such costs are approximated by a factor of the bare module cost, as shown in <xref ref-type="disp-formula" rid="e10">Eq. 10</xref>. The plant cost factors <inline-formula id="inf40">
<mml:math id="m49">
<mml:mrow>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mrow>
<mml:mi>P</mml:mi>
<mml:mi>C</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> considered in this work are summarised in <xref ref-type="sec" rid="s9">Supplementary Table S2</xref> (<xref ref-type="bibr" rid="B17">Garrett, 1989</xref>; <xref ref-type="bibr" rid="B52">Witte et al., 2018a</xref>).<disp-formula id="e10">
<mml:math id="m50">
<mml:mrow>
<mml:msubsup>
<mml:mi>C</mml:mi>
<mml:mi>k</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>n</mml:mi>
<mml:mi>v</mml:mi>
</mml:mrow>
</mml:msubsup>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mrow>
<mml:mi>B</mml:mi>
<mml:mi>M</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2b;</mml:mo>
<mml:mstyle displaystyle="true">
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mrow>
<mml:mi>P</mml:mi>
<mml:mi>C</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mstyle>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x22c5;</mml:mo>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mrow>
<mml:mi>C</mml:mi>
<mml:mi>E</mml:mi>
<mml:mi>P</mml:mi>
<mml:mi>C</mml:mi>
<mml:mi>I</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
<label>(10)</label>
</disp-formula>
</p>
<p>The investment cost also include a factor <italic>F<sub>CEPCI</sub>
</italic>, which quantifies the cost increase of equipment (<xref ref-type="bibr" rid="B48">Ulrich and Vasudevan, 2004</xref>) according to <xref ref-type="disp-formula" rid="e11">Eq. 11</xref>. It consists of the Chemical Engineering plant cost index (CEPCI, <italic>I</italic>
<sub>2019</sub>), listed in <xref ref-type="sec" rid="s9">Supplementary Table S3</xref>. The reference value <italic>I</italic>
<sub>
<italic>ref</italic>(2004)</sub> was 400 in 2004 (<xref ref-type="bibr" rid="B48">Ulrich and Vasudevan, 2004</xref>).<disp-formula id="e11">
<mml:math id="m54">
<mml:mrow>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mrow>
<mml:mi>C</mml:mi>
<mml:mi>E</mml:mi>
<mml:mi>P</mml:mi>
<mml:mi>C</mml:mi>
<mml:mi>I</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mrow>
<mml:mn>2019</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mn>2004</mml:mn>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
<label>(11)</label>
</disp-formula>
</p>
<p>The annualised capital expenses (CAPEX, <inline-formula id="inf46">
<mml:math id="m57">
<mml:mrow>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mrow>
<mml:mi>C</mml:mi>
<mml:mi>A</mml:mi>
<mml:mi>P</mml:mi>
<mml:mi>E</mml:mi>
<mml:mi>X</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>) are calculated, based on the annual interest rate <inline-formula id="inf47">
<mml:math id="m58">
<mml:mi>i</mml:mi>
</mml:math>
</inline-formula> and the plant lifetime <inline-formula id="inf48">
<mml:math id="m59">
<mml:mi>a</mml:mi>
</mml:math>
</inline-formula>, both shown in <xref ref-type="sec" rid="s9">Supplementary Table S3</xref>:<disp-formula id="e12">
<mml:math id="m60">
<mml:mrow>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mrow>
<mml:mi>C</mml:mi>
<mml:mi>A</mml:mi>
<mml:mi>P</mml:mi>
<mml:mi>E</mml:mi>
<mml:mi>X</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:msubsup>
<mml:mi>C</mml:mi>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>n</mml:mi>
<mml:mi>v</mml:mi>
</mml:mrow>
</mml:msubsup>
<mml:mfrac>
<mml:mrow>
<mml:msup>
<mml:mrow>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>i</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mi>a</mml:mi>
</mml:msup>
<mml:mo>&#x22c5;</mml:mo>
<mml:mi>i</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:msup>
<mml:mrow>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>i</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mi>a</mml:mi>
</mml:msup>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
<label>(12)</label>
</disp-formula>
</p>
<p>The cost for the nutrient solution for biological methanation was estimated based on a total nutrient concentration of 2&#xa0;g/L of slurry (<xref ref-type="bibr" rid="B33">Rachbauer et al., 2016</xref>; <xref ref-type="bibr" rid="B34">Rusmanis et al., 2019</xref>). In case of an artificial nutrient solution fed to the reactor, this growth medium has to be mixed from a variety of chemicals at an estimated cost of around 20&#xa0;&#x20ac;/kg. We therefore estimate the cost for growth media provision to around 40&#xa0;&#x20ac;/m<sup>3</sup>. Alternatively, the reactor can also be supplied with digester efflux, which provides all necessary nutrients (<xref ref-type="bibr" rid="B47">Ullrich and Lemmer, 2019</xref>). Nevertheless, also the integration of such a waste stream and its conditioning is associated with certain costs, which may be in a similar range as shown above.</p>
<p>For the membrane unit, a module geometry of 1&#xa0;m length and 4&#x2033; diameter is considered. This corresponds to a commercial biogas upgrading membrane. A 2&#x2033; version of such a module was tested in previous work for its capability for hydrogen recycle in Power-to-Gas applications (<xref ref-type="bibr" rid="B16">Gantenbein et al., 2021</xref>). The separation area of this module was estimated to be around 65.12&#xa0;m<sup>2</sup>, which corresponds to a fourfold area of the 2&#x201d; module (<xref ref-type="bibr" rid="B50">Witte, 2018</xref>). The module consists of a cartridge fixed into a tubular, pressure-resistant case, which itself provides all the connections to the plant. The cartridge contains the hollow fibres and is sealed with a resin. Due to ageing effects, the cartridges have to be replaced after a certain time interval. The replacement cost were obtained from Ulrich and Vasudevan (<xref ref-type="bibr" rid="B48">Ulrich and Vasudevan, 2004</xref>), which indicate 50&#xa0;$/m<sup>2</sup> in 2004. As the cartridges also maintain pressure and require specialised sealing, we assumed a twice as high cost for the membrane material. Considering the CEPCI, this cost is 126&#xa0;&#x20ac;/m<sup>2</sup>, which is in a similar range as a quote obtained by a supplier. This cost refers to the replacement cost of the membrane cartridge, which does not include further installations. We consider this replacement cost as operational expenses, analogously to sorbent and catalyst replacement cost. For the initial installation cost of the membrane modules, we consider the same amount, but include also the plant cost factors.</p>
<p>For the two methanation reactor types, a detailed cost analysis was performed, based on the reactor vessel itself, as well as the auxiliary units required for operation and heat utilisation. For the fluidised bed system, this includes the cooling cycle, the internal heat transfer tubes and a vessel for catalyst activation as already used by Witte et al. (<xref ref-type="bibr" rid="B52">Witte et al., 2018a</xref>). The biological system also includes an external cooling system, internal cooling coils, which also serve as baffles, as well as the multi-stage stirring system.</p>
</sec>
</sec>
<sec sec-type="results|discussion" id="s3">
<title>3 Results and Discussion</title>
<p>All process options evaluated produce biomethane which is ready for grid injection at a level of 10&#xa0;bar (absolute). The input parameters were taken from the publication of Inkeri et al. (<xref ref-type="bibr" rid="B23">Inkeri et al., 2018</xref>) in the case of biological methanation, or based on previous work of Witte et al. (<xref ref-type="bibr" rid="B53">Witte et al., 2018b</xref>; <xref ref-type="bibr" rid="B52">2018a</xref>) for BFB methanation. Compared to the earlier studies, the input parameters for the BFB cases were slightly modified to result in a process which is comparable to the biological one. These were namely the BFB reactor temperature (360&#xb0;C) and the system pressure (10&#xa0;bar). Nevertheless they are in a comparable range to the optimal conditions reported by Witte et al. A summary of the input parameters for the process simulations can be found in <xref ref-type="table" rid="T3">Table 3</xref>. The values obtained from the process simulation are given in <xref ref-type="table" rid="T4">Table 4</xref>; <xref ref-type="sec" rid="s9">Supplementary Table S4</xref>. Based on these parameters, further cost analysis was performed.</p>
<table-wrap id="T4" position="float">
<label>TABLE 4</label>
<caption>
<p>Technical data for the considered processes. The processes marked with &#x201c;bg&#x201d; are operated using biogas as feed.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="left">Technical parameters</th>
<th align="left"/>
<th align="center">B1</th>
<th align="center">B6</th>
<th align="center">C1</th>
<th align="center">C6</th>
<th align="center">C1S</th>
<th align="center">C6S</th>
<th align="center">C1bg</th>
<th align="center">C6bg</th>
<th align="center">C1Sbg</th>
<th align="center">C6Sbg</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="left">Input Flow</td>
<td align="center">m<sup>3</sup>/h</td>
<td align="center">54.77</td>
<td align="center">328.59</td>
<td align="center">54.47</td>
<td align="center">326.88</td>
<td align="center">54.47</td>
<td align="center">326.91</td>
<td align="center">136.67</td>
<td align="center">820.31</td>
<td align="char" char=".">136.67</td>
<td align="char" char=".">820.31</td>
</tr>
<tr>
<td align="left">Feed Compressor</td>
<td align="center">kW</td>
<td align="center">3.96</td>
<td align="center">23.76</td>
<td align="center">5.11</td>
<td align="center">30.30</td>
<td align="center">5.12</td>
<td align="center">30.38</td>
<td align="center">11.88</td>
<td align="center">70.80</td>
<td align="char" char=".">11.88</td>
<td align="char" char=".">70.80</td>
</tr>
<tr>
<td align="left">Evaporator area</td>
<td align="center">m<sup>2</sup>
</td>
<td align="center">n/a</td>
<td align="center">n/a</td>
<td align="center">n/a</td>
<td align="center">n/a</td>
<td align="center">n/a</td>
<td align="center">n/a</td>
<td align="center">2.3</td>
<td align="center">13.81</td>
<td align="char" char=".">2.3</td>
<td align="char" char=".">13.81</td>
</tr>
<tr>
<td align="left">H<sub>2</sub>/CO<sub>2</sub> (feed)</td>
<td align="center">mol/mol</td>
<td align="center">3.96</td>
<td align="center">3.96</td>
<td align="center">3.98</td>
<td align="center">3.98</td>
<td align="center">3.98</td>
<td align="center">3.98</td>
<td align="center">3.97</td>
<td align="center">3.97</td>
<td align="char" char=".">3.97</td>
<td align="char" char=".">3.97</td>
</tr>
<tr>
<td align="left">Electrolyser Power</td>
<td align="center">kW</td>
<td align="center">1000</td>
<td align="center">6000</td>
<td align="center">1000</td>
<td align="center">6000</td>
<td align="center">737</td>
<td align="center">4423</td>
<td align="center">1000</td>
<td align="center">6000</td>
<td align="char" char=".">737</td>
<td align="char" char=".">4423</td>
</tr>
<tr>
<td align="left">Reactor</td>
<td align="left"/>
<td align="left"/>
<td align="left"/>
<td align="left"/>
<td align="left"/>
<td align="left"/>
<td align="left"/>
<td align="left"/>
<td align="left"/>
<td align="left"/>
<td align="left"/>
</tr>
<tr>
<td align="left">&#xa0;&#xa0;bed diameter</td>
<td align="center">m</td>
<td align="center">1</td>
<td align="center">2.4</td>
<td align="center">0.38</td>
<td align="center">0.83</td>
<td align="center">0.42</td>
<td align="center">0.42</td>
<td align="center">0.48</td>
<td align="center">0.96</td>
<td align="char" char=".">0.48</td>
<td align="char" char=".">0.96</td>
</tr>
<tr>
<td align="left">&#xa0;&#xa0;bed height</td>
<td align="center">m</td>
<td align="center">12</td>
<td align="center">26.5</td>
<td align="center">1.77</td>
<td align="center">2.28</td>
<td align="center">1.50</td>
<td align="center">1.50</td>
<td align="center">1.16</td>
<td align="center">1.74</td>
<td align="char" char=".">1.16</td>
<td align="char" char=".">1.74</td>
</tr>
<tr>
<td align="left">&#xa0;&#xa0;catalyst mass</td>
<td align="center">kg</td>
<td align="center">n/a</td>
<td align="center">n/a</td>
<td align="center">71</td>
<td align="center">427</td>
<td align="center">71</td>
<td align="center">71</td>
<td align="center">72</td>
<td align="center">434</td>
<td align="char" char=".">72</td>
<td align="char" char=".">434</td>
</tr>
<tr>
<td align="left">&#xa0;&#xa0;Internal HEX area</td>
<td align="center">m<sup>2</sup>
</td>
<td align="center">22.9</td>
<td align="center">138.6</td>
<td align="center">11.5</td>
<td align="center">68.9</td>
<td align="center">11.5</td>
<td align="center">11.5</td>
<td align="center">11.7</td>
<td align="center">70.1</td>
<td align="char" char=".">11.7</td>
<td align="char" char=".">70.1</td>
</tr>
<tr>
<td align="left">&#xa0;&#xa0;Total condenser area</td>
<td align="center">m<sup>2</sup>
</td>
<td align="center">n/a</td>
<td align="center">n/a</td>
<td align="center">6.6</td>
<td align="center">39.4</td>
<td align="center">6.8</td>
<td align="center">41.0</td>
<td align="center">9.7</td>
<td align="center">58.1</td>
<td align="char" char=".">9.8</td>
<td align="char" char=".">58.4</td>
</tr>
<tr>
<td align="left">&#xa0;&#xa0;Total HEX area SOE</td>
<td align="center">m<sup>2</sup>
</td>
<td align="center">n/a</td>
<td align="center">n/a</td>
<td align="center">n/a</td>
<td align="center">n/a</td>
<td align="center">17.3</td>
<td align="center">103.9</td>
<td align="center">n/a</td>
<td align="center">n/a</td>
<td align="char" char=".">17.3</td>
<td align="char" char=".">103.9</td>
</tr>
<tr>
<td align="left">&#xa0;&#xa0;El. heater SOE</td>
<td align="center">kW</td>
<td align="center">n/a</td>
<td align="center">n/a</td>
<td align="center">n/a</td>
<td align="center">n/a</td>
<td align="center">6.9</td>
<td align="center">6.9</td>
<td align="center">n/a</td>
<td align="center">n/a</td>
<td align="char" char=".">6.9</td>
<td align="char" char=".">41.1</td>
</tr>
<tr>
<td align="left">Desulphurisation</td>
<td align="left"/>
<td align="left"/>
<td align="left"/>
<td align="left"/>
<td align="left"/>
<td align="left"/>
<td align="left"/>
<td align="left"/>
<td align="left"/>
<td align="left"/>
<td align="left"/>
</tr>
<tr>
<td align="left">&#xa0;&#xa0;Vessel diameter</td>
<td align="center">m</td>
<td align="center">0.48</td>
<td align="center">0.86</td>
<td align="center">0.47</td>
<td align="center">0.86</td>
<td align="center">0.47</td>
<td align="center">0.86</td>
<td align="center">0.64</td>
<td align="center">1.17</td>
<td align="char" char=".">0.64</td>
<td align="char" char=".">1.17</td>
</tr>
<tr>
<td align="left">&#xa0;&#xa0;Vessel height</td>
<td align="center">m</td>
<td align="center">1.43</td>
<td align="center">2.59</td>
<td align="center">1.42</td>
<td align="center">2.59</td>
<td align="center">1.42</td>
<td align="center">2.59</td>
<td align="center">1.93</td>
<td align="center">3.52</td>
<td align="char" char=".">1.93</td>
<td align="char" char=".">3.52</td>
</tr>
<tr>
<td align="left">&#xa0;&#xa0;Sorbent mass (total)</td>
<td align="center">kg</td>
<td align="center">513</td>
<td align="center">2792</td>
<td align="center">510</td>
<td align="center">3059</td>
<td align="center">510</td>
<td align="center">3059</td>
<td align="center">1279</td>
<td align="center">7677</td>
<td align="char" char=".">1279</td>
<td align="char" char=".">7677</td>
</tr>
<tr>
<td align="left">Excess heat</td>
<td align="left"/>
<td align="left"/>
<td align="left"/>
<td align="left"/>
<td align="left"/>
<td align="left"/>
<td align="left"/>
<td align="left"/>
<td align="left"/>
<td align="left"/>
<td align="left"/>
</tr>
<tr>
<td align="left">&#xa0;&#xa0;Reactor</td>
<td align="center">kW</td>
<td align="center">&#x2212;160</td>
<td align="center">&#x2212;970</td>
<td align="center">&#x2212;111</td>
<td align="center">&#x2212;667</td>
<td align="center">&#x2212;5</td>
<td align="center">&#x2212;32</td>
<td align="center">&#x2212;105</td>
<td align="center">&#x2212;627</td>
<td align="char" char=".">0</td>
<td align="char" char=".">0</td>
</tr>
<tr>
<td align="left">&#xa0;&#xa0;Condenser @165&#xb0;C</td>
<td align="center">kW</td>
<td align="center">n/a</td>
<td align="center">n/a</td>
<td align="center">&#x2212;54</td>
<td align="center">&#x2212;325</td>
<td align="center">&#x2212;56</td>
<td align="center">&#x2212;336</td>
<td align="center">&#x2212;80</td>
<td align="center">&#x2212;482</td>
<td align="char" char=".">-80</td>
<td align="char" char=".">-482</td>
</tr>
<tr>
<td align="left">Grid compliance/Product gas concentrations</td>
<td align="left"/>
<td align="left"/>
<td align="left"/>
<td align="left"/>
<td align="left"/>
<td align="left"/>
<td align="left"/>
<td align="left"/>
<td align="left"/>
<td align="left"/>
<td align="left"/>
</tr>
<tr>
<td align="left">&#xa0;&#xa0;x&#x005F;H<sub>2</sub>
</td>
<td align="left"/>
<td align="center">1.4%</td>
<td align="center">1.1%</td>
<td align="center">1.1%</td>
<td align="center">1.3%</td>
<td align="center">1.1%</td>
<td align="center">1.3%</td>
<td align="center">1.4%</td>
<td align="center">1.6%</td>
<td align="char" char=".">1.4%</td>
<td align="char" char=".">1.6%</td>
</tr>
<tr>
<td align="left">&#xa0;&#xa0;x&#x005F;CO<sub>2</sub>
</td>
<td align="left"/>
<td align="center">0.0%</td>
<td align="center">0.0%</td>
<td align="center">0.7%</td>
<td align="center">0.8%</td>
<td align="center">0.7%</td>
<td align="center">0.8%</td>
<td align="center">0.6%</td>
<td align="center">0.7%</td>
<td align="char" char=".">0.6%</td>
<td align="char" char=".">0.7%</td>
</tr>
<tr>
<td align="left">&#xa0;&#xa0;x&#x005F;CH<sub>4</sub>
</td>
<td align="left"/>
<td align="center">98.6%</td>
<td align="center">98.9%</td>
<td align="center">98.2%</td>
<td align="center">97.9%</td>
<td align="center">98.2%</td>
<td align="center">97.8%</td>
<td align="center">97.9%</td>
<td align="center">97.5%</td>
<td align="char" char=".">97.9%</td>
<td align="char" char=".">97.5%</td>
</tr>
<tr>
<td align="left">&#xa0;&#xa0;Total product gas flow</td>
<td align="center">mol/s</td>
<td align="center">0.679</td>
<td align="center">4.066</td>
<td align="center">0.683</td>
<td align="center">4.106</td>
<td align="center">0.683</td>
<td align="center">4.107</td>
<td align="center">1.719</td>
<td align="center">10.343</td>
<td align="char" char=".">1.719</td>
<td align="char" char=".">10.343</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>For each process, investment cost, operation cost and revenues from gas and heat sales are determined. The costs are grouped according to system compartments. The H<sub>2</sub>-path contains the electrolyser and required water feed. In the main reactor section, the main vessel is considered as well as all units and heat exchanger equipment, which are required to remove the reaction heat. In the case of biological methanation, this also includes the stirrer. Membrane upgrading is required to ensure grid injection compliance and is only needed in the case of BFB methanation. The membrane upgrading also removes water, which replaces TSA drying in the BFB cases. As supporting modules, the units required for grid injection are considered, i.e. process analytics and gas odorisation. These units are assumed to be independent of methanation technology and scale. The cost for external heat exchanger equipment considers the equipment required for heat integration between feed gas and product gas (condenser). In the case of high temperature electrolysis, this also includes the heat exchangers required to recover the off-heat of the SOE stack and pre-heat feed steam. For the biogas cases, a water evaporator is considered to add steam to the feed gas. This steam is required to reduce catalyst stress due to the high methane content.</p>
<sec id="s3-1">
<title>3.1 Heat Integration</title>
<p>The high-temperature BFB methanation processes presented here allow two ways of heat utilisation: The pre-heating of the feed gases using the off-heat of the water condenser and in the case of high-temperature electrolysis the steam generation via the use of off-heat from the reactor. The second option was justified via pinch point analysis, as shown in <xref ref-type="sec" rid="s9">Supplementary Figure S1</xref>. The cold streams (blue) correspond to the heat requirement of the feed water to the electrolyser. Water enters the plant at 20&#xb0;C and is heated to 193&#xb0;C, which corresponds to the vaporisation temperature at the operation pressure of 15 bara. The heat requirement in the evaporator is indicated by the horizontal part. The steam is then superheated to the operation temperature of the electrolyser. The last 50&#xa0;K of temperature increase have to be provided by an electrical heater, whose energy demand is indicated by a horizontal shift of the curves. This point also corresponds to the minimal temperature approach of the hot and cold streams, which was set to 50&#xa0;K. The hot streams correspond to the combined heat provision by the product gas streams from the electrolyser (hydrogen and oxygen). Specialised heat exchanger equipment is necessary to valorise the heat of both gases efficiently. Concepts for such equipment are under development (<xref ref-type="bibr" rid="B15">Flaviana et al., 2016</xref>).</p>
<p>The hot streams (red) correspond to the product gas streams cooled to 60&#xb0;C. In the first step, the hot product gases are cooled to 287&#xb0;C, where the freed enthalpy is used to superheat the incoming steam to 750&#xb0;C. The reactor provides heat at a temperature level of 287&#xb0;C, which is sufficient in quality and quantity to provide heat to the evaporator. Subsequently, the product gases are further cooled by the incoming water feed. Part of the steam remains in the hydrogen product gas, as it is not fully converted in the electrolyser. This steam is condensed at a temperature level of 193&#xb0;C, visible as small horizontal part in the hot composite curve (red). The remaining cooling duty can be provided by air cooling, which is not further discussed in this work.</p>
<p>The reactor analysed in this work is cooled by a thermo-oil system and internal heat exchangers. This cooling system is directly connected to the evaporator providing an efficient heat transfer. The remaining heat after the water evaporation can be further utilised at a temperature level of around 280&#xb0;C. Other concepts combining fixed bed methanation with high-temperature electrolysis make use of direct, internal steam generation in the reactor (<xref ref-type="bibr" rid="B18">Giglio et al., 2018</xref>; <xref ref-type="bibr" rid="B10">Dannesboe et al., 2020</xref>). In order to provide an efficient heat transfer, the water has to be kept in a boiling state. Depending on the reactor temperature, this requires high operation pressures, which increases equipment cost substantially and decreases cost efficiency.</p>
<p>As shown in this pinch point analysis, the heat provided by the off-gas streams and the reactor is sufficient to produce and superheat the steam required for the high-temperature electrolysis. Additional electrical heating is necessary to reach isothermal conditions in the stack. Nevertheless, the corresponding electrical energy demand represents only a small share of the total electrolyser power.</p>
</sec>
<sec id="s3-2">
<title>3.2 Investment Cost</title>
<p>As shown in previous work (<xref ref-type="bibr" rid="B52">Witte et al., 2018a</xref>), the investment cost of Power-to-Gas systems are highly dominated by the electrolyser cost. The cases considered in this work are shown in <xref ref-type="fig" rid="F3">Figures 3</xref>, <xref ref-type="fig" rid="F4">4</xref> and summarised in <xref ref-type="sec" rid="s9">Supplementary Table S5</xref>. The hydrogen path contributes to at least 50% of the total cost in case of 1&#xa0;MW biological methanation, and to up to 92% in the case of BFB methanation with 6&#xa0;MW<sub>eq</sub> SOE. This extreme originates mainly in the high assumed cost for the SOE stack in the base-case. Nevertheless, even in the case of a commercial PEM electrolyser, the cost share amounts to 69% in the 1&#xa0;MW and 86% in the 6&#xa0;MW scale.</p>
<fig id="F3" position="float">
<label>FIGURE 3</label>
<caption>
<p>Investment cost structure in the 1 MW<sub>el</sub> cases. The left hand side shows the total investment cost including the H<sub>2</sub>-path, the right side shows the cost structure of the methanation equipment only.</p>
</caption>
<graphic xlink:href="fenrg-09-775259-g003.tif"/>
</fig>
<fig id="F4" position="float">
<label>FIGURE 4</label>
<caption>
<p>Cost structure of the 6 MW<sub>eq</sub> cases, including the H<sub>2</sub>-path <bold>(left)</bold> and focused on the methanation equipment only <bold>(right)</bold>.</p>
</caption>
<graphic xlink:href="fenrg-09-775259-g004.tif"/>
</fig>
<p>When a Power-to-Gas plant is in planning, two main considerations can be identified: An overall cost analysis for a newly installed plant, or a retrofit project, where the decision between methanation technologies has to be taken. In the first case, the hydrogen path has to be included in the considerations. For the second case, as the hydrogen path is identical for all methanation technologies, it can be excluded from the consideration. In some of these cases, existing electrolyser units are retrofitted with PtG equipment. It is therefore important to compare technologies based on total investment cost, as well as only based on methanation equipment cost.</p>
<p>The second-largest share of investment cost is the reactor vessel itself. As shown in previous work (<xref ref-type="bibr" rid="B23">Inkeri et al., 2018</xref>) and pilot plant projects, biological methanation requires a large reactor volume to produce utility-scale biomethane. This is also represented in the data discussed in this work, as for the operation with pressurised hydrogen and potentially corrosive slurries, i.e. CO<sub>2</sub> dissolved in water, a pressure resistant stainless steel tank is required. Although BFB methanation needs a smaller reactor volume, it also requires pressure and corrosion resistance at high temperatures. Nevertheless, mainly due to the size difference, the investment cost for the BFB vessel is less than a third of the cost of the bioreactor.</p>
<p>The data provided by the simulations of Inkeri et al. indicate a methane evolution rate of 168 vvd in the 1&#xa0;MW case and 70 vvd in the 6&#xa0;MW case, based on the reactor sizes given in <xref ref-type="table" rid="T2">Table 2</xref>. Other biological reactor concepts are available and under investigation (<xref ref-type="bibr" rid="B34">Rusmanis et al., 2019</xref>). For fixed film reactors, which include trickle bed reactors, the projects reported by Rusmanis et al. show an average methane evolution rate of around 18 vvd (<xref ref-type="bibr" rid="B34">Rusmanis et al., 2019</xref>). Compared to the total methane formation in CSTR systems, such rates would result in large vessel sizes. Even though these reactors are operated at atmospheric pressure, such volumes will result in high investment cost for the main reactor, exceeding the system considered in this study. The demonstration project ORBIT is a trickle bed reactor operated at pressures up to 12 bar, but also here, a methane formation rate of only 8.4 vvd is reported (<xref ref-type="bibr" rid="B45">Thema et al., 2021</xref>). This indicates that further development towards higher methane formation rates is necessary.</p>
<p>For catalytic BFB methanation, a larger main compressor is required than for biological methanation, since the membrane upgrading unit creates a recycle stream, which increases compression duty.</p>
<p>A further cost factor in the BFB-based processes is the heat exchanger equipment required for utilisation of the product gas stream and subsequent integration with the feed gas stream. This cost is not considered in the bioreactor, since the latter is operated at a low temperature level, which results in a low water content in the product gas. In the case with SOE electrolysis, more exchangers are required for heat management in the electrolyser. In the 1&#xa0;MW-case, this doubles the cost for heat exchanger equipment compared to the PEM-BFB case, for the 6&#xa0;MW-case this results even in a three-fold increase.</p>
<p>In addition to the CO<sub>2</sub> base cases, the process was furthermore simulated using biogas as carbon source. The feed stream was adjusted to match the CO<sub>2</sub> content with the electrolyser power. This results in a 2.5 times higher total feed flow rate than in the pure CO<sub>2</sub>-cases due to the additional methane. This higher total flow rate causes a cost increase of the supporting units, namely the compressor and desulphurisation. The reactor cost on the other hand is almost not affected. As the absolute amount of CO<sub>2</sub> remains identical, the same amount of catalyst is required for the methanation unit. Furthermore, the internal heat exchanger area remains equal, as the same reaction heat has to be dissipated. The gas composition itself has therefore less influence on the reactor dimension and subsequently on the corresponding investment cost than on the supporting units.</p>
<p>Rough estimates for plant scaling are obtained by relating the capacity increase of equipment to the cost increase by an exponential law according to <xref ref-type="disp-formula" rid="e13">Eq. 13</xref> (<xref ref-type="bibr" rid="B46">Tribe and Alpine, 1986</xref>):<disp-formula id="e13">
<mml:math id="m61">
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x3d;</mml:mo>
<mml:msup>
<mml:mrow>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mn>6</mml:mn>
<mml:mi>M</mml:mi>
<mml:mi>W</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mi>M</mml:mi>
<mml:mi>W</mml:mi>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mi>&#x3b1;</mml:mi>
</mml:msup>
</mml:mrow>
</mml:math>
<label>(13)</label>
</disp-formula>
</p>
<p>Such scaling exponents are obtained for all three processes and listed in <xref ref-type="table" rid="T5">Table 5</xref>. For the electrolyser we assumed a linear relationship of cost to the rated power, implying that there is no economy of scale applicable. Again, due to the dominance of the expenses for the H<sub>2</sub>-path, this somewhat conservative assumption has a high influence on the scaling exponent. For the total plant, the exponent is close to unity. When only the methanation plant is considered, the exponent for the biological system is close to 0.6, which is a common rule of thumb value for the scaling of vessels and therefore reflects the dominance of the cost for the reactor vessel. The scaling exponents for the catalytic system are in the range of 0.4&#x2013;0.5, which implies that there is a strong decrease in specific expenses with increase of the scale.</p>
<table-wrap id="T5" position="float">
<label>TABLE 5</label>
<caption>
<p>Scaling exponents for the considered methanation processes.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="left"/>
<th align="left">PEM-bio</th>
<th align="left">PEM-Cat</th>
<th align="left">SOE-Cat</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="left">Full System w/o H<sub>2</sub>-Path</td>
<td align="char" char=".">0.63</td>
<td align="char" char=".">0.43</td>
<td align="char" char=".">0.46</td>
</tr>
<tr>
<td align="left">Full System incl. H<sub>2</sub>-Path</td>
<td align="char" char=".">0.85</td>
<td align="char" char=".">0.87</td>
<td align="char" char=".">0.93</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>The annualised capital cost is calculated according to <xref ref-type="disp-formula" rid="e12">Eq. 12</xref>.</p>
<p>The specific investment cost for the biological methanation without the hydrogen path amount to 2975&#xa0;&#x20ac;/kW in the 1&#xa0;MW<sub>el</sub>-case (B1) and 1538&#xa0;&#x20ac;/kW in the 6&#xa0;MW<sub>el</sub>-case (B6), related to the energy output of the produced biomethane. Schlautmann et al. (<xref ref-type="bibr" rid="B31">M&#xf6;rs et al., 2020</xref>; <xref ref-type="bibr" rid="B37">Schlautmann et al., 2021</xref>) estimate the specfic investment cost for the pilot plant in Solothurn, which corresponds to a biomethane output of 0.325&#xa0;MW, to 4320&#xa0;&#x20ac;/kW. They furthermore provide estimates for specific investment cost for plants scaled to 1&#xa0;MW<sub>biomethane</sub> and 5&#xa0;MW<sub>biomethane</sub>, which amount to 2172&#xa0;&#x20ac;/kW and 940&#xa0;&#x20ac;/kW. In <xref ref-type="sec" rid="s9">Supplementary Figure S2</xref>, the estimates of the Store&#x26;Go project are compared to the results provided in this study. Although the estimates provided in the current study fit well within the estimates reported by Schlautmann et al., they are slightly more conservative. Especially the value for the larger plant size indicates that the cost-related scaling factor provided in this study is lower than the one by Schlautmann et al. It is difficult to assess such differences in detail, as the full cost assessment and detailed methodology of the Store&#x26;Go project is not known to the authors.</p>
</sec>
<sec id="s3-3">
<title>3.3 Operational Expenses (OPEX)</title>
<p>The annual operational expenses for the processes are shown in <xref ref-type="fig" rid="F5">Figures 5</xref>, <xref ref-type="fig" rid="F6">6</xref> and summarised <xref ref-type="sec" rid="s9">Supplementary Table S6</xref>.</p>
<fig id="F5" position="float">
<label>FIGURE 5</label>
<caption>
<p>Operation cost (OPEX) of the 1 MW<sub>el</sub> processes in k&#x20ac; per year. The OPEX including the H<sub>2</sub>-path are shown on the left, the right side shows the cost structure without the H<sub>2</sub>-path and the corresponding O&#x0026;M cost.</p>
</caption>
<graphic xlink:href="fenrg-09-775259-g005.tif"/>
</fig>
<fig id="F6" position="float">
<label>FIGURE 6</label>
<caption>
<p>Operation cost (OPEX) of the 6 MW<sub>el</sub> processes in k&#x20ac; per year. The OPEX including the H<sub>2</sub>-path are shown in on the left, the right side shows the cost structure without the H<sub>2</sub>-path and the corresponding O&#x0026;M cost.</p>
</caption>
<graphic xlink:href="fenrg-09-775259-g006.tif"/>
</fig>
<p>The obvious main driver for operational cost of the Power-to-Gas process is the electricity cost. In both scales, the cases PEM-Bio (B1 and B6) and PEM-BFB (C1 and C6) are comparable in terms of electricity consumption. The SOE introduced in cases C1S and C6S has a higher electrical efficiency, which can be considered as a reduction of electricity use per hydrogen produced.</p>
<p>The annual operation cost also includes labour and maintenance cost (O&#x26;M), which are estimated by 1.5% of the investment cost for the electrolyser and 5% of the total investment cost for the rest of the plant. This differentiation takes into account that electrolysers do not have any moving parts and require therefore less servicing. Due to the dominance of electrolyser costs, the O&#x26;M costs are also separated in two parts.</p>
<p>When excluding the H<sub>2</sub>-path in this analysis, the main contribution to the OPEX is the O&#x26;M cost for the rest of the plant. The actual consumables and utilities used in the process amount only to around one half of the remaining OPEX. For the biological methanation, the largest share of its OPEX are related to the nutrient cost. In the BFB processes, the largest contribution comes from the bi-annual membrane replacement, followed by the catalyst replacement, which occurs on an annual basis. According to our assumptions, the nutrient provision cost are around 4 times higher than the annual catalyst replacement cost in the 1&#xa0;MW catalytic case. For the 6&#xa0;MW case the amount increases tenfold. The amount of nutrients required scales with the volume of the biological system, subsequently there is a disproportionate increase compared to the energy scale.</p>
<p>The stirring power required for the bioreactor amounts to 0.9&#x2013;1.1% of the electrolyser power, as shown in <xref ref-type="table" rid="T2">Table 2</xref>. This corresponds to less than 1% of the total OPEX, but remains in a similar range as the cost caused by other supporting units, such as the main compressor and sorbent replacement costs.</p>
<p>Nevertheless, when including the H<sub>2</sub>-path in the calculation of the pure CO<sub>2</sub> cases, the costs for utilities and consumables are marginal. The cost for the H<sub>2</sub>-path, including the corresponding O&#x26;M expenses make up 76&#x2013;92% of the total.</p>
<p>The charges for membrane replacement have a direct influence on the OPEX, however it represents only a minor part of the total OPEX. Therefore, the sensitivity towards the replacement interval is only low. An annual replacement would double the share of membrane replacement cost compared to the base-case represented in <xref ref-type="fig" rid="F5">Figures 5</xref>, <xref ref-type="fig" rid="F6">6</xref>.</p>
<p>The total OPEX for catalytic methanation operated with biogas are almost twice as high as in the corresponding base cases. The additional cost originates from the price paid for the raw biogas, which can be considered as an opportunity cost for biogas upgrading through an alternative method. Apart from the increase of the total OPEX, the influence on the OPEX related to the methanation equipment is also visible: Due to higher gas throughput, the electricity consumption of the compressor increases, as well as the amount of sorbent required each year. Additionally, the membrane replacement cost approximately doubles. Although roughly the same amount of CO<sub>2</sub> and H<sub>2</sub> has to be separated from the product gas stream a larger amount of membrane modules is required. The increased number of modules compensates for the higher total gas flow rate, as the permeation rate of CO<sub>2</sub> and H<sub>2</sub> remains constant.</p>
</sec>
<sec id="s3-4">
<title>3.4 Revenues and Production Cost</title>
<p>The energy content of the biogas injected to the grid is determined according to <xref ref-type="disp-formula" rid="e14">Eq. 14</xref> from the higher heating values of the species CH<sub>4</sub> and H<sub>2</sub>.<disp-formula id="e14">
<mml:math id="m62">
<mml:mrow>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mrow>
<mml:mi>b</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>m</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>t</mml:mi>
<mml:mi>h</mml:mi>
<mml:mi>a</mml:mi>
<mml:mi>n</mml:mi>
<mml:mi>e</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mover accent="true">
<mml:mi>V</mml:mi>
<mml:mo>&#x2d9;</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mrow>
<mml:mi>b</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>m</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>t</mml:mi>
<mml:mi>h</mml:mi>
<mml:mi>a</mml:mi>
<mml:mi>n</mml:mi>
<mml:mi>e</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>x</mml:mi>
<mml:mrow>
<mml:mi>C</mml:mi>
<mml:msub>
<mml:mi>H</mml:mi>
<mml:mn>4</mml:mn>
</mml:msub>
</mml:mrow>
</mml:msub>
<mml:mi>H</mml:mi>
<mml:mi>H</mml:mi>
<mml:msub>
<mml:mi>V</mml:mi>
<mml:mrow>
<mml:mi>C</mml:mi>
<mml:msub>
<mml:mi>H</mml:mi>
<mml:mn>4</mml:mn>
</mml:msub>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>x</mml:mi>
<mml:mrow>
<mml:msub>
<mml:mi>H</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:mrow>
</mml:msub>
<mml:mi>H</mml:mi>
<mml:mi>H</mml:mi>
<mml:msub>
<mml:mi>V</mml:mi>
<mml:mrow>
<mml:msub>
<mml:mi>H</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(14)</label>
</disp-formula>
</p>
<p>The total revenues from biogas sales are then obtained by multiplication with the biogas selling price. The heat revenues are calculated by considering the high and mid temperature level off-heat. In order to account for heat losses in the heat extraction equipment, we assumed that only 70% of the heat produced can be sold. The specific production cost per kWh of gas injected to the grid are calculated by dividing the total annual cost of the plant (OPEX &#x2b; CAPEX) by the biogas sales. The values obtained for the different processes is summarised in <xref ref-type="sec" rid="s9">Supplementary Table S7</xref>.</p>
<p>The absolute production cost per kWh of biomethane injected to the grid are indicated in <xref ref-type="fig" rid="F7">Figure 7</xref>. The values correspond to the base-case scenario represented by the input parameters in <xref ref-type="sec" rid="s9">Supplementary Table S3</xref>. For this case, neither of the pure CO<sub>2</sub> processes can be operated economically. Using biogas as feed reduces production cost by around 26&#x2013;39%. This causes the two PEM-based catalytic processes to reach a profitable range (cases C1bg and C6bg) at a production cost of 10.2 and 9.5 &#x20ac;-ct./kWh. The process including high-temperature electrolysis is right at the profit margin at a scale of 6&#xa0;MW.</p>
<fig id="F7" position="float">
<label>FIGURE 7</label>
<caption>
<p>Production cost of biomethane for each considered process under the base-case assumptions indicated in <xref ref-type="sec" rid="s9">Supplementary Table S3</xref>. The projected cost for biological methanation operated with biogas is indicated with shaded bars.</p>
</caption>
<graphic xlink:href="fenrg-09-775259-g007.tif"/>
</fig>
<p>As we have no full model for biological methanation available, the production cost for biological methanation (cases B1 and B2) using biogas were projected as shown in <xref ref-type="fig" rid="F7">Figure 7</xref>. The analysis of the investment and operation expenses showed that for the catalytic cases, the change of the feed gas to biogas does hardly affect the cost related to the reactor but increases the expenditures for auxiliary units. We therefore assume that approximately the same cost decrease occurs for the biological system as for the catalytic technology. Even though there is a significant step towards profitability possible, the margin is not yet reached under base-case assumptions.</p>
<p>The 1&#xa0;MW-scale SOE-BFB process has the highest production cost, mostly caused by high investment cost for the SOE stack. The cost decreases with scale, but remains the highest among the technologies. When switching to biogas methanation, the production cost decreases by 39%, which is the highest relative decrease among the technologies. This is a result of the higher electrical efficiency of the SOE stack, which reduces the relative OPEX due to lower expenses for electricity while still producing the same amount of biomethane.</p>
<p>For both feed gases, the 6&#xa0;MW PEM-BFB case indicates the lowest production cost, which is a result of scaling effects, lower-cost electrolyser technology, as well as a small reactor size. Therefore, switching the technology from biological methanation to BFB-based technology can lead to a cost reduction of 17% (1&#xa0;MW) and 19% (6&#xa0;MW) in the case of pure CO<sub>2</sub> methanation. When using biogas, this benefit even increases to 24% (1&#xa0;MW) and 26% (6&#xa0;MW).</p>
<p>A further decrease of production cost is achieved by valorizing the process heat. For the PEM-BFB cases a cost reduction of around 5%, and for the SOE-BFB cases a reduction of around 2% is possible. The latter cases do have a limited heat output due to the heat requirement of the SOE. Furthermore, when biogas is used only very limited heat sales are possible, as part of it is redirected to feed water evaporation. The valorization of the oxygen produced in the electrolysis would be another opportunity to decrease cost. However, this is challenging to realize at most sites and is therefore not considered in this work.</p>
</sec>
<sec id="s3-5">
<title>3.5 Sensitivity on Electricity Price and SOE Investment Cost</title>
<p>The electricity cost states the largest influence on the operation cost of the processes. The base-case used above is based on an electricity price of 5&#xa0;&#x20ac;-cts./kWh. This value is comparable to the production cost of wind and large-scale photovoltaic electricity (<xref ref-type="bibr" rid="B3">Bauer et al., 2017</xref>). <xref ref-type="fig" rid="F8">Figure 8</xref> indicates the production cost for each process as a function of the electricity price. The profit margin corresponds to the selling price of renewable gas indicated at 0.11 &#x20ac;/kWh<sub>biomethane</sub>. As shown in <xref ref-type="fig" rid="F7">Figure 7</xref> for the pure CO<sub>2</sub> cases, the 6&#xa0;MW PEM-BFB case shows the lowest production cost and therefore reaches profitability at the highest electricity cost of approximately 3.2&#xa0;&#x20ac;-cts./kWh<sub>el</sub>. The corresponding 1&#xa0;MW-case with the second-lowest production cost reaches profitability at 2.2&#xa0;&#x20ac;-cts./kWh<sub>el</sub>, whereas the larger-scale biological methanation process requires an electricity price as low as 1.2&#xa0;&#x20ac;-cts./kWh<sub>el</sub> to reach the profit margin. The 1&#xa0;MW-scale biological process only reaches a profitable range when electricity is provided at 0.3&#xa0;&#x20ac;-cts./kWh<sub>el</sub>. The performance of the SOE-based processes is inferior, as no profitability is reached, even if no electricity cost is considered. Nevertheless, it is notable, that the SOE-based processes are less sensitive towards electricity cost, i.e. the production cost do not rise as fast with the electricity price as in the PEM-based cases. This is a direct consequence of the higher efficiency of the SOE stack and the heat integration with the BFB reactor.</p>
<fig id="F8" position="float">
<label>FIGURE 8</label>
<caption>
<p>Production cost in cents per kWh of biomethane as a function of the electricity price for each considered process. The profit margin is indicated as line at 11 &#x20ac;-ct/kWh<sub>biomethane</sub>, which indicates the selling price for renewable gas.</p>
</caption>
<graphic xlink:href="fenrg-09-775259-g008.tif"/>
</fig>
<p>Biogas as feedstock increases the total amount of methane produced per kWh of electricity. Generally, the processes are therefore less sensitive towards the electricity price than in the case of pure CO<sub>2</sub> methanation. This can be seen in <xref ref-type="fig" rid="F8">Figure 8</xref> by a decreased slope of the corresponding lines. The 6&#xa0;MW-scale PEM-BFB process reaches profitability at 7.4&#xa0;&#x20ac;-cts./kWh<sub>el</sub> and at 6.2&#xa0;&#x20ac;-cts./kWh<sub>el</sub> in the 1&#xa0;MW case. Lower, but still above the pure CO<sub>2</sub> processes, the methanation with high temperature electrolysis reaches the margin: At 5.4&#xa0;&#x20ac;-cts./kWh<sub>el</sub> for the 6&#xa0;MW case and 3.5&#xa0;&#x20ac;-cts./kWh<sub>el</sub> for the 1&#xa0;MW scale, respectively.</p>
<p>In order to make the SOE-based processes profitable, the investment cost and subsequent CAPEX for the SOE stack have to decrease considerably. Schmidt et al. (<xref ref-type="bibr" rid="B39">Schmidt et al., 2017</xref>) estimate the cost for the SOE technology to be in a range of 3000&#x2013;1900&#xa0;&#x20ac;/kW<sub>el</sub> in the year 2030, depending on R&#x26;D funding and production scale. This is still around three times as high as their estimate for PEM electrolysers. <xref ref-type="fig" rid="F9">Figure 9</xref> shows the influence of the investment cost of the SOE stack on the production cost of biomethane. Electricity prices of 5&#xa0;&#x20ac;-cts./kWh<sub>el</sub> or 2.5&#xa0;&#x20ac;-cts./kWh<sub>el</sub> are taken into account.</p>
<fig id="F9" position="float">
<label>FIGURE 9</label>
<caption>
<p>Production cost in &#x20ac;-ct/kWh<sub>biomethane</sub> of the SOE-based processes as a function of investment cost for the SOE stack. Two different electricity prices are considered: 5 &#x20ac;-ct/kWh<sub>el</sub> and 2.5 &#x20ac;-ct/kWh<sub>el</sub>.</p>
</caption>
<graphic xlink:href="fenrg-09-775259-g009.tif"/>
</fig>
<p>Depending on the feed gas, two further cost ranges can be identified: The case for pure CO<sub>2</sub> methanation and the case for direct methanation of biogas. As in the latter a larger amount of methane is fed to the grid based on the same nominal power of the electrolyser, biogas methanation has a higher profitability range and allows for higher investment cost for the electrolyser. For the same reason, the production cost for direct methanation is less sensitive towards the electrolyser investment cost and the corresponding lines appear with a reduced slope in <xref ref-type="fig" rid="F9">Figure 9</xref>.</p>
<p>As shown in <xref ref-type="fig" rid="F7">Figure 7</xref>, the two SOE-based processes for CO<sub>2</sub> methanation exhibit production costs of 19.4&#xa0;&#x20ac;-cts./kWh<sub>el</sub> for the 1&#xa0;MW-scale and 17.4&#xa0;&#x20ac;-cts./kWh<sub>el</sub> for the 6&#xa0;MW-scale respectively. To reach a profitable range, the production cost has to be reduced to below 11 &#x20ac;-cts./kWh<sub>el</sub>. As shown in <xref ref-type="fig" rid="F9">Figure 9</xref>, for the higher electricity price this will be at investment costs lower than 500&#xa0;&#x20ac;/kW<sub>el</sub> for the 1&#xa0;MW case and lower than 960&#xa0;&#x20ac;/kW<sub>el</sub> for the 6&#xa0;MW-scale. This is a very ambitious cost decrease, nevertheless, under favourable circumstances, a capital cost lower than 1000&#xa0;&#x20ac;/kW<sub>el</sub> in future is considered possible (<xref ref-type="bibr" rid="B39">Schmidt et al., 2017</xref>; <xref ref-type="bibr" rid="B1">Adam Christensen, 2020</xref>). Literature indicates that capital cost for SOEs may adjust to a similar range as PEM or AEC technology (<xref ref-type="bibr" rid="B39">Schmidt et al., 2017</xref>; <xref ref-type="bibr" rid="B6">Buttler and Spliethoff, 2018</xref>). A halving of the electricity cost to 2.5&#xa0;&#x20ac;-cts./kWh<sub>el</sub> results in an increase of the allowed investment cost in the 1 and 6&#xa0;MW case of approximately 2.6-fold to 1260&#xa0;&#x20ac;/kW and 1.9-fold to 1840&#xa0;&#x20ac;/kW, respectively.</p>
<p>Due to the reduced production expenses for direct methanation, the maximum allowed electrolyser investment increases to 2260 &#x20ac;/kW<sub>el</sub> (1 MW<sub>el</sub>) and 2930&#xa0;&#x20ac;/kWh<sub>el</sub> (6&#xa0;MW<sub>el</sub>) for the two scales based on an electricity price of 5&#xa0;&#x20ac;-cts./kWh<sub>el</sub>. Using the lower electricity cost, this numbers increase to 3150&#xa0;&#x20ac;/kW<sub>el</sub> and 3820&#xa0;&#x20ac;/kW<sub>el</sub>. This indicates, that even with moderate improvements in the investment cost of the SOE, profitability of the process can be achieved.</p>
<p>Again, the profitability range can be further extended by also including heat utilisation and oxygen sales. For the heat utilisation, a production cost decrease of around 2% is possible in the pure CO<sub>2</sub> case.</p>
</sec>
</sec>
<sec id="s4">
<title>4 Conclusion</title>
<p>A techno-economic analysis of three CO<sub>2</sub> methanation process chains in two different scales (1&#xa0;MW<sub>el</sub> and 6&#xa0;MW<sub>el</sub>) was performed. The focus was set on the comparison of biological CO<sub>2</sub> methanation in a stirred bubble column with catalytic bubbling fluidised bed methanation. Furthermore, the benefits of including high-temperature electrolysis compared to industrially mature PEM electrolysis was investigated. For all investigated cases, fulfilling the grid injection requirements (mainly a H<sub>2</sub>-content of &#x3c;2%) was required. For the investigated catalytic processes, we also included the case of direct methanation of biogas, which showed a clear cost benefit over CO<sub>2</sub> methanation.</p>
<p>In general, the investment and operation cost of the process chains are highly dominated by the hydrogen path. On one hand, there is the high specific investment cost for the electrolyser unit, which is more than twice as high in the case of high-temperature electrolysis. On the other hand, electricity cost directly influences the operation cost of the process. Nevertheless, these costs are specific to the electrolyser technology applied and are independent from the methanation technology. Therefore, for cost comparison in between the technologies, a focussed view on the respective cost without the H<sub>2</sub>-path may be advisable.</p>
<p>The investment cost for biological methanation is highly dominated by the large reactor vessel required. The investment cost benefit for the BFB technology mainly originates from the significantly smaller reactor vessel required for catalytic methanation. The investment cost of biological methanation can only be reduced when smaller reactor volumes are applicable. The required reactor volume on the other hand is directly linked to the gas-to-liquid mass transfer. As long as this coefficient cannot be increased, the methane formation per volume is limited. A possibility of decreasing the vessel cost would be the reduction of the operation pressure. Nevertheless, this results in a reduced mass transfer rate and therefore again to a larger required reactor volume.</p>
<p>The nutrient provision for biological methanation states the largest influence on the OPEX of that process option. This cost causes the OPEX to exceed the annual cost of the PEM-BFB process. Concepts exist where nutrients are provided by feeding digester efflux to the reactor, which may cause a cost decrease. Still, even when no cost is associated with nutrient provision the production cost for the biological system remains inferior to the catalytic system.</p>
<p>The integration of high-temperature electrolysis leads to a higher electrical efficiency of the Power-to-Gas process chain. Nevertheless, the reduction in electricity cost does not compensate for added capital expensed due to the high stack investment cost and additional heat transfer equipment needed. Currently, SOE cells are still under investigation, which results in high investment cost in near future. Based on an electricity price of 5&#xa0;ct./kWh<sub>el</sub> and a biomethane selling price of 0.11&#xa0;&#x20ac;/kWh, the investment cost for SOE technology have to decrease to 950&#xa0;&#x20ac;/kW<sub>el</sub> for 6&#xa0;MW plant size to reach profitability. However, such a cost decrease is not unfeasible if the technology is further developed.</p>
<p>Direct methanation of biogas reduced the production cost in all investigated cases by 26&#x2013;39%, which resulted in biomethane production cost of less than 11&#xa0;ct/kWh for the 1&#xa0;MW<sub>el</sub> and 6&#xa0;MW<sub>el</sub> PEM-BFB cases, as well as the 6&#xa0;MW<sub>el</sub> SOE-BFB case. Even though direct methanation showed also great cost reduction potential in biological methanation, profitability could not be reached in the 1&#xa0;MW<sub>el</sub> case. In the 6&#xa0;MW<sub>el</sub>-scale, our calculations indicate a production cost slightly below the profitability limit, in the same range as the 6&#xa0;MW<sub>el</sub> SOE-BFB case.</p>
<p>Methanation of CO<sub>2</sub> and biogas are technologies on a high TRL (catalytic 6-7, biological 7-8) and allow for fast integration of biogenic gases into the future energy landscape. This study shows that the underlying technologies can play a cost-efficient role in the seasonal storage of renewable electricity.</p>
</sec>
</body>
<back>
<sec id="s5">
<title>Data Availability Statement</title>
<p>The original contribution presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.</p>
</sec>
<sec id="s6">
<title>Author Contributions</title>
<p>AG: Investigation, Methodology, Calculations and Software, Writing&#x2014;original draft. TS: Conceptualization, Writing&#x2014;review and editing, Funding acquisition, Supervision, Validation, Methodology. SB: Funding acquisition, Project administration, Resources. OK: Funding acquisition, Supervision, Writing&#x2014;editing.</p>
</sec>
<sec sec-type="COI-statement" id="s7">
<title>Conflict of Interest</title>
<p>The 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="s8">
<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>
<ack>
<p>This research project is financially supported by the Swiss Federal Office of Energy within the project &#x201c;CarbonATE&#x201d;. Wolfgang Merkle (ZHAW) is sincerely thanked for their contribution to discussion on nutrients.</p>
</ack>
<sec id="s9">
<title>Supplementary Material</title>
<p>The Supplementary Material for this article can be found online at: <ext-link ext-link-type="uri" xlink:href="https://www.frontiersin.org/articles/10.3389/fenrg.2021.775259/full#supplementary-material">https://www.frontiersin.org/articles/10.3389/fenrg.2021.775259/full&#x23;supplementary-material</ext-link>
</p>
<supplementary-material xlink:href="DataSheet2.pdf" id="SM1" mimetype="application/pdf" xmlns:xlink="http://www.w3.org/1999/xlink"/>
<supplementary-material xlink:href="DataSheet3.pdf" id="SM2" mimetype="application/pdf" xmlns:xlink="http://www.w3.org/1999/xlink"/>
<supplementary-material xlink:href="DataSheet1.pdf" id="SM3" mimetype="application/pdf" xmlns:xlink="http://www.w3.org/1999/xlink"/>
</sec>
<ref-list>
<title>References</title>
<ref id="B1">
<citation citation-type="web">
<collab>Adam Christensen</collab> (<year>2020</year>). <article-title>Assessment of Hydrogen Production Costs from Electrolysis: United States and Europe</article-title>. <comment>Available at: <ext-link ext-link-type="uri" xlink:href="https://theicct.org/sites/default/files/publications/final_icct2020_assessment_of_hydrogen_production_costsv2.pdf">https://theicct.org/sites/default/files/publications/final_icct2020_assessment_of_hydrogen_production_costsv2.pdf</ext-link>
</comment>. </citation>
</ref>
<ref id="B2">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Baker</surname>
<given-names>R. W.</given-names>
</name>
</person-group> (<year>2002</year>). <article-title>Future Directions of Membrane Gas Separation Technology</article-title>. <source>Ind. Eng. Chem. Res.</source> <volume>41</volume>, <fpage>1393</fpage>&#x2013;<lpage>1411</lpage>. <pub-id pub-id-type="doi">10.1021/ie0108088</pub-id> </citation>
</ref>
<ref id="B3">
<citation citation-type="web">
<person-group person-group-type="author">
<name>
<surname>Bauer</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Hirschberg</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>B&#xe4;uerle</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Biollaz</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Calbry-Muzyka</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Cox</surname>
<given-names>B.</given-names>
</name>
<etal/>
</person-group> (<year>2017</year>). <article-title>Potentials, Costs and Environmental Assessment of Electricity Generation Technologies. Villigen PSI</article-title>. <comment>Available at: <ext-link ext-link-type="uri" xlink:href="https://www.psi.ch/lea/HomeEN/Final-Report-BFE-Project.pdf">https://www.psi.ch/lea/HomeEN/Final-Report-BFE-Project.pdf</ext-link>
</comment>. </citation>
</ref>
<ref id="B4">
<citation citation-type="web">
<person-group person-group-type="author">
<name>
<surname>Bertuccioli</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Chan</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Hart</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Lehner</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Madden</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Standen</surname>
<given-names>E.</given-names>
</name>
</person-group> (<year>2014</year>). <article-title>Study on Development of Water Electrolysis in the EU</article-title>. <comment>Available at: <ext-link ext-link-type="uri" xlink:href="https://www.fch.europa.eu/sites/default/files/FCHJUElectrolysisStudy_FullReport_(ID_199214).pdf">https://www.fch.europa.eu/sites/default/files/FCHJUElectrolysisStudy_FullReport_(ID_199214).pdf</ext-link>
</comment>. </citation>
</ref>
<ref id="B5">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>B&#xf6;hm</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Zauner</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Rosenfeld</surname>
<given-names>D. C.</given-names>
</name>
<name>
<surname>Tichler</surname>
<given-names>R.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Projecting Cost Development for Future Large-Scale Power-To-Gas Implementations by Scaling Effects</article-title>. <source>Appl. Energ.</source> <volume>264</volume>, <fpage>114780</fpage>. <pub-id pub-id-type="doi">10.1016/j.apenergy.2020.114780</pub-id> </citation>
</ref>
<ref id="B6">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Buttler</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Spliethoff</surname>
<given-names>H.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>Current Status of Water Electrolysis for Energy Storage, Grid Balancing and Sector Coupling via Power-To-Gas and Power-To-Liquids: A Review</article-title>. <source>Renew. Sustain. Energ. Rev.</source> <volume>82</volume>, <fpage>2440</fpage>&#x2013;<lpage>2454</lpage>. <pub-id pub-id-type="doi">10.1016/j.rser.2017.09.003</pub-id> </citation>
</ref>
<ref id="B7">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Calbry-Muzyka</surname>
<given-names>A. S.</given-names>
</name>
<name>
<surname>Gantenbein</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Schneebeli</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Frei</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Knorpp</surname>
<given-names>A. J.</given-names>
</name>
<name>
<surname>Schildhauer</surname>
<given-names>T. J.</given-names>
</name>
<etal/>
</person-group> (<year>2019</year>). <article-title>Deep Removal of Sulfur and Trace Organic Compounds from Biogas to Protect a Catalytic Methanation Reactor</article-title>. <source>Chem. Eng. J.</source> <volume>360</volume>, <fpage>577</fpage>&#x2013;<lpage>590</lpage>. <pub-id pub-id-type="doi">10.1016/j.cej.2018.12.012</pub-id> </citation>
</ref>
<ref id="B8">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Calbry-Muzyka</surname>
<given-names>A. S.</given-names>
</name>
<name>
<surname>Schildhauer</surname>
<given-names>T. J.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Direct Methanation of Biogas-Technical Challenges and Recent Progress</article-title>. <source>Front. Energ. Res.</source> <volume>8</volume>. <pub-id pub-id-type="doi">10.3389/fenrg.2020.570887</pub-id> </citation>
</ref>
<ref id="B9">
<citation citation-type="journal">
<collab>Chemical Engineering Plant Cost Index (CEPCI)</collab> (<year>2020</year>). <article-title>Process Integration, Simulation and Economics</article-title>. <source>Chem. Eng.</source>. <comment>Available at: <ext-link ext-link-type="uri" xlink:href="https://www.chemengonline.com/pci-home">https://www.chemengonline.com/pci-home</ext-link>
</comment> <volume>127</volume>, <fpage>72</fpage>. </citation>
</ref>
<ref id="B10">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Dannesboe</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Hansen</surname>
<given-names>J. B.</given-names>
</name>
<name>
<surname>Johannsen</surname>
<given-names>I.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Catalytic Methanation of CO<sub>2</sub> in Biogas: Experimental Results from a Reactor at Full Scale</article-title>. <source>React. Chem. Eng.</source> <volume>5</volume>, <fpage>183</fpage>&#x2013;<lpage>189</lpage>. <pub-id pub-id-type="doi">10.1039/C9RE00351G</pub-id> </citation>
</ref>
<ref id="B11">
<citation citation-type="book">
<collab>DIN 1343:1990&#x2013;01</collab> (<year>1990</year>). <source>Reference Conditions, Normal Conditions, Normal Volume; Concepts and Values</source>. <pub-id pub-id-type="doi">10.31030/2333155</pub-id> </citation>
</ref>
<ref id="B12">
<citation citation-type="web">
<collab>Electrochaea.dk ApS</collab> (<year>2014</year>). <article-title>BioCatProject</article-title>. <comment>Available at: <ext-link ext-link-type="uri" xlink:href="http://biocat-project.com/">http://biocat-project.com/</ext-link>
</comment> (<comment>Accessed February 11, 2020</comment>). </citation>
</ref>
<ref id="B13">
<citation citation-type="web">
<collab>Energie360&#xb0;</collab> (<year>2021</year>). <article-title>Gaspreise</article-title>. <comment>Available at: <ext-link ext-link-type="uri" xlink:href="https://www.energie360.ch/de/energie-360/wissen/erdgas-biogas/gaspreise/">https://www.energie360.ch/de/energie-360/wissen/erdgas-biogas/gaspreise/</ext-link>
</comment>. </citation>
</ref>
<ref id="B14">
<citation citation-type="web">
<collab>Evonik Industries AG SEPURAN&#xae;</collab> (<year>2021</year>). <article-title>Efficient Gas Separation Using Hollow-Fiber Membranes</article-title>. <comment>Available at: <ext-link ext-link-type="uri" xlink:href="https://www.membrane-separation.com/en/efficient-gas-separation-with-sepuran-by-evonik">https://www.membrane-separation.com/en/efficient-gas-separation-with-sepuran-by-evonik</ext-link>
</comment>. </citation>
</ref>
<ref id="B15">
<citation citation-type="confproc">
<person-group person-group-type="author">
<name>
<surname>Flaviana</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Monaco</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Michele</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Marchese</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Ambrosio</surname>
<given-names>P. E.</given-names>
</name>
<name>
<surname>Lorusso</surname>
<given-names>M.</given-names>
</name>
<etal/>
</person-group> (<year>2016</year>). &#x201c;<article-title>Selective Laser Melting for Heat Exchanger</article-title>,&#x201d; in <conf-name>Proceeding 6th Int. Conf. Addit. Technol.</conf-name>, <fpage>383</fpage>&#x2013;<lpage>388</lpage>. </citation>
</ref>
<ref id="B16">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Gantenbein</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Witte</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Biollaz</surname>
<given-names>S. M. A.</given-names>
</name>
<name>
<surname>Kr&#xf6;cher</surname>
<given-names>O.</given-names>
</name>
<name>
<surname>Schildhauer</surname>
<given-names>T. J.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Flexible Application of Biogas Upgrading Membranes for Hydrogen Recycle in Power-To-Methane Processes</article-title>. <source>Chem. Eng. Sci.</source> <volume>229</volume>, <fpage>116012</fpage>. <pub-id pub-id-type="doi">10.1016/j.ces.2020.116012</pub-id> </citation>
</ref>
<ref id="B17">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Garrett</surname>
<given-names>D. E.</given-names>
</name>
</person-group> (<year>1989</year>). <source>Chemical Engineering Economics</source>. <publisher-loc>Dordrecht</publisher-loc>: <publisher-name>Springer Netherlands</publisher-name>. <pub-id pub-id-type="doi">10.1007/978-94-011-6544-0</pub-id> </citation>
</ref>
<ref id="B18">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Giglio</surname>
<given-names>E.</given-names>
</name>
<name>
<surname>Deorsola</surname>
<given-names>F. A.</given-names>
</name>
<name>
<surname>Gruber</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Harth</surname>
<given-names>S. R.</given-names>
</name>
<name>
<surname>Morosanu</surname>
<given-names>E. A.</given-names>
</name>
<name>
<surname>Trimis</surname>
<given-names>D.</given-names>
</name>
<etal/>
</person-group> (<year>2018</year>). <article-title>Power-to-Gas through High Temperature Electrolysis and Carbon Dioxide Methanation: Reactor Design and Process Modeling</article-title>. <source>Ind. Eng. Chem. Res.</source> <volume>57</volume>, <fpage>4007</fpage>&#x2013;<lpage>4018</lpage>. <pub-id pub-id-type="doi">10.1021/acs.iecr.8b00477</pub-id> </citation>
</ref>
<ref id="B19">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Giglio</surname>
<given-names>E.</given-names>
</name>
<name>
<surname>Lanzini</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Santarelli</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Leone</surname>
<given-names>P.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>Synthetic Natural Gas via Integrated High-Temperature Electrolysis and Methanation: Part I-Energy Performance</article-title>. <source>J. Energ. Storage</source> <volume>1</volume>, <fpage>22</fpage>&#x2013;<lpage>37</lpage>. <pub-id pub-id-type="doi">10.1016/j.est.2015.04.002</pub-id> </citation>
</ref>
<ref id="B20">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>G&#xf6;tz</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Lefebvre</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>M&#xf6;rs</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>McDaniel Koch</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Graf</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Bajohr</surname>
<given-names>S.</given-names>
</name>
<etal/>
</person-group> (<year>2016</year>). <article-title>Renewable Power-To-Gas: A Technological and Economic Review</article-title>. <source>Renew. Energ.</source> <volume>85</volume>, <fpage>1371</fpage>&#x2013;<lpage>1390</lpage>. <pub-id pub-id-type="doi">10.1016/j.renene.2015.07.066</pub-id> </citation>
</ref>
<ref id="B21">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Gruber</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Weinbrecht</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>Biffar</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Harth</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Trimis</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Brabandt</surname>
<given-names>J.</given-names>
</name>
<etal/>
</person-group> (<year>2018</year>). <article-title>Power-to-Gas through thermal Integration of High-Temperature Steam Electrolysis and Carbon Dioxide Methanation - Experimental Results</article-title>. <source>Fuel Process. Techn.</source> <volume>181</volume>, <fpage>61</fpage>&#x2013;<lpage>74</lpage>. <pub-id pub-id-type="doi">10.1016/j.fuproc.2018.09.003</pub-id> </citation>
</ref>
<ref id="B22">
<citation citation-type="web">
<person-group person-group-type="author">
<name>
<surname>Huettenrauch</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Hoffmann</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>F&#xf6;cker</surname>
<given-names>H.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Innovative Large-Scale Energy Storage Tech-Nologies and Power-To-Gas Concepts after Optimisation</article-title>. <comment>Report on Economic Analysis of Test Cases (D5.10)</comment>. <comment>Available at: <ext-link ext-link-type="uri" xlink:href="https://www.storeandgo.info/publications/deliverables/">https://www.storeandgo.info/publications/deliverables/</ext-link>
</comment>. </citation>
</ref>
<ref id="B23">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Inkeri</surname>
<given-names>E.</given-names>
</name>
<name>
<surname>Tynj&#xe4;l&#xe4;</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Laari</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Hypp&#xe4;nen</surname>
<given-names>T.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>Dynamic One-Dimensional Model for Biological Methanation in a Stirred Tank Reactor</article-title>. <source>Appl. Energ.</source> <volume>209</volume>, <fpage>95</fpage>&#x2013;<lpage>107</lpage>. <pub-id pub-id-type="doi">10.1016/j.apenergy.2017.10.073</pub-id> </citation>
</ref>
<ref id="B24">
<citation citation-type="web">
<collab>Helmeth project</collab> (<year>2021</year>). <article-title>Karlsruher Institut f&#xfc;r Technologie (KIT) Helmeth Project</article-title>. <comment>Available at: <ext-link ext-link-type="uri" xlink:href="http://helmeth.eu/">http://helmeth.eu/</ext-link>
</comment>. </citation>
</ref>
<ref id="B25">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Kober</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Bauer</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Bach</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Beuse</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Georges</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Held</surname>
<given-names>M.</given-names>
</name>
<etal/>
</person-group> (<year>2019</year>). <source>Perspectives of Power-To-X Technologies in Switzerland - A White Paper</source>. <pub-id pub-id-type="doi">10.3929/ethz-b-000352294</pub-id> </citation>
</ref>
<ref id="B26">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kopyscinski</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Schildhauer</surname>
<given-names>T. J.</given-names>
</name>
<name>
<surname>Biollaz</surname>
<given-names>S. M. A.</given-names>
</name>
</person-group> (<year>2011</year>). <article-title>Methanation in a Fluidized Bed Reactor with High Initial CO Partial Pressure: Part II- Modeling and Sensitivity Study</article-title>. <source>Chem. Eng. Sci.</source> <volume>66</volume>, <fpage>1612</fpage>&#x2013;<lpage>1621</lpage>. <pub-id pub-id-type="doi">10.1016/j.ces.2010.12.029</pub-id> </citation>
</ref>
<ref id="B27">
<citation citation-type="web">
<collab>Limeco</collab> (<year>2020</year>). <article-title>Project Web Page: Power-To-Gas Plant Dietikon</article-title>. <comment>Available at: <ext-link ext-link-type="uri" xlink:href="https://www.powertogas.ch/">https://www.powertogas.ch/</ext-link>
</comment>. </citation>
</ref>
<ref id="B28">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Makaruk</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Harasek</surname>
<given-names>M.</given-names>
</name>
</person-group> (<year>2009</year>). <article-title>Numerical Algorithm for Modelling Multicomponent Multipermeator Systems</article-title>. <source>J. Membr. Sci.</source> <volume>344</volume>, <fpage>258</fpage>&#x2013;<lpage>265</lpage>. <pub-id pub-id-type="doi">10.1016/j.memsci.2009.08.013</pub-id> </citation>
</ref>
<ref id="B29">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Maurer</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Gschwend</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Wagner</surname>
<given-names>E. C.</given-names>
</name>
<name>
<surname>Schildhauer</surname>
<given-names>T. J.</given-names>
</name>
<name>
<surname>Ruud van Ommen</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Biollaz</surname>
<given-names>S. M. A.</given-names>
</name>
<etal/>
</person-group> (<year>2016</year>). <article-title>Correlating Bubble Size and Velocity Distribution in Bubbling Fluidized Bed Based on X-ray Tomography</article-title>. <source>Chem. Eng. J.</source> <volume>298</volume>, <fpage>17</fpage>&#x2013;<lpage>25</lpage>. <pub-id pub-id-type="doi">10.1016/j.cej.2016.02.012</pub-id> </citation>
</ref>
<ref id="B30">
<citation citation-type="web">
<collab>microbEnergy GmbH</collab> (<year>2019</year>). <article-title>microbEnergy</article-title>. <comment>Available at: <ext-link ext-link-type="uri" xlink:href="https://www.microbenergy.de/">https://www.microbenergy.de/</ext-link>
</comment> (<comment>Accessed February 11 2020</comment>). </citation>
</ref>
<ref id="B31">
<citation citation-type="web">
<person-group person-group-type="author">
<name>
<surname>M&#xf6;rs</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Schlautmann</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Gorre</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Leonhard</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Graf</surname>
<given-names>F.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Innovative Large-Scale Energy Storage Technologies and Power-To-Gas Concepts after Optimisation</article-title>. <comment>Final Report on Evaluation of Technologies and Processes (D5.9)</comment>. <comment>Available at: <ext-link ext-link-type="uri" xlink:href="https://www.storeandgo.info/publications/deliverables/">https://www.storeandgo.info/publications/deliverables/</ext-link>
</comment>. </citation>
</ref>
<ref id="B32">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Posdziech</surname>
<given-names>O.</given-names>
</name>
<name>
<surname>Schwarze</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Brabandt</surname>
<given-names>J.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Efficient Hydrogen Production for Industry and Electricity Storage via High-Temperature Electrolysis</article-title>. <source>Int. J. Hydrogen Energ.</source> <volume>44</volume>, <fpage>19089</fpage>&#x2013;<lpage>19101</lpage>. <pub-id pub-id-type="doi">10.1016/j.ijhydene.2018.05.169</pub-id> </citation>
</ref>
<ref id="B33">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Rachbauer</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Voitl</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Bochmann</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Fuchs</surname>
<given-names>W.</given-names>
</name>
</person-group> (<year>2016</year>). <article-title>Biological Biogas Upgrading Capacity of a Hydrogenotrophic Community in a Trickle-Bed Reactor</article-title>. <source>Appl. Energ.</source> <volume>180</volume>, <fpage>483</fpage>&#x2013;<lpage>490</lpage>. <pub-id pub-id-type="doi">10.1016/j.apenergy.2016.07.109</pub-id> </citation>
</ref>
<ref id="B34">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Rusmanis</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>O&#x2019;Shea</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Wall</surname>
<given-names>D. M.</given-names>
</name>
<name>
<surname>Murphy</surname>
<given-names>J. D.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Biological Hydrogen Methanation Systems - an Overview of Design and Efficiency</article-title>. <source>Bioengineered</source> <volume>10</volume>, <fpage>604</fpage>&#x2013;<lpage>634</lpage>. <pub-id pub-id-type="doi">10.1080/21655979.2019.1684607</pub-id> </citation>
</ref>
<ref id="B35">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Schill</surname>
<given-names>N.</given-names>
</name>
<name>
<surname>van Gulik</surname>
<given-names>W. M.</given-names>
</name>
<name>
<surname>von Stockar</surname>
<given-names>U.</given-names>
</name>
<name>
<surname>von Stockar</surname>
<given-names>U.</given-names>
</name>
</person-group> (<year>1996</year>). <article-title>Continuous Cultures Limited by a Gaseous Substrate: Development of a Simple, Unstructured Mathematical Model and Experimental Verification with Methanobacterium Thermoautotrophicum</article-title>. <source>Biotechnol. Bioeng.</source> <volume>51</volume>, <fpage>645</fpage>&#x2013;<lpage>658</lpage>. <pub-id pub-id-type="doi">10.1002/(SICI)1097-0290(19960920)51:6&#x3c;645:AID-BIT4&#x3e;3.0.CO;2-H</pub-id> </citation>
</ref>
<ref id="B36">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Schillinger</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Maurer</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Wagner</surname>
<given-names>E. C.</given-names>
</name>
<name>
<surname>van Ommen</surname>
<given-names>J. R.</given-names>
</name>
<name>
<surname>Schildhauer</surname>
<given-names>T. J.</given-names>
</name>
<name>
<surname>Schildhauer</surname>
<given-names>T. J.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Influence of Vertical Heat Exchanger Tubes, Their Arrangement and the Column Diameter on the Hydrodynamics in a Gas-Solid Bubbling Fluidized Bed</article-title>. <source>Int. J. Multiphase Flow</source> <volume>97</volume>, <fpage>46</fpage>&#x2013;<lpage>59</lpage>. <pub-id pub-id-type="doi">10.1016/j.ijmultiphaseflow.2017.07.013</pub-id> </citation>
</ref>
<ref id="B37">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Schlautmann</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>B&#xf6;hm</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Zauner</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>M&#xf6;rs</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Tichler</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Graf</surname>
<given-names>F.</given-names>
</name>
<etal/>
</person-group> (<year>2021</year>). <article-title>Renewable Power-To-Gas: A Technical and Economic Evaluation of Three Demo Sites within the STORE&#x26;GO Project</article-title>. <source>Chem. Ingenieur Technik</source> <volume>93</volume>, <fpage>568</fpage>&#x2013;<lpage>579</lpage>. <pub-id pub-id-type="doi">10.1002/cite.202000187</pub-id> </citation>
</ref>
<ref id="B38">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Schmid</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Decurtins</surname>
<given-names>D.</given-names>
</name>
</person-group> (<year>2018</year>). <source>Energieversorgung der Zukunft</source>. </citation>
</ref>
<ref id="B39">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Schmidt</surname>
<given-names>O.</given-names>
</name>
<name>
<surname>Gambhir</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Staffell</surname>
<given-names>I.</given-names>
</name>
<name>
<surname>Hawkes</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Nelson</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Few</surname>
<given-names>S.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Future Cost and Performance of Water Electrolysis: An Expert Elicitation Study</article-title>. <source>Int. J. Hydrogen. Energ.</source> <volume>42</volume>, <fpage>30470</fpage>&#x2013;<lpage>30492</lpage>. <pub-id pub-id-type="doi">10.1016/j.ijhydene.2017.10.045</pub-id> </citation>
</ref>
<ref id="B56">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Specht</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Brellochs</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Frick</surname>
<given-names>V.</given-names>
</name>
<name>
<surname>St&#x00FC;rmer</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Zuberb&#x00FC;hler</surname>
<given-names>U.</given-names>
</name>
</person-group> (<year>2016</year>). <article-title>The Power to Gas Process</article-title>. <source>Synthetic Natural Gas from Coal, Dry Biomass, and Power&#x2010;to&#x2010;Gas Applications</source>. <publisher-name>John Wiley &#x0026; Sons, Ltd</publisher-name>, <fpage>191</fpage>&#x2013;<lpage>220</lpage>. <pub-id pub-id-type="doi">10.1002/9781119191339.ch7</pub-id> </citation>
</ref>
<ref id="B40">
<citation citation-type="web">
<collab>Store and Go Project</collab> (<year>2016b</year>). <comment>Available at: <ext-link ext-link-type="uri" xlink:href="https://www.storeandgo.info/">https://www.storeandgo.info/</ext-link>
</comment> (<comment>Accessed February 11, 2020</comment>).</citation>
</ref>
<ref id="B41">
<citation citation-type="web">
<collab>Store and Go Project</collab> (<year>2016a</year>). <article-title>Store&#x26;Go Demonstration Site in Switzerland</article-title>. <comment>Available at: <ext-link ext-link-type="uri" xlink:href="https://www.storeandgo.info/demonstration-sites/switzerland/">https://www.storeandgo.info/demonstration-sites/switzerland/</ext-link>
</comment>(<comment>Accessed February 11, 2020</comment>). </citation>
</ref>
<ref id="B42">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Teske</surname>
<given-names>L. S.</given-names>
</name>
</person-group> (<year>2014</year>). <source>Integrating Rate Based Models into a Multi-Objective Process Design &#x26; Optimisation Framework Using Surrogate Models</source>, <fpage>6302</fpage>. <pub-id pub-id-type="doi">10.5075/epfl-thesis-6302</pub-id> </citation>
</ref>
<ref id="B43">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Thees</surname>
<given-names>O.</given-names>
</name>
<name>
<surname>Burg</surname>
<given-names>V.</given-names>
</name>
<name>
<surname>Erni</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Bowman</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Lemm</surname>
<given-names>R.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Biomassepotenziale der Schweiz f&#xfc;r die energetische Nutzung. Ergebnisse des Schweizerischen Energiekompetenzzentrums SCCER BIOSWEET</article-title>. <source>WSL Berichte</source> <volume>57</volume>, <fpage>299</fpage>. <comment>Available at: <ext-link ext-link-type="uri" xlink:href="https://www.dora.lib4ri.ch/wsl/islandora/object/wsl:13277">https://www.dora.lib4ri.ch/wsl/islandora/object/wsl:13277</ext-link>
</comment>. </citation>
</ref>
<ref id="B44">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Thema</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Bauer</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Sterner</surname>
<given-names>M.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Power-to-Gas: Electrolysis and Methanation Status Review</article-title>. <source>Renew. Sustain. Energ. Rev.</source> <volume>112</volume>, <fpage>775</fpage>&#x2013;<lpage>787</lpage>. <pub-id pub-id-type="doi">10.1016/j.rser.2019.06.030</pub-id> </citation>
</ref>
<ref id="B45">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Thema</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Weidlich</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Kaul</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>B&#xf6;llmann</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Huber</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Bellack</surname>
<given-names>A.</given-names>
</name>
<etal/>
</person-group> (<year>2021</year>). <article-title>Optimized Biological CO<sub>2</sub>-methanation with a Pure Culture of Thermophilic Methanogenic Archaea in a Trickle-Bed Reactor</article-title>. <source>Bioresour. Techn.</source> <volume>333</volume>, <fpage>125135</fpage>. <pub-id pub-id-type="doi">10.1016/j.biortech.2021.125135</pub-id> </citation>
</ref>
<ref id="B46">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Tribe</surname>
<given-names>M. A.</given-names>
</name>
<name>
<surname>Alpine</surname>
<given-names>R. L. W.</given-names>
</name>
</person-group> (<year>1986</year>). <article-title>Scale Economies and the &#x201C;0.6 Rule&#x201d;</article-title>. <source>Eng. Costs Prod. Econ.</source> <volume>10</volume>, <fpage>271</fpage>&#x2013;<lpage>278</lpage>. <pub-id pub-id-type="doi">10.1016/0167-188X(86)90053-4</pub-id> </citation>
</ref>
<ref id="B47">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ullrich</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Lemmer</surname>
<given-names>A.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Performance Enhancement of Biological Methanation with Trickle Bed Reactors by Liquid Flow Modulation</article-title>. <source>GCB Bioenergy</source> <volume>11</volume>, <fpage>63</fpage>&#x2013;<lpage>71</lpage>. <pub-id pub-id-type="doi">10.1111/gcbb.12547</pub-id> </citation>
</ref>
<ref id="B48">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Ulrich</surname>
<given-names>G. D.</given-names>
</name>
<name>
<surname>Vasudevan</surname>
<given-names>P. T.</given-names>
</name>
</person-group> (<year>2004</year>). <source>Chemical Engineering Process Design and Economics: A Practical Guide</source>. <edition>2nd ed.</edition> <publisher-loc>Durham</publisher-loc>: <publisher-name>N.H.: Process Publishing</publisher-name>. </citation>
</ref>
<ref id="B49">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Werther</surname>
<given-names>J.</given-names>
</name>
</person-group> (<year>1976</year>). <article-title>Die Bedeutung der Blasenkoaleszenz f&#xfc;r die Auslegung von Gas/Feststoff-Wirbelschichten</article-title>. <source>Chem. Ingenieur Technik</source> <volume>48</volume>, <fpage>339</fpage>. <pub-id pub-id-type="doi">10.1002/cite.330480421</pub-id> </citation>
</ref>
<ref id="B50">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Witte</surname>
<given-names>J.</given-names>
</name>
</person-group> (<year>2018</year>). <source>Experimental and Techno-Economic Assessment of Catalytic Methanation of Biogas for Power-To-Gas Processes</source>. <pub-id pub-id-type="doi">10.3929/ETHZ-B-000304900</pub-id> </citation>
</ref>
<ref id="B51">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Witte</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Calbry-Muzyka</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Wieseler</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Hottinger</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>Biollaz</surname>
<given-names>S. M. A.</given-names>
</name>
<name>
<surname>Schildhauer</surname>
<given-names>T. J.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Demonstrating Direct Methanation of Real Biogas in a Fluidised Bed Reactor</article-title>. <source>Appl. Energ.</source> <volume>240</volume>, <fpage>359</fpage>&#x2013;<lpage>371</lpage>. <pub-id pub-id-type="doi">10.1016/j.apenergy.2019.01.230</pub-id> </citation>
</ref>
<ref id="B52">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Witte</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Kunz</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Biollaz</surname>
<given-names>S. M. A.</given-names>
</name>
<name>
<surname>Schildhauer</surname>
<given-names>T. J.</given-names>
</name>
</person-group> (<year>2018a</year>). <article-title>Direct Catalytic Methanation of Biogas - Part II: Techno-Economic Process Assessment and Feasibility Reflections</article-title>. <source>Energ. Convers. Manag.</source> <volume>178</volume>, <fpage>26</fpage>&#x2013;<lpage>43</lpage>. <pub-id pub-id-type="doi">10.1016/j.enconman.2018.09.079</pub-id> </citation>
</ref>
<ref id="B53">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Witte</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Settino</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Biollaz</surname>
<given-names>S. M. A.</given-names>
</name>
<name>
<surname>Schildhauer</surname>
<given-names>T. J.</given-names>
</name>
</person-group> (<year>2018b</year>). <article-title>Direct Catalytic Methanation of Biogas - Part I: New Insights into Biomethane Production Using Rate-Based Modelling and Detailed Process Analysis</article-title>. <source>Energ. Convers. Manag.</source> <volume>171</volume>, <fpage>750</fpage>&#x2013;<lpage>768</lpage>. <pub-id pub-id-type="doi">10.1016/j.enconman.2018.05.056</pub-id> </citation>
</ref>
<ref id="B54">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhang</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Musselman</surname>
<given-names>I. H.</given-names>
</name>
<name>
<surname>Ferraris</surname>
<given-names>K. J.</given-names>
</name>
<name>
<surname>Balkus</surname>
<given-names>K. J.</given-names>
</name>
</person-group> (<year>2008</year>). <article-title>Gas Permeability Properties of Matrimid Membranes Containing the Metal-Organic Framework Cu-BPY-HFSS</article-title>. <source>J. Membr. Sci.</source> <volume>313</volume>, <fpage>170</fpage>&#x2013;<lpage>181</lpage>. <pub-id pub-id-type="doi">10.1016/j.memsci.2008.01.005</pub-id> </citation>
</ref>
</ref-list>
</back>
</article>