<!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>
      <journal-title-group>
        <journal-title>Khabarovsk, Russia</journal-title>
      </journal-title-group>
    </journal-meta>
    <article-meta>
      <title-group>
        <article-title>The Calculation Parameters for the Effective Seismic Sensors Placements to Monitor Burst-Hazard Rock Massif</article-title>
      </title-group>
      <contrib-group>
        <contrib contrib-type="author">
          <string-name>Andrey V. Gladyr</string-name>
        </contrib>
        <contrib contrib-type="author">
          <string-name>Igor Ju. Rasskazov</string-name>
        </contrib>
        <contrib contrib-type="author">
          <string-name>Alexander V. Konstantinov</string-name>
        </contrib>
        <aff id="aff0">
          <label>0</label>
          <institution>Institute of Mining</institution>
          ,
          <addr-line>Far East Branch</addr-line>
          ,
          <institution>Russian Academy of Sciences</institution>
          ,
          <addr-line>Khabarovsk, Russian Federation</addr-line>
        </aff>
      </contrib-group>
      <pub-date>
        <year>2019</year>
      </pub-date>
      <volume>1</volume>
      <fpage>6</fpage>
      <lpage>19</lpage>
      <abstract>
        <p>An algorithm for constructing the optimal configuration of an observational network of a seismo-acoustic system for monitoring rock pressure has been developed. Software module was implemented based on the developed algorithm. It allows the construction of maps of the observation network sensitivity at current time and after the addition of seismic receivers. Two operating modes of the software module are implemented: in the first, a map is constructed based on gradient levels, in the second, based on control levels reflecting location confidence of seismic acoustic signals. Recommendations for working with the software module are given. Main shortcomings of the developed algorithm for constructing the optimal configuration of the observational network are indicated and directions for further improvement of the algorithm are formulated.</p>
      </abstract>
    </article-meta>
  </front>
  <body>
    <sec id="sec-1">
      <title>Introduction</title>
      <p>
        As long-term practice shows, underground mining at a number of mining enterprises is complicated by shock
hazardous conditions. An increase in the number and intensity of various geodynamic phenomena is predicted with
a further decrease in the horizons of mining [
        <xref ref-type="bibr" rid="ref1 ref2">1, 2</xref>
        ]. In these conditions, more in-depth research is needed. It
includes the study of geodynamic fields laws and technogenic influence processes of mining operations using
seismic, seismo acoustic, geodesic, seismic-deformation and other methods [
        <xref ref-type="bibr" rid="ref3 ref4">3, 4</xref>
        ].
      </p>
      <p>
        Monitoring of the geomechanical state of a shock-hazardous rock mass using automated systems is one of the
promising methods for preventing the occurrence of mountain and mountain tectonic impacts. These systems
provide real-time registration and operational processing of the parameters of seismo acoustic events that contain
information about geomechanical and geodynamic processes occurring in a massif [
        <xref ref-type="bibr" rid="ref5 ref6 ref7">5 -7</xref>
        ].
      </p>
      <p>
        The seismo acoustic monitoring system of rock pressure “Prognoz ADS” [
        <xref ref-type="bibr" rid="ref8 ref9">8, 9</xref>
        ] is a measuring and computing
complex. An innovative approach was taken in developing this system. It provides the conversion of analog
information into digital in the immediate neighborhood of the primary transducer (geophone). Further processing
of monitoring data takes place in digital form.
      </p>
      <p>This approach has the following advantages:
– wide dynamic range of transmitted signals;
– high noise immunity;
– the ability to remotely control, configure and remotely update the digital receiver firmware.</p>
      <p>The system is divided into surface and underground parts. The underground part is represented by a distributed
network of digital receivers-converters with connected primary converters and an underground equipment room in
which repeaters and power supply are located. They also perform functions of a time synchronization signal
generator.</p>
      <p>The surface part is a center for receiving and processing seismic pulses. It is organized on the basis of a
