GEO-ACOUSTIC EVENT LOCATION METHOD AND INSTABILITY DISASTER EARLY WARNING METHOD BASED ON SAME, GEO-ACOUSTIC SENSOR, MONITORING SYSTEM, AND READABLE STORAGE MEDIUM

Abstract

Sound signal when a wave generated by a geo-acoustic event source reaches any monitoring point (S1), constructing a theoretical propagation difference model and an observed propagation difference model of the waveform characterization quantity between monitoring points, to calculate a waveform characterization quantity difference value between two monitoring points (S2); and constructing an objective function based on the theoretical propagation difference model and the observed propagation difference model, and obtaining the location of the geo-acoustic event by means of inversion based on the objective function (S3). According to the geo-acoustic event location method, the arrival time, time domain parameters, spectral information, and waveform shape of the geo-acoustic signal when the wave generated by the geo-acoustic event source reaches any monitoring point are considered, then the non-uniformity of a propagation medium is comprehensively reflected, and the inversion precision of geo-acoustic event location is finally improved.

Claims

1. A geo-acoustic event location method, comprising the following steps: S1: acquiring an observed value of a waveform characterization quantity of a geo-acoustic event, the waveform characterization quantity comprising: any combination of an arrival time, time domain parameters, spectral information, and waveform shape of a geo-acoustic signal when a wave generated by a geo-acoustic event source reaches any monitoring point; S2: constructing a theoretical propagation difference model and an observed propagation difference model of the waveform characterization quantity between monitoring points, wherein a waveform characterization quantity difference value between two monitoring points is calculated; and S3: constructing an objective function based on the theoretical propagation difference model and the observed propagation difference model, and obtaining a location of the geo-acoustic event by inversion based on the objective function.

2. The geo-acoustic event location method according to claim 1, wherein the objective function is as follows:
G(x.sub.p,y.sub.p,z.sub.p, . . . )=minΣ(DA.sub.xyz+DF.sub.xyz+DS.sub.xyz+DT.sub.xyz) wherein G(x.sub.p,y.sub.p,z.sub.p, . . . ) indicates the objective function of the geo-acoustic event (x.sub.p, y.sub.p, z.sub.p), x.sub.p, y.sub.p, z.sub.p is a position coordinate of the geo-acoustic event, and DF.sub.xyz, DA.sub.xyz, DS.sub.xyz and DT.sub.xyz indicate a degree of deviation of an observed quantity and the a theoretical value of a propagation difference of the waveform characterization quantity, respectively.

3. The geo-acoustic event location method according to claim 2, wherein formulas for calculating DF.sub.xyz, DA.sub.xyz, DS.sub.xyz and DT.sub.xyz are respectively as follows:
DT.sub.xyz=Σ(Δt.sub.xyx.sup.nm−Δt.sub.0.sup.nm).sup.2
DA.sub.xyz=Σ(ΔA.sub.xyx.sup.nm−ΔA.sub.0.sup.nm).sup.2
DF.sub.xyz=Σ(Δf.sub.xyx.sup.nm−Δf.sub.0.sup.nm).sup.2
DS.sub.xyz=Σ(ΔS.sub.xyx.sup.nm−ΔS.sub.0.sup.nm).sup.2 wherein Δt.sub.xyx.sup.nm and Δt.sub.0.sup.nm respectively indicate a theoretical arrival time difference and the an observed quantity of the arrival time difference between two different monitoring points m and n; ΔA.sub.xyx.sup.nm and ΔA.sub.0.sup.nm respectively indicate a theoretical attenuation difference and an observed value of an attenuation difference of the time domain parameters a.sub.0 between two different monitoring points m and n; Δf.sub.xyx.sup.nm and Δf.sub.0.sup.nm respectively indicate a theoretical difference value and an observed difference value of the spectral information f between two different monitoring points m and n; and ΔS.sub.xyx.sup.nm and ΔS.sub.0.sup.nm respectively indicate a theoretical value of a waveform shape evolution difference and an observed value of the waveform shape evolution difference between two different monitoring points m and n.

4. The geo-acoustic event location method according to claim 1, wherein the time domain parameters comprise any one or any combination of a rise time, a duration, an amplitude and energy time domain parameters of a waveform, and the spectral information is any combination of all frequency information of the waveform after Fourier decomposition; and the waveform shape is a function, wherein the function defines a shape of a geo-acoustic signal waveform in a time domain.

5. An instability disaster early warning method based on the geo-acoustic event location method according to claim 1, comprising the following steps: step 1: forming self-clustering clusters according to a location result of the geo-acoustic event, the self-clustering clusters being obtained by clustering a spatial distribution of the geo-acoustic event; step 2: counting precursor index information of each of the self-clustering clusters according to a pre-set time window, wherein the precursor index information is any combination of indexes of the a relative energy, a moment magnitude, an apparent volume, b value of each of the self-clustering clusters, a time series change of a rate of geo-acoustic events, a dominant frequency and a corner frequency of the geo-acoustic event in each of the self-clustering clusters; and if a sum of a category of precursor indexes of all the geo-acoustic events in each of the self-clustering clusters under one time window becomes larger compared with a sum of a same category of precursor indexes under a previous time window, a precursor factor of the corresponding precursor index is 1; and step 3: performing instability disaster early warning analysis based on precursor factors of precursor indexes of each of the self-clustering clusters, wherein a risk increases as a sum of the precursor factors of all the precursor indexes of each of the self-clustering clusters increases.

