CN108693561B - Coal mining seismic acquisition system and method based on wave detector subdivision array - Google Patents

Abstract

The present invention relates to a kind of coal mining seismic acquisition systems and method based on wave detector subdivision array.Working face is divided into multiple divided areas using the two neighboring seismic detector string being set in deep hole by the present invention, effective seismic signal will be converted into from cutting the continuous shaking noise intraocularly transmitted using multidimensional deconvolution interferometry, " void " big gun collection of high s/n ratio, high spatial resolution is obtained in divided area.

Description

Coal mining seismic acquisition system and method based on wave detector subdivision array

Technical field

The present invention relates to a kind of coal mining seismic acquisition system and methods, belong to coal mine intelligence production technique field, tool Body is to be related to a kind of coal mining seismic acquisition system and method based on wave detector subdivision array.

Background technique

The key technology that coal mine is intelligently exploited is intelligent detecting, intelligent navigation and intelligent control, wherein " intelligent detecting " belongs to In geophysical exploration field, it is most important technology in intelligence exploitation, is also most difficult to realize.

Intelligence exploitation is a kind of unmanned operation, recovery scheme and various operating parameters all judged automatically by artificial intelligence and It completes, many influences little craven fault, the variation of coal thickness etc. in people-machine mining system can all make in automatic intelligent exploitation At bigger influence.

The tradition method of seismic prospecting such as Three Dimensional Ground earthquake, underground slot wave can provide figure for the construction in working face Picture, however be limited to seismic observation system unfavorable factors, these images such as distant at a distance from detection target and intelligence is also far not achieved Precision required for capable of exploiting.

Intelligence exploitation not only needs the detailed configuration information in front of working face, the information such as stress evolution of top plate and bottom plate It requires to obtain in time, to adjust mining operations scheme in time, achievees the purpose that safe and efficient coal mining, common seismic prospecting Method is all static instrumentation method, is unable to complete this real-time detection task.

Summary of the invention

The technical problems to be solved by the invention

Although there are some in the prior art is able to carry out real-time detection with adopting detection method, it is all along working face two Side roadway layout geophone, using coalcutter as passive source, the seismic ray of this method does not cover entire work Face, and the angle of cut is little each other between seismic ray, is unfavorable for inversion imaging.In addition, the geophone offset of this method is bigger, it is unfavorable for visiting Micrometer little structure.

In view of the above problems, the present invention is visited using the subdivision receiving array that geophone group in multiple groups lateral aperture synthesizes It surveys, so that coal working face is subdivided into multiple lesser search coverages, reaches the coal mine intelligently fine detection purpose of exploitation.

The present invention utilizes seismic interference side using the intraocular coal-cutting machinery noise of working face open-off cut, microseism etc. as passive source Method converts these passive source energy to " void " the big gun collection received in network, realizes that the in-built piecemeal of working face is finely visited It surveys.

Solve scheme used by technical problem

Above-mentioned technical problem of the invention is mainly to be addressed by following technical proposals:

A kind of coal mining seismic acquisition system based on wave detector subdivision array, comprising:

Seismic detector string array, including multiple seismic detector strings, each seismic detector string are located at work Make in the parallel deep hole arranged in face, the parallel deep hole both ends are separately connected tunnel；

Data collector connects several detector strings, and logical shielded cable accesses optical transmitter and receiver；

Wherein, working face is divided into a divided area by two neighboring seismic detector string, dry using multidimensional deconvolution Effective seismic signal will be converted into from cutting the continuous shaking noise intraocularly transmitted by relating to method, obtained in divided area high s/n ratio, " void " big gun collection of high spatial resolution.

In at least one embodiment of the present invention, further includes: the shallow bore hole wave detector in the tunnel.

In at least one embodiment of the present invention, the wave detector is placed in pushing device, and the pushing device includes For placing the shell of wave detector, a hinged backup arm on the shell, the hinged backup screw rod in the middle part of the backup arm, institute The other end for stating backup screw rod is connected on the stroke lever of a backup motor；Also, institute is driven when the stroke lever transverse shifting Stating backup screw rod drives at least one end of the backup arm to deviate the shell.