personal computer with installed specialized software. The connection between the undergr ound and surface parts
is realized through the main communication line.</p>
      <p>The task of optimal configuration of a distributed observational network.</p>
      <p>A key element of the developed system is a distributed observational network of underground digital receiving
transducers. Each of them is made in the borehole version and consists of a piezoceramic-type primary transducer
(geophone) and an acoustic signal preprocessing unit (digital receiver of the transducer).</p>
      <p>Digital receivers and primary transducers are installed in horizontal wells 2.5 m long and with a diameter of 76 to
105 mm drilled in the sides of underground mines with a bottom elevation of 2 ÷ 3 °.</p>
      <p>The long-term practical application of the mountain pressure monitoring system “Prognoz ADS” in a number of
enterprises has shown that one of the most actual tasks requiring solution at the projection stage is the task of optimal
configuration of the underground observation network of receivers of acoustic signals.</p>
      <p>The increase in the scale of the work performed in the controlled areas and the displacement of their front requires
a gradual expansion of the controlled spatial zone, and, consequently, an increase in the number of seismic acoustic
sensors to cover the required volume of the monitored array with the observation network. Therefore, as the number
of geophones increases, the problem of choosing the optimal place for their installation is most acute.
2</p>
    </sec>
    <sec id="sec-2">
      <title>Development of Methods for the Optimal Geophone Location</title>
      <p>The construction of the optimal configuration of the observational network is a nontrivial task and requires
consideration of many factors. However, the solution is greatly simplified with the introduction of a number of
conventions and restrictions. In this case, it allows formulating methodological approaches to the optimal location of
geophones in the first approach and to highlight directions for further improvement of developed approach.</p>
      <p>Constructing the optimal configuration of the observational network provides for the choice of a location of
geophones that would ensure the maximum possible volume of the zone of confident seismo acoustic signal
recording. The choice of locations is proposed to be carried out on the basis of horizon plans since it is intended to
deploy geophones at different horizons of the mine field and sensitivity maps of the observation network calculated
for each horizon. They should take into account the existing mining and geological situation.</p>
      <p>The algorithm for calculating the sensitivity map of the observational network is the main analytical algorithm for
the optimal location method of geophones and provides for the calculation of the sensitivity coefficient for a spatial
zone represented by a volume domain set of cubic form. Their size is allowed to customize.</p>
      <p>The input information of the algorithm provides for the acquisition of the following data in a table form and
abstract objects constructed on their basis (Fig. 1):</p>
      <p>– sensitivity limit of a geophone is a scalar quantity. It characterizes the limiting value of the distance between the
signal source and the geophone, at which the average useful seismo acoustic signal becomes indistinguishable at
noise level;
– table of the location of the actually installed geophones;
– table of virtual geophones location;
– spatial zone characteristic table for building sensitivity map;
– graphic representations of plans for horizons of underground mines.</p>
      <p>а)</p>
      <p>б)</p>
      <p>The Use of Modern Information Technology for Research of Technical and Natural Systems
________________________________________________________________________________________________</p>
      <p>Each of domains in spatial zone is an abstract object and is described by a set of properties that characterize its
position in space, etc. (Fig. 2).
CenterX, CenterY – center domain coordinates; GeofonesNumber – sensor number fixing this domain; Sensitivity –
sensitivity; Xmax, Xmin, Ymax, Ymin, Z – domain boundary points coordinate</p>
      <p>The algorithm for calculating sensitivity for each of domains of selected spatial zone is represented by the
following steps:
geophone sensitivity limit;
1. The division of the spatial zone into domain regions of cubic form and their initialization;
2. Execution of a double cycle: external performs passage through the domains, internal - through geophones.
Following operations are performed at each iteration:
– calculation of the distance between domain center and sensor coordinates;
– increasing geophone counter, which can location event in a specific domain;
– increase in advance value of the relative sensitivity, if calculated in the previous operation distance is less than
– zeroing domain sensitivity if geophone counter is insufficient (less than 4), or adjusting domain sensitivity
parameter taking into account primary transducer number;
3. Removal of domains with zero sensitivity.</p>
      <p>The value characterizing domain sensitivity is calculated using following formula:</p>
      <p>0 =  ∙ ∑ = (1 − √ ⁄ д),