6. The instability disaster early warning method according to claim 5, wherein the precursor indexes comprise the indexes of the relative energy, the moment magnitude, the apparent volume, the b value of each of the self-clustering clusters, the time series change of the rate of geo-acoustic events, the dominant frequency and the corner frequency of the geo-acoustic event in each of the self-clustering clusters and the precursor indexes are divided into three main categories; and wherein the relative energy, the moment magnitude, the apparent volume and the b value of each of the self-clustering clusters are categorized into a first main category, the time series change of the rate of geo-acoustic events is categorized into a second main category, the dominant frequency and the corner frequency are categorized into a third main category, the precursor factors respectively corresponding to the first main category, the second main category and the third main category are expressed as P.sub.1, P.sub.2 and P.sub.3, and an early warning rule based on the precursor factors is as follows: if P.sub.1+P.sub.2+P.sub.3=3, it indicates that operation in an area is prohibited; if P.sub.1+P.sub.2+P.sub.3=2, it indicates that a hidden danger of a potential source of danger requires to be removed; if P.sub.1+P.sub.2+P.sub.3=1, it indicates that operation is allowed to be carried out under real-time monitoring; and if P.sub.1+P.sub.2+P.sub.3=0, it indicates safe operation.

7. An intelligent geo-acoustic sensor, comprising a geo-acoustic sensing unit, an automatic gain and acoustic-electric conversion unit, an RC filtering network unit, an analog-to-digital conversion unit, an intelligent sensing filter unit and a photoelectric conversion unit, wherein the geo-acoustic sensing unit, the automatic gain and acoustic-electric conversion unit, the RC filtering network unit, the analog-to-digital conversion unit, the intelligent sensing filter unit and the photoelectric conversion unit are connected to each other in a communication mode; wherein the geo-acoustic sensing unit acquires multi-band geo-acoustic signals; the automatic gain and acoustic-electric conversion unit is configured for converting a geo-acoustic signal into an analog signal, and performing initial adjustment on the analog signal to realize a fast response of a weak signal gain; the RC filtering network unit is configured for truncating the analog signal in a frequency domain to complete preliminary filtering; the analog-to-digital conversion unit is configured for converting the analog signal into a digital signal; the intelligent sensing filter unit is configured for performing denoising and filtering processing on the digital signal; and the photoelectric conversion unit is configured for converting the digital signal into an optical signal, wherein a signal acquired by the intelligent geo-acoustic sensor is used for locating a geo-acoustic event according to the steps of the geo-acoustic event location method according to claim 1.

8. The intelligent geo-acoustic sensor according to claim 7, wherein a metal full-shielding frame is disposed outside the intelligent geo-acoustic sensor.

9. A geo-acoustic event monitoring system, comprising a geo-acoustic signal sensing module, a geo-acoustic event location module and a geo-acoustic event instability disaster early warning module, wherein the geo-acoustic signal sensing module, the geo-acoustic event location module and the geo-acoustic event instability disaster early warning module are connected in a communication mode, wherein the geo-acoustic signal sensing module is configured for sensing geo-acoustic signals; the geo-acoustic event location module is configured for locating geo-acoustic events by using the geo-acoustic event location method according to claim 1; and the geo-acoustic event instability disaster early warning module is configured for performing instability disaster early warning analysis by using an instability disaster early warning method based on the geo-acoustic event location method according to claim 1, wherein the instability disaster early warning method comprises the following steps: step 1: forming self-clustering clusters according to a location result of the geo-acoustic event, the self-clustering clusters being obtained by clustering a spatial distribution of the geo-acoustic event step 2: counting precursor index information of each of the self-clustering clusters according to a pre-set time window, wherein the precursor index information is any combination of indexes of the a relative energy, a moment magnitude, an apparent volume, b value of each of the self-clustering clusters, a time series change of a rate of geo-acoustic events, a dominant frequency and a corner frequency of the geo-acoustic event in each of the self-clustering clusters; and if a sum of a category of precursor indexes of all the geo-acoustic events in each of the self-clustering clusters under one time window becomes larger compared with a sum of a same category of precursor indexes under a previous time window, a precursor factor of the corresponding precursor index is 1; and step 3: performing instability disaster early warning analysis based on precursor factors of precursor indexes of each of the self-clustering clusters, wherein a risk increases as a sum of the precursor factors of all the precursor indexes of each of the self-clustering clusters increases.

