<!DOCTYPE article PUBLIC "-//NLM//DTD JATS (Z39.96) Journal Archiving and Interchange DTD v1.0 20120330//EN" "JATS-archivearticle1.dtd">
<article xmlns:xlink="http://www.w3.org/1999/xlink">
  <front>
    <journal-meta />
    <article-meta>
      <title-group>
        <article-title>Investigation of Indoor Positioning Technologies for Underground Mine Environments</article-title>
      </title-group>
      <contrib-group>
        <contrib contrib-type="author">
          <string-name>Binghao Li</string-name>
          <email>binghao.li@unsw.edu.au</email>
          <xref ref-type="aff" rid="aff2">2</xref>
        </contrib>
        <contrib contrib-type="author">
          <string-name>Kai Zhao</string-name>
          <xref ref-type="aff" rid="aff2">2</xref>
        </contrib>
        <contrib contrib-type="author">
          <string-name>Serkan Saydam</string-name>
          <xref ref-type="aff" rid="aff2">2</xref>
        </contrib>
        <contrib contrib-type="author">
          <string-name>Chris Rizos</string-name>
          <xref ref-type="aff" rid="aff2">2</xref>
        </contrib>
        <contrib contrib-type="author">
          <string-name>Qiang Wang</string-name>
          <xref ref-type="aff" rid="aff1">1</xref>
        </contrib>
        <contrib contrib-type="author">
          <string-name>Jian Wang</string-name>
          <xref ref-type="aff" rid="aff0">0</xref>
        </contrib>
        <aff id="aff0">
          <label>0</label>
          <institution>Beijing University of Civil Engineering and Architecture</institution>
          ,
          <addr-line>Beijing</addr-line>
          ,
          <country country="CN">China</country>
        </aff>
        <aff id="aff1">
          <label>1</label>
          <institution>China University of Mining and Technology</institution>
          ,
          <addr-line>Xuzhou</addr-line>
          ,
          <country country="CN">China</country>
        </aff>
        <aff id="aff2">
          <label>2</label>
          <institution>University of New South Wales</institution>
          ,
          <addr-line>Sydney</addr-line>
          ,
          <country country="AU">Australia</country>
        </aff>
      </contrib-group>
      <abstract>
        <p>Positioning in underground mining environments is a key requirement for ensuring the safety of mine workers. It is also a critical technological capability in resolving mine productivity bottlenecks, which has a great economic impact in Australia. To support the growth of the mining sector, innovative technologies need to be developed, with underground positioning an important though significant engineering challenge. This paper aims to demonstrate a robust high accuracy positioning system for underground mining environments to meet the requirements of worker safety and mine efficiency improvement. Several technologies which could be part of the “mix” in the solution to the challenges have been identified. Tests have been carried out in a metal mine. Radio Frequency signal plus Inertial Measurement Unit, multisensor integration, geomagnetic field positioning and ultra-wide band have been tested. The data are analysed and the results are reported.</p>
      </abstract>
      <kwd-group>
        <kwd>Indoor positioning</kwd>
        <kwd>Underground mining</kwd>
      </kwd-group>
    </article-meta>
  </front>
  <body>
    <sec id="sec-1">
      <title>-</title>
      <p>
        The safety of mine workers is one of the highest priorities of Australia’s mining
industry. When an accident occurs, the immediate initiation of a search and rescue response
is vital, because the survival rate decreases rapidly as time passes - the so-called
“golden 72 hours” [
        <xref ref-type="bibr" rid="ref1">1</xref>
        ]. The greatest unknown for the rescue team is to know how many
people are trapped and where they are located. If the locations of the victims are known,
time could be saved and the chances of success markedly improved. Consequently,
knowledge of mine workers’ position is highly correlated with their safety. In the worst
case, when an underground positioning and communication system is completely
disabled, the last known position of the miners is extremely useful.
      </p>
      <p>
        A key objective of the mining industry is to achieve zero-harm in every work place
through continuous improvement, intensive training, introducing advanced work
practices and implementing new technologies. Considerable effort has been made to
develop safety-related technologies, via new training systems, legislation and regulations
[
        <xref ref-type="bibr" rid="ref2">2</xref>
        ][
        <xref ref-type="bibr" rid="ref3">3</xref>
        ]. However, mine disasters still occur. The ordeal of the two survivors of the