where  – geophone number detecting seismo acoustic event in the domain;</p>
      <p>– distance between domain center and sensor coordinates, m;
 д – sensor sensitivity limit, m.</p>
      <p>The developed algorithm is implemented in C # and integrated into the GeoAcoustics-ADS. It is seismic acoustic
data processing and analysis software package, which is one of the most important components of the Prognoz ADS
system.</p>
      <p>Window of this program module is called from main menu of GeoAcoustics-ADS (Fig. 3).</p>
      <p>The possibility of cooperative graphic display of sensitivity zone calculation results and mine plans is implemented
for visual presentation.</p>
      <p>Golden Software Surfer software is used as auxiliary visualization tools. It provides advanced features for
visualizing data sets, including by running a script written in VisualBasic. Setting the script parameters is provided in
developed module interface.</p>
      <p>The sensitivity map display of spatial domain provides for two main modes.</p>
      <p>Mode 1. This mode allows you to select a certain maximum value that characterizes limiting value of domain
sensitivity parameter of mine or its area. Next, a user-preset color scheme is loaded that contains level set of
percentages and their corresponding colors, represented in RGB model. Interpolation is performed on discrete set of
sensitivity values based on specified parameters, during which intermediate color values are determined.</p>
      <p>This mode is intended for:
– analysis of geophone distribution within zone boundaries;
– identify maximum sensitivity level, as well as the sensitivity levels of particular areas.</p>
      <p>Current situation in one of controlled areas of mining enterprise with active mining operations, which was formed
using gradient transitions in the first mode, is shown in Fig. 4.
Mode 2. This mode is designed to divide entire geometric zone into 3 characteristic areas:
1. Uncontrolled area. Acoustic events can be located after improving the characteristics of geophones, or after
increasing their number at a sufficient distance from source in this control zone;
2. Non-guaranteed control area. Acoustic events can be located with accuracy that is on limit of permissible
technical characteristics of monitoring system in this control zone;</p>
      <p>3. Zone of guaranteed control. Acoustic events can be located with accuracy according to technical characteristics
of monitoring system and with sufficient veracity of their parametric description in this zone.</p>
      <p>The plan of site control zone considered earlier is combined with a sensitivity map and is shown in figure 5
(mode 2). Gradient fill level values are selected empirically based on test data analysis and are 0, 5 and 15 units for
the uncontrolled zone, the non-guaranteed control zone and the guaranteed control zone, accordingly.
This mode allows you to:
– to split the research area into three separate zones;
– to make qualitative evaluation of efficiency of installed geophones;
– to check prone to shock hazard mine field areas for inclusion in the zone of confident location.</p>
      <p>Following method is proposed for choosing number and placements of geophones and expanding controlled area
in order to increase shock hazard control efficiency:</p>
      <p>1. If the analysis is performed for the first time, it is necessary to use the first mode and define maximum
sensitivity level. This limit will characterize the research area zone. Next, it is necessary to estimate the useful
information loss about seismic events by conducting industrial test series, or by similar research on an adjacent region
with similar mining and geological conditions. The maximum sensitivity reaches about 40 units for situation in Fig. 4.
As an example, select following sensitivity levels of control zones:
– below 5 - uncontrolled zone;
– 5 - 10 - non-guaranteed control zone;
– above 10 - zone of guaranteed control
2. Monitored region coordinates are determined sequentially. Then they are checked for entry into the guaranteed
control zone using the second mode. Take for testing point with coordinates 850, 350 from WTO 32 region (BDO).
This region belongs to the non-guaranteed control zone, as follows from situation presented in Fig. 5.</p>
      <p>3. The situation of sensitivity domain zone redistribution is simulated for each of the areas that do not meet the
