<?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. Earth Sci.</journal-id>
<journal-title>Frontiers in Earth Science</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Earth Sci.</abbrev-journal-title>
<issn pub-type="epub">2296-6463</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="publisher-id">1219544</article-id>
<article-id pub-id-type="doi">10.3389/feart.2023.1219544</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Earth Science</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Prediction for the surface settlement of double-track subway tunnels for shallow buried loess based on peck formula</article-title>
<alt-title alt-title-type="left-running-head">Wang et al.</alt-title>
<alt-title alt-title-type="right-running-head">
<ext-link ext-link-type="uri" xlink:href="https://doi.org/10.3389/feart.2023.1219544">10.3389/feart.2023.1219544</ext-link>
</alt-title>
</title-group>
<contrib-group>
<contrib contrib-type="author">
<name>
<surname>Wang</surname>
<given-names>Shuangyuan</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/2304634/overview"/>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Song</surname>
<given-names>Zhanping</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<xref ref-type="corresp" rid="c001">&#x2a;</xref>
<uri xlink:href="https://loop.frontiersin.org/people/1354881/overview"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Tian</surname>
<given-names>Xiaoxu</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/1756204/overview"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Sun</surname>
<given-names>Yinhao</given-names>
</name>
<xref ref-type="aff" rid="aff3">
<sup>3</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/2351996/overview"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Zhang</surname>
<given-names>Yuwei</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/1645587/overview"/>
</contrib>
</contrib-group>
<aff id="aff1">
<sup>1</sup>
<institution>School of Civil Engineering</institution>, <institution>Xi&#x2019;an University of Architecture and Technology</institution>, <addr-line>Xi&#x2019;an</addr-line>, <country>China</country>
</aff>
<aff id="aff2">
<sup>2</sup>
<institution>Key Laboratory of Geotechnical and Underground Space Engineering of Shaanxi Province</institution>, <addr-line>Xi&#x2019;an</addr-line>, <country>China</country>
</aff>
<aff id="aff3">
<sup>3</sup>
<institution>China Railway 20th Bureau Group First Engineering Co., Ltd.</institution>, <addr-line>Suzhou</addr-line>, <country>China</country>
</aff>
<author-notes>
<fn fn-type="edited-by">
<p>
<bold>Edited by:</bold> <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/1900689/overview">Yi Xue</ext-link>, Xi&#x2019;an University of Technology, China</p>
</fn>
<fn fn-type="edited-by">
<p>
<bold>Reviewed by:</bold> <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/1379526/overview">Pengpeng Ni</ext-link>, Sun Yat-Sen University, China</p>
<p>
<ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/2037276/overview">Kun Fang</ext-link>, China University of Geosciences Wuhan, China</p>
<p>
<ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/1871743/overview">Zhengzheng Cao</ext-link>, Henan Polytechnic University, China</p>
</fn>
<corresp id="c001">&#x2a;Correspondence: Zhanping Song, <email>songzhpyt@xauat.edu.cn</email>
</corresp>
</author-notes>
<pub-date pub-type="epub">
<day>26</day>
<month>06</month>
<year>2023</year>
</pub-date>
<pub-date pub-type="collection">
<year>2023</year>
</pub-date>
<volume>11</volume>
<elocation-id>1219544</elocation-id>
<history>
<date date-type="received">
<day>09</day>
<month>05</month>
<year>2023</year>
</date>
<date date-type="accepted">
<day>15</day>
<month>06</month>
<year>2023</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#xa9; 2023 Wang, Song, Tian, Sun and Zhang.</copyright-statement>
<copyright-year>2023</copyright-year>
<copyright-holder>Wang, Song, Tian, Sun and Zhang</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/">
<p>This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.</p>
</license>
</permissions>
<abstract>
<p>In the process of constructing double-track subway tunnels in shallow buried loess areas, the interaction of double-track tunnels is significantly influenced by the net distance and the cross-section size, which is challenging to control the surface settlement. Therefore, the surface settlement prediction is essential while constructing double-track subway tunnels in shallow buried loess areas. The paper analyzed the surface settlement law of shallow buried double-track tunnels in loess areas through theoretical research and numerical simulation. The research results show that with the decrease of the net distance, the surface settlement superimposed curve was double <italic>V</italic> shape -<italic>W</italic> shape - single <italic>V</italic> shape. When the superimposed curve is double <italic>V</italic> shape and <italic>W</italic> shape, the Peck formula was used to calculate the surface settlement curve of the single-track tunnel, then superimposed to obtain the final surface settlement curve. When the superimposed surface settlement curve was <italic>V</italic> shape, based on the Peck formula, the formula for predicting the surface settlement suitable for symmetry and asymmetry was established. The net distance ratio and the area ratio were defined, and considering the tunnel&#x2019;s interaction, the value and position of the maximum were corrected. Then numerical tests were carried out 16 times with different net distance ratios and area ratios, to determine the parameters of increments and position offsets of the maximum regarding the net distance ratio and the area ratio. Finally, two engineering were conducted for verifying the rationality and applicability exhibited by the prediction formula. The prediction formula served for predicting the surface settlement of double-track subway tunnels in shallow buried loess areas. Which can reduce construction risks and assure the safety of buildings above the ground.</p>
</abstract>
<kwd-group>
<kwd>prediction</kwd>
<kwd>surface settlement</kwd>
<kwd>double-track</kwd>
<kwd>subway tunnels</kwd>
<kwd>shallow buried</kwd>
<kwd>peck formula</kwd>
</kwd-group>
<custom-meta-wrap>
<custom-meta>
<meta-name>section-at-acceptance</meta-name>
<meta-value>Geohazards and Georisks</meta-value>
</custom-meta>
</custom-meta-wrap>
</article-meta>
</front>
<body>
<sec id="s1">
<title>1 Introduction</title>
<p>Loess can be broadly found in western China (<xref ref-type="bibr" rid="B11">He et al., 2013</xref>; <xref ref-type="bibr" rid="B25">Song et al., 2021</xref>; <xref ref-type="bibr" rid="B36">Xue et al., 2023a</xref>), during the implementation of the Western Development Strategy, the number of subway tunnels has increased in loess areas. Subway tunnels in loess areas are generally in the form of double-track (<xref ref-type="bibr" rid="B5">Ding, 2019</xref>; <xref ref-type="bibr" rid="B27">Song et al., 2020</xref>; <xref ref-type="bibr" rid="B29">Su et al., 2020</xref>; <xref ref-type="bibr" rid="B35">Xue et al., 2023b</xref>), affected by the occupied space and adjacent buildings, the forms of loess subway tunnels are gradually increasing. The double-track subway tunnel in shallow buried loess areas is affected by the tunnel burial depth, tunnel clearance and surrounding rock, and the construction is difficult to carry out (<xref ref-type="bibr" rid="B2">Chapman et al., 2004</xref>; <xref ref-type="bibr" rid="B7">Feng and Yu, 2019</xref>; <xref ref-type="bibr" rid="B6">Fan et al., 2021</xref>). The original equilibrium state of the stratum will be changed during the excavation, which brings risks to the safety of buildings above the ground (<xref ref-type="bibr" rid="B46">Zhou et al., 2021</xref>; <xref ref-type="bibr" rid="B1">Cao et al., 2023</xref>; <xref ref-type="bibr" rid="B4">Chen et al., 2023</xref>; <xref ref-type="bibr" rid="B37">Xue et al., 2023c</xref>).</p>
<p>A lot of previous studies focus on the laws of surface settlement resulting from the excavation of double-track tunnels (DTT) in loess areas. Based on the modified random medium theory, <xref ref-type="bibr" rid="B8">Han and Su (2022)</xref> used the standard transformation principle to obtain the polar coordinate expression of the analytical formula, and the research results could serve for predicting the surface settlement resulting from the construction of loess tunnels with arbitrary cross-section. <xref ref-type="bibr" rid="B24">She et al. (2015)</xref> summarized the characteristics exhibited by the surface longitudinal settlement trough resulting from the subway tunnel excavation, analyzed the surface displacement release rate and the largest slope of the surface longitudinal settlement, and proposed a prediction method for surface longitudinal settlement. <xref ref-type="bibr" rid="B9">Han and Wang (2012)</xref> analyzed the monitoring data of the Chaoyang Gate&#x2014;KangFu Road tunnel of Xi&#x2019;an subway Line 1, and studied the longitudinal and transverse surface settlement laws of tunnel excavation in loess areas, then the Peck formula served for establishing the surface settlement prediction model of the excavation of subway tunnels. <xref ref-type="bibr" rid="B41">Zhang et al. (2017)</xref> built a 3D elastoplastic model and studied the surface settlement law of DTT in loess areas under different net distances, the research results will serve for optimizing the net distance of double-track subway tunnels in loess areas. <xref ref-type="bibr" rid="B12">He et al. (2012)</xref> summarized the monitoring data of the surface settlement in the process of constructing DTT in saturated soft loess areas, and studied the laws regarding the maximum settlement as well as the ground loss rate. <xref ref-type="bibr" rid="B30">Suwansawat and Einstein (2007)</xref> calculated the surface settlement resulting from the single-track tunnel (STT) excavation using the Gaussian formula, and superimposed two curves to obtain the final surface settlement curve. Since the interaction between tunnels is not considered, the calculation results are largely different from the monitoring data. <xref ref-type="bibr" rid="B21">Oack (2014)</xref> calculated the ground settlement of the first excavation tunnel by using the Peck formula, and calculated the surface settlement of the second excavation tunnel by using the modified Gaussian curve, followed by superimposing the surface settlement of the above two tunnels for obtaining the final surface settlement. <xref ref-type="bibr" rid="B20">Mahmuto&#x11f;lu (2011)</xref> analyzed factors affecting the surface settlement of double-track symmetric tunnels in shallow buried areas, and found the association of the settlement trough change with the surrounding rock disturbance. <xref ref-type="bibr" rid="B15">Hunt (2005)</xref> studied the law of surface settlement resulting from DTT excavation through numerical simulation, considering the impact of the second excavation tunnel construction, and corrected the surface settlement curve. <xref ref-type="bibr" rid="B46">Zhou et al. (2021)</xref> studied the influencing mechanism regarding DTT excavation by using the Peck formula, for the disturbance resulting from the first excavation tunnel, the disturbance correction factor was input for deriving the prediction formula for the surface settlement resulting from DTT construction. <xref ref-type="bibr" rid="B18">Li (2014)</xref> analyzed the factors that caused the surface settlement by tunnel excavation in loess areas, and established a calculation model for surface settlement by using random medium theory. <xref ref-type="bibr" rid="B28">Song et al. (2022)</xref> conducted a study on the construction sequence of shallow buried tunnels with small spacing asymmetric section in loess areas, to control the surface settlement, the large-section tunnel should be excavated first. <xref ref-type="bibr" rid="B16">Jiang (2022)</xref> optimized the construction parameters of the double side heading method in shallow buried loess subway tunnels, and the result shows that increasing the spacing between the tunnel faces can reduce surface settlement, and the surface settlement curve is a single <italic>V</italic>-shape after excavation.</p>
<p>Taken together, there are more studies on the surface settlement resulting from double-track subway tunnel excavation in loess areas, but there is a lack of comprehensive research on the symmetry and net distance on the effect of the surface settlement. The paper first analyzes the mechanism of the surface settlement in double-track subway tunnels in shallow buried loess areas, then establishes a surface settlement prediction formula suitable for symmetric and asymmetric double-track subway tunnels in shallow buried loess areas, determines the parameters of increment and position offset of the maximum on surface settlement regarding the net distance ratio and the area ratio. Finally, the reliability of the prediction formula is proved through the application and monitoring comparison of actual projects, and the application sphere of the prediction formula is discussed.</p>
</sec>
<sec id="s2">
<title>2 Surface settlement mechanisms and prediction formula for double-track subway tunnels</title>
<sec id="s2-1">
<title>2.1 Peck formula</title>
<p>The Peck formula and the random medium theory are commonly used methods for predicting surface settlement, the Peck formula is simpler and easier to calculate than the random medium theory, in order to increase the practicality of the prediction formula, the paper studied the surface settlement based on the Peck formula, the Peck formula is Eqs <xref ref-type="disp-formula" rid="e1">1</xref>, <xref ref-type="disp-formula" rid="e2">2</xref>.<disp-formula id="e1">
<mml:math id="m1">
<mml:mrow>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>x</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mi>max</mml:mi>
</mml:msub>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>exp</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mfrac>
<mml:msup>
<mml:mi>x</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:msup>
<mml:mi>i</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(1)</label>
</disp-formula>
<disp-formula id="e2">
<mml:math id="m2">
<mml:mrow>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mi>max</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mo>&#x2212;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi mathvariant="normal">&#x3c0;</mml:mi>
<mml:msup>
<mml:mi>R</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mi>&#x3b7;</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:msqrt>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mi mathvariant="normal">&#x3c0;</mml:mi>
</mml:mrow>
</mml:msqrt>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
<label>(2)</label>
</disp-formula>
</p>
<p>Where <italic>S</italic>
<sub>(x)</sub> is the surface settlement value at any point of the surface settlement curve, <italic>S</italic>
<sub>max</sub> is the maximum of the surface settlement, <italic>x</italic> is the horizontal distance between the calculated point and the origin, <italic>i</italic> is the settlement trough width, <italic>R</italic> is the tunnel diameter, <italic>&#x3b7;</italic> is the ground loss rate.</p>
</sec>
<sec id="s2-2">
<title>2.2 Analysis of surface settlement mechanism</title>
<p>
<xref ref-type="fig" rid="F1">Figure 1</xref> displays the surface settlement curves with the net distance of double-track subway tunnels in shallow buried loess areas, where <italic>D</italic>
<sub>i</sub> is the tunnel diameter (the value of <italic>i</italic> is 1, 2) and <italic>L</italic>
<sub>j</sub> is the net distance of the tunnel (the value of <italic>j</italic> is 1, 2, 3, 4, 5, and 6). According to <xref ref-type="fig" rid="F1">Figure 1A</xref>, when the tunnel section is symmetric, two settlement curves do not affect each other when the tunnel net distance is <italic>L</italic>
<sub>1</sub>, the surface settlement curve is a double <italic>V</italic> shape, and the Peck formula can serve for calculating the curve and obtaining the maximum of settlement.</p>
<fig id="F1" position="float">
<label>FIGURE 1</label>
<caption>
<p>Surface settlement mechanism: <bold>(A)</bold> The net distance of the tunnel is <italic>L</italic>
<sub>1</sub>, <bold>(B)</bold> The net distance of the tunnel is <italic>L</italic>
<sub>2</sub>, <bold>(C)</bold> The net distance of the tunnel is <italic>L</italic>
<sub>3</sub>, <bold>(D)</bold> The net distance of the tunnel is <italic>L</italic>
<sub>4</sub>, <bold>(E)</bold> The net distance of the tunnel is <italic>L</italic>
<sub>5</sub>, <bold>(F)</bold> The net distance of the tunnel is <italic>L</italic>
<sub>6</sub>.</p>
</caption>
<graphic xlink:href="feart-11-1219544-g001.tif"/>
</fig>
<p>According to the research on the law of the surface settlement regarding DTT influenced by the net distance (<xref ref-type="bibr" rid="B3">Chen et al., 2014</xref>; <xref ref-type="bibr" rid="B19">Lu et al., 2019</xref>), as <italic>L</italic>
<sub>j</sub> decreases, the superimposed surface settlement curve is a <italic>W</italic> shape and the tunnel interaction is small, which is displayed in <xref ref-type="fig" rid="F1">Figure 1B</xref>. The superimposed part has no influence on the value and the position of the maximum, so the surface settlement curve can be obtained by the superposition of the Peck formula. In <xref ref-type="fig" rid="F1">Figure 1C</xref>, when <italic>L</italic>
<sub>j</sub> continues to decrease, the tunnel exhibits a large interaction, which changes the value and position of the maximum, the superimposed surface settlement curve is a single <italic>V</italic> shape, and the surface settlement curve conforms to the normal distribution. Therefore, for double-track symmetric subway tunnels in shallow buried loess areas, with the decrease of <italic>L</italic>
<sub>j</sub>, the surface settlement superimposed curve is double <italic>V</italic> shape -<italic>W</italic> shape - single <italic>V</italic> shape. From <xref ref-type="fig" rid="F1">Figures 1D&#x2013;F</xref>, the surface settlement mechanism of double-track asymmetric subway tunnels is the same as that of double-track symmetric subway tunnels in shallow buried loess areas. When the surface settlement curve is double <italic>V</italic> shape and <italic>W</italic> shape, the Peck formula can be used to calculate the single-track tunnel settlement curve and superimposed to obtain the final surface settlement curve. When the surface settlement curve is a single <italic>V</italic> shape, the surface settlement is normally distributed.</p>
</sec>
<sec id="s2-3">
<title>2.3 Establishment of the prediction formula</title>
<p>From the results in <xref ref-type="sec" rid="s2-2">Section 2.2</xref>, the surface settlement curve can be superimposed by the Peck formula when it is double <italic>V</italic> shape and <italic>W</italic> shape. The tunnel center is the coordinate origin when calculating the surface settlement, so the surface settlement curve from the STT cannot be superimposed directly (<xref ref-type="bibr" rid="B26">Song et al., 2018</xref>; <xref ref-type="bibr" rid="B39">Yin et al., 2022</xref>). In <xref ref-type="fig" rid="F2">Figure 2A</xref>, the formula for the surface settlement in single-track tunnels is Eqs <xref ref-type="disp-formula" rid="e3">3</xref>, <xref ref-type="disp-formula" rid="e4">4</xref>.<disp-formula id="e3">
<mml:math id="m3">
<mml:mrow>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mrow>
<mml:msub>
<mml:mrow>
</mml:mrow>
<mml:mi mathvariant="normal">L</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>x</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mrow>
<mml:mi>max</mml:mi>
<mml:mo>&#x2061;</mml:mo>
<mml:mi mathvariant="normal">L</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>exp</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mfrac>
<mml:msubsup>
<mml:mi>x</mml:mi>
<mml:mn>1</mml:mn>
<mml:mn>2</mml:mn>
</mml:msubsup>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:msubsup>
<mml:mi>i</mml:mi>
<mml:mi mathvariant="normal">L</mml:mi>
<mml:mn>2</mml:mn>
</mml:msubsup>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(3)</label>
</disp-formula>
</p>
<fig id="F2" position="float">
<label>FIGURE 2</label>
<caption>
<p>Translation of coordinate axis: <bold>(A)</bold> Before translation, <bold>(B)</bold> After translation.</p>
</caption>
<graphic xlink:href="feart-11-1219544-g002.tif"/>
</fig>
<p>Where <italic>S</italic>
<sub>L(x)</sub> is the value at any point of the surface settlement curve of the left tunnel, <italic>S</italic>
<sub>maxl</sub> is the maximum of the surface settlement result of the left tunnel, <italic>x</italic>
<sub>1</sub> is the horizontal distance between the calculated point and the origin in the left tunnel coordinate system, <italic>i</italic>
<sub>L</sub> is the settlement trough width of the left tunnel.<disp-formula id="e4">
<mml:math id="m4">
<mml:mrow>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">R</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>x</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mrow>
<mml:mi>max</mml:mi>
<mml:mo>&#x2061;</mml:mo>
<mml:mi mathvariant="normal">R</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>exp</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mfrac>
<mml:msup>
<mml:msub>
<mml:mi>x</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:msup>
<mml:msub>
<mml:mi>i</mml:mi>
<mml:mi mathvariant="normal">R</mml:mi>
</mml:msub>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(4)</label>
</disp-formula>
</p>
<p>Where <italic>S</italic>
<sub>R(x)</sub> is the value at any point of the surface settlement curve of the right tunnel, <italic>S</italic>
<sub>maxR</sub> is the maximum of the surface settlement result of the right tunnel, <italic>x</italic>
<sub>2</sub> is the horizontal distance between the calculated point and the origin in the right tunnel coordinate system, and <italic>i</italic>
<sub>R</sub> is the settlement trough width of the right tunnel.</p>
<p>In <xref ref-type="fig" rid="F2">Figure 2B</xref>, taking the intersection of the tunnel circle center line and the central axis as the origin, a new coordinate system is established. Translate and superimpose the coordinate axis of Eqs <xref ref-type="disp-formula" rid="e3">3</xref>, <xref ref-type="disp-formula" rid="e4">4</xref> to obtain Eq. <xref ref-type="disp-formula" rid="e5">5</xref>. When both <italic>a</italic> and <italic>b</italic> are 0.5, Eq. <xref ref-type="disp-formula" rid="e5">5</xref> can superimpose the surface settlement curve for the double-track symmetric tunnel, when <italic>a</italic> &#x2260; <italic>b</italic>, Eq. <xref ref-type="disp-formula" rid="e5">5</xref> can superimpose the surface settlement curve for the double-track asymmetric tunnel.<disp-formula id="e5">
<mml:math id="m5">
<mml:mrow>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>x</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mrow>
<mml:mi>max</mml:mi>
<mml:mo>&#x2061;</mml:mo>
<mml:mi mathvariant="normal">L</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>exp</mml:mi>
<mml:mrow>
<mml:mfenced open="[" close="]" separators="|">
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mfrac>
<mml:msup>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>x</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>a</mml:mi>
<mml:mi>L</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:msup>
<mml:msub>
<mml:mi>i</mml:mi>
<mml:mi mathvariant="normal">L</mml:mi>
</mml:msub>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mrow>
<mml:mi>max</mml:mi>
<mml:mo>&#x2061;</mml:mo>
<mml:mi mathvariant="normal">R</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>exp</mml:mi>
<mml:mrow>
<mml:mfenced open="[" close="]" separators="|">
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mfrac>
<mml:msup>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>x</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>b</mml:mi>
<mml:mi>L</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:msup>
<mml:msub>
<mml:mi>i</mml:mi>
<mml:mi mathvariant="normal">R</mml:mi>
</mml:msub>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(5)</label>
</disp-formula>
</p>
<p>Where <italic>S</italic>
<sub>(x)</sub> is the value at any point of the superimposed surface settlement curve, <italic>a</italic> is the <italic>L</italic>
<sub>L</sub>/<italic>L</italic>, <italic>b</italic> is the <italic>L</italic>
<sub>R</sub>/<italic>L</italic>, <italic>L</italic> denotes the horizontal distance between the center of the left and right tunnels, <italic>L</italic>
<sub>L</sub> denotes the horizontal distance between the central axis of the new coordinate system and the central axis of the left tunnel, <italic>L</italic>
<sub>R</sub> denotes the horizontal distance between the central axis of the new coordinate system and the central axis of the right tunnel.</p>
<p>When the superimposed surface settlement curve is a single <italic>V</italic> shape, the maximum of the surface settlement is resulting from the superimposed settlement curve, and the interaction cannot be neglected. Since the surface settlement curve exhibits a normal distribution, it is assumed that the superimposed surface settlement curve is a modified form of Peck formula, and Eq. <xref ref-type="disp-formula" rid="e5">5</xref> is rewritten as Eq. <xref ref-type="disp-formula" rid="e6">6</xref>.<disp-formula id="e6">
<mml:math id="m6">
<mml:mrow>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>x</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msubsup>
<mml:mi>S</mml:mi>
<mml:mi>max</mml:mi>
<mml:mo>&#x2032;</mml:mo>
</mml:msubsup>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>exp</mml:mi>
<mml:mrow>
<mml:mfenced open="[" close="]" separators="|">
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mfrac>
<mml:msup>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>x</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:msup>
<mml:mi>l</mml:mi>
<mml:mo>&#x2032;</mml:mo>
</mml:msup>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:msup>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msup>
<mml:mi>i</mml:mi>
<mml:mo>&#x2032;</mml:mo>
</mml:msup>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(6)</label>
</disp-formula>
</p>
<p>Where <inline-formula id="inf1">
<mml:math id="m7">
<mml:mrow>
<mml:msubsup>
<mml:mi>S</mml:mi>
<mml:mi>max</mml:mi>
<mml:mo>&#x2032;</mml:mo>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> is the modified maximum of the surface settlement curve, <italic>x</italic> is the horizontal distance between the calculation point and the origin, <inline-formula id="inf2">
<mml:math id="m8">
<mml:mrow>
<mml:msup>
<mml:mi>l</mml:mi>
<mml:mo>&#x2032;</mml:mo>
</mml:msup>
</mml:mrow>
</mml:math>
</inline-formula> is the horizontal distance between the maximum position and the origin, <inline-formula id="inf3">
<mml:math id="m9">
<mml:mrow>
<mml:msup>
<mml:mi>i</mml:mi>
<mml:mo>&#x2032;</mml:mo>
</mml:msup>
</mml:mrow>
</mml:math>
</inline-formula> is the modified settlement trough width.</p>
<p>When the surface settlement superimposed curve is a single V shape, the value and position of the maximum will be changed by the tunnel interaction. Specific to DTT, the area difference and the net distance are important factors affecting the surface settlement, so <italic>m</italic> is defined as the net distance of the tunnel/small section tunnel diameter, and <italic>n</italic> is the area of the large section tunnel/small section tunnel area. Using <inline-formula id="inf4">
<mml:math id="m10">
<mml:mrow>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>m</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>n</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> to correct the maximum value and <inline-formula id="inf5">
<mml:math id="m11">
<mml:mrow>
<mml:msub>
<mml:mi>w</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>m</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>n</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> to correct the maximum position, the final surface settlement expression is Eq. <xref ref-type="disp-formula" rid="e7">(7)</xref>.<disp-formula id="e7">
<mml:math id="m12">
<mml:mrow>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>x</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>m</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>n</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:msub>
<mml:msubsup>
<mml:mi>S</mml:mi>
<mml:mi>max</mml:mi>
<mml:mo>&#x2032;</mml:mo>
</mml:msubsup>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>exp</mml:mi>
<mml:mrow>
<mml:mfenced open="[" close="]" separators="|">
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mfrac>
<mml:msup>
<mml:mrow>
<mml:mfenced open="[" close="]" separators="|">
<mml:mrow>
<mml:mi>x</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>w</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>m</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>n</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msup>
<mml:mi>l</mml:mi>
<mml:mo>&#x2032;</mml:mo>
</mml:msup>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:msup>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msup>
<mml:mi>i</mml:mi>
<mml:mo>&#x2032;</mml:mo>
</mml:msup>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(7)</label>
</disp-formula>
</p>
<p>Where <inline-formula id="inf6">
<mml:math id="m13">
<mml:mrow>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>m</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>n</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the modified expression of the maximum value, <inline-formula id="inf7">
<mml:math id="m14">
<mml:mrow>
<mml:msub>
<mml:mi>w</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>m</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>n</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the modified expression of the maximum position.</p>
<p>It is necessary to obtain the surface settlement value resulting from the STT excavation before the surface settlement curve is superimposed (<xref ref-type="bibr" rid="B22">Ou et al., 2019</xref>). According to <xref ref-type="fig" rid="F3">Figure 3</xref>, the surface settlement during the excavation of DTT falls into five stages: the initial stage of excavation (I)&#x2014;the stage of the first excavation near the characteristic surface (II)&#x2014;the stage of the first excavation away from the characteristic surface, the second excavation near the characteristic surface (III)&#x2014;the stage of the second excavation near the characteristic surface (&#x2163;)&#x2014;the stage of the second excavation away from the characteristic surface (&#x2164;). Since the second tunnel is not yet excavated upon the excavation of the first tunnel, the surface settlement value at the upper extent of the first tunnel at the end of Stage II is taken as the value of surface settlement resulting from the separate excavation regarding the first tunnel. As the excavation of the first tunnel affected the surface settlement in the second tunnel&#x2019;s upper part, the variation in the surface settlement in the part at the end of Stage IV and Stage II is taken as the value of surface settlement resulting from the separate excavation regarding the second tunnel.</p>
<fig id="F3" position="float">
<label>FIGURE 3</label>
<caption>
<p>Surface settlement in the upper part of the tunnel.</p>
</caption>
<graphic xlink:href="feart-11-1219544-g003.tif"/>
</fig>
</sec>
</sec>
<sec id="s3">
<title>3 Parameters analysis of prediction formula</title>
<sec id="s3-1">
<title>3.1 Scheme of numerical simulation</title>
<p>In order to determine <italic>f</italic>
<sub>(m,n)</sub> and <italic>w</italic>
<sub>(m,n)</sub> in Eq. <xref ref-type="disp-formula" rid="e7">(7)</xref>, the paper compares the surface settlement curve after superposition and completion of excavation through the numerical test, and establishes the equation for increment and position offset of the surface settlement maximum at different <italic>m</italic> and <italic>n</italic>, the control variables method is used for the 16 tests, with only <italic>m</italic> and <italic>n</italic> changed during the simulation. <italic>m</italic> is taken as 1, 2, 3, and 4, and <italic>n</italic> is taken as 1, 2, 3, and 4 for each value of <italic>m</italic>. In order to control the surface settlement (<xref ref-type="bibr" rid="B31">Wang, 2017</xref>; <xref ref-type="bibr" rid="B19">Lu et al., 2019</xref>; <xref ref-type="bibr" rid="B44">Zheng and Wu, 2019</xref>), the large section tunnel was excavated first in the numerical test, and the results of the numerical test are valid when the superimposed surface settlement curve is a single <italic>V</italic> shape.</p>
</sec>
<sec id="s3-2">
<title>3.2 Finite element mesh and model parameters</title>
<p>The paper used MIDAS/GTS software to establish a 3D numerical model, the majority of single-track subway tunnels are around 6&#xa0;m in diameter (<xref ref-type="bibr" rid="B38">Yang, 2017</xref>), so the diameter of the small cross-section tunnel is taken as 6&#xa0;m. <xref ref-type="fig" rid="F4">Figure 4</xref> illustrates the model (with <italic>m</italic> is 2 and <italic>n</italic> is 2), the large section tunnel has a diameter of 8.80&#xa0;m. Given to which degree the tunnel excavation influences the surrounding rock (<xref ref-type="bibr" rid="B38">Yang, 2017</xref>; <xref ref-type="bibr" rid="B19">Lu et al., 2019</xref>; <xref ref-type="bibr" rid="B40">Yu et al., 2019</xref>), the lengths of the X, Y, and Z direction are 96&#xa0;m, 120&#xa0;m, and 60&#xa0;m respectively, where the Y direction stands for the tunnel excavation direction. The small section tunnel has a burial depth of 20&#xa0;m. The distance between the excavation faces is 36&#xa0;m and the step method is adopted for tunnel excavation. The surrounding boundaries of the model are set as horizontal constraints, and the bottom boundary is set as horizontal and vertical constraints.</p>
<fig id="F4" position="float">
<label>FIGURE 4</label>
<caption>
<p>Three-Dimensional model.</p>
</caption>
<graphic xlink:href="feart-11-1219544-g004.tif"/>
</fig>
<p>For eliminating the way the model boundary influences the simulation results (<xref ref-type="bibr" rid="B17">Li et al., 2020</xref>; <xref ref-type="bibr" rid="B33">Wei et al., 2023</xref>), 100&#xa0;m away from the initial excavation face is selected as the characteristic face. The soil parameters, the burial depth of the small section tunnel, the construction method, the support measures, and the distance between the excavation faces are kept consistent during the 16 tests.</p>
<p>The surrounding rock adopts the Mohr-Coulomb model and is simulated with 3 D solid element (<xref ref-type="bibr" rid="B34">Xiang et al., 2020</xref>; <xref ref-type="bibr" rid="B42">Zhang et al., 2021</xref>; <xref ref-type="bibr" rid="B14">Hu et al., 2023</xref>), the initial support adopts elastic model and is simulated with 2 D plate element. <xref ref-type="table" rid="T1">Table 1</xref> lists the material parameters.</p>
<table-wrap id="T1" position="float">
<label>TABLE 1</label>
<caption>
<p>Model material parameters.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="center">Type</th>
<th align="center">Weight (KN/m<sup>3</sup>)</th>
<th align="center">Modulus (KPa)</th>
<th align="center">Poisson ratio (&#x3bc;)</th>
<th align="center">Friction angle (&#xb0;)</th>
<th align="center">Cohesive force (KPa)</th>
<th align="center">Thickness (m)</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="center">Loss</td>
<td align="center">18.50</td>
<td align="center">5.00&#xd7;10<sup>4</sup>
</td>
<td align="center">0.30</td>
<td align="center">22.50</td>
<td align="center">25.00</td>
<td align="center">&#x2014;</td>
</tr>
<tr>
<td align="center">Initial support</td>
<td align="center">23.00</td>
<td align="center">2.10&#xd7;10<sup>7</sup>
</td>
<td align="center">0.20</td>
<td align="center">&#x2014;</td>
<td align="center">&#x2014;</td>
<td align="center">0.30</td>
</tr>
<tr>
<td align="center">Temporary support</td>
<td align="center">25.00</td>
<td align="center">2.70&#xd7;10<sup>7</sup>
</td>
<td align="center">0.20</td>
<td align="center">&#x2014;</td>
<td align="center">&#x2014;</td>
<td align="center">0.50</td>
</tr>
</tbody>
</table>
</table-wrap>
</sec>
<sec id="s3-3">
<title>3.3 Simulation process</title>
<p>The excavation footage in the simulation process is 2&#xa0;m, the upper step of the right tunnel is excavated first, then applying the initial support for the upper step. The distance between the upper and lower steps is 6&#xa0;m, and then excavating the lower step and applying the initial support for the lower step, the left tunnel undergoes the same excavation process as the right tunnel, and the distance between the upper steps of the left and right tunnels is always 36&#xa0;m.</p>
</sec>
<sec id="s3-4">
<title>3.4 Analysis of simulation results</title>
<p>
<xref ref-type="fig" rid="F5">Figure 5</xref> shows the surface settlement change curve in the process of excavation (with <italic>m</italic> is 2 and <italic>n</italic> is 1). From <xref ref-type="fig" rid="F5">Figure 5</xref>, there is a symmetric surface settlement curve regarding STT excavation at the tunnel center. The surface settlement maximum after superposition is 12.91&#xa0;mm, located 4.4&#xa0;m to the right of the central axis. After excavation, the maximum value is 21.58 mm, which is 3.8&#xa0;m on the right of the central axis. The increment is 8.67&#xa0;mm and the position offset is 0.5&#xa0;m of maximum on the surface settlement, which shows the interaction of tunnels has a greater influence on the value of the maximum and a smaller influence on the position of the maximum.</p>
<fig id="F5" position="float">
<label>FIGURE 5</label>
<caption>
<p>The variation on the surface settlement at the characteristic face.</p>
</caption>
<graphic xlink:href="feart-11-1219544-g005.tif"/>
</fig>
<p>The simulation is continued according to the method in <xref ref-type="fig" rid="F5">Figure 5</xref>, the final settlement curve and the superimposed settlement curve are compared to obtain all the maximum increments and position offsets, the simulation results are shown in <xref ref-type="table" rid="T2">Tables 2</xref>, <xref ref-type="table" rid="T3">3</xref>.</p>
<table-wrap id="T2" position="float">
<label>TABLE 2</label>
<caption>
<p>Value increments of the maximum (mm).</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="center">
<italic>n</italic>
</th>
<th rowspan="2" align="center">1</th>
<th rowspan="2" align="center">2</th>
<th rowspan="2" align="center">3</th>
<th rowspan="2" align="center">4</th>
</tr>
<tr>
<th align="center">
<italic>m</italic>
</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="center">1</td>
<td align="center">4.44</td>
<td align="center">8.67</td>
<td align="center">9.96</td>
<td align="center">14.28</td>
</tr>
<tr>
<td align="center">2</td>
<td align="center">3.78</td>
<td align="center">6.51</td>
<td align="center">8.82</td>
<td align="center">12.39</td>
</tr>
<tr>
<td align="center">3</td>
<td align="center">1.50</td>
<td align="center">5.25</td>
<td align="center">7.86</td>
<td align="center">10.98</td>
</tr>
<tr>
<td align="center">4</td>
<td align="center">0.03</td>
<td align="center">4.98</td>
<td align="center">7.08</td>
<td align="center">10.38</td>
</tr>
</tbody>
</table>
</table-wrap>
<table-wrap id="T3" position="float">
<label>TABLE 3</label>
<caption>
<p>Position offsets of the maximum (m).</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="center">
<italic>n</italic>
</th>
<th rowspan="2" align="center">1</th>
<th rowspan="2" align="center">2</th>
<th rowspan="2" align="center">3</th>
<th rowspan="2" align="center">4</th>
</tr>
<tr>
<th align="center">
<italic>m</italic>
</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="center">1</td>
<td align="center">&#x2212;0.7</td>
<td align="center">&#x2212;0.6</td>
<td align="center">&#x2212;0.3</td>
<td align="center">0</td>
</tr>
<tr>
<td align="center">2</td>
<td align="center">&#x2212;0.5</td>
<td align="center">&#x2212;0.4</td>
<td align="center">0</td>
<td align="center">0</td>
</tr>
<tr>
<td align="center">3</td>
<td align="center">&#x2212;0.3</td>
<td align="center">&#x2212;0.2</td>
<td align="center">0</td>
<td align="center">0</td>
</tr>
<tr>
<td align="center">4</td>
<td align="center">&#x2212;0.2</td>
<td align="center">&#x2212;0.1</td>
<td align="center">0</td>
<td align="center">0</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>
<xref ref-type="fig" rid="F6">Figure 6</xref> shows the results of data fitting in <xref ref-type="table" rid="T2">Tables 2</xref>, <xref ref-type="table" rid="T3">3</xref>, the value increments and position offsets of surface settlement maximum regarding different <italic>m</italic> and <italic>n</italic> are fitted. <xref ref-type="fig" rid="F6">Figure 6A</xref> shows the surface after fitting the 16 maximum value increments, the value increment of the maximum increases with the increase in <italic>n</italic> for a certain <italic>m</italic>, as <italic>n</italic> becomes larger, the effect of large section tunnel on small section tunnel becomes more pronounced. For a given <italic>n</italic>, the weakened interaction between the left and right tunnels, the increment of the maximum decreases with the increase in <italic>m</italic>. When <italic>n</italic> is fixed, the largest change in the maximum increment is 14.28&#xa0;mm, when <italic>m</italic> is fixed, the Largest change in the maximum increment is 10.38&#xa0;mm, which shows that the maximum increment is greatly influenced by <italic>m</italic> and <italic>n</italic>. The increments are fitted and the fitted expression is <italic>f</italic>
<sub>(m,n)</sub>&#x3d;&#x2212;0.1706<italic>m</italic>
<sup>2</sup>&#x2212;0.0844<italic>n</italic>
<sup>2</sup>&#x2212;0.0993 <italic>m n</italic>&#x2b;2.3651 <italic>m</italic>&#x2212;3.2524 <italic>n</italic>&#x2212;3.8213, the fitting degree is 0.9721 and the fit is good.</p>
<fig id="F6" position="float">
<label>FIGURE 6</label>
<caption>
<p>Fitting of the surface settlement maximum: <bold>(A)</bold> Increment of the value, <bold>(B)</bold> Offset of the position.</p>
</caption>
<graphic xlink:href="feart-11-1219544-g006.tif"/>
</fig>
<p>
<xref ref-type="fig" rid="F6">Figure 6B</xref> shows the surface after fitting the 16 maximum position offsets. From <xref ref-type="fig" rid="F6">Figure 6B</xref>, when <italic>n</italic> is 1 and 2, the maximum position offset can be influenced by <italic>m</italic> in the same way as the maximum increment change law, when <italic>n</italic> belongs to 3<italic>&#x2013;</italic>4 and <italic>m</italic> belongs to 1<italic>&#x2013;</italic>4, the position offset is 0. When <italic>n</italic> is 1 and <italic>m</italic> is 1, the position offset is 0.7&#xa0;m, which indicates that the maximum position offset is slightly influenced by <italic>m</italic> and <italic>n</italic>. The position offsets of the maximum are fitted and the fitted expression is <italic>w</italic>
<sub>(m,n)</sub>&#x3d;&#x2212;0.0313<italic>m</italic>
<sup>2</sup>&#x2212;0.0063<italic>n</italic>
<sup>2</sup>&#x2212;0.0590 <italic>m n</italic>&#x2b;0.4113 m&#x2b;0.3313 <italic>n</italic>&#x2212;1.4125, the fitting degree is 0.9120 and the fit is good.</p>
</sec>
</sec>
<sec id="s4">
<title>4 Verification of engineering</title>
<sec id="s4-1">
<title>4.1 Double-track symmetric subway tunnel</title>
<p>A double-track symmetric subway tunnel of Xi&#x2019;an Subway Line 3 is located underground on Jinhua North Road, Xi&#x2019;an, The mileage of the left tunnel is ZDK31&#x2b;443.908 to ZDK31&#x2b;503.908, and the mileage of the right tunnel is YDK31&#x2b;443.965 to YDK31&#x2b;503.965. The dimension is 8.08&#xa0;m <inline-formula id="inf8">
<mml:math id="m15">
<mml:mrow>
<mml:mo>&#x2219;</mml:mo>
</mml:mrow>
</mml:math>
</inline-formula> 8.67&#xa0;m, the burial depth is 9.8&#xa0;m and the net distance is 2.05&#x2013;4.75&#xa0;m. The net distance is simplified to 2.05&#xa0;m (<xref ref-type="bibr" rid="B32">Wei, 2018</xref>), where <italic>m</italic> is 0.25 and <italic>n</italic> is 1. The geological conditions are miscellaneous fill, plain fill, new loess, and old loess from top to bottom. <xref ref-type="fig" rid="F7">Figure 7</xref> displays the tunnel cross-section.</p>
<fig id="F7" position="float">
<label>FIGURE 7</label>
<caption>
<p>Cross-Section of the double-track symmetric tunnel.</p>
</caption>
<graphic xlink:href="feart-11-1219544-g007.tif"/>
</fig>
<p>The tunnel radius, internal friction angle, and tunnel burial depth are substituted into the Peck formula. The internal friction angle of various geology is calculated by the weighted average method (<xref ref-type="bibr" rid="B13">Hu et al., 2012</xref>), and the surface loss rate is taken as 4% for acquiring the surface settlement curve regarding the STT excavation (<xref ref-type="bibr" rid="B10">Han et al., 2007</xref>). The two surface settlement curves are superimposed and fitted, then the value of the maximum is corrected by using <italic>f</italic>
<sub>(0.25,1)</sub> &#x3d; &#x2212;12.41&#xa0;mm, and the position offset of the maximum is corrected by using <italic>w</italic>
<sub>(0.25,1)</sub> &#x3d; &#x2212;0.53&#xa0;m.</p>
<p>The distance between adjacent monitoring sections along the excavation direction is 5&#xa0;m, and the lateral arrangement of the monitoring sections is shown in <xref ref-type="fig" rid="F8">Figure 8</xref>, the distance between the monitoring points of the main influence area is 5&#xa0;m, and the distance between the monitoring points of the secondary influence area is 3.35&#xa0;m.</p>
<fig id="F8" position="float">
<label>FIGURE 8</label>
<caption>
<p>Monitoring points layout of Double-track symmetric subway tunnel.</p>
</caption>
<graphic xlink:href="feart-11-1219544-g008.tif"/>
</fig>
<p>Comparing the superposed surface settlement curve with the monitoring data (DK31 &#x2b; 480), the results displayed in <xref ref-type="fig" rid="F9">Figure 9</xref>, the surface settlement maximum after superposition is 11.94&#xa0;mm, which is located at 2.05&#xa0;m to the right of the central axis. The corrected surface settlement curve and the monitoring data are both single <italic>V</italic> shape, and the largest values are 24.35&#xa0;mm and 21.35&#xa0;mm respectively, with a difference of 3&#xa0;mm. The maximum is located at 1.52&#xa0;m and 3&#xa0;m on the right side of the central axis, and the difference is 1.48&#xa0;m, which proves the reliability exhibited by the formula in predicting the surface settlement regarding double-track symmetric subway tunnels in shallow buried loess areas.</p>
<fig id="F9" position="float">
<label>FIGURE 9</label>
<caption>
<p>Surface settlement of the shallow-buried symmetric subway tunnel.</p>
</caption>
<graphic xlink:href="feart-11-1219544-g009.tif"/>
</fig>
</sec>
<sec id="s4-2">
<title>4.2 Double-track asymmetric subway tunnel</title>
<p>Xi&#x2019;an Subway Line 4 Yannan Fourth Road Station &#x223c; Datang Furongyuan Fourth Road Station is mainly along the Furong West Road in Yanta District. A 246-meter double-track asymmetric subway tunnel in shallow buried loess areas is built at the southern end of Datang Furongyuan. The mileage is CK9 &#x2b; 358.100 to CK9 &#x2b; 619.325, the geological conditions are miscellaneous fill, new loess, and old loess from top to bottom. The sizes of the right and left tunnels are 12.2&#xa0;m <inline-formula id="inf9">
<mml:math id="m16">
<mml:mrow>
<mml:mo>&#x2219;</mml:mo>
</mml:mrow>
</mml:math>
</inline-formula> 10.0&#xa0;m and 6.4&#xa0;m <inline-formula id="inf10">
<mml:math id="m17">
<mml:mrow>
<mml:mo>&#x2219;</mml:mo>
</mml:mrow>
</mml:math>
</inline-formula> 6.7&#xa0;m respectively, and the maximum buried depth is 22.1&#xa0;m. The right tunnel is 4.1&#xa0;m away from the left tunnel, then <italic>m</italic> is 0.68, and <italic>n</italic> is 2.72. <xref ref-type="fig" rid="F10">Figure 10</xref> displays the tunnel cross-section.</p>
<fig id="F10" position="float">
<label>FIGURE 10</label>
<caption>
<p>Cross-Section of the double-track asymmetric tunnel.</p>
</caption>
<graphic xlink:href="feart-11-1219544-g010.tif"/>
</fig>
<p>The surface settlement curve regarding the STT can be obtained according to the method in <xref ref-type="sec" rid="s4-1">Section 4.1</xref>. After fitting the superimposed settlement curve, <italic>f</italic>
<sub>(0.68,2.72)</sub>&#x3d; &#x2212;10.70&#xa0;mm is used to correct the fitted maximum value, and <italic>w</italic>
<sub>(0.68,2.72)</sub>&#x3d; &#x2212;0.40&#xa0;m is used to correct the fitted maximum position. Surface settlement monitoring points are arranged as shown in <xref ref-type="fig" rid="F11">Figure 11</xref>, the distance between adjacent monitoring sections along the excavation direction is 10&#xa0;m. Most of the monitoring sections are set up with 2 points and a few with 11 points. In order to fully reflect the final surface settlement caused by the tunnel excavation, the monitoring sections CK9&#x2b;430 and CK9&#x2b;530 are selected for analysis.</p>
<fig id="F11" position="float">
<label>FIGURE 11</label>
<caption>
<p>Monitoring points layout of double-track asymmetric subway tunnel.</p>
</caption>
<graphic xlink:href="feart-11-1219544-g011.tif"/>
</fig>
<p>Compare the final surface settlement curve and the monitoring data to get <xref ref-type="fig" rid="F12">Figure 12</xref>. In <xref ref-type="fig" rid="F12">Figure 12</xref>, the maximum value of the superimposed surface settlement is 14.52&#xa0;mm, located at 5.19&#xa0;m to the right of the central axis. The corrected surface settlement prediction curve and the monitoring data are single <italic>V</italic> shape, and the largest corrected surface settlement is 25.22&#xa0;mm, with an increment of 73.69%, the maximum is located at 4.73&#xa0;m to the right of the tunnel central axis. The largest surface settlements are 20.66&#xa0;mm and 20.23&#xa0;mm in the two monitoring sections respectively, and the locations are 2.05&#xa0;m and 5.05&#xa0;m. The corrected surface settlement prediction curve is slightly different from the maximum monitoring data, and the maximum location is in the middle of the maximum locations in the two monitoring sections, which proves the reliability of the prediction formula for predicting surface settlement in double-track asymmetric subway tunnels in shallow buried loess areas.</p>
<fig id="F12" position="float">
<label>FIGURE 12</label>
<caption>
<p>Surface settlement of the shallow-buried asymmetrical subway tunnel.</p>
</caption>
<graphic xlink:href="feart-11-1219544-g012.tif"/>
</fig>
</sec>
</sec>
<sec sec-type="discussion" id="s5">
<title>5 Discussion</title>
<p>Based on the International Tunneling Association (ITA) classification criteria for the tunnel section area, when the left tunnel has a diameter of 6&#xa0;m and <italic>n</italic> of 4, the right tunnel is an extra large section tunnel. Thus, in this section, only <italic>m</italic> is analyzed and studied to determine the applicable range of <italic>m</italic> when <italic>n</italic> is 1&#x2013;4. From the discussion in <xref ref-type="sec" rid="s2-3">Section 2.3</xref>, the prediction formula is applicable to double-track subway tunnels in shallow buried loess areas where the superimposed surface settlement curve is a single <italic>V</italic> shape, Therefore, the applicable standard of the formula is that the surface settlement superimposed curve is a single <italic>V</italic> shape.</p>
<p>
<xref ref-type="fig" rid="F13">Figure 13</xref> shows the variation of the superimposed surface settlement curve with increasing <italic>m</italic> for different <italic>n</italic>. For double-track symmetric subway tunnels in shallow buried loess areas, it can be seen from <xref ref-type="fig" rid="F13">Figure 13A</xref> that the superimposed curve is a single <italic>V</italic> shape when <italic>m</italic> is 1, 2, 3, 4, <italic>W</italic> shape when <italic>m</italic> is 5, and double <italic>V</italic> shape when <italic>m</italic> is 6. This indicates that for double-track symmetric subway tunnels in shallow buried loess areas, the applicable range of <italic>m</italic> is 1&#x2013;4.</p>
<fig id="F13" position="float">
<label>FIGURE 13</label>
<caption>
<p>Superimposed curves of surface settlement: <bold>(A)</bold> <italic>n</italic> &#x3d; 1, <bold>(B)</bold> <italic>n</italic> &#x3d; 2, <bold>(C)</bold> <italic>n</italic> &#x3d; 3, <bold>(D)</bold> <italic>n</italic> &#x3d; 4.</p>
</caption>
<graphic xlink:href="feart-11-1219544-g013.tif"/>
</fig>
<p>For double-track asymmetric subway tunnels in shallow buried loess areas, <xref ref-type="fig" rid="F13">Figures 13B&#x2013;D</xref> show that when <italic>n</italic> is 2, 3 and <italic>m</italic> is 1,2,3,4, the superimposed curve is a single <italic>V</italic> shape, the superimposed curve is a <italic>W</italic> shape when <italic>m</italic> is 5, and the superimposed curve is a double <italic>V</italic> shape when <italic>m</italic> is 6. When <italic>n</italic> is 4, <italic>m</italic> is 1&#x2013;5, the superimposed curve is single a <italic>V</italic> shape, when <italic>m</italic> is 6, the superimposed curve is a <italic>W</italic> shape, when <italic>m</italic> is 7, the superimposed curve is a double <italic>V</italic> shape. The application range of the formula is that when <italic>n</italic> is 2&#x2013;3, <italic>m</italic> is 1&#x2013;4, and when <italic>n</italic> is 4, <italic>m</italic> is 1&#x2013;5.</p>
<p>The results of the analysis in <xref ref-type="fig" rid="F13">Figure 13</xref> show that for a certain <italic>n</italic>, as <italic>m</italic> increases, the surface settlement superimposed curve changes in a single <italic>V</italic> shape&#x2014;<italic>W</italic> shape&#x2014;double <italic>V</italic> shape trend, this indicates the surface settlement mechanism in loess areas is the same for double-track symmetric subway tunnels and double-track asymmetric subway tunnels, which validates the conclusions in <xref ref-type="sec" rid="s2-2">Section 2.2</xref>.</p>
<p>The prediction formula proposed in the paper is simple in form and easy to take values, and can predict the surface settlement for shallow buried double-track tunnels in loess areas, which is more practical. However, the diameter, the burial depth, and the section shape of the small cross-section tunnel are fixed, which limits the applicability of the prediction formula. Therefore, in the future, the tunnel diameter, the burial depth, and the cross-section shape should be considered to increase the applicability of the prediction formula.</p>
</sec>
<sec sec-type="conclusion" id="s6">
<title>6 Conclusion</title>
<p>The paper examines the surface settlement mechanism regarding double-track subway tunnels in shallow buried loess areas, and divides the superimposed curve of surface settlement into double <italic>V</italic> shape, <italic>W</italic> shape, and single <italic>V</italic> shape. Then the prediction formula is built and relevant parameters are determined by using numerical tests. Finally, the prediction formula is applied to shallow buried double-track subway tunnels project in shallow loess areas. Data monitoring serves for verifying the prediction formula&#x2019;s reliability and the application range of the prediction formula is discussed. The study has four main four conclusions.<list list-type="simple">
<list-item>
<p>(1) The surface settlement mechanism in loess areas is the same for double-track symmetric subway tunnels and double-track asymmetric subway tunnels. When the superimposed curve of the surface settlement regarding double-track subway tunnels in shallow buried loess areas is double <italic>V</italic> shape and <italic>W</italic> shape, the interaction between the left and right tunnels can be ignored, and the Peck formula serves for calculating the surface settlement curve resulted from the single-track tunnel excavation, and then superimposing the surface settlement curve to obtain the final surface settlement curve.</p>
</list-item>
<list-item>
<p>(2) When the superimposed surface settlement curve of double-track subway tunnels in shallow buried loess areas is a single <italic>V</italic> shape, the net distance ratio and the area ratio are defined. According to the Peck formula, the value and position of the maximum are corrected by considering the interaction of the tunnel, and the final surface settlement expression of double-track subway tunnels in shallow buried loess areas is obtained. Through 16 numerical tests, the parameters on the surface settlement regarding the net distance ratio and the area ratio in the expression are determined.</p>
</list-item>
<list-item>
<p>(3) The formula is applied in symmetric and asymmetric subway tunnels in shallow buried loess areas, which verifies the formula&#x2019;s rationality and the accuracy of the parameters. Then the scope of application of the formula is discussed. The research results show that when the area ratio is 1 and the net distance ratio is 1&#x2013;3, the prediction formula is suitable for double-track symmetric subway tunnels in shallow buried loess areas. When the area ratio is 2&#x2013;3, the net distance ratio is 1&#x2013;4, the area ratio is 4 and the net distance ratio is 1&#x2013;5, the prediction formula is suitable for double-track asymmetric subway tunnels in shallow buried loess areas.</p>
</list-item>
<list-item>
<p>(4) The prediction formula served for predicting the surface settlement of double-track subway tunnels in shallow buried loess areas, which can reduce construction risks and ensure the safety of buildings above ground. However, the diameter, the burial depth, and the cross-section shape of the small-section tunnel are fixed. In actual engineering, there has a major difference in terms of the diameter, the burial depth, and the cross-section shape of the tunnel<italic>,</italic> the predictions in the paper have some limitations. Therefore, the variation of the diameter, the burial depth, and the cross-section shape of the tunnel should be considered, which will be the focus of the prediction modification.</p>
</list-item>
</list>
</p>
</sec>
</body>
<back>
<sec sec-type="data-availability" id="s7">
<title>Data availability statement</title>
<p>The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author.</p>
</sec>
<sec id="s8">
<title>Author contributions</title>
<p>SW conceived and designed the methods, completed the numerical model. ZS provided the guide of monitoring and theoretical research for the study, and economically supported the project. XT and YS completed data analysis, and the writing of paper manuscript. YZ provided the monitoring materials and completed the data collection. All authors contributed to the article and approved the submitted version.</p>
</sec>
<sec id="s9">
<title>Funding</title>
<p>This study was supported by the Provincial Natural Science Foundation of Shaanxi (No. 2021JM-373) and the Science and Technology Innovation Team of Shaanxi Innovation Capability Support Plan (No. 2020TD005).</p>
</sec>
<sec sec-type="COI-statement" id="s10">
<title>Conflict of interest</title>
<p>Author YS was employed by China Railway 20th Bureau Group First Engineering Co., Ltd.</p>
<p>The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.</p>
</sec>
<sec sec-type="disclaimer" id="s11">
<title>Publisher&#x2019;s note</title>
<p>All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.</p>
</sec>
<ref-list>
<title>References</title>
<ref id="B1">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Cao</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>Hu</surname>
<given-names>X. B.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>E. L.</given-names>
</name>
<name>
<surname>Guo</surname>
<given-names>H. Y.</given-names>
</name>
</person-group> (<year>2023</year>). <article-title>Mechanical properties and numerical simulation analysis on large deformation of Jiangluling Carbonaceous Shale Tunnel</article-title>. <source>Front. Earth Sci.</source> <volume>2023</volume> (<issue>11</issue>), <fpage>1125410</fpage>. <pub-id pub-id-type="doi">10.3389/feart.2023.1125410</pub-id>
</citation>
</ref>
<ref id="B2">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Chapman</surname>
<given-names>D. N.</given-names>
</name>
<name>
<surname>Rogers</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Hunt</surname>
<given-names>D.</given-names>
</name>
</person-group> (<year>2004</year>). <article-title>Predicting the settlements above twin tunnels constructed in soft ground</article-title>. <source>Tunn. Undergr. Sp. Tech.</source> <volume>19</volume> (<issue>4-5</issue>), <fpage>78</fpage>&#x2013;<lpage>380</lpage>. <pub-id pub-id-type="doi">10.1016/j.tust.2004.02.008</pub-id>
</citation>
</ref>
<ref id="B3">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Chen</surname>
<given-names>C. L.</given-names>
</name>
<name>
<surname>Zhao</surname>
<given-names>C. L.</given-names>
</name>
<name>
<surname>Wei</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Ding</surname>
<given-names>Z.</given-names>
</name>
</person-group> (<year>2014</year>). <article-title>Prediction of soil settlement induced by double-line shield tunnel based on Peck formula</article-title>. <source>Rock Soil Mech.</source> <volume>35</volume> (<issue>08</issue>), <fpage>2212</fpage>&#x2013;<lpage>2218</lpage>. <pub-id pub-id-type="doi">10.16285/j.rsm.2014.08.024</pub-id>
</citation>
</ref>
<ref id="B4">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Chen</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Tan</surname>
<given-names>L. X.</given-names>
</name>
<name>
<surname>Xiao</surname>
<given-names>N.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>K. W.</given-names>
</name>
<name>
<surname>Jia</surname>
<given-names>P. J.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>W.</given-names>
</name>
</person-group> (<year>2023</year>). <article-title>The hydro-mechanical characteristics and micro-structure of loess enhanced by microbially induced carbonate precipitation</article-title>. <source>Geomech. ENERGY ENVIR</source> <volume>34</volume>, <fpage>100469</fpage>. <pub-id pub-id-type="doi">10.1016/j.gete.2023.100469</pub-id>
</citation>
</ref>
<ref id="B5">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Ding</surname>
<given-names>Z. Y.</given-names>
</name>
</person-group> (<year>2019</year>). <source>Research on ground surface settlement and control induced by double-line tunnels construction</source>. <publisher-loc>Xi&#x2019;an</publisher-loc>: <publisher-name>Xi&#x2019;an University of Science and Technology</publisher-name>.</citation>
</ref>
<ref id="B6">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Fan</surname>
<given-names>S. Y.</given-names>
</name>
<name>
<surname>Song</surname>
<given-names>Z. P.</given-names>
</name>
<name>
<surname>Tian</surname>
<given-names>X. X.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>K. M.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>Y. W.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Tunnel deformation and stress response under the bilateral foundation pit construction - a case study</article-title>. <source>Arch. Civ. Mech. Eng.</source> <volume>21</volume> (<issue>3</issue>), <fpage>1</fpage>&#x2013;<lpage>19</lpage>. <pub-id pub-id-type="doi">10.1007/s43452-021-00259-7</pub-id>
</citation>
</ref>
<ref id="B7">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Feng</surname>
<given-names>H. J.</given-names>
</name>
<name>
<surname>Yu</surname>
<given-names>R. G.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Numerical analysis of influence of double - line tunnel shield tunneling on ground surface settlement</article-title>. <source>J. Railw. Eng. Soc.</source> <volume>36</volume> (<issue>3</issue>), <fpage>78</fpage>&#x2013;<lpage>83</lpage>. <pub-id pub-id-type="doi">10.3969/j.issn.1006-2106.2019.03.014</pub-id>
</citation>
</ref>
<ref id="B8">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Han</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Su</surname>
<given-names>S, Q.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Prediction of surface settlement caused by construction of shallow-buried and under-excavated subway tunnel in loess stratum</article-title>. <source>Urb Rd Bdg F C</source> <volume>5</volume> (<issue>11</issue>), <fpage>168</fpage>&#x2013;<lpage>170&#x2b;177&#x2b;21</lpage>. <pub-id pub-id-type="doi">10.16799/j.cnki.csdqyfh.2022.11.045</pub-id>
</citation>
</ref>
<ref id="B9">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Han</surname>
<given-names>R. M.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>Y. Z.</given-names>
</name>
</person-group> (<year>2012</year>). <article-title>Analysis and prediction for ground subsidence of subway tunnel built in saturated soft loess</article-title>. <source>Urb Rpd Rail Transit</source> <volume>25</volume> (<issue>2</issue>), <fpage>74</fpage>&#x2013;<lpage>78</lpage>. <pub-id pub-id-type="doi">10.3969/j.issn.1672-6073.2012.02.019</pub-id>
</citation>
</ref>
<ref id="B10">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Han</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>N.</given-names>
</name>
<name>
<surname>Jamie</surname>
<given-names>R. S.</given-names>
</name>
</person-group> (<year>2007</year>). <article-title>An adaptability study of Gaussian equation applied to predicting ground settlements induced by tunneling in China</article-title>. <source>Rock Soil Mech.</source> <volume>28</volume> (<issue>01</issue>), <fpage>23</fpage>&#x2013;<lpage>28&#x2b;35</lpage>. <pub-id pub-id-type="doi">10.3969/j.issn.1000-7598.2007.01.005</pub-id>
</citation>
</ref>
<ref id="B11">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>He</surname>
<given-names>L. C.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>M. S.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>Y. J.</given-names>
</name>
</person-group> (<year>2013</year>). <article-title>Supporting technology for shallow large-span neighborhood loess tunnel</article-title>. <source>Rock Soil Mech.</source> <volume>34</volume> (<issue>S2</issue>), <fpage>306</fpage>&#x2013;<lpage>310</lpage>. <pub-id pub-id-type="doi">10.16285/j.rsm.2013.s2.048</pub-id>
</citation>
</ref>
<ref id="B12">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>He</surname>
<given-names>N. N.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>Shao</surname>
<given-names>S. J.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>J. K.</given-names>
</name>
<name>
<surname>Jiao</surname>
<given-names>Y. Y.</given-names>
</name>
</person-group> (<year>2012</year>). <article-title>Ground settlement monitoring above Xi&#x27;an Metro tunnel through the saturated soft loess</article-title>. <source>J. Earth Sci. Environ.</source> <volume>34</volume> (<issue>1</issue>), <fpage>96</fpage>&#x2013;<lpage>103</lpage>. <pub-id pub-id-type="doi">10.3969/j.issn.1672-6561.2012.01.013</pub-id>
</citation>
</ref>
<ref id="B13">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Hu</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>Y. L.</given-names>
</name>
<name>
<surname>Tang</surname>
<given-names>H. M.</given-names>
</name>
<name>
<surname>Xu</surname>
<given-names>H. Q.</given-names>
</name>
</person-group> (<year>2012</year>). <article-title>Research on ground subsidence due to tunnel excavation in huquan-mingdu section of wuhan subway</article-title>. <source>Chin. J. Rock Mech. Eng.</source> <volume>31</volume> (<issue>05</issue>), <fpage>908</fpage>&#x2013;<lpage>913</lpage>. <pub-id pub-id-type="doi">10.3969/j.issn.1000-6915.2012.05.006</pub-id>
</citation>
</ref>
<ref id="B14">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Hu</surname>
<given-names>Z. N.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>Y. G.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>Y. F.</given-names>
</name>
</person-group> (<year>2023</year>). <article-title>Mechanical behaviors of surrounding rock and supporting structure of shallowburied unsymmetrical pressure tunnel crossing soil&#x2013;rock interface</article-title>. <source>Front. Earth Sci.</source> <volume>10</volume>, <fpage>1113430</fpage>. <pub-id pub-id-type="doi">10.3389/feart.2022.1113430</pub-id>
</citation>
</ref>
<ref id="B15">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Hunt</surname>
<given-names>D.</given-names>
</name>
</person-group> (<year>2005</year>). <source>Predicting the ground movements above twin tunnels constructed in London Clay</source>. <publisher-loc>Birmingham</publisher-loc>: <publisher-name>University of Birmingham</publisher-name>.</citation>
</ref>
<ref id="B16">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Jiang</surname>
<given-names>Z. L.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Optimization of construction parameters for the double sidewall guide trench method in shallow buried metro tunnels in loess areas</article-title>. <source>Railw. Eng.</source> <volume>62</volume> (<issue>03</issue>), <fpage>126</fpage>&#x2013;<lpage>131</lpage>. <pub-id pub-id-type="doi">10.3969/j.issn.1003&#x2010;1995.2022.03.28</pub-id>
</citation>
</ref>
<ref id="B17">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Li</surname>
<given-names>C. L.</given-names>
</name>
<name>
<surname>Wu</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Qin</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Shao</surname>
<given-names>Z. S.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Theoretical investigation to the effect of bolt reinforcement on tunnel viscoelastic behavior</article-title>. <source>Arab. J. Sci. Eng.</source> <volume>45</volume> (<issue>5</issue>), <fpage>3707</fpage>&#x2013;<lpage>3718</lpage>. <pub-id pub-id-type="doi">10.1007/s13369-019-04215-9</pub-id>
</citation>
</ref>
<ref id="B18">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Li</surname>
<given-names>M.</given-names>
</name>
</person-group> (<year>2014</year>). <source>Study and prediction the surface settlement in the loess strata during shield construction</source>. <publisher-loc>Xi&#x2019;an</publisher-loc>: <publisher-name>Xi&#x2019;an University of Science and Technology</publisher-name>.</citation>
</ref>
<ref id="B19">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lu</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Yao</surname>
<given-names>A. J.</given-names>
</name>
<name>
<surname>Zheng</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>J. T.</given-names>
</name>
<name>
<surname>Tian</surname>
<given-names>T.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Study on the law and computational method of ground surface settlement induced by double-line tunnel excavation</article-title>. <source>C<italic>hin. J. Rock Mech. Eng</italic>.</source> <volume>38</volume> (<issue>S2</issue>), <fpage>3735</fpage>&#x2013;<lpage>3747</lpage>. <pub-id pub-id-type="doi">10.13722/j.cnki.jrme.2019.0630</pub-id>
</citation>
</ref>
<ref id="B20">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Mahmuto&#x11f;lu</surname>
<given-names>Y.</given-names>
</name>
</person-group> (<year>2011</year>). <article-title>Surface subsidence induced by twin subway tunnelling in soft ground conditions in Istanbul</article-title>. <source>B Eng. Geol. Environ.</source> <volume>70</volume> (<issue>1</issue>), <fpage>115</fpage>&#x2013;<lpage>131</lpage>. <pub-id pub-id-type="doi">10.1007/s10064-010-0289-8</pub-id>
</citation>
</ref>
<ref id="B21">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ocak</surname>
<given-names>I.</given-names>
</name>
</person-group> (<year>2014</year>). <article-title>A new approach for estimating the transverse surface settlement curve for twin tunnels in shallow and soft soils</article-title>. <source>Environ. Earth Sci.</source> <volume>72</volume> (<issue>7</issue>), <fpage>2357</fpage>&#x2013;<lpage>2367</lpage>. <pub-id pub-id-type="doi">10.1007/s12665-014-3145-5</pub-id>
</citation>
</ref>
<ref id="B22">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ou</surname>
<given-names>X. D.</given-names>
</name>
<name>
<surname>Luo</surname>
<given-names>B. X.</given-names>
</name>
<name>
<surname>Jiang</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Qin</surname>
<given-names>Jinxi.</given-names>
</name>
<name>
<surname>Yang</surname>
<given-names>S. N.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Study on prediction model of surface settlement of shield double-ine tunnel in nanning metro</article-title>. <source>
<italic>J Guangxi Univ (Nat Sci Ed</italic>).</source> <volume>44</volume> (<issue>04</issue>), <fpage>1010</fpage>&#x2013;<lpage>1020</lpage>. <pub-id pub-id-type="doi">10.13624/j.cnki.issn.1001-7445.2019.1010</pub-id>
</citation>
</ref>
<ref id="B23">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Qi</surname>
<given-names>W. Q.</given-names>
</name>
<name>
<surname>Yang</surname>
<given-names>Z. Y.</given-names>
</name>
<name>
<surname>Jiang</surname>
<given-names>Y. S.</given-names>
</name>
<name>
<surname>Shao</surname>
<given-names>X. K.</given-names>
</name>
<name>
<surname>Yang</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>He</surname>
<given-names>Q.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Structural deformation of existing horseshoe-shaped tunnels by shield overcrossing</article-title>. <source>Ksce J. Civ. Eng.</source> <volume>25</volume> (<issue>2</issue>), <fpage>735</fpage>&#x2013;<lpage>749</lpage>. <pub-id pub-id-type="doi">10.1007/s12205-020-0599-8</pub-id>
</citation>
</ref>
<ref id="B24">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>She</surname>
<given-names>F. T.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>Y. X.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>Y.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>Surface longitudinal settlement induced by subway tunneling in loess and its prediction analysis method</article-title>. <source>Rock Soil Mech.</source> <volume>36</volume> (<issue>S1</issue>), <fpage>287</fpage>&#x2013;<lpage>292</lpage>. <pub-id pub-id-type="doi">10.16285/j.rsm.2015.S1.049</pub-id>
</citation>
</ref>
<ref id="B25">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Song</surname>
<given-names>Z. P.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>Y, L.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>Y. W.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Research and prospect of demarcation of deep and shallow buried in loess tunnel</article-title>. <source>Hazard Control Tunn Undergr Eng</source> <volume>3</volume> (<issue>2</issue>), <fpage>1</fpage>&#x2013;<lpage>15</lpage>. <pub-id pub-id-type="doi">10.19952/j.cnki.2096-5052.2021.02.01</pub-id>
</citation>
</ref>
<ref id="B26">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Song</surname>
<given-names>Z. P.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>S. H.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>X. G.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>J. B.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>Tao.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>Study on strata settlement regular pattern inducedby tunnel construction based on Peck formula</article-title>. <source>Xi&#x27; Univ. Arch. Tech. (Nat Sci Ed).</source> <volume>50</volume> (<issue>02</issue>), <fpage>190</fpage>&#x2013;<lpage>195</lpage>. <pub-id pub-id-type="doi">10.15986/j.1006-7930.2018.02.006</pub-id>
</citation>
</ref>
<ref id="B27">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Song</surname>
<given-names>Z. P.</given-names>
</name>
<name>
<surname>Tian</surname>
<given-names>X. X.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>Y., W.</given-names>
</name>
<name>
<surname>Zhou</surname>
<given-names>G. N.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Numerical analysis and application of the construction method for the small interval tunnel in the turn line of metro</article-title>. <source>Sci. Progress-Uk</source> <volume>103</volume> (<issue>3</issue>), <fpage>0036850420932067</fpage>. <pub-id pub-id-type="doi">10.1177/0036850420932067</pub-id>
</citation>
</ref>
<ref id="B28">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Song</surname>
<given-names>Z. P.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>S. Y.</given-names>
</name>
<name>
<surname>Su</surname>
<given-names>C. S.</given-names>
</name>
<name>
<surname>Pan</surname>
<given-names>H. W.</given-names>
</name>
<name>
<surname>Wei</surname>
<given-names>S. F.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Optimization of construction scheme for shallow buried tunnel with small spacing asymmetric section in loess stratum</article-title>. <source>J. Xi&#x27; Univ. Arch. Tech. (Nat Sci Ed).</source> <volume>54</volume> (<issue>5</issue>), <fpage>646</fpage>&#x2013;<lpage>656</lpage>. <pub-id pub-id-type="doi">10.15986/j.1006-7930.2022.05.002</pub-id>
</citation>
</ref>
<ref id="B29">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Su</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>Y. J.</given-names>
</name>
<name>
<surname>Cui</surname>
<given-names>J. P.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>C. A.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Installation time of ground support during tunnel excavation: A novel graph methodology</article-title>. <source>Ksce J. Civ. Eeg</source> <volume>24</volume> (<issue>12</issue>), <fpage>3866</fpage>&#x2013;<lpage>3874</lpage>. <pub-id pub-id-type="doi">10.1007/s12205-020-1079-x</pub-id>
</citation>
</ref>
<ref id="B30">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Suwansawat</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Einstein</surname>
<given-names>H. H.</given-names>
</name>
</person-group> (<year>2007</year>). <article-title>Describing settlement troughs over twin tunnels using a superposition technique</article-title>. <source>J. Geotech. Geoenviron</source> <volume>133</volume> (<issue>4</issue>), <fpage>445</fpage>&#x2013;<lpage>468</lpage>. <pub-id pub-id-type="doi">10.1061/(asce)1090-0241(2007)133:4(445)</pub-id>
</citation>
</ref>
<ref id="B31">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Wang</surname>
<given-names>T.</given-names>
</name>
</person-group> (<year>2017</year>). <source>Study on construction optimization and deformation control of shallow buried large span small spacing tunnel in loess areas</source>. <publisher-loc>Xi&#x2019;an</publisher-loc>: <publisher-name>Xi&#x2019;an University of Architecture and Technology</publisher-name>.</citation>
</ref>
<ref id="B32">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wei</surname>
<given-names>J. J.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>The law of surrounding rock and ground surface deformation in loess tunnel construction</article-title>. <source>Sci Tech Eng.</source> <volume>18</volume> (<issue>11</issue>), <fpage>287</fpage>&#x2013;<lpage>292</lpage>. <pub-id pub-id-type="doi">10.3969/j.issn.1671-1815.2018.11.045</pub-id>
</citation>
</ref>
<ref id="B33">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wei</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Hong</surname>
<given-names>Z. L.</given-names>
</name>
<name>
<surname>Yang</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>N.</given-names>
</name>
<name>
<surname>Tan</surname>
<given-names>T. X.</given-names>
</name>
</person-group> (<year>2023</year>). <article-title>Impact of subway shield tunnel construction on deformation of existing utility tunnel</article-title>. <source>Front. Earth Sci.</source> <volume>11</volume>, <fpage>1104865</fpage>. <pub-id pub-id-type="doi">10.3389/feart.2023.1104865</pub-id>
</citation>
</ref>
<ref id="B34">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Xiang</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Jin</surname>
<given-names>B. C.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Tian</surname>
<given-names>C. M.</given-names>
</name>
<name>
<surname>Ou</surname>
<given-names>Y.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Surrounding rock pressure of loess deep large-diameter shield tunnel</article-title>. <source>Chin. Civ. Eng. J.</source> <volume>53</volume> (<issue>S1</issue>), <fpage>112</fpage>&#x2013;<lpage>118</lpage>. <pub-id pub-id-type="doi">10.15951/j.tmgcxb.2020.s1.019</pub-id>
</citation>
</ref>
<ref id="B35">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Xue</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Liang</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>S. H.</given-names>
</name>
<name>
<surname>Ma</surname>
<given-names>Z. Y.</given-names>
</name>
<etal/>
</person-group>(<year>2023b</year>). <article-title>Influence mechanism of brine-gas two-phase flow on sealing property of anisotropic caprock for hydrogen and carbon energy underground storage</article-title>. <source>Int. J. HYDROGEN ENERG</source> <volume>48</volume>(<issue>30</issue>): <fpage>11287</fpage>&#x2013;<lpage>11302</lpage>. <pub-id pub-id-type="doi">10.1016/j.ijhydene.2022.05.173</pub-id>
</citation>
</ref>
<ref id="B36">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Xue</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Chai</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Ranjith</surname>
<given-names>P. G.</given-names>
</name>
<name>
<surname>Cai</surname>
<given-names>C. Z.</given-names>
</name>
<etal/>
</person-group> (<year>2023a</year>). <article-title>Effect of water-cooling shock on fracture initiation and morphology of high-temperature granite: Application of hydraulic fracturing to enhanced geothermal systems</article-title>. <source>Appl. ENERG</source> <volume>337</volume>, <fpage>120858</fpage>. <pub-id pub-id-type="doi">10.1016/j.apenergy.2023.120858</pub-id>
</citation>
</ref>
<ref id="B37">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Xue</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Ranjith</surname>
<given-names>P. G.</given-names>
</name>
<name>
<surname>Gao</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>Z. Z.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>S. H.</given-names>
</name>
</person-group> (<year>2023c</year>). <article-title>Experimental investigations on effects of gas pressure on mechanical behaviors and failure characteristic of coals</article-title>. <source>J. ROCK Mech. Geotech.</source> <volume>15</volume> (<issue>2</issue>), <fpage>412</fpage>&#x2013;<lpage>428</lpage>. <pub-id pub-id-type="doi">10.1016/j.jrmge.2022.05.013</pub-id>
</citation>
</ref>
<ref id="B38">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Yang</surname>
<given-names>Z. D.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Study on shield tunnel section size for Tianjin metro new line</article-title>. <source>Mod. Urt</source> (<issue>10</issue>), <fpage>47</fpage>&#x2013;<lpage>50</lpage>.</citation>
</ref>
<ref id="B39">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Yin</surname>
<given-names>G. M.</given-names>
</name>
<name>
<surname>Fu</surname>
<given-names>H. L.</given-names>
</name>
<name>
<surname>Hou</surname>
<given-names>W. Z.</given-names>
</name>
<name>
<surname>Yang</surname>
<given-names>S. S.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Introduction of several methods of value of Peck formula parameters</article-title>. <source>J. Railw. Sci. Eng.</source> <volume>19</volume> (<issue>07</issue>), <fpage>2015</fpage>&#x2013;<lpage>2022</lpage>. <pub-id pub-id-type="doi">10.19713/j.cnki.43-1423/u.T20210810</pub-id>
</citation>
</ref>
<ref id="B40">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Yu</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Jiang</surname>
<given-names>X. H.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>D. L.</given-names>
</name>
<name>
<surname>Xu</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Kong</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Cui</surname>
<given-names>C.</given-names>
</name>
<etal/>
</person-group> (<year>2019</year>). <article-title>RBCK1 promotes p53 degradation via ubiquitination in renal cell carcinoma</article-title>. <source>Mod. Tunn. Technol.</source> <volume>56</volume> (<issue>S2</issue>), <fpage>254</fpage>&#x2013;<lpage>262</lpage>. <pub-id pub-id-type="doi">10.1038/s41419-019-1488-2</pub-id>
</citation>
</ref>
<ref id="B41">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhang</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Zhou</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Lai</surname>
<given-names>J. X.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Feng</surname>
<given-names>Z. H,</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Optimization of clearance of tunnel in xi&#x27;an metro constructed with double shield TBM</article-title>. <source>Road. Mach. Constr. Mech.</source> <volume>34</volume> (<issue>10</issue>), <fpage>78</fpage>&#x2013;<lpage>83</lpage>. <pub-id pub-id-type="doi">10.3969/j.issn.1000-033X.2017.10.011</pub-id>
</citation>
</ref>
<ref id="B42">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhang</surname>
<given-names>Y. T.</given-names>
</name>
<name>
<surname>Sun</surname>
<given-names>X. L.</given-names>
</name>
<name>
<surname>Yang</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Sun</surname>
<given-names>Y. J.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Influence of local expansion on safety of tunnel support</article-title>. <source>Tunn. Constr.</source> <volume>41</volume> (<issue>S2</issue>), <fpage>306</fpage>&#x2013;<lpage>312</lpage>. <pub-id pub-id-type="doi">10.3973/j.issn.2096-44982021.S2.039</pub-id>
</citation>
</ref>
<ref id="B43">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhang</surname>
<given-names>Z. Y.</given-names>
</name>
<name>
<surname>Zhao</surname>
<given-names>C. Y.</given-names>
</name>
<name>
<surname>Long</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>X. H.</given-names>
</name>
<name>
<surname>Lei</surname>
<given-names>M. F.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Research on the stability of shallow-buried large cross-section tunnel by construction method conversion</article-title>. <source>Front. Earth Sci.</source> <volume>10</volume>, <fpage>37</fpage>. <pub-id pub-id-type="doi">10.3389/feart.2022.831169</pub-id>
</citation>
</ref>
<ref id="B44">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zheng</surname>
<given-names>X. X.</given-names>
</name>
<name>
<surname>Wu</surname>
<given-names>S. S.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Research on deformation control of asymmetric neighborhood large-span shallow buried city tunnel</article-title>. <source>Tech Hwy Trans.</source> <volume>35</volume> (<issue>03</issue>), <fpage>94</fpage>&#x2013;<lpage>99</lpage>. <pub-id pub-id-type="doi">10.13607/j.cnki.Gljt.2019.03.016</pub-id>
</citation>
</ref>
<ref id="B45">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhou</surname>
<given-names>P. Y.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>J. B.</given-names>
</name>
<name>
<surname>Song</surname>
<given-names>Z. P.</given-names>
</name>
<name>
<surname>Cao</surname>
<given-names>Z. L.</given-names>
</name>
<name>
<surname>Pei</surname>
<given-names>Z. M.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Construction method optimization for transfer section between cross passage and main tunnel of metro station</article-title>. <source>Front. Earth Sci.</source> <volume>10</volume>, <fpage>183</fpage>. <pub-id pub-id-type="doi">10.3389/feart.2022.770888</pub-id>
</citation>
</ref>
<ref id="B46">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhou</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Ding</surname>
<given-names>H. H.</given-names>
</name>
<name>
<surname>Miao</surname>
<given-names>L. W.</given-names>
</name>
<name>
<surname>Gong</surname>
<given-names>C. J.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Predictive model for the surface settlement caused by the excavation of twin tunnels</article-title>. <source>Tunn. Undergr. Sp. Tech.</source> <volume>114</volume> (<issue>3</issue>), <fpage>104014</fpage>. <pub-id pub-id-type="doi">10.1016/j.tust.2021.104014</pub-id>
</citation>
</ref>
</ref-list>
</back>
</article>