A kind of coal mining seismic detection method based on wave detector subdivision array, comprising:

Working face is divided by region division step using the seismic detector string being arranged in working face in parallel deep hole Multiple divided areas；

Data collection steps are obtained using the intraocular coal-cutting machinery noise of working face open-off cut, microseism as passive source from cutting eye The continuous shaking noise inside transmitted；It is converted into effective seismic signal；

Data processing step will be converted into effectively from cutting the continuous shaking noise intraocularly transmitted as unit of divided area Seismic signal obtains " void " big gun collection of high s/n ratio, high spatial resolution in divided area；

Data interpretation procedure carries out geologic interpretation based on described " void " big gun collection.

In at least one embodiment of the present invention, in the data collection steps, seismic detection after coal work stops Device string delays scheduled duration closing, and the subsequent microseismic signals induced by coal work are recorded in delay periods.

In at least one embodiment of the present invention, it in the data processing step, is executed in each divided area Following sub-step:

Observation system modeling sub-step, the coordinate in input service face and each wave detector, to establish observation system, in real time more The data that fresh cut eye position promotes provide foundation for data explanation；

Working condition determines sub-step, when the kurtosis of analysing energy level and cross-correlation function is to filter out coal mining activity The data of section, interception coal mining period carry out subsequent processing；

The anti-pleat of multidimensional interferes sub-step, will be turned from cutting the continuous shaking noise intraocularly transmitted using multidimensional deconvolution interferometry Effective seismic signal is turned to, " void " big gun collection of high s/n ratio, high spatial resolution is obtained in divided area；

Coal exploration constructs sub-step, is handled the CT obtained segmenting in array, slot wave frequency dispersion song to empty big gun collection is transmitted One of line, comprehensive coal thickness variation are a variety of；

Stress variation detect sub-step, between former and later two coal cuttings single time " void " big gun collection carry out coda wave interference and handle with Obtain stress variation figure.

In at least one embodiment of the present invention, the working condition determines in sub-step, carries out to the amplitude of data Analysis is judged as the non-coal cutting of coalcutter when amplitude is less than threshold value；When amplitude is greater than threshold value, data are grouped, In The kurtosis that cross-correlation function is sought in group is determined as coalcutter just in coal cutting, in data screening file when kurtosis is greater than a threshold value The middle record coal cutting period.

In at least one embodiment of the present invention, in the anti-pleat interference sub-step of the multidimensional,

It will be defined as one group of image shotpoint near the string for cutting eye in divided area, will be far from the string for cutting eye It is defined as receiving point, the seismic data of the seismic data of each image shotpoint and all receiving points carries out cross-correlation and obtains initial empty big gun Initial empty big gun collection is arranged successively by image shotpoint as a Multidimensional numerical by collection；

The seismic data of the seismic data of each image shotpoint and other image shotpoints is done into cross-correlation and obtains the point expansion of image shotpoint Function is dissipated, point spread function is arranged successively by image shotpoint as another Multidimensional numerical；Utilize point spread function Multidimensional numerical pair Initial void big gun collection Multidimensional numerical does multidimensional deconvolution processing, obtains final empty big gun collection Multidimensional numerical.

In at least one embodiment of the present invention, in the coal exploration construction sub-step, ART or SIRT method is utilized Inverting obtains the CT in subdivision array；

And/or using reflection seismic Prestack Migration Technology handle reflected waveform data, obtain search coverage in tomography at Picture；

And/or the slot wave dispersion curve of every empty big gun collection seismic ray is extracted, inverting obtains the thickness of the coal on every ray, will The coal thickness informix of every ray gets up, and is changed using the synthesis coal thickness that chromatography imaging method obtains each region.

In at least one embodiment of the present invention, in the data interpretation procedure,

By analysis CT image, reflected image, slot wave dispersion curve to be quantified to tomography, karst collapse col umn, the variation of coal thickness Mark；

And/or it is based on stress variation dynamic image, combined structure weakness band and safety monitoring information, predict geological disaster；

And/or tomography, karst collapse col umn, the variation of coal thickness are drawn a circle to approve according to structural exploration result and presented with three-dimensional visualization.