10. A readable storage medium, storing a computer program, the computer program being called by a processor to execute the geo-acoustic event location method according to claim 1 or an instability disaster early warning method based on the geo-acoustic event location method according to claim 1, wherein the instability disaster early warning method comprises the following steps: step 1: forming self-clustering clusters according to a location result of the geo-acoustic event, the self-clustering clusters being obtained by clustering a spatial distribution of the geo-acoustic event; step 2: counting precursor index information of each of the self-clustering clusters according to a pre-set time window, wherein the precursor index information is any combination of indexes of the a relative energy, a moment magnitude, an apparent volume, b value of each of the self-clustering clusters, a time series change of a rate of geo-acoustic events, a dominant frequency and a corner frequency of the geo-acoustic event in each of the self-clustering clusters; and if a sum of a category of precursor indexes of all the geo-acoustic events in each of the self-clustering clusters under one time window becomes larger compared with a sum of a same category of precursor indexes under a previous time window, a precursor factor of the corresponding precursor index is 1; and step 3: performing instability disaster early warning analysis based on precursor factors of precursor indexes of each of the self-clustering clusters, wherein a risk increases as a sum of the precursor factors of all the precursor indexes of each of the self-clustering clusters increases.

11. The instability disaster early warning method according to claim 5, wherein in the geo-acoustic event location method, the objective function is as follows:
G(x.sub.p, y.sub.p, z.sub.p, . . . )=minΣ(DA.sub.xyz+DF.sub.xyz+DS.sub.xyz+DT.sub.xyz) wherein G(x.sub.p, y.sub.p, z.sub.p, . . . ) indicates the objective function of the geo-acoustic event (x.sub.p, y.sub.p, z.sub.p), x.sub.p, y.sub.p, z.sub.p is a position coordinate of the geo-acoustic event, and DF.sub.xyz, DA.sub.xyz, DS.sub.xyz and DT.sub.xyz indicate a degree of deviation of an observed quantity and a theoretical value of a propagation difference of the waveform characterization quantity, respectively.

12. The instability disaster early warning method according to claim 11, wherein in the geo-acoustic event location method, formulas for calculating DF, DA.sub.xyz, DS.sub.xyz and DT.sub.xyz are respectively as follows:
DT.sub.xyz=Σ(Δt.sub.xyx.sup.nm−Δt.sub.0.sup.nm).sup.2
DA.sub.xyz=Σ(ΔA.sub.xyx.sup.nm−ΔA.sub.0.sup.nm).sup.2
DF.sub.xyz=Σ(Δf.sub.xyx.sup.nm−Δf.sub.0.sup.nm).sup.2
DS.sub.xyz=Σ(ΔS.sub.xyx.sup.nm−ΔS.sub.0.sup.nm).sup.2 wherein Δt.sub.xyx.sup.nm and Δt.sub.0.sup.nm respectively indicate a theoretical arrival time difference and an observed quantity of the arrival time difference between two different monitoring points m and n; ΔA.sub.xyx.sup.nm and ΔA.sub.0.sup.nm respectively indicate a theoretical attenuation difference and an observed value of an attenuation difference of the time domain parameters a.sub.0 between two different monitoring points m and n; Δf.sub.xyx.sup.nm and Δf.sub.0.sup.nm respectively indicate a theoretical difference value and an observed difference value of the spectral information f between two different monitoring points m and n; and ΔS.sub.xyx.sup.nm and ΔS.sub.0.sup.nm respectively indicate a theoretical value of a waveform shape evolution difference and an observed value of the waveform shape evolution difference between two different monitoring points m and n.

13. The instability disaster early warning method according to claim 5, wherein in the geo-acoustic event location method, the time domain parameters comprise any one or any combination of a rise time, a duration, an amplitude and energy time domain parameters of a waveform, and the spectral information is any combination of all frequency information of the waveform after Fourier decomposition; and the waveform shape is a function, wherein the function defines a shape of a geo-acoustic signal waveform in a time domain.

14. The intelligent geo-acoustic sensor according to claim 7, wherein in the geo-acoustic event location method, the objective function is as follows:
G(x.sub.p, y.sub.p,z.sub.p, . . . )=minΣ(DA.sub.xyz+DF.sub.xyz+DS.sub.xyz+DT.sub.xyz) wherein G(x.sub.p, y.sub.p,z.sub.p, . . . ) indicates the objective function of the geo-acoustic event (x.sub.p, y.sub.p, z.sub.p), x.sub.p, y.sub.p, z.sub.p is a position coordinate of the geo-acoustic event, and DF.sub.xyz, DA.sub.xyz, DS.sub.xyz and DT.sub.xyz indicate a degree of deviation of an observed quantity and a theoretical value of a propagation difference of the waveform characterization quantity, respectively.