Beaconsfield Mine in Australia in April 2006 still resonates. On 19 November 2010, the
Pike River Mine accident killed 29 people and injured 2 in New Zealand. Just a few
months earlier in the United States 29 people were killed in the Upper Big Branch
Mine disaster. In developing countries, mine disasters occur more frequently and are
more serious. For instance, in 2014 the Soma coal mine disaster in Turkey killed 301
workers [
        <xref ref-type="bibr" rid="ref4">4</xref>
        ]. In China, although the number of deaths has decreased over the past 10
years, it was still almost 1500 in 2012 [
        <xref ref-type="bibr" rid="ref5">5</xref>
        ].
      </p>
      <p>
        In addition, according to [
        <xref ref-type="bibr" rid="ref6">6</xref>
        ], position determination has been identified by the mining
industry as one of the technologies critical to resolving major mine productivity
bottlenecks. Positioning is associated with process optimisation and control, operation and
maintenance. For instance, if any equipment needs repair, the location of the
equipment is needed before a technician can be sent underground. Knowing the location of
people and equipment, and the patterns of their movement can also improve
productivity.
      </p>
      <p>
        For open-pit mines or other surface activities, the Global Navigation Satellite System
(GNSS) is the preferred technology. However, positioning in the underground mine
environment is a challenge due to the lack of a GNSS-like technology. The positioning
systems on the market can be classified as belonging to one of two generations. The
first-generation products use radio-frequency identification scanners to monitor the
tags carried by the workers passing a scanner. It can only report proximity to a scanner
if a worker moves into an area [
        <xref ref-type="bibr" rid="ref7">7</xref>
        ]. The second-generation products typically use
Zigbee or WiFi to trilaterate the location of the workers using the communications
signals. This kind of system requires deployment of many Zigbee nodes or WiFi access
points, which typically makes the implementation of such systems more expensive [
        <xref ref-type="bibr" rid="ref8">8</xref>
        ].
A new generation system requires quick deployment, positioning on server and client
sides, high accuracy, and large coverage area, as well as cost effectiveness. We are
seeking an elegant solution: utilising innovation in a variety of fields - including radio
frequency (RF) signal, geomagnetic field, low-cost inertial measurement unit, and
others - deploy as few as possible of the signal transmitters, develop a novel algorithm
to integrate these technologies to achieve global (i.e. whole underground mine area)
positioning.
2
2.1
      </p>
    </sec>
    <sec id="sec-2">
      <title>Our Approaches</title>
      <sec id="sec-2-1">
        <title>BLE Signal Propagation Model</title>
        <p>Bluetooth Low Energy (BLE) can be used to transmit the RF signal which can give the
distance from the transmitter to the receiver. It is a low-cost technology that has low
power consumption – an AA size Lithium-ion battery (2700mAh) can power it for
over 3 years.</p>
        <p>Radio waves are affected by many phenomena such as reflection, refraction,
diffraction, absorption and scattering etc. The propagation of radio waves is characterized by
several factors. When all the factors are considered, the loss is
2
n  4 
L  d     La
(1)
where d is the transmitter-receiver separated distance, La is the attenuation caused by
obstacles, n is the path loss exponent (2 in the case of free space) and λ is the
wavelength (2.45 GHz, the centre frequency of 2.4 GHz band). Expressed in decibels</p>
        <p>L ( dB )  40.23  10  n  log( d )  La ( dB )
When a signal is propagating along a corridor, normally the wave guiding effect might
make the exponent lower than 2. However, since the BLE node was installed 3 meters
above the ground and when the range is calculated, the height is not considered, the
value of 2.1 was used. If the propagation is LOS, La can be ignored. The d can be
estimated based on Equation (2).</p>
        <p>Obviously, using model to convert signal strength to range is not very accurate,
however, based on the signal strength, one can detect if the receiver is close to a BLE node.
To estimate more accurate position of the receiver, other technologies are required.
2.2</p>
      </sec>
      <sec id="sec-2-2">
        <title>Using IMU to Estimate the Distance</title>
        <p>A low-cost IMU is used to detect the worker’s steps, and then used to estimate the