Invention effect

Therefore, the present invention has the advantage that

1, a series of a kind of method that working face is integrally detected and is converted into divided areas and finely detects, detection essence are provided Degree is much higher than conventional method, is very suitable to the fine intelligent detecting for solving the problems, such as intelligently to exploit.

2, using be parallel to the drilling for cutting eye lay string, and usually in tunnel receive seismic signal with fief Shake method is compared, and what is obtained is " void " big gun collection between wave detector, and distance measurement target is closer, and " void " shot point and receiving point are all It immobilizes, the precision of detection pressure stress variation is higher.

3, coda wave interference detection essence is compared using the detection of coda wave interference technique stress variation and bulk wave, slot wave etc. between hole Du Genggao, it is very sensitive to the velocity of wave variation of medium because coda wave is a kind of wave of Multiple Scattering in medium；Coda wave is dry between hole Relating to is subdivision section detection, can further improve spatial resolution.

Detailed description of the invention

Coal working face and wave detector arrangement schematic diagram Fig. 1 of the invention；

Piecemeal void big gun collection ray Fig. 2 of the invention covers schematic diagram；

Deep hole string schematic diagram Fig. 3 of the invention；

Keyway hanging-connecting means for jointing rods Fig. 4 of the invention.

In figure, wave detector 1, backup motor 2, backup back 3, backup screw rod 4, key 5, slot 6.

Specific embodiment

Below with reference to the embodiments and with reference to the accompanying drawing the technical solutions of the present invention will be further described.

Embodiment:

Geophone group in multiple groups lateral aperture is combined into subdivision receiving array by the present embodiment, and coal working face is subdivided into one A small search coverage reaches the coal mine intelligently fine detection purpose of exploitation.

The present embodiment utilizes seismic interference using the intraocular coal-cutting machinery noise of working face open-off cut, microseism etc. as passive source Method converts these passive source energy to " void " the big gun collection received in network, realizes that the in-built piecemeal of working face is fine Detection.

With reference to the accompanying drawing, this implementation is further illustrated.

As shown in Figure 1, the present embodiment arranges a series of parallel deep hole in working face, hole and working face to cut eye mutually equal Row, about 100 meters of pitch of holes (detection accuracy requires higher distance between borehole smaller).Seismic detector string, tunnel spacing are arranged in hole 10 meters, the string in multiple holes collectively constitutes wave detector subdivision array, and working face is divided into lesser detection one by one Region.Cut intraocular coalcutter, scratch board conveyor mine when can issue seismic wave, Sub-Level Caving, microseism etc. can also issue earthquake Wave, these seismic energies constitute the passive source of this method.All these seismic waves in coal mining activity continuously to Working face front and wave detector Internet communication.It is converted into using the seismic energy that seismic interference technology continuously mixes these similar " void " big gun collection of dynamite source can be used as wave detector shape in " void " focus, with adjacent holes away from wave detector in the close hole of eye is cut It is right at " big gun-inspection ".The detection result of divided area between a series of this available holes of method, because of the big gun-of this " empty big gun " Away from closer, the working face between device to hole of having ready conditions finely is detected for inspection.

The present embodiment has studied the seismic data processing technique based on wave detector subdivision array in hole, including multidimensional deconvolution Chromatography imaging technique, coda wave interference detection method of coal rock layer stress variation etc. between interference technique, hole.

Seismic interference is a kind of data re-establishing method, can by geophone array record environmental noise, mechanicalness noise, The energy of the passive sources such as a large amount of microseisms occurred at random is converted into the Green's function between different wave detectors, i.e. some of them is examined Wave device is as shot point, and other wave detectors still act as wave detector, thus this reconstructed Green's function is also referred to as " empty big gun Collection ", it is not the earthquake big gun collection really blown out, but essentially identical with true big gun collection.

Multidimensional deconvolution interference is a kind of new method of seismic interference, more practical than traditional cross-correlation interference technique, Traditional cross-correlation interference technique requires passive source uniform ring around detector array, and multidimensional deconvolution interference adapts to more one As hypocenter distributing, for example in the present invention, hypocenter distributing is intraocular in cutting, and cannot surround detector array, and hypocenter distributing is excessively Rule is easy to produce glitch using traditional seismic interference method.