15. The intelligent geo-acoustic sensor according to claim 14, wherein in the geo-acoustic event location method, formulas for calculating DF.sub.xyz, DA.sub.xyz, DS.sub.xyz and DT.sub.xyz are respectively as follows:
DT.sub.xyz=Σ(Δt.sub.xyx.sup.nm−Δt.sub.0.sup.nm).sup.2
DA.sub.xyz=Σ(ΔA.sub.xyx.sup.nm−ΔA.sub.0.sup.nm).sup.2
DF.sub.xyz=Σ(Δf.sub.xyx.sup.nm−Δf.sub.0.sup.nm).sup.2
DS.sub.xyz=Σ(ΔS.sub.xyx.sup.nm−ΔS.sub.0.sup.nm).sup.2 wherein Δt.sub.xyx.sup.nm and Δt.sub.0.sup.nm respectively indicate a theoretical arrival time difference and an observed quantity of the arrival time difference between two different monitoring points m and n; ΔA.sub.xyx.sup.nm and ΔA.sub.0.sup.nm respectively indicate a theoretical attenuation difference and an observed value of an attenuation difference of the time domain parameters a.sub.0 between two different monitoring points m and n; Δf.sub.xyx.sup.nm and Δf.sub.0.sup.nm respectively indicate a theoretical difference value and an observed difference value of the spectral information f between two different monitoring points m and n; and ΔS.sub.xyx.sup.nm and ΔS.sub.0.sup.nm respectively indicate a theoretical value of a waveform shape evolution difference and an observed value of the waveform shape evolution difference between two different monitoring points m and n.

16. The intelligent geo-acoustic sensor according to claim 7, wherein in the geo-acoustic event location method, the time domain parameters comprise any one or any combination of a rise time, a duration, an amplitude and energy time domain parameters of a waveform, and the spectral information is any combination of all frequency information of the waveform after Fourier decomposition; and the waveform shape is a function, wherein the function defines a shape of a geo-acoustic signal waveform in a time domain.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0061] FIG. 1 is a schematic view of an application process provided by the present disclosure; and

[0062] FIG. 2 is a structural schematic view of an intelligent geo-acoustic sensor provided by the present disclosure.

DETAILED DESCRIPTION OF THE EMBODIMENTS

[0063] The present disclosure is further described below with reference to the embodiments.

Example 1

[0064] Scientific research technicians have been trying to explore effective ways to predict geo-acoustic disaster instability, and a large number of traditional monitoring methods such as displacement, stress, vibration, strain and microseism have achieved some results at a certain level. However, most of the existing monitoring methods are point monitoring, and microseismic monitoring is regional coverage monitoring, it is difficult to capture tiny precursor characteristic information in the instability and sub-instability stages of the geotechnical project, there is an urgent need to invent a high-sensitivity sensing apparatus capturing the occurrence of small deformation, small rupture and small collapse of rock and soil, and an intelligent sensor capable of eliminating the superimposed noise in the high-sensitivity monitoring process so as to acquire the change information in the instability evolution process of the rock as a whole. For this purpose, the present example provides an intelligent multi-band geo-acoustic sensor as shown in FIG. 1.

[0065] As shown in FIG. 1, the intelligent geo-acoustic sensor provided by the present example includes a geo-acoustic sensing unit, an automatic gain and acoustic-electric conversion unit, an RC filtering network unit, an analog-to-digital conversion unit, an intelligent sensing filter unit, a photoelectric conversion unit and a power supply module which are connected to each other in a communication mode.

[0066] The geo-acoustic sensing unit is equipped with high-sensitivity, high-bandwidth and high-specification vibration sensing elements, such as: piezoelectric PZT ceramics to cover and capture geo-acoustic signals at different scales and different frequency bands in rock.

[0067] The automatic gain and acoustic-electric conversion unit is equipped with AGC automatic gain control, which is specifically configured for dynamically adjusting the gain according to the intensity of vibration signals and controlling the output signal level. Firstly, the geo-acoustic signal is converted into an analog signal, and then initial adjustment is performed on the analog signal to realize a fast response of a weak signal gain so that a wide gain range is provided for the intelligent sensor.

[0068] The RC filtering network unit sets the truncating frequency for acquiring the geo-acoustic signal and completes preliminary filtering by truncating the analog signal in the frequency domain. The analog-to-digital conversion unit is configured for converting the analog signal into a digital signal.

[0069] The intelligent sensing filter unit can be equipped with FPGA or DSP according to engineering practice and apparatus production conditions; and can also be configured with modules of low pass, high pass, band pass, band stop, Kalman filtering, wavelet filtering and the like for filtering and denoising. The type and number of specifically configured filtering modules are set according to actual requirements, and are not specifically defined in the present disclosure.

[0070] In general, the intelligent sensing filter unit is configured for eliminating the noise of field apparatuses such as fans, drilling rigs and crushers, mingled in the signal, magnetic field interference as well as the temperature drift and zero offset generated inside the geo-acoustic sensing device under certain conditions. In the case that an effective signal and a noise signal are not superimposed in the frequency domain, specific frequency bands of some noise are filtered out directly by using the low-pass, high-pass, band-pass and band-stop filtering modules. In the superimposition of Gaussian white noise caused by the noise signal and pulse interference with the effective signal in the frequency domain, Kalman filtering and wavelet filtering are used to remove peaks and abrupt changes in the acquired signal.