distance the person has travelled (using a vehicle is another case to be investigated).
RF signal and magnetic field sensors can be used to (re)initialise the IMU.
In order for the step detection algorithm to work independently of the device’s
orientation, only the amplitude of the acceleration is retained.</p>
        <p>
          a 
k
a x2k
 a 2
yk
 a 2
zk
Then the Earth’s gravity is removed from the values by applying a high-pass filter to
the acceleration amplitudes. General speaking, there are two types of algorithm can be
applied for step detection, one is operating in time domain; the other is operating in
frequency domain [
          <xref ref-type="bibr" rid="ref10">10</xref>
          ]. A Fast Fourier Transform (FFT) based algorithm and a
windowed peak detection algorithm were tested. The FFT based algorithm was chosen as
it was more reliable in most of the scenarios.
        </p>
        <p>
          Once the step frequency has been estimated, the step size is still needed in order to
compute the estimated speed of the user. In our context, the step size is defined as the
distance from the heel of one foot to the heel of the other foot. There are many
different models to estimate the step size. For instance, in [
          <xref ref-type="bibr" rid="ref11">11</xref>
          ], the step size is assumed to
increase linearly with the step frequency. However, the step size is not only a function
of step frequency but also depends on the height of the user. In [
          <xref ref-type="bibr" rid="ref12">12</xref>
          ], the step size is
modelled as a linear function of the height ls=k*h where k only depends on the sex of
the user. In our application k varies linearly with the step frequency following this
model:
ls  (a  b  fs   ) * h

  ~ N (0, s )
The parameters of the model (a,b,σs) are estimated from empirical data.
(2)
(3)
(4)
        </p>
      </sec>
      <sec id="sec-2-3">
        <title>Magnetic Field Positioning</title>
        <p>
          The magnetic field can be utilised to improve the positioning accuracy between two
RF signal transmitters. When some of the RF transmitters are switch off, due to an
accident or a flat battery, magnetic field positioning can be used instead. However, a
smart way must be developed to create and maintain the magnetic field database.
Previous research [
          <xref ref-type="bibr" rid="ref9">9</xref>
          ] has revealed that the stability of the magnetic field implies that
travelling along the same path will generate the same curve and applying fingerprinting
technology is not reliable as in reality only two elements (the vertical component and
horizontal component of the magnetic field) can be used.
        </p>
        <p>As the tunnel restrains a worker’s movement, it is more reliable to apply pattern
matching for positioning. The magnetic curve matching algorithms is required using
magnetic field for positioning. A few matching algorithm have been investigated.
2.4</p>
      </sec>
      <sec id="sec-2-4">
        <title>Other Technologies</title>
        <p>
          Barometers can detect the change of air pressure, and the change can indicate the
movement of the sensor in the vertical direction [
          <xref ref-type="bibr" rid="ref13">13</xref>
          ]. Ultra-wideband (UWB) systems
can provide 10cm level positioning accuracy.
3
3.1
        </p>
      </sec>
    </sec>
    <sec id="sec-3">
      <title>Testing and Analysing</title>
      <sec id="sec-3-1">
        <title>Testing Bed</title>
        <p>The testing bed was selected in a copper and gold mine located in New South
Wales, Australia which is a block cave underground mine. The area chosen for our
test is about 5km away from the entrance of the mine and 500m under the surface
level. The part highlighted with green is chosen to deploy our BLE tags and carry out
most of the tests (see Figure 1). The drill drives along the tunnel have been mined out;
it is quiet and ideal for testing. The straight line part of the highlighted tunnel is about
100m and the curve part is about 70m. The floor of the first 90m is relatively flat, the
altitude only decreases slightly, however, the change of altitude is significant although
the exact value is unknown. The tunnel is about 5m wide and 5.5m high (refer to
Figure 1).</p>
      </sec>
      <sec id="sec-3-2">
        <title>Testing setup</title>
        <p>For this test, two data logging system has been developed. One is based on a STM32