Tomography is a kind of Detection Techniques of maturation between hole, usually utilizes 2 parallel holes, installs detection in a hole Device string gradually puts a column big gun in another hole, forms sufficient ray alternate covering to two interpore regions, finally utilizes layer It analyses imaging technique and inversion imaging is carried out to ray overlay area.The present invention carries out hole interbed using the empty big gun collection that seismic interference obtains Analysis imaging, thus string is all arranged in two holes, it is not only easy for construction, moreover it is possible to repeated multiple times detection.

Coda wave interference detection is a kind of new method detected using Seismic Coda to the faint variation of propagation medium, tail Wave and direct wave, back wave etc. are different, it is a kind of seismic wave wave that Multiple Scattering is formed in medium, since Multiple Scattering passes Path length is broadcast, is occurred the latest on seismic chart, thus is visually referred to as coda wave.Because being formed by Multiple Scattering, coda wave Apparent velocity is significantly less than direct wave, is changed using the very indistinguishable media properties of direct wave, is then readily recognized using coda wave. Coda wave interference, which refers to, carries out interference processing to the coda wave that front and back records twice, obtains the variable quantity of coda wave then, these variable quantities With regard to having reacted variation of the medium during front and back detects twice.

The technical solution adopted in the present invention include in lateral aperture wave detector subdivision array design installation, data acquisition with Four transmission, data processing, data interpretation parts, are introduced separately below.

1, the design of wave detector subdivision array is installed in lateral aperture

Since it is desired that with detection is adopted, it is desirable that wave detector pipe nipple can fit closely for a long time with hole wall in hole, and the present embodiment is set Mechanical arm pushing device is counted, backup powers off at once after stablizing, and backup arm still can close backup, and this design for a long time It can also easily recycle, equipment can be drawn back from hole as long as connecting electricity recycling mechanical arm, as shown in Figure 3；Detection in lateral aperture Device can only generally lean on manpower to push into hole under mine, and wave detector pipe nipple can be sent into deep hole using push rod push In, since lane space is limited, bar is designed to sectionally assembled, and the junction of bar is designed as keyway hanging-connecting, such as Fig. 4 institute Show, the advantages of this design is that scene assembly and disassembly speed of application is fast, moreover it is possible to increase the freedom degree between each bar, not only avoid casing The interference of wave can also make the more preferable by wall effect of wave detector pipe nipple.This string is all laid in multiple parallel holes, jointly Collectively form subdivision detector array.

2, data acquisition and transmission

Working face both sides tunnel be parallel to cut eye one group of deep hole in arrange wave detector (Fig. 1): pacify in parallel deep hole Installation tool backup string, wave detector are bonded close (Fig. 3) using mechanical backup arm and hole wall, are to push away between wave detector pipe nipple Connecting rod is sent, every bar is 1.2 meters long, and bar directly utilizes key-slot hanging device to connect (Fig. 4) with bar, and not only installation and dismantling are very convenient, This hanging device can also effectively eliminate casing wave, and bar and the loosely connected of bar can also allow the backup of wave detector adherent more Closely, the road spacing of deep hole wave detector is 10 meters.Cable is through in the cavity of connecting rod, uses thunder with the connector of wave detector pipe nipple Not water-proof connector plugs convenient and waterproof infiltration.Every 6 wave detector accesses a data collector, and more number of units are according to acquisition Device accesses an optical transmitter and receiver by shielded cable, and more optical transmitter and receivers transmit data by optical fiber to the ground in real time.Coal work is opened System for acquiring seismic data opens work simultaneously when the beginning, and seismic acquisition configuration is delayed 20 minutes and closed after coal work stops, In The recordable subsequent microseismic signals induced by coal work of delay periods.