[0071] The photoelectric conversion unit is configured for converting the digital signal into an optical signal. The photoelectric conversion unit can resist electromagnetic interference in special environment and long-distance transmission, which realizes high-performance fidelity.

[0072] In the present example, the intelligent sensing filter unit can be configured for effectively sensing a useful signal, removing interference noise and ensuring the fidelity effect of the signal. On the other hand, a weak signal can be amplified by using AGC automatic gain control of the automatic gain and acoustic-electric conversion unit, and thus the intelligent geo-acoustic sensor provided by the present disclosure has the effects of high sensitivity, wide sensing frequency and intelligent filtering.

[0073] The present example also preferably disposes an outer metal full-shielding frame outside the intelligent geo-acoustic sensor to further block external interference. The power supply module provided in the present example is configured with a high transient absorption or lightning protection device which is controlled to mitigate the influence caused by the interference of an internal transient electric pulse. In other feasible examples, this is not required in the present disclosure.

Example 2

[0074] The present example provides a geo-acoustic event location method, including the following steps:

[0075] S1: acquiring an observed value of a waveform characterization quantity of a geo-acoustic event.

[0076] The waveform characterization quantity selected in the present example includes the arrival time, time domain parameters, spectral information, and waveform shape of a geo-acoustic signal when a wave generated by a geo-acoustic event source P.sub.xyz reaches any monitoring point. In other feasible examples, characteristics may be combined as the waveform characterization quantity according to actual requirements, which is not specifically required in the present disclosure.

[0077] For example, observed quantities in the present example when the wave generated by the geo-acoustic event source P.sub.xyz reaches any monitoring point m include: observed arrival time t.sub.0.sup.m, observed time domain parameters a.sub.0.sup.m, observed spectral information f.sub.0.sup.m, and an observed waveform shape S.sub.0.sup.m.

[0078] Preferably, the time domain parameters include characteristic parameters of the rise time, duration, amplitude and energy of the waveform, the spectral information is a combination of all frequency information of the waveform after Fourier decomposition, and the waveform shape is a function that defines the shape of a geo-acoustic signal waveform in a time domain. In other feasible examples, the selection of parameters may be adjusted adaptively according to locating precision requirements.

[0079] S2: constructing a theoretical propagation difference model and an observed propagation difference model of the waveform characterization quantity between monitoring points. A waveform characterization quantity difference value between two monitoring points is calculated.

Theoretical Propagation Difference Model:

[0080] Based on the theoretical geo-acoustic travel time model, geo-acoustic time domain parameter attenuation model, spectral information dissipation model, and waveform shape evolution model, the propagation time when an acoustic wave generated by a geo-acoustic event source P.sub.xyz reaches any monitoring point m is t.sub.xyz.sup.m, the geo-acoustic time domain parameter attenuation quantity is a.sub.xyz.sup.m, the spectral information dissipation quantity is f.sub.xyz.sup.m, and the waveform shape evolution quantity is S.sub.xyz.sup.m. The theoretical values of the above parameters can all be determined using the prior art, e.g. geophysical model inversion.

[0081] Then, the theoretical arrival time difference Δt.sub.xyz.sup.m, the time domain parameter attenuation difference Δa.sub.xyz.sup.nm, the spectral information dissipation difference Δf.sub.xyz.sup.nm and the waveform shape evolution difference Δs.sub.xyz.sup.nm between two different monitoring points m and n are respectively noted as:

Δ.sub.xyz.sup.nm=t.sub.xyz.sup.m−t.sub.xyz.sup.n

Δa.sub.xyz.sup.nm=a.sub.xyz.sup.m−a.sub.xyz.sup.n

Δf.sub.xyz.sup.nm=f.sub.xyz.sup.m−f.sub.xyz.sup.n

ΔS.sub.xyz.sup.nm=S.sub.xyz.sup.m=S.sub.xyz.sup.n

Observed Propagation Difference Model:

[0082] Similarly, the observed propagation time difference Δt.sub.0.sup.nm between two different monitoring points m and n, the observed time domain parameters Δa.sub.0.sup.nm, the observed spectral information Δf.sub.0.sup.nm, and the waveform shape ΔS.sub.0.sup.nm are respectively noted as:

Δt.sub.0.sup.nm=t.sub.0.sup.m−t.sub.0.sup.n

Δa.sub.0.sup.nm=a.sub.0.sup.m−a.sub.0.sup.n

Δf.sub.0.sup.nm=f.sub.0.sup.m−f.sub.0.sup.n

ΔS.sub.0.sup.nm=S.sub.0.sup.m−S.sub.0.sup.n

[0083] S3: constructing an objective function based on the theoretical propagation difference model and the observed propagation difference model, and obtaining the location of the geo-acoustic event by means of inversion based on the objective function.