conditions of guaranteed control by addition, removal or movement of geophones. The final sensitivity map for the
considered example after introduction of new geophones is shown in Fig. 6.</p>
      <p>It should be noted that possibility of installing geophones must be confirmed in the rock massif selected location
on basis of enterprise technical conditions in each specific case.
3</p>
    </sec>
    <sec id="sec-3">
      <title>Developed Model Limitations and Directions for Further Improvement</title>
      <p>The proposed algorithm for constructing sensitivity maps for considered controlled zone and method for choosing
the optimal installation geophone placements in mine field provide important practical problem to be solved.
However, further improvement and modernization is required due to accepted limitations of design model. Let us
consider some possible directions for further development of this approach.</p>
      <p>1. Anisotropic rock massif properties.</p>
      <p>
        Environment anisotropic properties, caused by physic mechanical properties of rocks and existence of empty
spaces in seismic acoustic waves propagation analysis from the source to the geophone, should be taken into account
[
        <xref ref-type="bibr" rid="ref10 ref11">10, 11</xref>
        ]. It is assumed that entire research area is composed of rocks and its environment is isotropic. Such statement
greatly simplifies the design model. It requires creation of more advanced mathematical description of domains to
solve this problem. For example, by using the velocity model created using results of experimental measurements.
2. Volumetric characteristics of receiving antennas.
      </p>
      <p>
        The decrease in event location accuracy should be taken into account if it was recorded by geophones lying in
plane (flat antenna) [
        <xref ref-type="bibr" rid="ref12">12</xref>
        ]. In such case, it is necessary to take into account flat antenna situation and prevent such
sensor placements, or introduce additional sensitivity correction and characterizing antenna volume coefficient.
3. Automate optimal geophone placements.
      </p>
      <p>It is necessary to change from manual selection of geophone installation places to automatic one, due to the
subjective judgment of designing person. At initial stage, developed algorithm should allow the most efficient
geophone installation in automatic mode. Then it is should consider possibility of selecting sensitivity levels for
guaranteed and non-guaranteed control zones without user influence.</p>
      <p>4. Use of three-dimensional domain sensitivity mathematical model.</p>
      <p>It is proposed to use a two-dimensional geometric model of structure of the domain when calculating sensitivity
maps in the developed method. This restriction is due to the existence of only flat plans for the mine field horizons.
However, use of sufficient number of mine horizons plans for different elevations will allow to calculate
pseudothree-dimensional models with possibility of obtaining output information in volumetric representation. It is
necessary to consider possibility of creating module that allows to interact with existing three-dimensional geometric
models on controlled objects in order to obtain more complete and easy-to-understand information about controlled
object.
4</p>
    </sec>
    <sec id="sec-4">
      <title>Conclusion</title>
      <p>The developed algorithm for constructing sensitivity maps allows one to get represent about geophone placement
effectiveness in conditions of working seismic acoustic monitoring system, or to modeling zone parameters change in
new configuration.</p>
      <p>The developed method of selecting optimal location for installing geophones in mine field allows to formulate
recommendations for optimal monitoring system monitoring zone extension based on model stages results of updated
underground observation network configurations.</p>
      <p>The directions for further improvement of proposed approach are formulated and options for a possible solution of
model limitations are proposed.</p>
      <p>The algorithm for constructing sensitivity maps for controlled zone and method of choosing optimal geophone