It in addition to parallel deep hole geophone group, is also arranged auxiliary geophone array (Fig. 1) in working face both sides tunnel, tunnel inspection Wave device is mounted on coal wall in the hole of 2 meters of depths, using spring backup or cementitious；Reduce the interference of tunnel noise, road to the greatest extent Spacing is 10 meters；When coal work starts, system for acquiring seismic data also starts simultaneously；The seismic data of each acquisition station exists Underground converges to optical transmitter and receiver, then by optical fiber real-time Transmission to the ground, handles in central station storage.

3, data processing:

1) observation system, input service face, the coordinate for cutting eye and each wave detector are established；Observation system is initially set up, it is defeated The coordinate for entering working face and each wave detector, is ready for data processing, and subsequent processing wants the measurement accuracy of space coordinate Ask higher；It cuts the data that a position promotes and is also required to real-time update, provide important evidence for data explanation；As coal work pushes away Into some distances are cut wave detector of the eye less than 20 meters and will be removed；

2) it determines working condition, can directly be judged using seismic data instantly whether in coal mining activity, analysing energy water The kurtosis of gentle cross-correlation function can relatively accurately judge, seismic energy levels are high when coal mining activity, and each track data is mutual The kurtosis of pass is also big；Hash is rejected, only the data of interception coal mining period carry out subsequent processing, and data can be greatly decreased in this The workload of processing.It is specific good, 4 can be extracted from the data that hole 1 records, the amplitude of data be analyzed, when amplitude is small It is judged as coalcutter non-coal cutting when empirical value, carries out next step analysis when being greater than threshold value, is divided to two groups for 4 track datas, often Group 2 does cross-correlation, seeks the kurtosis of cross-correlation function, be determined as when kurtosis is greater than empirical value coalcutter just in coal cutting, The coal cutting period is recorded in data screening file；

3) multidimensional deconvolution is interfered:

Coal cutting period data is chosen to be handled；As unit of each deep hole divided area, data screening is divided into more A unit, the wave detector of small size deep hole is defined as image shotpoint in each unit, remaining wave detector is defined as receiving point；Utilize multidimensional Deconvolution interference technique is effective seismic signal is converted into from cutting the continuous shaking noise intraocularly transmitted, in divided area one by one Interior " void " big gun collection for obtaining high s/n ratio, high spatial resolution；

Specifically, one group of image shotpoint will be defined as near No. 1 hole string for cutting eye, No. 2 holes and two sides tunnel Wave detector is defined as receiving point, and the seismic data of the seismic data of each image shotpoint and all receiving points does cross-correlation and obtains initially Initial empty big gun collection is arranged successively by image shotpoint as a Multidimensional numerical by empty big gun collection；The seismic data of each image shotpoint and other The seismic data of image shotpoint does cross-correlation and obtains the point spread function of the image shotpoint, and point spread function is arranged successively by image shotpoint For another Multidimensional numerical；Initial empty big gun collection Multidimensional numerical is done at multidimensional deconvolution using this point spread function Multidimensional numerical Reason, obtains final empty big gun collection Multidimensional numerical；

4) detection construction: seismic interference has the characteristics that energy is stackable, coal seam rank of the coalcutter before exploiting No. 1 hole Section, all energy can be applied, and the signal-to-noise ratio of empty big gun collection can further increase after superposition every time；

With exploitation, continuous stack power obtains the empty big gun collection of the highest transmission of signal-to-noise ratio as far as possible, utilizes ART or the side SIRT Method inverting obtains the CT in subdivision array；Reflected waveform data is handled using the Prestack Migration Technology of reflection seismic, obtains detecting area Tomographic imaging in domain；The slot wave dispersion curve of every empty big gun collection seismic ray is extracted, inverting obtains the thickness of the coal on every ray, The coal thickness informix of every ray is got up, is changed using the synthesis coal thickness that chromatography imaging method obtains each region；

It is specific good, it, can using the covering ray between adjacent holes (such as between No. 1 hole and No. 2 holes or between No. 2 holes and No. 3 holes) To do three kinds of imagings: a, doing tomography to the interpore region of segmentation, common inversion algorithm has ART or SIRT；B, reflection is utilized Seismic method does catoptric imaging to the tomography in region, determines position of fault；C, the dispersion curve for extracting slot wave utilizes coal seam thickness The relationship of degree and slot wave dispersion curve, does thickness inversion imaging to coal seam；