[0084] The difference values DT.sub.xyz, DA.sub.xyz, DF.sub.xyz and DS.sub.xyz between each observed value and theoretical value are respectively used for describing the degree of deviation between the spatial location of the geo-acoustic events under different theoretical propagation models and the actual spatial location of the geo-acoustic events, and are respectively noted as:

DT.sub.xyz=Σ(Δt.sub.xyx.sup.nm−Δt.sub.0.sup.nm).sup.2

DA.sub.xyz=Σ(ΔA.sub.xyx.sup.nm−ΔA.sub.0.sup.nm).sup.2

DF.sub.xyz=Σ(Δf.sub.xyx.sup.nm−Δf.sub.0.sup.nm).sup.2

DS.sub.xyz=Σ(ΔS.sub.xyx.sup.nm−ΔS.sub.0.sup.nm).sup.2

[0085] the objective function is constructed as:

G(x.sub.p, y.sub.p, z.sub.p, . . . )=minΣ(DA.sub.xyz+DF.sub.xyz+DS.sub.xyz+DT.sub.xyz)

[0086] the above objective function is a non-linear fitting problem, and the position (x.sub.p, y.sub.p, z.sub.p) of the geo-acoustic event source P.sub.xyz can be found by solving the above system of equations by the existing computational means.

[0087] It should be understood that the objective function provided in the present example is set based on a selected waveform characterization quantity, and in other feasible examples, different combinations of waveform characterization quantities are chosen for use and then the objective function is adaptively modified.

[0088] It should be noted that, the present disclosure preferably uses the intelligent geo-acoustic sensor described in Example 1 to acquire data, and uses the geo-acoustic event location method to perform space location inversion.

Example 3

[0089] Instability disaster early warning of the rock and soil structure has been a world-wide problem for a long time, and is limited by the one-sided early warning criterion, narrow scale coverage of index information and the like. It is difficult for the existing technical means to adapt to the engineering application scenario of initial fracture and extensional instability in complex multi-field environment. It can be seen that based on the multi-type fused monitoring device of the high-sensitivity sensing technology, by relying on multi-physical field and multi-scale comprehensive monitoring indexes, seeking suitable and proper entry points for early warning indexes, and constructing an intelligent discrimination system for the geo-acoustic disaster risk are of great engineering value. For this purpose, the present example provides an instability disaster early warning method based on the geo-acoustic event location method, including the following steps:

[0090] step 1: forming self-clustering clusters according to the location result of a geo-acoustic event, the self-clustering cluster being obtained by clustering the spatial distribution of the geo-acoustic event.

[0091] In general, one self-clustering cluster corresponds to a potential instability disaster, and a high-density region of the self-clustering cluster corresponds to the source of potential instability disasters.

[0092] step 2: counting precursor index information of each self-clustering cluster according to a pre-set time window.

[0093] In the present example, the precursor indexes include the indexes of the relative energy, moment magnitude, apparent volume, b value of the self-clustering cluster, time series change of the rate of geo-acoustic events, dominant frequency and corner frequency of the geo-acoustic event in the self-clustering cluster; and in other feasible examples, the precursor indexes may be other combinations of the above indexes, and the present disclosure is not specifically limited in this respect.

[0094] Specifically, the relative energy E.sub.r is calculated as follows:

E.sub.r=Σ.sub.1.sup.ndE/n=Σ.sub.1.sup.nd∫f(t).sup.2dt/n

[0095] in the formula, n indicates the number of stations; f(t) is the speed amplitude of a geo-acoustic signal; t is duration of the signal; d is a distance attenuation factor, given by d=L.sup.qm/L.sup.pm, L.sup.qm is a theoretical propagation path distance from the center q of the cluster to the monitoring point m, and L.sup.qm is a theoretical propagation path distance from the geo-acoustic event q to the monitoring point m.

[0096] The moment magnitude M.sub.w characterizes the strength level of a geo-acoustic event.

[0097] The apparent volume V characterizes the total volume of rock that is plastically deformed by a geo-acoustic event.

[0098] The b value characterizes the relationship of the frequency distribution of the moment magnitude of the geo-acoustic event in the self-clustering cluster over one time window.

[0099] The event rate R characterizes the rate of geo-acoustic events generated in the self-clustering cluster over one time window.

[0100] The dominant frequency f.sub.d characterizes the dominant frequency of the geo-acoustic event signal in one self-clustering cluster.

[0101] The corner frequency f.sub.c value characterizes the corner frequency of the geo-acoustic event signal in one self-clustering cluster.

[0102] In order to effectively identify whether a hidden danger exists, the present disclosure determines the precursor factor corresponding to the precursor index according to the following rules.

[0103] P.sub.M indicates the sum of the moment magnitudes M.sub.w of all geo-acoustic event signals generated in the self-clustering cluster over one time window, and when the sum of the moment magnitudes increases compared with that over a previous time window, a precursor factor P.sub.M corresponding to the moment magnitude M.sub.w is noted as 1.