installation placements in mine field allow analyzing sensitivity zones of seismo-acoustic monitoring system, provide
detailed represent of influence of geophone location on monitoring zone spatial dimensions and suggest mathematical
justification for optimal geophone location.</p>
    </sec>
  </body>
  <back>
    <ref-list>
      <ref id="ref1">
        <mixed-citation>
          1.
          <string-name>
            <surname>Rasskazov</surname>
            ,
            <given-names>I. Ju.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Saksin</surname>
            ,
            <given-names>B. G.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Petrov</surname>
            ,
            <given-names>V. A.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Prosekin</surname>
            ,
            <given-names>B. A.</given-names>
          </string-name>
          :
          <article-title>Geomechanical conditions and features of the dynamic manifestations of rock pressure at the Antey deposits, Physical and technical problems in the development of mineral resources 3</article-title>
          :
          <fpage>3</fpage>
          -
          <lpage>13</lpage>
          (
          <year>2012</year>
          )
        </mixed-citation>
      </ref>
      <ref id="ref2">
        <mixed-citation>
          2.
          <string-name>
            <surname>Manchao</surname>
            ,
            <given-names>H.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Fuqiang</surname>
            ,
            <given-names>R.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Dongqiao</surname>
            ,
            <given-names>L.</given-names>
          </string-name>
          :
          <article-title>Rockburst mechanism research and its control</article-title>
          .
          <source>In: International Journal of Mining Science and Technology 28</source>
          .5:
          <fpage>829</fpage>
          -
          <lpage>837</lpage>
          (
          <year>2018</year>
          )
        </mixed-citation>
      </ref>
      <ref id="ref3">
        <mixed-citation>
          3.
          <string-name>
            <surname>Rasskazov</surname>
            ,
            <given-names>I. Ju.</given-names>
          </string-name>
          :
          <article-title>Control and management of mining pressure in the mines of the Far Eastern region</article-title>
          , M .: Gornaya Kniga Publishing House,
          <volume>329</volume>
          (
          <year>2008</year>
          )
        </mixed-citation>
      </ref>
      <ref id="ref4">
        <mixed-citation>
          4.
          <string-name>
            <surname>Hu</surname>
          </string-name>
          ,
          <string-name>
            <surname>Sh</surname>
          </string-name>
          .-Ch.,
          <string-name>
            <surname>Tan</surname>
            ,
            <given-names>Y.-L.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Ning</surname>
          </string-name>
          , J.-G.,
          <string-name>
            <surname>Guo</surname>
          </string-name>
          , W.-Y.,
          <string-name>
            <surname>Liu</surname>
          </string-name>
          , X.
          <article-title>-Sh.: Multiparameter Monitoring and Prevention of FaultSlip Rock Burst</article-title>
          .
          <source>In: Shock and Vibration</source>
          <year>2017</year>
          .
          <volume>7580109</volume>
          :
          <issue>8</issue>
          (
          <issue>2017</issue>
          ), https://doi.org/10.1155/
          <year>2017</year>
          /7580109
        </mixed-citation>
      </ref>
      <ref id="ref5">
        <mixed-citation>
          5.
          <string-name>
            <surname>Gladyr</surname>
            ,
            <given-names>A. V.</given-names>
          </string-name>
          :
          <article-title>System of integration of microseismic and geoacoustic data of geomechanical control</article-title>
          ,
          <source>Mining information and analytical bulletin</source>
          <volume>6</volume>
          :
          <fpage>220</fpage>
          -
          <lpage>234</lpage>
          (
          <year>2017</year>
          )
        </mixed-citation>
      </ref>
      <ref id="ref6">
        <mixed-citation>
          6.
          <string-name>
            <surname>Shabarov</surname>
            ,
            <given-names>A. N.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Cirel</surname>
            ,
            <given-names>S. V.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Morozov</surname>
            ,
            <given-names>K. V.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Rasskazov</surname>
            ,
            <given-names>I. Ju.:</given-names>
          </string-name>
          <article-title>The concept of integrated geodynamic monitoring of underground mining</article-title>
          ,
          <source>Physical and technical problems of mining</source>
          <volume>9</volume>
          :
          <fpage>59</fpage>
          -
          <lpage>64</lpage>
          (
          <year>2017</year>
          )
        </mixed-citation>
      </ref>
      <ref id="ref7">
        <mixed-citation>
          7.
          <string-name>
            <surname>Konstantinov</surname>
            ,
            <given-names>A.</given-names>
          </string-name>
          :
          <article-title>Development of Multi-Channel Portable Impact Control Device for Local Assessment of The State of The Edge Parts of The Rock Massif</article-title>
          , VII International Scientific Conference Problems of Complex Development of Georesources. Khabarovsk,
          <year>02024</year>
          (
          <year>2018</year>
          )
        </mixed-citation>
      </ref>
      <ref id="ref8">
        <mixed-citation>
          8.
          <string-name>
            <surname>Rasskazov</surname>
            ,
            <given-names>M. I.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Anikin</surname>
            ,
            <given-names>P. A.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Migunov</surname>
            ,
            <given-names>D. S.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Terjoshkin</surname>
            ,
            <given-names>A. A.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Gladyr</surname>
            ,
            <given-names>A. V.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Coj</surname>
            ,
            <given-names>D. I.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Lomov</surname>
            ,
            <given-names>M. A.</given-names>
          </string-name>
          :
          <article-title>Automated Seismic Acoustic Pressure Monitoring System at the Underground Mine World, VIII International Scientific and Practical Conference Innovative Directions in the Design of Mining Enterprises: Geomechanical Support for Mining Design and Maintenance</article-title>
          .
          <source>St. Petersburg</source>
          ,
          <volume>95</volume>
          -
          <fpage>101</fpage>
          (
          <year>2017</year>
          )
        </mixed-citation>
      </ref>
      <ref id="ref9">
        <mixed-citation>
          9.
          <string-name>
            <surname>Rasskazov</surname>
            ,
            <given-names>M. I.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Gladyr</surname>
            ,
            <given-names>A. V.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Katella</surname>
            ,
            <given-names>E. A.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Miroshnikov</surname>
            ,
            <given-names>V. I.</given-names>
          </string-name>
          :
          <article-title>Development of a multi-level system for complex geodynamic monitoring of the impact hazard of a rock mass</article-title>
          ,
          <source>Mining Informational Analytical Bulletin</source>
          <volume>9</volume>
          :
          <fpage>113</fpage>
          -
          <lpage>120</lpage>
          (
          <year>2014</year>
          )
        </mixed-citation>
      </ref>
      <ref id="ref10">
        <mixed-citation>
          10.
          <string-name>
            <surname>Vernik</surname>
            ,
            <given-names>L.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Nur</surname>
            ,
            <given-names>A.</given-names>
          </string-name>
          :
          <article-title>Ultrasonic velocity and anisotropy of hydrocarbon source rocks</article-title>
          ,
          <source>Geophysics</source>
          <volume>57</volume>
          (
          <year>1992</year>
          ), doi: 10.1190/1.1443286
        </mixed-citation>
      </ref>
      <ref id="ref11">
        <mixed-citation>
          11.
          <string-name>
            <surname>Vernik</surname>
            ,
            <given-names>L.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Liu</surname>
            ,
            <given-names>X.</given-names>
          </string-name>
          :
          <article-title>Velocity anisotropy in shales: A petrophysical study</article-title>
          ,
          <source>Geophysics</source>
          <volume>62</volume>
          (
          <year>1997</year>
          ), doi: 10.1190/1.1444162
        </mixed-citation>
      </ref>
      <ref id="ref12">
        <mixed-citation>
          12.
          <string-name>
            <surname>Rasskazov</surname>
            ,
            <given-names>I. Ju.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Iskra</surname>
            ,
            <given-names>A. Ju.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Kjanno</surname>
            ,
            <given-names>A. A.</given-names>
          </string-name>
          :
          <article-title>Algorithms and software for location of sources of acoustic emission in the system of geomechanical monitoring</article-title>
          ,
          <source>Mining information and analytical bulletin</source>
          <volume>15</volume>
          :
          <fpage>130</fpage>
          -
          <lpage>142</lpage>
          (
          <year>2007</year>
          )
        </mixed-citation>
      </ref>
    </ref-list>
  </back>
</article>