5) detect stress variation: stress detection concerns that stress changes with time, and is former and later two in embodiment Stress variation (coal cutting single time refer to that coalcutter is cut from cutting one, a face to other one) between coal cutting single time；

" void " big gun collection that one singly plows will often be cut to extract, by the front and back of each big gun inspection pair, " void " big gun is drawn in one twice It rises, the seismic phase of front is almost consistent, and coda wave part was not only similar but also variant, and coda wave interference processing is carried out to it, obtains medium change Change, be mainly stress variation figure；It connects and plays to a series of stress variation figures, obtain stress variation cloud atlas.

It does multidimensional deconvolution interference for each coal cutting single time, obtains empty big gun collection, the coda wave of front and back twice is interfered Processing, obtains coda wave arrival time difference, and the coda wave arrival time difference that multiple coal cuttings are single time is concatenated, and obtains variation cloud atlas and has just reacted medium Stress variation；

6) between all carried out the adjacent parallel hole of every two it is above-mentioned 3)~5) processing of step, obtain the detection of entire working face As a result；Also it can simplify as only with two adjacent parallel holes, with mining, the propulsion rolling for cutting eye is detected forward, although so not Entire working face can be detected simultaneously, but can detect and cut within 100 meters of side at the moment, and advantage is can to save detecting devices.

4, data interpretation

By analyzing CT image, reflected image, slot wave dispersion curve, tomography, karst collapse col umn, the variation of coal thickness etc. are quantified Mark, is sent to intelligent mining system, and the further planning intelligently to exploit provides geological information；For utilizing coda wave interference side The stress variation dynamic image that method obtains, combined structure weakness band and other safety monitoring information judge the disaster that may occur, Anticipation is made to position and risk, anticipation result is sent to intelligent mining system.

Under conditions of abundant collection work face geologic information, it is thick that tomography, karst collapse col umn, coal are drawn a circle to approve according to structural exploration result Variation etc., is presented, data are transmitted to intelligent mining system with three-dimensional visualization；Stress variation " cloud atlas " is presented in the form of animation, number According to intelligent mining system and ground safety monitoring department is transmitted to, early warning judgement is made in due course.

Specific embodiment described herein is only an example for the spirit of the invention.The neck of technology belonging to the present invention The technical staff in domain can make various modifications or additions to the described embodiments or replace by a similar method In generation, however, it does not deviate from the spirit of the invention or beyond the scope of the appended claims.

Claims (10)

1. a kind of coal mining seismic acquisition system based on wave detector subdivision array characterized by comprising

Seismic detector string array, including multiple seismic detector strings, each seismic detector string are located at working face In the parallel deep hole of interior arrangement, the parallel deep hole both ends are separately connected tunnel；

Data collector connects several strings, and accesses optical transmitter and receiver by shielded cable；

Wherein, working face is divided into a divided area by two neighboring seismic detector string, utilizes multidimensional deconvolution interferometry Effective seismic signal will be converted into from cutting the continuous shaking noise intraocularly transmitted, high s/n ratio and high-altitude are obtained in divided area Between resolution ratio " void " big gun collection.

2. a kind of coal mining seismic acquisition system based on wave detector subdivision array according to claim 1, feature It is, further includes: the shallow bore hole wave detector in the tunnel.

3. a kind of coal mining seismic acquisition system based on wave detector subdivision array according to claim 1, feature It is, the wave detector is placed in pushing device, and the pushing device includes the shell for placing wave detector, on the shell A hinged backup arm, the hinged backup screw rod in the middle part of the backup arm, the other end of the backup screw rod are connected to a backup On the stroke lever of motor；Also, the backup screw rod is driven to drive the backup arm at least when the stroke lever transverse shifting Deviate the shell in one end.

4. a kind of coal mining seismic detection method based on wave detector subdivision array characterized by comprising

Working face is divided into multiple by region division step using the seismic detector string being arranged in working face in parallel deep hole Divided area；

Data collection steps obtain intraocular from cutting using the intraocular coal-cutting machinery noise of working face open-off cut or microseism as passive source The continuous shaking noise transmitted；It is converted into effective seismic signal；