[0104] P.sub.E indicates the sum of the relative energy E.sub.r of all geo-acoustic event signals generated in the self-clustering cluster over one time window, and when the sum of the relative energy increases compared with that over a previous time window, a precursor factor P.sub.E corresponding to the relative energy E.sub.r is noted as 1.

[0105] P.sub.A indicates the sum of the apparent volume V of all geo-acoustic event signals generated in the self-clustering cluster over one time window, and when the sum of the apparent volume increases compared with that over a previous time window, a precursor factor P.sub.A corresponding to the apparent volume V is noted as 1.

[0106] P.sub.b indicates a decrease in the b value over one time window compared with the b value over a previous time window, and is noted as 1.

[0107] P.sub.R indicates the change in the event rate R over one time window compared with that over a previous time window, and is noted as 1.

[0108] P.sub.f.sub.d indicates a frequency band formed by counting the dominant frequency f.sub.d of all geo-acoustic event signals generated in the self-clustering cluster over one time window, and is noted as 1 when there is a change in the number of dominant frequency bands compared with dominant frequency bands over a previous time window.

[0109] P.sub.f.sub.c indicates the frequency band formed by counting the corner frequencies f.sub.c of all the geo-acoustic event signals generated in the self-clustering cluster in a time window, and is noted as 1 when there is a change in the number of corner frequency bands compared to corner frequency bands over a previous time window.

[0110] Further, the above-mentioned indexes are divided into three main categories in the present example, as follows:

[0111] the relative energy, moment magnitude, apparent volume and b value of the self-clustering cluster are categorized into the first main category, the time series change of the rate of geo-acoustic events is categorized into the second main category, the dominant frequency and corner frequency are categorized into the third main category, and the corresponding precursor factors thereof are expressed as P.sub.1, P.sub.2 and P.sub.3, they are:

P.sub.1=P.sub.M+P.sub.E+P.sub.b+P.sub.A

P.sub.2=P.sub.R

P.sub.3=P.sub.f.sub.d+P.sub.f.sub.c

[0112] step 3: performing instability disaster early warning analysis based on precursor factors of precursor indexes of each self-clustering cluster. The greater the sum of the precursor factors of all the precursor indexes of the self-clustering cluster is, the greater the risk is.

[0113] It should be understood that the greater the sum of P.sub.1+P.sub.2+P.sub.3 is, the higher the risk is; and the less the sum is, the higher the safety is. Thus, the present disclosure sets rule indexes for actual operation based on experience research, as follows:

[0114] if P.sub.1+P.sub.2+P.sub.3=3, a red level, indicating that operation in the area is prohibited;

[0115] if P.sub.1+P.sub.2+P.sub.3=2, a yellow level, indicating that the hidden danger of a potential source of danger requires to be removed;

[0116] P.sub.1+P.sub.2+P.sub.3=1, a blue level, indicating that operation can be carried out under real-time monitoring; and

[0117] if P.sub.1+P.sub.2+P.sub.3=0, a green level, indicating safe operation.

[0118] The needs of the present disclosure can also be met without classification in other feasible examples, and the present disclosure is not specifically limited in this respect.

[0119] In other feasible examples, the precursor index also includes an S value, and the S value characterizes the sum of the frequency, average energy level, and maximum energy level composition of a geo-acoustic event in a self-clustering cluster over one time window. P.sub.S indicates that the S value over one time window increases compared with the S value over the previous single time window, and is noted as 1. In the present example, it is preferable that P.sub.S indicates that the S value over a single time window in the time series shows a continuous increasing trend and the value is at a high level, and noted as 1.

[0120] In other feasible examples, the relative energy, moment magnitude and apparent volume are categorized into the first main-category precursor factor P.sub.1 (seismic source parameter), the b value, the S value and the rate of geo-acoustic events is the second main-category precursor factor P.sub.2 (time series), the dominant frequency and the corner frequency are the third main-category precursor factor P.sub.3 (spectral characteristics), and they have a Boolean operation relationship. Then, in step 3, instability disaster early warning analysis is performed based on the precursor factor formed by the precursor indexes of each self-clustering cluster according to the Boolean operation relationship.

[0121] The time window is set based on the experience or need, for example, when historical training sample data is available, the time window with the best early warning efficiency is taken, and the precursor factors are counted based on this. When no historical training sample data is available, the default time window length is set to 1 month, and the interval time counted in a mobile mode is 7 days, with no specific limitations in other feasible examples.

[0122] It should be noted that, preferably, in the present example, the early warning method described in Example 3 is used to guide operations based on the location of the sound source determined in Example 2. As shown in FIG. 1, Examples 1-3 of the present disclosure constitute a monitoring system for geo-acoustic events, integrating geo-acoustic signal acquisition, geo-acoustic event location and geo-acoustic event instability disaster early warning. The monitoring system includes a geo-acoustic signal sensing module, a geo-acoustic event location module and a geo-acoustic event instability disaster early warning module connected in a communication mode, and the geo-acoustic signal sensing module is configured for sensing geo-acoustic signals and acquiring relevant data for locating subsequent geo-acoustic events. The geo-acoustic event location module is configured for locating the geo-acoustic event, and reference may be made to the description of Example 2 described above for the specific implementation process. The geo-acoustic event instability disaster early warning module is configured for performing instability disaster early warning analysis, and reference may be made to the description of Example 3 described above for the specific implementation process.