development board, the other is based on Raspberry Pi 3B, details can be found in
Figure 2. The reasons that two loggers were developed are to evaluate more IMU
unites and to make sure at least one set of data can be collected successfully. 170m
was measured using tap measure and every ten meters a mark was painted on the
wall. A BLE tag was installed on the square plates attached to the ground support
(bolt) on the side-wall about 3.5m above the floor. The tag was installed align with
the mark as good as possible. In total, 17 BLE tags were installed.
When there are a few BLE transmitters deployed in a line (as in the case of a tunnel), a
moving receiver (carried by a person) can detect the change of the nearby BLE
transmitters’ signal strength, and these changes should be useful for locating the receiver.
Two testers carried data logger A and B respectively (the rectangular box was attached
to the waist of a tester), walked along the testing path started from 0m mark to 170m
mark and then turned around, back to 0m with a speed that was judged as constant as
possible. The tests have been done four times and then the testers swapped the devices,
carried out another three rounds of tests. Figure 4 (left) shows the variation in received
signal strength of one of the test. Clearly, there are 33 peaks of signal strength, which
indicate the location when the BLE receiver passed by a transmitter. These peaks,
together with other signal strength measurements, can be used to determine the receiver’s
range from a reference location. Figure 4 (right) compares the estimated range based
on the received signal strength and the true (known) range.
“Dead reckoning” is a well-known relative navigation technique. Starting from a
known position, successive position displacements are estimated and then
accumulated. Pedestrian Dead-Reckoning (PDR) solutions integrate step lengths and orientation
estimations at each detected step, so as to compute the final position and orientation of
a person. Inertial Measurement Units (IMU) typically comprises several
accelerometers, gyroscopes, and perhaps magnetometers. A low-cost IMU can be used to estimate
a person’s relative position by implementing a PDR-type solution.
Figure 5 (left) is a plot of the angle data collected during the test. It can be seen that the
tester walked along the straight line first, and then started to turn in the curve part. The
significant change indicates that the tester turned around and moved back to the start
point along the same path. After applying the step detection and step length estimation,
the trajectory of the IMU (the tester) is shown in Figure 5 (right) (the green line).
Obviously, the trajectory gradually drifts away from the real path. After applying
landmark (where the turning around was detected) correction, the result can be improved
significantly (the blue line). The further improved red one is the final trajectory after
applying map matching. BLE and magnetic field positioning can be used as landmarks
to correct the drift of the PDR solution.
3.5</p>
      </sec>
      <sec id="sec-3-3">
        <title>Magnetic Field Based Positioning</title>
        <p>
          The main issues of magnetic field positioning are the database creation and the
positioning algorithm. The proposed database creation will be separated into two phases. In
the first phase, when the system is deployed, the magnetic field data will be collected.
This is the basic database. In the second phase, while the positioning system is running,
the mine workers will collect the data which will be used to refine the database.
The test was carried out in the same testbed. An Xsens Mti-100 was used to collect the
magnetic field data. A tester carried the IMU, walked from start point to end point and
back to start point, twice. The data collected in the first traverses were used as the
reference dataset, and that of the second traverse was used as the testing data. After
applying a low-pass filter to remove the noise and using an averaging window to reduce the
sampling frequency (from 100Hz to 20Hz), the dynamic time warping (DTW) [
          <xref ref-type="bibr" rid="ref14">14</xref>
          ]