Data processing step will be converted into effective earthquake from cutting the continuous shaking noise intraocularly transmitted as unit of divided area Signal obtains " void " big gun collection of high s/n ratio and high spatial resolution in divided area；

Data interpretation procedure carries out geologic interpretation based on described " void " big gun collection.

5. a kind of coal mining seismic detection method based on wave detector subdivision array according to claim 4, feature It is, in the data collection steps, seismic detector string delays scheduled duration closing after coal work stops, in delay periods Record the subsequent microseismic signals induced by coal work.

6. a kind of coal mining seismic detection method based on wave detector subdivision array according to claim 4, feature It is, in the data processing step, following sub-step is executed in each divided area:

Observation system models sub-step, and the coordinate in input service face and each wave detector, to establish observation system, real-time update is cut The data that eye position promotes provide foundation for data explanation；

Working condition determines sub-step, and the kurtosis of analysing energy level and cross-correlation function is cut with filtering out the coal mining activity period The data for period of mining are taken to carry out subsequent processing；

The anti-pleat of multidimensional interferes sub-step, will be converted into from cutting the continuous shaking noise intraocularly transmitted using multidimensional deconvolution interferometry Effective seismic signal obtains " void " big gun collection of high s/n ratio and high spatial resolution in divided area；

Coal exploration constructs sub-step, is handled to obtain CT in subdivision array, slot wave dispersion curve, comprehensive to empty big gun collection is transmitted Close one of coal thickness variation or a variety of；

Stress variation detects sub-step, carries out coda wave interference and processing to " void " big gun collection between former and later two coal cuttings single time to obtain Stress variation figure.

7. a kind of coal mining seismic detection method based on wave detector subdivision array according to claim 6, feature It is, the working condition determines in sub-step, analyzes the amplitude of data, is judged as adopting when amplitude is less than threshold value The non-coal cutting of coal machine；When amplitude is greater than threshold value, data are grouped, the kurtosis of cross-correlation function is sought in group, when kurtosis is big It is determined as that coalcutter just in coal cutting, records the coal cutting period in data screening file when a threshold value.

8. a kind of coal mining seismic detection method based on wave detector subdivision array according to claim 6, feature It is, in the anti-pleat interference sub-step of multidimensional,

It will be defined as one group of image shotpoint near the string for cutting eye in divided area, will be far from the string definition for cutting eye For receiving point, the seismic data of the seismic data of each image shotpoint and all receiving points carries out cross-correlation and obtains initial empty big gun collection, Initial empty big gun collection is arranged successively by image shotpoint as a Multidimensional numerical；

The seismic data of the seismic data of each image shotpoint and other image shotpoints is done into cross-correlation and obtains the point spread function of image shotpoint Number, point spread function is arranged successively by image shotpoint as another Multidimensional numerical；Using point spread function Multidimensional numerical to initial Empty big gun collection Multidimensional numerical does multidimensional deconvolution processing, obtains final empty big gun collection Multidimensional numerical.

9. a kind of coal mining seismic detection method based on wave detector subdivision array according to claim 6, feature It is, in the coal exploration construction sub-step, obtains the CT in subdivision array using ART or SIRT inversion method；

And/or reflected waveform data is handled using the Prestack Migration Technology of reflection seismic, obtain the tomographic imaging in search coverage；

And/or the slot wave dispersion curve of every empty big gun collection seismic ray is extracted, inverting obtains the thickness of the coal on every ray, will be every The coal thickness informix of ray gets up, and is changed using the synthesis coal thickness that chromatography imaging method obtains each region.

10. a kind of coal mining seismic detection method based on wave detector subdivision array according to claim 4, feature It is, in the data interpretation procedure,

By analysis CT image, reflected image and slot wave dispersion curve to carry out quantitative mark to the variation of tomography, karst collapse col umn and coal thickness Note；

And/or it is based on stress variation dynamic image, combined structure weakness band and safety monitoring information, predict geological disaster；

And/or tomography, karst collapse col umn and the variation of coal thickness are drawn a circle to approve according to structural exploration result and presented with three-dimensional visualization.