<?xml version="1.0" encoding="utf-8"?><!DOCTYPE article PUBLIC "-//NLM//DTD JATS (Z39.96) Journal Publishing DTD v1.0 20120330//EN" "JATS-journalpublishing1.dtd"><article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" article-type="research-article">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">TRAN</journal-id>
<journal-title-group><journal-title>Transport</journal-title></journal-title-group>
<issn pub-type="epub">1648-3480</issn><issn pub-type="ppub">1648-4142</issn><issn-l>1648-4142</issn-l>
</journal-meta>
<article-meta>
<article-id pub-id-type="publisher-id">TRAN000001</article-id>
<article-id pub-id-type="doi">10.3846/transport.2018.145 </article-id>
<article-categories><subj-group subj-group-type="heading">
<subject>Research Article</subject></subj-group></article-categories>
<title-group>
<article-title>Optimization of natural gas transport pipeline network layout: A new methodology based on dominance degree model</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author">
<name><surname>Zhu</surname><given-names>Zhenjun</given-names></name><xref ref-type="aff" rid="aff1">1</xref><xref ref-type="corresp" rid="cor1">*</xref>
</contrib>
<contrib contrib-type="author">
<name><surname>Sun</surname><given-names>Chaoxu</given-names></name><xref ref-type="aff" rid="aff2">2</xref>
</contrib>
<contrib contrib-type="author">
<name><surname>Zeng</surname><given-names>Jun</given-names></name><xref ref-type="aff" rid="aff1">1</xref>
</contrib>
<contrib contrib-type="author">
<name><surname>Chen</surname><given-names>Guowei</given-names></name><xref ref-type="aff" rid="aff3">3</xref>
</contrib>
<aff id="aff1"><label>1</label>Jiangsu Key Laboratory of Urban ITS, School of Transportation, Southeast University, Nanjing, China</aff>
<aff id="aff2"><label>2</label>Zhejiang Provincial Natural Gas Development Co. Ltd, Hangzhou, China</aff>
<aff id="aff3"><label>3</label>Zhejiang Electric power construction Co. Ltd, Ningbo, China</aff>
</contrib-group>
<author-notes>
<corresp id="cor1"><label>*</label>Corresponding author. E-mail: <email>zhuzhenjun@seu.edu.cn</email>.</corresp>
</author-notes>
<pub-date pub-type="ppub"><year>2017</year></pub-date><volume>00</volume><issue>0</issue><fpage>1</fpage><lpage>8</lpage><history><date date-type="received"><day>23</day> <month>01</month> <year>2017</year></date><date date-type="rev-recd"><day>13</day> <month>04</month> <year>2017</year></date><date date-type="accepted"><day>20</day> <month>04</month> <year>2017</year></date></history>
<permissions><copyright-statement>Copyright © 2017 The Author(s) Published by VGTU Press</copyright-statement><copyright-year>2017</copyright-year>
<license xlink:href="https://creativecommons.org/licenses/by/4.0/" license-type="open-access" xlink:type="simple">
<license-p>This is an Open Access article distributed under the terms of the Creative Commons Attribution License (<ext-link ext-link-type="uri" xlink:href="https://creativecommons.org/licenses/by/4.0/" xlink:type="simple">https://creativecommons.org/licenses/by/4.0/</ext-link>, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.</license-p></license></permissions>
<abstract>
<p>At the phase of 13th five-year plan in China, natural gas will play an important role in energy revolution. With the growth of consumption, natural gas infrastructures will become hot spots of future investment and pipeline network construction will also usher in a period of rapid development. Therefore, it is of great theoretical and practical significance to study layout methods of transport pipeline network. This paper takes natural gas transport pipeline network as a research object, introduces dominance degree to analyse benefits of pipeline projects. Then, this paper proposes Dominance Degree Model (DDM) of transport pipeline projects based on Potential Model (PM) and Economic Potential Theory (EPT). According to DDM of gas transport pipeline projects, layout methods of pipeline network are put forward, which is simple and easy to obtain the overall optimal solution and ensure maximum comprehensive benefits. What’s more, construction sequences of gas transport pipeline projects can be also determined. Finally, the model is applied to a real case of natural gas transport pipeline projects in Zhejiang Province, China. The calculation results suggest that the model should deal with the transport pipeline network layout problem well, which have important implications for other potential pipeline networks not only in the Zhejiang Province but also throughout China and beyond.</p>
</abstract>
<kwd-group>
<label>Keywords</label>
<kwd>natural gas</kwd>
<kwd>transport pipeline network</kwd>
<kwd>dominance degree model</kwd>
<kwd>potential model</kwd>
<kwd>economic potential theory</kwd>
<kwd>layout method</kwd>
</kwd-group>
</article-meta>
</front>
<body>
<sec id="x1-10001">
<label>1.</label>
<title>Introduction</title>
<p>As an efficient and cleaner energy, natural gas is an important bridge to make energy consumption transit to low-carbon (Gillingham et al., <xref ref-type="bibr" rid="ref007">2009</xref>). According to 13th five-year plan and National Energy Development Strategy Plan of China, with respect to the proportion of primary energy consumption, natural gas demand will continue to increase, its share will reach 10% by the year of 2020. Continued growth in consumption is bound to promote the construction of natural gas infrastructures, coverage area of trunk pipelines will in further expansion. Thus, gas transport and distribution network will be improved and facilities in different pipe network will also achieve interoperability gradually.</p>
<p>Natural gas transport pipeline network shoulders the task of gas supply, which plays a significant role in improving the overall socio-economic benefits for the region. However, pipeline projects’ investment is large and payback period is long, which cannot be easily reconstructed or expanded. So when we make the pipeline network layout, market supply and demand, benefits and construction costs of pipeline projects’ investments, etc. – should be considered to determine a reasonable layout scheme.</p>
<p>A number of studies have been conducted on layout optimization with respect to natural gas transport pipeline network, where most of this research has focused on mathematical optimization. For example, Edgar et al. (<xref ref-type="bibr" rid="ref004">1978</xref>) used Generalized Reduced Gradient Method to optimize natural gas transport pipeline network system. Bhaskaran and Salzborn (<xref ref-type="bibr" rid="ref003">1979</xref>) studied the optimal design problem of gas gathering pipeline network, which was divided into three sub-problem: system layout, node locations, and the diameter distribution. Pedrycz et al. (<xref ref-type="bibr" rid="ref012">1992</xref>) examined the design and application of neural networks in the context of a specific decision-making problem, the selection of an appropriate layout for natural gas distribution system. Singh and Nain (<xref ref-type="bibr" rid="ref016">2012</xref>), Sanaye and Mahmoudimehr (<xref ref-type="bibr" rid="ref015">2013</xref>) applied genetic algorithm to pipeline system optimization and simulated natural gas transport pipeline network. Üster and Dilaveroğlu (<xref ref-type="bibr" rid="ref018">2014</xref>) developed an integrated large-scale mixed-integer nonlinear optimization model to determine pipelines in the network. Schmidt et al. (<xref ref-type="bibr" rid="ref017">2015</xref>) presented stationary neuro-linguistic programming type models of gas networks that are primarily designed to include detailed nonlinear physics in the final optimization steps for mid-term planning problems. An and Peng (<xref ref-type="bibr" rid="ref001">2016</xref>) used risk cost functions to study the synchronization of the minimum risk loss and total cost of natural gas pipeline networks at the planning stage. However, these methods are essentially solving network structure to satisfy a given criterion, a local optimal solution can be only obtained.</p>
<table-wrap id="x1-10011">
<label>Table 1</label>
<caption>
<p>The methodology analysis</p>
</caption>
<table frame="hsides" rules="groups">
<thead>
<tr>
<td valign="top" align="left">Methods</td>
<td valign="top" align="left">Advantages</td>
<td valign="top" align="left">Limitations</td>
</tr>
</thead>
<tbody>
<tr>
<td valign="top" align="left">Graph theory</td>
<td valign="top" align="left">∙ easy access to computer program processing;<break/>∙ higher computing efficiency;<break/>∙ solving the shortest tree problem within a large scale network effectively</td>
<td valign="top" align="left">∙ only solving the network layout among known fixed points;<break/>∙ project investment costs are not considered;<break/>∙ solutions can only be used as the initial network layout</td>
</tr>
</tbody><tbody>
<tr>
<td valign="top" align="left">Dynamic programming</td>
<td valign="top" align="left">∙ solving optimization problems with multiple decision-making variables</td>
<td valign="top" align="left">∙ not suitable for dealing with large-scale network systems, ‘dimension obstacles’ exist in solving</td>
</tr>
</tbody><tbody>
<tr>
<td valign="top" align="left">Neural network method</td>
<td valign="top" align="left">∙ solving optimization problems with multiple decision-making variables</td>
<td valign="top" align="left">∙ only obtaining the local optimal solution</td>
</tr>
</tbody><tbody>
<tr>
<td valign="top" align="left">Genetic algorithm</td>
<td valign="top" align="left">∙ higher adaptability, which can overcome the difficulties in solving nonlinear optimization</td>
<td valign="top" align="left">∙ lower computing efficiency;<break/>∙ no effective quantitative analysis concerning algorithm precision, feasibility and computational complexity</td>
</tr>
</tbody><tbody>
<tr>
<td valign="top" align="left">Complex method</td>
<td valign="top" align="left">∙ the algorithm is simple and suitable for dealing with constrained optimization problems</td>
<td valign="top" align="left">∙ unable to deal with multi-variable, multi-constraint optimization problems</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>Minimum Spanning Tree Method (MSTM) and Dynamic Programming (DP) are the most commonly used solutions. MSTM abstracts pipeline network into undirected network, including classic solutions of Dijkstra Algorithms, Kruskal Algorithms, Steiner Algorithms (Graham and Hell, <xref ref-type="bibr" rid="ref006">1985</xref>). Compared to solutions of traditional graph theory, these algorithms adopt computer programs for processing, whose operation efficiency is relatively high, Steiner Algorithms effectively solve to the shortest path problem of large-scale network (Han and Lim, <xref ref-type="bibr" rid="ref008">2010</xref>; Sniedovich, <xref ref-type="bibr" rid="ref014">2006</xref>; Zhou, <xref ref-type="bibr" rid="ref019">2004</xref>). However, the above three algorithms do not consider investment costs of pipeline projects, whose results can be only regarded as the initial pipeline network layout. DP can deal the optimization problem of multiple decision-making variables, while ‘dimension obstacles’ exist in solving, that is computation will increase exponentially with growth of the number of variables. When the dimension of this problem increases to a certain degree, it cannot be solved (Bellman and Dreyfus, <xref ref-type="bibr" rid="ref002">2016</xref>). Through the above analysis, we can see that current commonly used methods have certain advantages and limitations, as shown in Table <xref rid="x1-10011">1</xref>.</p>
<p>In summary, layout optimization of regional natural gas transmission pipeline network is a multi-objective nonlinear programming problem; the above methods are difficult to solve the problem effectively in the practical application. Compared to the existing methods, a new methodology considering practical application is proposed to simplify computing and obtain the overall optimal solution in this paper, which can not only determine the optimal layout scheme, but also obtain construction sequences of gas transport pipeline projects, which ensure maximum comprehensive benefits of pipeline network. Specifically, we combine Potential Model (PM) and Economic Potential Theory (EPT) to establish the Dominance Degree Model (DDM) for pipeline projects, and use dominance degree to measure comprehensive benefits of pipeline projects. Then, a layout method based on the DDM for natural gas transport pipeline network is proposed, which can search out the most advantageous pipeline project successively through cycle calculations. The calculation process is simple and easy to obtain the overall optimal solution. To verify validity of the method, natural gas transport pipeline projects in Zhejiang Province are taken as a case study.</p>
</sec>
<sec id="x1-20002">
<label>2.</label>
<title>Dominance degree of pipeline projects</title>
<p>Dominance generally refers to the advantageous form that one party can overcome or overwhelm its counterpart. Dominance degree is an index, which can be used to reflect and compare the degree of the pros and cons among system elements (Shamsie, <xref ref-type="bibr" rid="ref013">2003</xref>). Due to differences in location, economy and geographical conditions, transport pipeline projects have made a tremendous impact on social and economic benefits and influenced dominance degree of network layout. Therefore, this paper uses dominance degree to reflect comprehensive benefits of transport pipeline projects and compare their comparative dominance.</p>
<p>As for regional natural gas transport pipeline network, dominance degree of transport pipeline projects are closely related to influential factors, such as functional orientation, scale, economic development level of cities and towns, construction costs of transport pipeline projects, etc.</p>
<sec id="x1-30002.1">
<label>2.1.</label>
<title>Functional orientation of cities and towns</title>
<p>Functional orientation of cities and towns is the core of their development and competition, cities and towns that play a dominant role in socio-economic development are important strategic nodes of regional natural gas transport pipeline network layout. They have different development advantages, such as location, resource and policy, etc. They can also maximize the optimal resource configuration, which offers good conditions for expanding natural gas consumption market.</p>
</sec>
<sec id="x1-40002.2">
<label>2.2.</label>
<title>Scale and economic development level of cities and towns</title>
<p>Natural gas consumption is closely related to scale and economic development level of cities and towns, which determine consumption capacity and growth speed of natural gas. In order to make natural gas consumption market have more potential, transport pipeline projects require cities and towns with a certain concentration of consumer groups and economic strength. Therefore, scale and economic development level are the main factors that affect natural gas consumption.</p>
</sec>
<sec id="x1-50002.3">
<label>2.3.</label>
<title>Construction costs of pipeline projects</title>
<p>Construction costs of pipeline projects are related to its length and costs of crossing obstacles, which directly affect the investment and benefits of projects. Natural gas transport pipeline projects have a certain service life, the more construction costs we input, the longer payback period and the weaker profitability it will be. Therefore, construction costs should be considered to ensure profitability when making natural gas transport pipeline network layout.</p>
</sec>
</sec>
<sec id="x1-60003">
<label>3.</label>
<title>Related theories and models</title>
<sec id="x1-70003.1">
<label>3.1.</label>
<title>Potential model</title>
<p>PM is a commonly used model in the empirical study of regional economy, which is mainly used to measure a point influenced by a set of given points in the space. Firstly, interaction of a node with each other node is calculated. Then, comprehensive influence is obtained by summing up interaction of a node with each other point (Mátyás, <xref ref-type="bibr" rid="ref010">1998</xref>; Lu and Taur, <xref ref-type="bibr" rid="ref009">2006</xref>). Total potential of node <italic>i</italic> can be expressed as follows: 
<disp-formula>
<mml:math display="block" id="math001">
<mml:mtable displaystyle="true"><mml:mlabeledtr>
<mml:mtd id="x1-70011">
<mml:mtext>(1)</mml:mtext>
</mml:mtd>
<mml:mtd>
<mml:msub>
<mml:mrow>
<mml:mi mathvariant="italic">V</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">i</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>=</mml:mo>
<mml:mi mathvariant="italic">K</mml:mi>
<mml:mo>·</mml:mo>
<mml:munder>
<mml:mrow>
<mml:mstyle displaystyle="true">
<mml:mo largeop="true" movablelimits="false">∑</mml:mo></mml:mstyle>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">j</mml:mi>
</mml:mrow>
</mml:munder><mml:mstyle displaystyle="true">
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi mathvariant="italic">m</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">j</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi mathvariant="italic">d</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">i</mml:mi>
<mml:mi mathvariant="italic">j</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfrac>
</mml:mstyle>
<mml:mo mathvariant="normal">,</mml:mo>
</mml:mtd>
</mml:mlabeledtr></mml:mtable></mml:math>
</disp-formula> 
where: <italic>K</italic> is the dielectric constant; <inline-formula><mml:math id="math002">
<mml:msub>
<mml:mrow>
<mml:mi mathvariant="italic">m</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">j</mml:mi>
</mml:mrow>
</mml:msub></mml:math></inline-formula> is the scale of point <italic>j</italic>, according to different research questions, population and economy scale can be selected; <inline-formula><mml:math id="math003">
<mml:msub>
<mml:mrow>
<mml:mi mathvariant="italic">d</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">i</mml:mi>
<mml:mi mathvariant="italic">j</mml:mi>
</mml:mrow>
</mml:msub></mml:math></inline-formula> is the impedance between <italic>i</italic> and <italic>j</italic>.</p>
<p>In terms of cities and towns within the study area, PM can be used to analyse potential social and economic benefits generated by transport pipeline projects and compare their advantages and disadvantages. Therefore, this paper introduces PM to analyse dominance degree of natural gas transport pipeline projects.</p>
</sec>
<sec id="x1-80003.2">
<label>3.2.</label>
<title>Economic potential theory</title>
<p>EPT is a manifestation trend of economic energy (Chang, <xref ref-type="bibr" rid="ref005">2008</xref>). Natural gas transport pipelines undertake the important task of transporting natural gas, which plays a significant role in optimizing regional energy consumption structure, promoting economic development and improving environment quality. There, we consider natural gas transport pipelines will generate strong economic potential in its gas supply region.</p>
<p>According to the electromagnetic theory, current-carrying conductors will generate a magnetic field in its surroundings. Thus, each point in vicinity of current-carrying conductors owns the potential condition (Novoselov et al., <xref ref-type="bibr" rid="ref011">2004</xref>). Pipeline projects are just like current-carrying conductors, the impact on cities and towns is related to gas transport flow, the distance between cities and off-take stations, development level. In this paper, pipeline projects are simulated into current-carrying conductors, gas transport flow is simulated into electric current, the city <italic>P</italic> and off-take station <italic>A</italic> are simulated into nodes. Field intensity can be used to describe economic potential of node <italic>P</italic>, which can be repressed as follows: 
<disp-formula>
<mml:math display="block" id="math004">
<mml:mtable displaystyle="true"><mml:mlabeledtr>
<mml:mtd id="x1-80012">
<mml:mtext>(2)</mml:mtext>
</mml:mtd>
<mml:mtd>
<mml:mi mathvariant="italic">B</mml:mi>
<mml:mo>=</mml:mo><mml:mstyle displaystyle="true">
<mml:mfrac>
<mml:mrow>
<mml:mi mathvariant="italic">μ</mml:mi>
<mml:mo>·</mml:mo>
<mml:mi mathvariant="italic">I</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">a</mml:mi>
</mml:mrow>
</mml:mfrac>
</mml:mstyle>
<mml:mo mathvariant="normal">,</mml:mo>
</mml:mtd>
</mml:mlabeledtr></mml:mtable></mml:math>
</disp-formula> 
where: <italic>μ</italic> is undetermined coefficient; <italic>I</italic> is natural gas transport flow from node <italic>A</italic> to <italic>P</italic>; <italic>a</italic> is the distance between node <italic>P</italic> (cities or towns) and node <italic>A</italic> (off-take stations).</p>
</sec>
</sec>
<sec id="x1-90004" sec-type="methods">
<label>4.</label>
<title>Methodology</title>
<p>Comprehensive benefits of pipeline projects are related to functional orientation, scale, economic development level of cities or towns. Through the analysis of these influential factors, we find that PM can be used to analyse the dominance degree of pipeline projects. The specific construction method is shown in Fig. <xref rid="x1-90011">1</xref>.</p>
<fig id="x1-90011">
<label>Fig. 1.</label>
<caption>
<p>Flow chart of model construction.</p>
</caption>
<graphic xlink:href="f01.jpg"/>
</fig>
<sec id="x1-100004.1">
<label>4.1.</label>
<title>Model construction</title>
<p>Based on PM and EPT, this paper proposes following assumptions and methods when constructing DDM:</p>
<list list-type="bullet">
<list-item>
<p>pipeline projects will generate economic potential, which can be regarded as an index reflecting city scale;</p>
</list-item>
<list-item>
<p>construction costs of pipeline projects will affect their economical benefits directly, therefore, construction costs are taken as parameter reflecting impedance in PM;</p>
</list-item>
<list-item>
<p>the city’s economic potential influenced by pipeline projects is closely related to its development level, pipeline projects will directly affect the transformation of economic potential.</p>
</list-item>
</list>
<p>Therefore, some specific indicators can be selected to reflect these influential factors. Above all, the DDM of gas transport pipeline projects is established as follows: 
<disp-formula>
<mml:math display="block" id="math005">
<mml:mtable displaystyle="true"><mml:mlabeledtr>
<mml:mtd id="x1-100013">
<mml:mtext>(3)</mml:mtext>
</mml:mtd>
<mml:mtd>
<mml:msub>
<mml:mrow>
<mml:mi mathvariant="italic">A</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">i</mml:mi>
<mml:mi mathvariant="italic">j</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>=</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mi mathvariant="italic">K</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">i</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>·</mml:mo><mml:mstyle displaystyle="true">
<mml:mfrac>
<mml:mrow>
<mml:mo mathvariant="normal" fence="true" stretchy="false">(</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mi mathvariant="italic">μ</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">i</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>·</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mi mathvariant="italic">I</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">i</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo mathvariant="normal" fence="true" stretchy="false">)</mml:mo>
<mml:mo mathvariant="normal" stretchy="false">/</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mi mathvariant="italic">a</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">i</mml:mi>
<mml:mi mathvariant="italic">j</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi mathvariant="italic">f</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">i</mml:mi>
<mml:mi mathvariant="italic">j</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfrac>
</mml:mstyle>
<mml:mo mathvariant="normal">,</mml:mo>
</mml:mtd>
</mml:mlabeledtr></mml:mtable></mml:math>
</disp-formula> 
where: <inline-formula><mml:math id="math006">
<mml:msub>
<mml:mrow>
<mml:mi mathvariant="italic">A</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">i</mml:mi>
<mml:mi mathvariant="italic">j</mml:mi>
</mml:mrow>
</mml:msub></mml:math></inline-formula> is the dominance degree of pipeline projects between city <italic>i</italic> and gas transport station <italic>j</italic>; <inline-formula><mml:math id="math007">
<mml:msub>
<mml:mrow>
<mml:mi mathvariant="italic">K</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">i</mml:mi>
</mml:mrow>
</mml:msub></mml:math></inline-formula> is the dielectric constant, which is a constant within a specific period for specific cities; <inline-formula><mml:math id="math008">
<mml:msub>
<mml:mrow>
<mml:mi mathvariant="italic">μ</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">i</mml:mi>
</mml:mrow>
</mml:msub></mml:math></inline-formula> is the sensitivity coefficient and reflects sensitive degree that cities put on pipeline projects; <inline-formula><mml:math id="math009">
<mml:msub>
<mml:mrow>
<mml:mi mathvariant="italic">I</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">i</mml:mi>
</mml:mrow>
</mml:msub></mml:math></inline-formula> is the gas demand of city <italic>I</italic>; <inline-formula><mml:math id="math010">
<mml:msub>
<mml:mrow>
<mml:mi mathvariant="italic">a</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">i</mml:mi>
<mml:mi mathvariant="italic">j</mml:mi>
</mml:mrow>
</mml:msub></mml:math></inline-formula> is the distance between city <italic>i</italic> and gas transport station <italic>j</italic>, which can be expressed by straight-line distance; <inline-formula><mml:math id="math011">
<mml:msub>
<mml:mrow>
<mml:mi mathvariant="italic">f</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">i</mml:mi>
<mml:mi mathvariant="italic">j</mml:mi>
</mml:mrow>
</mml:msub></mml:math></inline-formula> is construction costs of pipeline projects between city <italic>i</italic> and natural gas transport station <italic>j</italic>.</p>
</sec>
<sec id="x1-110004.2">
<label>4.2.</label>
<title>Parameter calibration</title>
<p>Parameter calibration is critical for model application, this paper proposes following calibration methods for model parameters.</p>
<p><inline-formula><mml:math id="math012">
<mml:msub>
<mml:mrow>
<mml:mi mathvariant="italic">K</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">i</mml:mi>
</mml:mrow>
</mml:msub></mml:math></inline-formula> is the dielectric constant, As for DDM of natural gas transport pipeline projects, which is closely related to urban functional orientation of cities and towns. Functional orientation is a manifestation of government functions, which guides many planning system, such as urban policies, industry structure, etc. City scale is of great significance for functional orientation and relatively stable in a long period of time. Therefore, in order to simplify this question, urban scale <italic>P</italic> is selected to reflect this parameter. With the reference to notification about adjusting classification criteria for city scale of 2014, <inline-formula><mml:math id="math013">
<mml:msub>
<mml:mrow>
<mml:mi mathvariant="italic">K</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">i</mml:mi>
</mml:mrow>
</mml:msub></mml:math></inline-formula> can be calibrated as shown in Table <xref rid="x1-110012">2</xref>.</p>
<table-wrap id="x1-110012">
<label>Table 2</label>
<caption>
<p>Calibration for the value of <inline-formula><mml:math id="math014">
<mml:msub>
<mml:mrow>
<mml:mi mathvariant="italic">K</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">i</mml:mi>
</mml:mrow>
</mml:msub></mml:math></inline-formula></p>
</caption>
<table frame="hsides" rules="groups">
<thead>
<tr>
<td valign="top" align="left"><italic>P</italic> [10000 people]</td>
<td valign="top" align="left"><inline-formula><mml:math id="math015">
<mml:msub>
<mml:mrow>
<mml:mi mathvariant="italic">K</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">i</mml:mi>
</mml:mrow>
</mml:msub></mml:math></inline-formula></td>
</tr>
</thead>
<tbody>
<tr>
<td valign="top" align="left">&lt;20</td>
<td valign="top" align="left">1</td>
</tr>
<tr>
<td valign="top" align="left">20–50</td>
<td valign="top" align="left">1.25</td>
</tr>
<tr>
<td valign="top" align="left">50–100</td>
<td valign="top" align="left">1.5</td>
</tr>
<tr>
<td valign="top" align="left">100–300</td>
<td valign="top" align="left">2.0</td>
</tr>
<tr>
<td valign="top" align="left">300–500</td>
<td valign="top" align="left">2.5</td>
</tr>
<tr>
<td valign="top" align="left">500–1000</td>
<td valign="top" align="left">3.0</td>
</tr>
</tbody>
</table>
</table-wrap>
<p><inline-formula><mml:math id="math016">
<mml:msub>
<mml:mrow>
<mml:mi mathvariant="italic">μ</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">i</mml:mi>
</mml:mrow>
</mml:msub></mml:math></inline-formula> is related to the richness of urban resources, market development level, living standards of residents, the existing state of industry. In this paper, Gross Domestic Product (GDP) per capital of macroeconomic indicators is used to represent this parameter.</p>
<p><inline-formula><mml:math id="math017">
<mml:msub>
<mml:mrow>
<mml:mi mathvariant="italic">I</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">i</mml:mi>
</mml:mrow>
</mml:msub></mml:math></inline-formula> is the gas demand of city <italic>i</italic>, which is related to local energy consumption structure and policies. In order to facilitate analysis of this question, energy consumption per unit GDP is selected to reflect energy consumption and saving status. To some extent, it reflects current gas demand and development potential of future natural gas.</p>
<p><inline-formula><mml:math id="math018">
<mml:msub>
<mml:mrow>
<mml:mi mathvariant="italic">f</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">i</mml:mi>
<mml:mi mathvariant="italic">j</mml:mi>
</mml:mrow>
</mml:msub></mml:math></inline-formula> is construction costs of pipeline projects between city <italic>i</italic> and gas transport station <italic>j</italic>, which will directly affect investment income. Natural gas transport pipeline projects will cross rivers, roads, railways, etc. Thus, construction costs are mainly determined by pipeline length and investment for crossing obstacles. Therefore, the calculation method of <inline-formula><mml:math id="math019">
<mml:msub>
<mml:mrow>
<mml:mi mathvariant="italic">f</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">i</mml:mi>
<mml:mi mathvariant="italic">j</mml:mi>
</mml:mrow>
</mml:msub></mml:math></inline-formula> is shown as follows: 
<disp-formula>
<mml:math display="block" id="math020">
<mml:mtable displaystyle="true"><mml:mlabeledtr>
<mml:mtd id="x1-110024">
<mml:mtext>(4)</mml:mtext>
</mml:mtd>
<mml:mtd>
<mml:msub>
<mml:mrow>
<mml:mi mathvariant="italic">f</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">i</mml:mi>
<mml:mi mathvariant="italic">j</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>=</mml:mo>
<mml:mi mathvariant="italic">V</mml:mi>
<mml:mo>·</mml:mo>
<mml:mi mathvariant="italic">G</mml:mi>
<mml:mo>·</mml:mo>
<mml:mi mathvariant="italic">β</mml:mi>
<mml:mo>·</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mi mathvariant="italic">a</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">i</mml:mi>
<mml:mi mathvariant="italic">j</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>+</mml:mo>
<mml:munderover accentunder="false" accent="false">
<mml:mrow>
<mml:mstyle displaystyle="true">
<mml:mo largeop="true" movablelimits="false">∑</mml:mo></mml:mstyle>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">k</mml:mi>
<mml:mo>=</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mn>4</mml:mn>
</mml:mrow>
</mml:munderover>
<mml:msub>
<mml:mrow>
<mml:mi mathvariant="italic">m</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">k</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>·</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mi mathvariant="italic">w</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">k</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo mathvariant="normal">,</mml:mo>
</mml:mtd>
</mml:mlabeledtr></mml:mtable></mml:math>
</disp-formula> 
where: <italic>V</italic> is pipeline price per unit weight; <italic>G</italic> is pipeline weight per unit length (the value of <italic>V</italic> and <italic>G</italic> can take pipelines of the same diameter, material, thickness as standard); <italic>β</italic> is terrain correction coefficient, which is used to modify pipeline route length (different topography conditions are endowed with different values, which can be determined by referring to engineering data within the study area); <inline-formula><mml:math id="math021">
<mml:msub>
<mml:mrow>
<mml:mi mathvariant="italic">m</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">k</mml:mi>
</mml:mrow>
</mml:msub></mml:math></inline-formula> is crossing times of the <italic>k</italic>th type’ <inline-formula><mml:math id="math022">
<mml:msub>
<mml:mrow>
<mml:mi mathvariant="italic">w</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">k</mml:mi>
</mml:mrow>
</mml:msub></mml:math></inline-formula> is the corresponding construction costs (<inline-formula><mml:math id="math023">
<mml:mi mathvariant="italic">k</mml:mi>
<mml:mo>=</mml:mo>
<mml:mn>1</mml:mn>
<mml:mo mathvariant="normal">,</mml:mo>
<mml:mn>2</mml:mn>
<mml:mo mathvariant="normal">,</mml:mo>
<mml:mn>3</mml:mn>
<mml:mo mathvariant="normal">,</mml:mo>
<mml:mn>4</mml:mn></mml:math></inline-formula> respectively represents the river, grade highway, substandard highway and railway).</p>
</sec>
<sec id="x1-120004.3">
<label>4.3.</label>
<title>The layout method based on DDM</title>
<p>DDM of pipeline projects fully embody comprehensive socio-economic benefits of the projects. Apply DDM to natural gas transport pipeline network layout will maximum socio-economic benefits. Through the analysis of DDM of pipeline projects, the optimal layout scheme and construction sequences are determined, which will give full play to advantages of pipeline projects and promote regional socio-economic development. The basic application process is shown in Fig. <xref rid="x1-120012">2</xref>.</p>
<fig id="x1-120012">
<label>Fig. 2.</label>
<caption>
<p>Natural gas transport pipeline network layout process.</p>
</caption>
<graphic xlink:href="f02.jpg"/>
</fig>
<p>According to the layout process, layout methods of natural gas transport network based on DDM have four steps, which are shown as follows:</p>
<list>
<list-item>
<label>–</label>
<p><italic>Step 1</italic>: Abstract pipeline network in the regional area. First, abstract gas transport stations and valve chests into spatial nodes, which record a set of nodes <inline-formula><mml:math id="math024">
<mml:mi mathvariant="italic">j</mml:mi>
<mml:mo>=</mml:mo>
<mml:mo fence="true" stretchy="false">{</mml:mo>
<mml:mn>1</mml:mn>
<mml:mo mathvariant="normal">,</mml:mo>
<mml:mn>2</mml:mn>
<mml:mo mathvariant="normal">,</mml:mo>
<mml:mo>…</mml:mo>
<mml:mo mathvariant="normal">,</mml:mo>
<mml:mi mathvariant="italic">m</mml:mi>
<mml:mo fence="true" stretchy="false">}</mml:mo></mml:math></inline-formula>. Using Geographic Information System (GIS), their locations in the map can be determined. Then, connect these nodes based on pipeline alignments and obtain spatial structure diagram of regional pipeline network;</p>
</list-item>
<list-item>
<label>–</label>
<p><italic>Step 2</italic>: Calculate the dominance degree of pipeline projects. Similarly, abstract cities and towns that uncovered by gas transmission pipeline network into a set of nodes <inline-formula><mml:math id="math025">
<mml:mi mathvariant="italic">i</mml:mi>
<mml:mo>=</mml:mo>
<mml:mo fence="true" stretchy="false">{</mml:mo>
<mml:mn>1</mml:mn>
<mml:mo mathvariant="normal">,</mml:mo>
<mml:mn>2</mml:mn>
<mml:mo mathvariant="normal">,</mml:mo>
<mml:mo>…</mml:mo>
<mml:mo mathvariant="normal">,</mml:mo>
<mml:mi mathvariant="italic">n</mml:mi>
<mml:mo fence="true" stretchy="false">}</mml:mo></mml:math></inline-formula>. Calibrate model parameters and use DDM to calculate dominance degree between node <italic>i</italic> and node <italic>j</italic> sequentially. Then, sort calculation results in descending order.</p>
</list-item>
<list-item>
<label>–</label>
<p><italic>Step 3</italic>: Update pipeline network in the regional area. According to calculation results, take the pipeline project of maximum dominance degree and connect its nodes at both ends to form a new pipeline route. Then, corresponding nodes <italic>i</italic> is incorporated into a set of node <italic>j</italic>, that is <inline-formula><mml:math id="math026">
<mml:mi mathvariant="italic">j</mml:mi>
<mml:mo>=</mml:mo>
<mml:mo fence="true" stretchy="false">{</mml:mo>
<mml:mn>1</mml:mn>
<mml:mo mathvariant="normal">,</mml:mo>
<mml:mn>2</mml:mn>
<mml:mo mathvariant="normal">,</mml:mo>
<mml:mo>…</mml:mo>
<mml:mo mathvariant="normal">,</mml:mo>
<mml:mi mathvariant="italic">m</mml:mi>
<mml:mo>+</mml:mo>
<mml:mn>1</mml:mn>
<mml:mo fence="true" stretchy="false">}</mml:mo></mml:math></inline-formula>. For example, assume that <inline-formula><mml:math id="math027">
<mml:mi mathvariant="italic">i</mml:mi>
<mml:mo>=</mml:mo>
<mml:mo fence="true" stretchy="false">{</mml:mo>
<mml:mn>1</mml:mn>
<mml:mo mathvariant="normal">,</mml:mo>
<mml:mn>2</mml:mn>
<mml:mo mathvariant="normal">,</mml:mo>
<mml:mn>3</mml:mn>
<mml:mo mathvariant="normal">,</mml:mo>
<mml:mn>4</mml:mn>
<mml:mo fence="true" stretchy="false">}</mml:mo></mml:math></inline-formula> and <inline-formula><mml:math id="math028">
<mml:mi mathvariant="italic">j</mml:mi>
<mml:mo>=</mml:mo>
<mml:mo fence="true" stretchy="false">{</mml:mo>
<mml:mn>1</mml:mn>
<mml:mo mathvariant="normal">,</mml:mo>
<mml:mn>2</mml:mn>
<mml:mo fence="true" stretchy="false">}</mml:mo></mml:math></inline-formula>, the calculated maximum dominance degree <inline-formula><mml:math id="math029">
<mml:msub>
<mml:mrow>
<mml:mi mathvariant="italic">A</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>21</mml:mn>
</mml:mrow>
</mml:msub></mml:math></inline-formula>, which means <inline-formula><mml:math id="math030">
<mml:mi mathvariant="italic">i</mml:mi>
<mml:mo>=</mml:mo>
<mml:mn>2</mml:mn></mml:math></inline-formula> is removed from a set of node <italic>i</italic>, subsumed into a set of node <italic>j</italic> and denoted as <inline-formula><mml:math id="math031">
<mml:mi mathvariant="italic">j</mml:mi>
<mml:mo>=</mml:mo>
<mml:mn>3</mml:mn></mml:math></inline-formula>, then <inline-formula><mml:math id="math032">
<mml:mi mathvariant="italic">i</mml:mi>
<mml:mo>=</mml:mo>
<mml:mo fence="true" stretchy="false">{</mml:mo>
<mml:mn>1</mml:mn>
<mml:mo mathvariant="normal">,</mml:mo>
<mml:mn>3</mml:mn>
<mml:mo mathvariant="normal">,</mml:mo>
<mml:mn>4</mml:mn>
<mml:mo fence="true" stretchy="false">}</mml:mo></mml:math></inline-formula> and <inline-formula><mml:math id="math033">
<mml:mi mathvariant="italic">j</mml:mi>
<mml:mo>=</mml:mo>
<mml:mo fence="true" stretchy="false">{</mml:mo>
<mml:mn>1</mml:mn>
<mml:mo mathvariant="normal">,</mml:mo>
<mml:mn>2</mml:mn>
<mml:mo mathvariant="normal">,</mml:mo>
<mml:mn>3</mml:mn>
<mml:mo fence="true" stretchy="false">}</mml:mo></mml:math></inline-formula>. Finally, updated spatial structure diagram of regional pipeline network is obtained.</p>
</list-item>
<list-item>
<label>–</label>
<p><italic>Step 4</italic>: Calculate cyclically and output the results. Repeat Step 2 and Step 3 until the set of nodes <inline-formula><mml:math id="math034">
<mml:mi mathvariant="italic">i</mml:mi>
<mml:mo>=</mml:mo>
<mml:mo fence="true" stretchy="false">{</mml:mo>
<mml:mn>1</mml:mn>
<mml:mo mathvariant="normal">,</mml:mo>
<mml:mn>2</mml:mn>
<mml:mo mathvariant="normal">,</mml:mo>
<mml:mo>…</mml:mo>
<mml:mo mathvariant="normal">,</mml:mo>
<mml:mi mathvariant="italic">n</mml:mi>
<mml:mo fence="true" stretchy="false">}</mml:mo></mml:math></inline-formula> is an empty set, while the set of nodes <inline-formula><mml:math id="math035">
<mml:mi mathvariant="italic">j</mml:mi>
<mml:mo>=</mml:mo>
<mml:mo fence="true" stretchy="false">{</mml:mo>
<mml:mn>1</mml:mn>
<mml:mo mathvariant="normal">,</mml:mo>
<mml:mn>2</mml:mn>
<mml:mo mathvariant="normal">,</mml:mo>
<mml:mo>…</mml:mo>
<mml:mo mathvariant="normal">,</mml:mo>
<mml:mi mathvariant="italic">m</mml:mi>
<mml:mo>+</mml:mo>
<mml:mi mathvariant="italic">n</mml:mi>
<mml:mo fence="true" stretchy="false">}</mml:mo></mml:math></inline-formula>. At this time, all cities and towns have been covered by pipeline network and layout scheme of pipeline network is obtained.</p>
</list-item>
</list>
</sec>
</sec>
<sec id="x1-130005">
<label>5.</label>
<title>Case study</title>
<p>In this paper, Lishui region along Jinhua–Lishui–Wenzhou gas transport trunk pipeline engineering in Zhejiang Province, China is taken as a case study. DDM is used to determine the optimal layout scheme and construction sequences.</p>
<sec id="x1-140005.1">
<label>5.1.</label>
<title>Pipeline network and model parameters</title>
<p>There are three natural gas transport stations and two valve chests, GIS is used to determine their locations. Then, abstract pipeline network in the regional area as shown in Fig. <xref rid="x1-130013">3</xref>.</p>
<fig id="x1-130013">
<label>Fig. 3.</label>
<caption>
<p>Natural gas transport pipeline network.</p>
</caption>
<graphic xlink:href="f03.jpg"/>
</fig>
<table-wrap id="x1-130023">
<label>Table 3</label>
<caption>
<p>Input parameters of cities and towns</p>
</caption>
<table frame="hsides" rules="groups">
<thead>
<tr>
<td valign="top" align="left">Name</td>
<td valign="top" align="left"><italic>K</italic></td>
<td valign="top" align="left"><italic>M</italic> [10000 yuan/person]</td>
<td valign="top" align="left"><italic>I</italic> [tons of standard coal/10000 yuan]</td>
</tr>
</thead>
<tbody>
<tr>
<td valign="top" align="left">Suiyang County</td>
<td valign="top" align="left">1.25</td>
<td valign="top" align="left">46.596</td>
<td valign="top" align="left">0.470</td>
</tr>
<tr>
<td valign="top" align="left">Songyang County</td>
<td valign="top" align="left">1.25</td>
<td valign="top" align="left">44.037</td>
<td valign="top" align="left">0.456</td>
</tr>
<tr>
<td valign="top" align="left">Yunhe County</td>
<td valign="top" align="left">1</td>
<td valign="top" align="left">45.627</td>
<td valign="top" align="left">0.532</td>
</tr>
<tr>
<td valign="top" align="left">Jingning County</td>
<td valign="top" align="left">1</td>
<td valign="top" align="left">38.541</td>
<td valign="top" align="left">0.458</td>
</tr>
<tr>
<td valign="top" align="left">Longquan City</td>
<td valign="top" align="left">1.25</td>
<td valign="top" align="left">39.546</td>
<td valign="top" align="left">0.491</td>
</tr>
<tr>
<td valign="top" align="left">Qinyuan County</td>
<td valign="top" align="left">1.25</td>
<td valign="top" align="left">43.943</td>
<td valign="top" align="left">0.490</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>According to Statistical Yearbook of Lishui City of 2015 and energy consumption statistics of each county of Lishui city (Mátyás, <xref ref-type="bibr" rid="ref010">1998</xref>), input parameters are shown in Table <xref rid="x1-130023">3</xref>.</p>
<p>Take related data of Lishui section in the Jinhua–Lishui–Wenzhou pipeline engineering as reference, use L450M steel for standard pipelines. The value of <italic>β</italic> is 1.032; <italic>V</italic> is 0.77 ten thousand yuan/tons; <italic>G</italic> is 222.3 tons/km; <inline-formula><mml:math id="math036">
<mml:msub>
<mml:mrow>
<mml:mi mathvariant="italic">w</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub></mml:math></inline-formula>, <inline-formula><mml:math id="math037">
<mml:msub>
<mml:mrow>
<mml:mi mathvariant="italic">w</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub></mml:math></inline-formula>, <inline-formula><mml:math id="math038">
<mml:msub>
<mml:mrow>
<mml:mi mathvariant="italic">w</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>3</mml:mn>
</mml:mrow>
</mml:msub></mml:math></inline-formula>, <inline-formula><mml:math id="math039">
<mml:msub>
<mml:mrow>
<mml:mi mathvariant="italic">w</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>4</mml:mn>
</mml:mrow>
</mml:msub></mml:math></inline-formula> are respectively 250, 15, 4, 25 thousand yuan/time.</p>
</sec>
<sec id="x1-150005.2">
<label>5.2.</label>
<title>Model calculation results</title>
<p>Take relevant parameters into the DDM, through cycle calculation, the results are obtained as shown in Table <xref rid="x1-150014">4</xref>.</p>
<table-wrap id="x1-150014">
<label>Table 4</label>
<caption>
<p>Model calculation results</p>
</caption>
<table frame="hsides" rules="groups">
<thead>
<tr>
<td valign="top" align="left">Cycle times</td>
<td valign="top" align="left"><italic>i</italic></td>
<td valign="top" align="left"><italic>j</italic></td>
<td valign="top" align="left"><inline-formula><mml:math id="math040">
<mml:msub>
<mml:mrow>
<mml:mi mathvariant="italic">A</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">i</mml:mi>
<mml:mi mathvariant="italic">j</mml:mi>
</mml:mrow>
</mml:msub></mml:math></inline-formula></td>
</tr>
</thead>
<tbody>
<tr>
<td valign="top" align="left">1</td>
<td valign="top" align="left">Songyang County</td>
<td valign="top" align="left">Terminal station in Lishui City</td>
<td valign="top" align="char" char=".">6.63</td>
</tr>
<tr>
<td valign="top" align="left">2</td>
<td valign="top" align="left">Songyang County</td>
<td valign="top" align="left">Suichang</td>
<td valign="top" align="char" char=".">18.85</td>
</tr>
<tr>
<td valign="top" align="left">3</td>
<td valign="top" align="left">Songyang County</td>
<td valign="top" align="left">Yunhe</td>
<td valign="top" align="char" char=".">8.69</td>
</tr>
<tr>
<td valign="top" align="left">4</td>
<td valign="top" align="left">Yunhe County</td>
<td valign="top" align="left">Jingning</td>
<td valign="top" align="char" char=".">28.52</td>
</tr>
<tr>
<td valign="top" align="left">5</td>
<td valign="top" align="left">Longquan City</td>
<td valign="top" align="left">Yunhe County</td>
<td valign="top" align="char" char=".">6.14</td>
</tr>
<tr>
<td valign="top" align="left">6</td>
<td valign="top" align="left">Qinyuan County</td>
<td valign="top" align="left">Longquan City</td>
<td valign="top" align="char" char=".">5.22</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>According to calculation results, update natural gas transport pipeline network within the study area successively. Finally, the layout scheme of natural gas transport pipeline network is obtained, as shown in Fig. <xref rid="x1-150024">4</xref>.</p>
<fig id="x1-150024">
<label>Fig. 4.</label>
<caption>
<p>Layout scheme of natural gas transport pipeline network.</p>
</caption>
<graphic xlink:href="f04.jpg"/>
</fig>
<p>Based on considering economic benefits and costs, the layout scheme calculated by DDM is the overall optimal solution. Meanwhile, according to cycle calculation results, construction sequences of natural gas transport pipeline projects can be determined and maximum comprehensive benefits of natural gas transport pipeline network are ensured. Under the condition of funding constraints, the optimal construction sequence in Lishui region is arranged successively: from teminal station in Lishui City to Songyang County, Songyang County to Suichang County, Songyang County to Yunhe County, Yunhe County to Jingning County, Yunhe County to Longquan City, Longquan City to Qinyuan County. Gas transport stations can be also obtained, which is arranged at Songyang County, Yunhe County and Longquan City.</p>
</sec>
</sec>
<sec id="x1-160006">
<label>6.</label>
<title>Conclusions</title>
<p>The layout optimization of natural gas transport pipeline is a subject with social and economic benefits, but it is also a complex system problem. Based on influential factors analysis of natural gas transport pipeline projects’ dominance degree, this paper proposes layout method of natural gas transport pipeline network based on DDM. The model can not only determine the optimal layout scheme, but also obtain construction sequences of gas transport pipeline projects, which ensure maximum comprehensive benefits of pipeline network. Natural gas transport pipeline projects in Zhejiang Province, China is taken as a case study, which suggest that the model should deal with the transport pipeline network layout problem well. The study finding has important implications for other potential pipeline networks not only in the Zhejiang Province but also throughout China and beyond.</p>
<p>With respect to calculating construction costs of pipeline projects, this paper mainly considers pipeline material consumption and costs of crossing obstacles, etc. However, during pipeline engineering design and construction process, construction costs of pipeline projects are also influenced by proportion of tunnels, topography and policy, etc. Therefore, how to quantify construction costs of natural gas transport pipeline projects comprehensively need to be further researched. In addition, the layout scheme of natural gas transport pipeline network can be obtained by DDM. However, distribution of natural gas, pipeline network technology and changes of airflow direction under different conditions should also be researched to determine the location and quantity of natural gas stations and valve chests reasonably and implement rational utilization of natural gas resources.</p>
</sec>
<sec id="x1-17000x">
<title>Contribution</title>
<p>Zhenjun Zhu and Chaoxu Sun conceived and designed the research; Zhenjun Zhu and Jun Zeng analyzed the data; Zhenjun Zhu wrote the paper; Chaoxu Sun and Guowei Chen helped improve the figures and manuscript.</p>
<p>All authors have read and approved the final manuscript.</p>
</sec>
<sec id="x1-18000x">
<title>Disclosure statement</title>
<p>The authors declare that there are no conflict competing financial, professional and personal interests from other parties.</p>
</sec>
</body>
<back>
<ref-list>
<title>References</title>
<ref id="ref001">
<mixed-citation publication-type="journal"><string-name><surname>An</surname>, <given-names>J.</given-names></string-name>, &amp; <string-name><surname>Peng</surname>, <given-names>S.</given-names></string-name> (<year>2016</year>). <article-title>Layout optimization of natural gas network planning: Synchronizing minimum risk loss with total cost</article-title>. <source>Journal of Natural Gas Science and Engineering</source>, <volume>33</volume>, <fpage>255</fpage>–<lpage>263</lpage>. <ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1016/j.jngse.2016.05.017" xlink:type="simple">http://dx.doi.org/10.1016/j.jngse.2016.05.017</ext-link>.</mixed-citation>
</ref>
<ref id="ref002">
<mixed-citation publication-type="other"><string-name><surname>Bellman</surname>, <given-names>R. E.</given-names></string-name>, &amp; <string-name><surname>Dreyfus</surname>, <given-names>S. E.</given-names></string-name> (<year>2016</year>). <source>Applied Dynamic Programming</source>. <publisher-loc>Princeton</publisher-loc>: <publisher-name>Princeton University Press</publisher-name>. <comment>390 p.</comment></mixed-citation>
</ref>
<ref id="ref003">
<mixed-citation publication-type="journal"><string-name><surname>Bhaskaran</surname>, <given-names>S.</given-names></string-name>, &amp; <string-name><surname>Salzborn</surname>, <given-names>F.</given-names></string-name> (<year>1979</year>). <article-title>Optimal diameter assignment for gas pipeline networks</article-title>. <source>The Journal of the Australian Mathematical Society. Series B. Applied Mathematics</source>, <volume>21</volume>(<issue>2</issue>), <fpage>129</fpage>–<lpage>144</lpage>. <ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1017/S0334270000002009" xlink:type="simple">http://dx.doi.org/10.1017/S0334270000002009</ext-link>.</mixed-citation>
</ref>
<ref id="ref004">
<mixed-citation publication-type="other"><string-name><surname>Edgar</surname>, <given-names>T. F.</given-names></string-name>, <string-name><surname>Himmelblau</surname>, <given-names>D. M.</given-names></string-name>, &amp; <string-name><surname>Bickel</surname>, <given-names>T. C.</given-names></string-name> (1978). Optimal design of gas transmission networks. <italic>Society of Petroleum Engineers Journal</italic>, <italic>18</italic>(2). <ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.2118/6034-PA" xlink:type="simple">http://dx.doi.org/10.2118/6034-PA</ext-link>.</mixed-citation>
</ref>
<ref id="ref005">
<mixed-citation publication-type="journal"><string-name><surname>Chang</surname>, <given-names>S. J.</given-names></string-name> (<year>2008</year>). <article-title>A cultural and philosophical perspective on Korea’s education reform: A critical way to maintain Korea’s economic momentum</article-title>. <source>Academic Paper Series</source>, <volume>3</volume>(<issue>2</issue>), <fpage>1</fpage>–<lpage>11</lpage>.</mixed-citation>
</ref>
<ref id="ref006">
<mixed-citation publication-type="journal"><string-name><surname>Graham</surname>, <given-names>R. L.</given-names></string-name>, &amp; <string-name><surname>Hell</surname>, <given-names>P.</given-names></string-name> (<year>1985</year>). <article-title>On the history of the minimum spanning tree problem</article-title>. <source>Annals of the History of Computing</source>, <volume>7</volume>(<issue>1</issue>), <fpage>43</fpage>–<lpage>57</lpage>. <ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1109/MAHC.1985.10011" xlink:type="simple">http://dx.doi.org/10.1109/MAHC.1985.10011</ext-link>.</mixed-citation>
</ref>
<ref id="ref007">
<mixed-citation publication-type="journal"><string-name><surname>Gillingham</surname>, <given-names>K.</given-names></string-name>, <string-name><surname>Newell</surname>, <given-names>R. G.</given-names></string-name>, &amp; <string-name><surname>Palmer</surname>, <given-names>K.</given-names></string-name> (<year>2009</year>). <article-title>Energy efficiency economics and policy</article-title>. <source>Annual Review of Resource Economics</source>, <volume>1</volume>, <fpage>597</fpage>–<lpage>620</lpage>. <ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1146/annurev.resource.102308.124234" xlink:type="simple">http://dx.doi.org/10.1146/annurev.resource.102308.124234</ext-link>.</mixed-citation>
</ref>
<ref id="ref008">
<mixed-citation publication-type="journal"><string-name><surname>Han</surname>, <given-names>D.-M.</given-names></string-name>, &amp; <string-name><surname>Lim</surname>, <given-names>J.-H.</given-names></string-name> (<year>2010</year>). <article-title>Design and implementation of smart home energy management systems based on ZigBee</article-title>. <source>IEEE Transactions on Consumer Electronics</source>, <volume>56</volume>(<issue>3</issue>), <fpage>1417</fpage>–<lpage>1425</lpage>. <ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1109/TCE.2010.5606278" xlink:type="simple">http://dx.doi.org/10.1109/TCE.2010.5606278</ext-link>.</mixed-citation>
</ref>
<ref id="ref009">
<mixed-citation publication-type="journal"><string-name><surname>Lu</surname>, <given-names>H.</given-names></string-name>, &amp; <string-name><surname>Taur</surname>, <given-names>Y.</given-names></string-name> (<year>2006</year>). <article-title>An analytic potential model for symmetric and asymmetric DG MOSFETs</article-title>. <source>IEEE Transactions on Electron Devices</source>, <volume>53</volume>(<issue>5</issue>), <fpage>1161</fpage>–<lpage>1168</lpage>. <ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1109/TED.2006.872093" xlink:type="simple">http://dx.doi.org/10.1109/TED.2006.872093</ext-link>.</mixed-citation>
</ref>
<ref id="ref010">
<mixed-citation publication-type="journal"><string-name><surname>Mátyás</surname>, <given-names>L.</given-names></string-name> (<year>1998</year>). <article-title>The gravity model: Some econometric considerations</article-title>. <source>The World Economy</source>, <volume>21</volume>(<issue>3</issue>), <fpage>397</fpage>–<lpage>401</lpage>. <ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1111/1467-9701.00136" xlink:type="simple">http://dx.doi.org/10.1111/1467-9701.00136</ext-link>.</mixed-citation>
</ref>
<ref id="ref011">
<mixed-citation publication-type="journal"><string-name><surname>Novoselov</surname>, <given-names>K. S.</given-names></string-name>, <string-name><surname>Geim</surname>, <given-names>A. K.</given-names></string-name>, <string-name><surname>Morozov</surname>, <given-names>S. V.</given-names></string-name>, <string-name><surname>Jiang</surname>, <given-names>D.</given-names></string-name>, <string-name><surname>Zhang</surname>, <given-names>Y.</given-names></string-name>, <string-name><surname>Dubonos</surname>, <given-names>S. V.</given-names></string-name>, <string-name><surname>Grigorieva</surname>, <given-names>I. V.</given-names></string-name>, &amp; <string-name><surname>Firsov</surname>, <given-names>A. A.</given-names></string-name> (<year>2004</year>). <article-title>Electric field effect in atomically thin carbon films</article-title>. <source>Science</source>, <volume>306</volume>(<issue>5696</issue>), <fpage>666</fpage>–<lpage>669</lpage>. <ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1126/science.1102896" xlink:type="simple">http://dx.doi.org/10.1126/science.1102896</ext-link>.</mixed-citation>
</ref>
<ref id="ref012">
<mixed-citation publication-type="chapter"><string-name><surname>Pedrycz</surname>, <given-names>W.</given-names></string-name>, <string-name><surname>Davidson</surname>, <given-names>J.</given-names></string-name>, &amp; <string-name><surname>Goulter</surname>, <given-names>I.</given-names></string-name> (<year>1992</year>). <chapter-title>Neural network based decision model, used for design of rural natural gas systems</chapter-title>. In <source>IEEE International Conference on Fuzzy Systems 1992</source>, <conf-loc>San Diego, CA, US</conf-loc>, <conf-date>8–12 March 1992</conf-date> (pp. <fpage>1219</fpage>–<lpage>1226</lpage>). <ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1109/FUZZY.1992.258651" xlink:type="simple">http://dx.doi.org/10.1109/FUZZY.1992.258651</ext-link>.</mixed-citation>
</ref>
<ref id="ref013">
<mixed-citation publication-type="journal"><string-name><surname>Shamsie</surname>, <given-names>J.</given-names></string-name> (<year>2003</year>). <article-title>The context of dominance: An industry-driven framework for exploiting reputation</article-title>. <source>Strategic Management Journal</source>, <volume>24</volume>(<issue>3</issue>), <fpage>199</fpage>–<lpage>215</lpage>. <ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1002/smj.291" xlink:type="simple">http://dx.doi.org/10.1002/smj.291</ext-link>.</mixed-citation>
</ref>
<ref id="ref014">
<mixed-citation publication-type="journal"><string-name><surname>Sniedovich</surname>, <given-names>M.</given-names></string-name> (<year>2006</year>). <article-title>Dijkstra’s algorithm revisited: The dynamic programming connexion</article-title>. <source>Control and Cybernetics</source>, <volume>35</volume>(<issue>3</issue>), <fpage>599</fpage>–<lpage>620</lpage>.</mixed-citation>
</ref>
<ref id="ref015">
<mixed-citation publication-type="journal"><string-name><surname>Sanaye</surname>, <given-names>S.</given-names></string-name>, &amp; <string-name><surname>Mahmoudimehr</surname>, <given-names>J.</given-names></string-name> (<year>2013</year>). <article-title>Optimal design of a natural gas transmission network layout</article-title>. <source>Chemical Engineering Research and Design</source>, <volume>91</volume>(<issue>12</issue>), <fpage>2465</fpage>–<lpage>2476</lpage>. <ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1016/j.cherd.2013.04.005" xlink:type="simple">http://dx.doi.org/10.1016/j.cherd.2013.04.005</ext-link>.</mixed-citation>
</ref>
<ref id="ref016">
<mixed-citation publication-type="journal"><string-name><surname>Singh</surname>, <given-names>R. R.</given-names></string-name>, &amp; <string-name><surname>Nain</surname>, <given-names>P. K. S.</given-names></string-name> (<year>2012</year>). <article-title>Optimization of natural gas pipeline design and its total cost using GA</article-title>. <source>International Journal of Scientific and Research Publications</source>, <volume>2</volume>(<issue>8</issue>), <fpage>1</fpage>–<lpage>10</lpage>.</mixed-citation>
</ref>
<ref id="ref017">
<mixed-citation publication-type="journal"><string-name><surname>Schmidt</surname>, <given-names>M.</given-names></string-name>, <string-name><surname>Steinbach</surname>, <given-names>M. C.</given-names></string-name>, &amp; <string-name><surname>Willert</surname>, <given-names>B. M.</given-names></string-name> (<year>2015</year>). <article-title>High detail stationary optimization models for gas networks</article-title>. <source>Optimization and Engineering</source>, <volume>16</volume>(<issue>1</issue>), <fpage>131</fpage>–<lpage>164</lpage>. <ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1007/s11081-014-9246-x" xlink:type="simple">http://dx.doi.org/10.1007/s11081-014-9246-x</ext-link>.</mixed-citation>
</ref>
<ref id="ref018">
<mixed-citation publication-type="journal"><string-name><surname>Üster</surname>, <given-names>H.</given-names></string-name>, &amp; <string-name><surname>Dilaveroğlu</surname>, <given-names>S.</given-names></string-name> (<year>2014</year>). <article-title>Optimization for design and operation of natural gas transmission networks</article-title>. <source>Applied Energy</source>, <volume>133</volume>, <fpage>56</fpage>–<lpage>59</lpage>. <ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1016/j.apenergy.2014.06.042" xlink:type="simple">http://dx.doi.org/10.1016/j.apenergy.2014.06.042</ext-link>.</mixed-citation>
</ref>
<ref id="ref019">
<mixed-citation publication-type="journal"><string-name><surname>Zhou</surname>, <given-names>H.</given-names></string-name> (<year>2004</year>). <article-title>Efficient Steiner tree construction based on spanning graphs</article-title>. <source>IEEE Transactions on Computer-Aided Design of Integrated Circuits and Systems</source>, <volume>23</volume>(<issue>5</issue>), <fpage>704</fpage>–<lpage>710</lpage>. <ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1109/TCAD.2004.826557" xlink:type="simple">http://dx.doi.org/10.1109/TCAD.2004.826557</ext-link>.</mixed-citation>
</ref>
</ref-list>
</back>
</article>