pattern matching algorithm was applied. Only the magnitude (total intensity) was
calculated. Figure 6 illustrates the results. The blue, green curves in figure denote the
reference data and testing data respectively. The testing data were divided into small
fractions (every 10s) to apply pattern matching separately. Most of the matching of the
fractions is successful; however the matching of the first and the last fraction is failed.
There could be two reasons: At the start point or end point, the tester changed his status
from static to dynamic or vice versa, it is much harder to control the speed; also, the
tester’s gesture might be different at each test.
The testing in a typical metalliferous mine has demonstrated the feasibility of a
costeffective positioning system for underground mining environment based on multiple
technologies. The received signal strength of BLE beacons alone can be used for
positioning if there are densely deployed beacons. Applying the step detection and step
length estimation based on IMU can generally provide an accurate estimation of the
position if re-initialisation based on landmarks and map matching can be applied. A
sparsely deployed BLE beacons can be used as a type of landmark as well as the
magnetic field. Collecting magnetic field data for underground mine is not an easy task;
however, crowd sourcing is a possible solution. After several rounds of iteration, it is
possible to build a reliable database. Other sensors such as a barometer which can
obtain the change of the height can also be used to provide the landmarks. Although
UWB can provide accurate range estimation, the requirement of external power supply
for the base stations and high power consumption of the tags makes this technology not
cost-effective.
        </p>
      </sec>
    </sec>
  </body>
  <back>
    <ref-list>
      <ref id="ref1">
        <mixed-citation>
          1.
          <string-name>
            <given-names>S.</given-names>
            ,
            <surname>Tadocoro</surname>
          </string-name>
          ,
          <string-name>
            <surname>S.</surname>
          </string-name>
          , “
          <article-title>Rescue Robotics: DDT Project on Robots and Systems for Urban Search</article-title>
          and Rescue”, Springer,
          <source>ISBN 978-1-84882-474-4</source>
          ,
          <year>2009</year>
          .
        </mixed-citation>
      </ref>
      <ref id="ref2">
        <mixed-citation>
          2.
          <string-name>
            <surname>R.</surname>
          </string-name>
          , Gurtunca, “
          <article-title>Possible Impact of New Safety Technology Developments on the Future of the United States Mining Industry</article-title>
          ,” First International Future Mining Conference,
          <volume>19</volume>
          -
          <fpage>21</fpage>
          November,
          <year>2008</year>
          , Sydney, Australia,
          <fpage>3</fpage>
          -
          <lpage>9</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref3">
        <mixed-citation>
          3.
          <string-name>
            <given-names>S.</given-names>
            ,
            <surname>Saydam</surname>
          </string-name>
          ,
          <string-name>
            <given-names>R.</given-names>
            ,
            <surname>Mitra</surname>
          </string-name>
          , and
          <string-name>
            <given-names>B.</given-names>
            ,
            <surname>Hebblewhite</surname>
          </string-name>
          , “
          <article-title>Implementation of Virtual Reality Technology in Mine Safety Training and Mining Engineering Education in Australia,” 19th Coal Congress</article-title>
          ,
          <fpage>21</fpage>
          -23 May,
          <year>2014</year>
          , Zonguldak, Turkey.
        </mixed-citation>
      </ref>
      <ref id="ref4">
        <mixed-citation>4. http://www.bbc.com/news/world-europe-27459912</mixed-citation>
      </ref>
      <ref id="ref5">
        <mixed-citation>5. http://www.fawan.com/Article/gn/sh/2013/01/18/142349183926.html</mixed-citation>
      </ref>
      <ref id="ref6">
        <mixed-citation>
          6.
          <string-name>
            <given-names>D.</given-names>
            ,
            <surname>Peterson</surname>
          </string-name>
          ,
          <string-name>
            <given-names>T.</given-names>
            ,
            <surname>LaTourrette</surname>
          </string-name>
          , and
          <string-name>
            <surname>J.T.</surname>
          </string-name>
          , Bartis, “New Forces at Work in Mining: Industry Views of Critical Technologies,”
          <source>Rand Corp, ISBN 0-8330-2967-3</source>
          ,
          <year>2001</year>
          .
        </mixed-citation>
      </ref>
      <ref id="ref7">
        <mixed-citation>
          7.
          <string-name>
            <given-names>W.</given-names>
            ,
            <surname>Zhao</surname>
          </string-name>
          , and R.,
          <source>Hu, “Research on 3D Positioning and Navigation Technologies in Underground Mine,” ISECS'10</source>
          ,
          <fpage>29</fpage>
          -
          <issue>31</issue>
          <year>July</year>
          ,
          <year>2010</year>
          , Guangzhou, China,
          <fpage>165</fpage>
          -
          <lpage>168</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref8">
        <mixed-citation>
          8.
          <string-name>
            <surname>J.</surname>
          </string-name>
          , Marshall, “
          <article-title>Navigating the Advances in Underground Navigation”</article-title>
          ,
          <source>CIM Magazine</source>
          , vol.
          <volume>5</volume>
          ,
          <fpage>20</fpage>
          -
          <lpage>21</lpage>
          ,
          <year>2010</year>
          .
        </mixed-citation>
      </ref>
      <ref id="ref9">
        <mixed-citation>
          9.
          <string-name>
            <given-names>B.</given-names>
            ,
            <surname>Li</surname>
          </string-name>
          ,
          <string-name>
            <given-names>T.</given-names>
            ,
            <surname>Gallagher</surname>
          </string-name>
          ,
          <string-name>
            <given-names>C.</given-names>
            ,
            <surname>Rizos</surname>
          </string-name>
          and
          <string-name>
            <surname>A.G.</surname>
          </string-name>
          , Dempster, “
          <article-title>Using Geomagnetic Field for Indoor Positioning</article-title>
          ,
          <source>” Journal of Applied Geodesy</source>
          ,
          <volume>7</volume>
          (
          <issue>4</issue>
          ),
          <fpage>299</fpage>
          -
          <lpage>308</lpage>
          ,
          <year>2013</year>
          .
        </mixed-citation>
      </ref>
      <ref id="ref10">
        <mixed-citation>
          10.
          <string-name>
            <given-names>A.</given-names>
            ,
            <surname>Brajdic</surname>
          </string-name>
          , and
          <string-name>
            <given-names>R.</given-names>
            <surname>Harle</surname>
          </string-name>
          , “September.
          <article-title>Walk detection and step counting on unconstrained smartphones,”</article-title>
          <source>In Proceedings of the 2013 ACM international joint conference on Pervasive and ubiquitous computing</source>
          (pp.
          <fpage>225</fpage>
          -
          <lpage>234</lpage>
          ). ACM.
        </mixed-citation>
      </ref>
      <ref id="ref11">
        <mixed-citation>
          11. H.
          <string-name>
            <surname>Lemelson</surname>
            ,
            <given-names>M. B.</given-names>
          </string-name>
          <string-name>
            <surname>Kjaergaard</surname>
            ,
            <given-names>R.</given-names>
          </string-name>
          <string-name>
            <surname>Hansen</surname>
            , and
            <given-names>T.</given-names>
          </string-name>
          <string-name>
            <surname>King</surname>
          </string-name>
          , “
          <source>Error Estimation for Indoor</source>
          <volume>802</volume>
          .11 Location Fingerprinting,”
          <source>in Proceedings of the 4th International Symposium on Location and Context Awareness</source>
          , Berlin, Heidelberg,
          <year>2009</year>
          , pp.
          <fpage>138</fpage>
          -
          <lpage>155</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref12">
        <mixed-citation>
          12. I. Bylemans,
          <string-name>
            <given-names>M.</given-names>
            <surname>Weyn</surname>
          </string-name>
          , and
          <string-name>
            <given-names>M.</given-names>
            <surname>Klepal</surname>
          </string-name>
          , “
          <article-title>Mobile Phone-Based Displacement Estimation for Opportunistic Localisation Systems,” in Mobile Ubiquitous Computing, Systems</article-title>
          , Services and Technologies, UBICOMM '
          <fpage>09</fpage>
          . Third International Conference on,
          <year>2009</year>
          , pp.
          <fpage>113</fpage>
          -
          <lpage>118</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref13">
        <mixed-citation>
          13.
          <string-name>
            <given-names>B.</given-names>
            ,
            <surname>Li</surname>
          </string-name>
          ,
          <string-name>
            <given-names>B.</given-names>
            ,
            <surname>Harvey</surname>
          </string-name>
          , and
          <string-name>
            <given-names>T.</given-names>
            ,
            <surname>Gallagher</surname>
          </string-name>
          ,
          <string-name>
            <surname>T.</surname>
          </string-name>
          ,
          <year>2013</year>
          ,
          <string-name>
            <surname>October.</surname>
          </string-name>
          <article-title>Using barometers to determine the height for indoor positioning</article-title>
          .
          <source>In Indoor Positioning and Indoor Navigation (IPIN)</source>
          , 2013 International Conference on (pp.
          <fpage>1</fpage>
          -
          <lpage>7</lpage>
          ). IEEE.
        </mixed-citation>
      </ref>
      <ref id="ref14">
        <mixed-citation>
          14. D.J.,
          <string-name>
            <surname>Berndt</surname>
          </string-name>
          , and J., Clifford, “
          <article-title>Using dynamic time warping to find patterns in time series,”</article-title>
          <source>In KDD workshop</source>
          (Vol.
          <volume>10</volume>
          , No.
          <volume>16</volume>
          , pp.
          <fpage>359</fpage>
          -
          <lpage>370</lpage>
          ),
          <year>July</year>
          ,
          <year>1994</year>
          .
        </mixed-citation>
      </ref>
    </ref-list>
  </back>
</article>