[0123] It should be understood that the division of the above-mentioned functional modules or units is merely division of logical function, and other division modes may exist during actual implementation. For example, a plurality of units or components may be combined or integrated into another system, or some characteristics may be ignored or not performed. At the same time, the above integrated unit can be implemented in the form of hardware, e.g. when the geo-acoustic signal sensing module is implemented in the form of hardware, it corresponds to the intelligent geo-acoustic sensor described in Example 1 of the present disclosure. The above integrated unit can also be implemented in the form of a software function unit.

[0124] In some feasible modes, the present disclosure further provides a readable storage medium storing a computer program, the computer program being called by a processor to execute:

[0125] S1: acquiring an observed value of a waveform characterization quantity of a geo-acoustic event;

[0126] S2: constructing a theoretical propagation difference model and an observed propagation difference model of the waveform characterization quantity between monitoring points; and

[0127] S3: constructing an objective function based on the theoretical propagation difference model and the observed propagation difference model, and obtaining the location of the geo-acoustic event by means of inversion based on the objective function;

[0128] or to execute:

[0129] step 1: forming self-clustering clusters according to the location result of a geo-acoustic event, the self-clustering cluster being obtained by clustering the spatial distribution of the geo-acoustic event;

[0130] step 2: counting precursor index information of each self-clustering cluster according to a pre-set time window; and

[0131] step 3: performing instability disaster early warning analysis based on precursor factors of precursor indexes of each self-clustering cluster. The greater the sum of the precursor factors of all the precursor indexes of the self-clustering cluster is, the greater the risk is.

[0132] Reference may be made to the description of the method for the specific implementation process of each step, which will not be repeated here in the present disclosure.

[0133] In some feasible modes, the present disclosure further provides a system, including a processor and a memory, the memory storing a computer program, the processor calling the computer program to execute:

[0134] S1: acquiring an observed value of a waveform characterization quantity of a geo-acoustic event;

[0135] S2: constructing a theoretical propagation difference model and an observed propagation difference model of the waveform characterization quantity between monitoring points; and

[0136] S3: constructing an objective function based on the theoretical propagation difference model and the observed propagation difference model, and obtaining the location of the geo-acoustic event by means of inversion based on the objective function;

[0137] or to execute:

[0138] step 1: forming self-clustering clusters according to the location result of a geo-acoustic event, the self-clustering cluster being obtained by clustering the spatial distribution of the geo-acoustic event;

[0139] step 2: counting precursor index information of each self-clustering cluster according to a pre-set time window; and

[0140] step 3: performing instability disaster early warning analysis based on precursor factors of precursor indexes of each self-clustering cluster. The greater the sum of the precursor factors of all the precursor indexes of the self-clustering cluster is, the greater the risk is.

[0141] For specific implementation processes of each step, reference may be made to the description of the foregoing method, which will not be repeated in the present disclosure.

[0142] It should be understood that, in this embodiment of the present disclosure, the mentioned processor may be a central processing unit (Central Processing Unit, CPU), or the processor may be another general purpose processor, a digital signal processor (Digital Signal Processor, DSP), an application-specific integrated circuit (Application-Specific Integrated Circuit, ASIC), a field-programmable gate array (Field-Programmable Gate Array, FPGA) or another programmable logic device, a discrete gate or a transistor logic device, a discrete hardware component, or the like. The general-purpose processor may be a microprocessor, or the processor may be any conventional processor and the like. The memory may include a read-only memory (ROM) and a random access memory (RAM), and provide instructions and data to the processor. A part of the memory may further include a non-volatile random access memory. For example, the memory may further store information of a device type.

[0143] The readable storage medium is a computer-readable storage medium, which may be an internal storage unit of the controller described in any of the foregoing embodiments, such as a hard disk or a memory of the controller. The readable storage medium may also be an external storage device of the controller, such as a plug-in hard disk, a smart media card (SMC), a secure digital (SD) card, or a flash card that is equipped on the controller. Further, the readable storage medium may further include an internal storage unit of the controller and an external storage device. the readable storage medium is configured to store the computer program and another program and data that are required by the controller. The readable storage medium may be further configured to temporarily store data that has been outputted or data to be outputted.

[0144] It should be emphasized that, examples described in the present disclosure are illustrative rather than restrictive, and therefore the present disclosure is not limited to the examples described in the specific implementations. All other implementations obtained by a person skilled in the art according to the technical solutions of the present disclosure do not depart from the spirit and scope of the present disclosure, and whether modified or replaced, also belong to the protection scope of the present disclosure.