CN110531415B - Three-dimensional small fault advanced detection method utilizing influence of surrounding rock loosening ring - Google Patents

Three-dimensional small fault advanced detection method utilizing influence of surrounding rock loosening ring Download PDF

Info

Publication number
CN110531415B
CN110531415B CN201910773665.9A CN201910773665A CN110531415B CN 110531415 B CN110531415 B CN 110531415B CN 201910773665 A CN201910773665 A CN 201910773665A CN 110531415 B CN110531415 B CN 110531415B
Authority
CN
China
Prior art keywords
seismic
surrounding rock
small fault
wave signal
head
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201910773665.9A
Other languages
Chinese (zh)
Other versions
CN110531415A (en
Inventor
黄兰英
姜慧
王圣程
宋雪娟
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xuzhou University of Technology
Original Assignee
Xuzhou University of Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Xuzhou University of Technology filed Critical Xuzhou University of Technology
Priority to CN201910773665.9A priority Critical patent/CN110531415B/en
Publication of CN110531415A publication Critical patent/CN110531415A/en
Application granted granted Critical
Publication of CN110531415B publication Critical patent/CN110531415B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/30Analysis
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/36Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/50Corrections or adjustments related to wave propagation
    • G01V2210/51Migration
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/60Analysis
    • G01V2210/64Geostructures, e.g. in 3D data cubes
    • G01V2210/642Faults

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Acoustics & Sound (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

The invention discloses a three-dimensional small fault advanced detection method utilizing the influence of a surrounding rock loose circle, wherein a three-component detector is placed in the middle of the surrounding rock loose circle, reflected Love type groove waves can be converted into Love type surface waves by utilizing the surrounding rock loose circle, the Love type surface waves with strong amplitude can be received by the three-component detector and then fed back to a seismograph for recording, the seismograph performs offset imaging on seismic signals of the Love type surface waves by adopting the existing pre-stack Fresnel beam offset method according to the recorded Love type surface waves, and finally the position of a small fault in front of the head is obtained through an imaging graph. The invention has simple construction and convenient use, can accurately detect whether a small fault exists in front of the head-on and can determine the position of the small fault, thereby providing guidance for the safe tunneling of the subsequent coal roadway.

Description

Three-dimensional small fault advanced detection method utilizing influence of surrounding rock loosening ring
Technical Field
The invention relates to a coal mine advanced detection method, in particular to a three-dimensional small fault advanced detection method utilizing the influence of a surrounding rock loose ring.
Background
According to the statistics of China coal industry Association, 2525 outbursts of 3082 coal and gas outbursts are related to geological structures, and the proportion is as high as 81.9%. And the geological formations that cause the outburst are most typically small faults, such as small faults. At present, geophysical advanced detection methods are more, and include seismic wave methods, electromagnetic methods and other methods. Seismic methods are most suitable for advanced fault prediction, with trough wave exploration being the most typical of downhole coal mines.
The feasibility numerical simulation test of the advanced detection of the front small fault by the channel wave obtains remarkable results, which are mainly reflected in that: for example, the advanced detection method of Rayleigh type trough wave coal roadway small faults, which is excited and received by a top plate and a bottom plate, is proposed by Yangtze and the like (numerical simulation of Rayleigh type trough wave advanced detection of coal roadway small faults, 2012; advanced detection method of coal roadway reflection trough waves, invented patent application No. 201711007670.6, for a thick coal seam, which is tunneled along the bottom plate or the top plate); a tangcapone and the like (TVSP advanced detection method based on a trough wave, 2018) perform a trough wave advanced detection study and the like by using a Rayleigh type trough wave perpendicular to a top plate direction through three-dimensional numerical simulation. However, the prior art still has certain disadvantages: the traditional reflection advanced detection technology has obvious effect on fault with large fall and obvious fault reflection interface, but has no obvious effect on advanced detection of small fault; secondly, the channel waves can be divided into Love channel waves and Rayleigh channel waves, the excitation condition of the Rayleigh channel waves is harsh, and when the parameters of surrounding rocks and the coal bed do not meet the excitation condition of the Rayleigh channel waves, effective signals cannot be acquired; and thirdly, the actual coal roadway tunneling geological condition is complex, loose circle influence exists around the coal roadway, the surrounding rock loose circle influence is considered by the road development (the research of mine roadway surface wave-converted transverse wave advanced detection technology, 2016), and the accuracy of the trough wave advanced detection is improved by avoiding the interference of the loose circle on the reflector wave. In addition, a two-dimensional section is considered, the actual exploration condition is a three-dimensional condition, the section direction of the top plate and the bottom plate is surrounding rock, the plane direction is a coal bed, and the loose circles in the section direction and the plane direction have difference, so that the two-dimensional section direction has limitation. In order to avoid the influence of the surrounding rock loose ring on the reflected groove wave advanced detection, the detectors are arranged in the region outside the surrounding rock loose ring by means of drilling or anchor rod construction on the coal roadway wall, so that the influence of the surrounding rock loose ring on the groove wave exploration can be effectively reduced without the surrounding rock loose ring by the received seismic waves, but the method is very difficult to implement, and needs deep drilling, deep hole installation of the detectors and other steps, so that the forward small fault of the groove wave advanced detection cannot be conveniently and quickly realized.
Disclosure of Invention
Aiming at the problems in the prior art, the invention provides a three-dimensional small fault advanced detection method utilizing the influence of a surrounding rock loose ring, which is simple in construction and convenient and fast to use, can accurately detect whether a small fault exists in front of a head-on part or not and can determine the position of the small fault, and provides guidance for the safe tunneling of a subsequent coal roadway.
In order to achieve the purpose, the invention adopts the technical scheme that: a three-dimensional small fault advanced detection method utilizing the influence of a surrounding rock loosening ring comprises the following specific steps:
A. arranging a shock excitation point on the central axis of the left wall or the right wall of the roadway and close to the head-on position, and then acquiring the range of the known surrounding rock loosening circle;
B. arranging a plurality of three-component detectors in rows at the middle position of a surrounding rock loosening ring of the left side or the right side of the roadway in a drilling mode, and then connecting each three-component detector with a seismometer to enable the seismometer to record seismic wave signals fed back by each three-component detector;
C. b, exciting SH type seismic waves by adopting a transverse wave seismic source at the shock excitation point determined in the step A, and detecting a small fault in front of the head of the coal roadway;
D. and (3) observing and analyzing seismic wave signals fed back by each three-component geophone in real time through a seismometer:
firstly, after observing a feedback seismic direct wave signal, performing inter-channel balance and intra-channel balance pretreatment on a subsequently received seismic wave signal, and if a seismic wave signal with an amplitude not lower than 1/3 of the amplitude of the direct wave and being a same-phase axis is not found, determining that no minor fault exists in a head-on front detection range;
after observing the fed-back seismic direct wave signal, performing inter-channel balance and intra-channel balance pretreatment on a subsequently received seismic wave signal, if the seismic wave signal with the amplitude not lower than 1/3 of the amplitude of the direct wave and the same-phase axis is found, analyzing the seismic wave signal, and if the seismic wave signal is linearly polarized and has a frequency dispersion characteristic, determining that the seismic wave signal is a Love type surface wave, thereby determining that a small fault exists in a detection range in front of the head;
the principle of determining the existence of the small fault according to Love surface waves is as follows: the inventor researches and discovers that the contact relationship between the coal seam loose circle and the coal seam is similar to the contact relationship between the surface covering layer and the underlying stratum, and the surface wave development is caused by the existence of the coal seam loose circle. When the Love trough waves are transmitted to the interface surface of the coal seam loose ring and the coal seam, Love surface waves are formed in the coal seam loose ring and on the surface of the complete coal seam, the energy of the Love surface waves is mainly concentrated in the coal seam loose ring, and the energy of the middle part of the loose ring is strongest.
Therefore, after SH type seismic waves are transmitted, Love type groove waves are formed in the coal seam and are transmitted to the front; if a small fault exists in front of the head-on, the Love-type trough wave can form a reflected Love-type trough wave after contacting the small fault, and the reflected Love-type trough wave returns to enter a surrounding rock loosening ring on the lateral side of the roadway and is converted into a Love-type surface wave; if no small fault exists in front of the head, the reflected Love-type groove wave and the converted Love-type surface wave cannot be formed; thus, Love-type surface waves are characteristic waves of whether small faults exist or not. In addition, after being received by the three-component detector, the Love-type surface wave is shown as a same-phase axis in seismograph recording and has a frequency dispersion characteristic, and the seismic wave signals of the three components are linearly polarized.
E. And D, performing migration imaging on the seismic wave signals determined as Love type surface waves in the step D by utilizing a known pre-stack Fresnel beam migration method, and then obtaining the position of the small fault in front of the head according to an imaging graph.
In the prior art, estimation is mainly carried out on the influence of a surrounding rock loosening ring on seismic waves so as to remove the influence of the surrounding rock loosening ring on detection, or the geophone is kept away from the surrounding rock loosening ring as far as possible to receive the seismic waves. The invention has simple construction and convenient use, can accurately detect whether a small fault exists in front of the head-on and can determine the position of the small fault, thereby providing guidance for the safe tunneling of the subsequent coal roadway.
Drawings
FIG. 1 is a schematic diagram of the position layout of a three-component detector according to the present invention;
fig. 2 is a top sectional view of fig. 1.
In the figure: 1. roadway, 2, three-component detector, 3, surrounding rock loosening ring, 4, shock excitation point, 5, head-on, 6 and small fault.
Detailed Description
The present invention will be further explained below.
As shown in the figure, the right side of fig. 1 is taken as the front side for description, and the specific steps of the invention are as follows:
A. arranging a shock excitation point on the central axis of the left upper or the right upper of the roadway 1 and close to the head-on position 5, and then acquiring the range of the known surrounding rock loosening ring 3;
B. arranging a plurality of three-component detectors 2 (namely, the detectors with X, Y, Z component receiving functions) in rows at the middle positions of surrounding rock loosening rings 3 on the left side or the right side of a roadway 1 in a drilling mode, and then connecting each three-component detector 2 with a seismometer to enable the seismometer to record seismic wave signals fed back by each three-component detector 2;
C. b, exciting SH type seismic waves by adopting a transverse wave seismic source at the shock excitation point 4 determined in the step A, and detecting a small fault 6 in front of the head of the coal roadway; the range of the front of each detection is 80-150 m;
D. and (3) observing and analyzing seismic wave signals fed back by each three-component geophone in real time through a seismometer:
firstly, after observing a feedback seismic direct wave signal, performing inter-channel balance and intra-channel balance pretreatment on a subsequently received seismic wave signal, and if a seismic wave signal with an amplitude not lower than 1/3 of the amplitude of the direct wave and being a same-phase axis is not found, determining that no small fault 6 exists in a detection range in front of a head 5;
after observing the fed-back seismic direct wave signal, performing inter-channel balance and intra-channel balance pretreatment on a subsequently received seismic wave signal, if the seismic wave signal with the amplitude not lower than 1/3 of the amplitude of the direct wave and the same-phase axis is found, analyzing the seismic wave signal, and if the seismic wave signal is linearly polarized and has a frequency dispersion characteristic, determining that the seismic wave signal is a Love type surface wave, thereby determining that a small fault 6 exists in a detection range in front of the head;
the principle of determining the existence of the small fault 6 according to Love type surface waves is as follows: the inventor researches and discovers that the contact relationship between the coal seam loose circle and the coal seam is similar to the contact relationship between the surface covering layer and the underlying stratum, and the surface wave development is caused by the existence of the coal seam loose circle. When the Love trough waves are transmitted to the interface surface of the coal seam loose ring and the coal seam, Love surface waves are formed in the coal seam loose ring and on the surface of the complete coal seam, the energy of the Love surface waves is mainly concentrated in the coal seam loose ring, and the energy of the middle part of the loose ring is strongest.
Therefore, after SH type seismic waves are transmitted, Love type groove waves are formed in the coal seam and are transmitted to the front; if a small fault exists in front of the head-on, the Love-type trough wave can form a reflected Love-type trough wave after contacting the small fault, and the reflected Love-type trough wave returns to enter a surrounding rock loosening ring on the lateral side of the roadway and is converted into a Love-type surface wave; if no small fault exists in front of the head, the reflected Love-type groove wave and the converted Love-type surface wave cannot be formed; thus, Love-type surface waves are characteristic waves of whether small faults exist or not. In addition, after being received by the three-component detector, the Love-type surface wave is shown as a same-phase axis in seismograph recording and has a frequency dispersion characteristic, and the three-component Love-type surface wave is linearly polarized.
E. And D, performing migration imaging on the seismic wave signals determined as Love type surface waves in the step D by utilizing a known pre-stack Fresnel beam migration method, and then obtaining the position of the small fault 6 in front of the head 5 according to an imaging graph.

Claims (1)

1. A three-dimensional small fault advanced detection method utilizing the influence of a surrounding rock loosening ring is characterized by comprising the following specific steps:
A. arranging a shock excitation point on the central axis of the left side or the right side of the roadway and close to the head-on position; then obtaining the range of the known surrounding rock loosening ring;
B. arranging a plurality of three-component detectors in rows at the middle position of a surrounding rock loosening ring on the left side or the right side of the roadway in a drilling mode, and then connecting each three-component detector with a seismometer to enable the seismometer to record seismic wave signals fed back by each three-component detector;
C. b, exciting SH type seismic waves by adopting a transverse wave seismic source at the shock excitation point determined in the step A, and detecting a small fault in front of the head of the coal roadway;
D. and (3) observing and analyzing seismic wave signals fed back by each three-component geophone in real time through a seismometer:
firstly, after observing a feedback seismic direct wave signal, performing inter-channel balance and intra-channel balance pretreatment on a subsequently received seismic wave signal, and if a seismic wave signal with an amplitude not lower than 1/3 of the amplitude of the direct wave and being a same-phase axis is not found, determining that no minor fault exists in a head-on front detection range;
after observing the fed-back seismic direct wave signal, performing inter-channel balance and intra-channel balance pretreatment on a subsequently received seismic wave signal, if the seismic wave signal with the amplitude not lower than 1/3 of the amplitude of the direct wave and the same-phase axis is found, analyzing the seismic wave signal, and if the seismic wave signal is linearly polarized and has a frequency dispersion characteristic, determining that the seismic wave signal is a Love type surface wave, thereby determining that a small fault exists in a detection range in front of the head;
E. and D, performing migration imaging on the seismic wave signals determined as Love type surface waves in the step D by using a pre-stack Fresnel beam migration method, and then obtaining the position of the small fault in front of the head according to an imaging graph.
CN201910773665.9A 2019-08-21 2019-08-21 Three-dimensional small fault advanced detection method utilizing influence of surrounding rock loosening ring Active CN110531415B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910773665.9A CN110531415B (en) 2019-08-21 2019-08-21 Three-dimensional small fault advanced detection method utilizing influence of surrounding rock loosening ring

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910773665.9A CN110531415B (en) 2019-08-21 2019-08-21 Three-dimensional small fault advanced detection method utilizing influence of surrounding rock loosening ring

Publications (2)

Publication Number Publication Date
CN110531415A CN110531415A (en) 2019-12-03
CN110531415B true CN110531415B (en) 2020-10-30

Family

ID=68664057

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910773665.9A Active CN110531415B (en) 2019-08-21 2019-08-21 Three-dimensional small fault advanced detection method utilizing influence of surrounding rock loosening ring

Country Status (1)

Country Link
CN (1) CN110531415B (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113138418B (en) * 2021-04-26 2022-06-07 中国石油大学(华东) Elastic medium undulating surface Fresnel beam offset imaging method
CN113960669A (en) * 2021-10-21 2022-01-21 中国石油化工股份有限公司 Reflection information enhancement method and device based on common imaging point gather combination calculation

Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6147929A (en) * 1998-03-03 2000-11-14 Gas Research Institute Method for predicting continuous and discontinuous waveguide targets using interwell seismic signature characteristics
CN102788991A (en) * 2012-07-16 2012-11-21 中煤科工集团西安研究院 High-density rapid detection method based on z-component transmitting channel waves
CN103217719A (en) * 2013-04-11 2013-07-24 中国矿业大学 Method of advanced detection of breaking-loss wing coal seam of coal road based on single offset pair observation system
CN103235333A (en) * 2013-04-11 2013-08-07 中国矿业大学 Coal roadway fault advanced-detection method based on reflecting channel wave signals
CN104166164A (en) * 2014-08-08 2014-11-26 山东科技大学 Three-component multiwave reflection three-dimensional earthquake advanced detection method for coal roadway tunneling geologic structure
CN105572745A (en) * 2015-12-10 2016-05-11 北京中矿大地地球探测工程技术有限公司 Seismic prospecting method of three-component slot waves under coalmine well
CN105911588A (en) * 2016-05-16 2016-08-31 中国矿业大学 Stope face small collapse column detection method based on diffracted wave
CN106019371A (en) * 2016-05-13 2016-10-12 中国矿业大学 Outburst coal seam roadway minor fault advanced qualitative forecast method
CN106680868A (en) * 2017-01-05 2017-05-17 南京君诺鑫成网络科技有限公司 In-seam seismic multi-attribute parameter tomographic method
CN107831530A (en) * 2017-10-25 2018-03-23 山东科技大学 High seam reflects slot wave forward probe method along bottom plate or top plate coal driving
CN208506263U (en) * 2018-07-25 2019-02-15 北京思凯维科地球物理信息技术有限公司 A kind of wave detector and advanced forecasting system for tunnel geological prediction

Patent Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6147929A (en) * 1998-03-03 2000-11-14 Gas Research Institute Method for predicting continuous and discontinuous waveguide targets using interwell seismic signature characteristics
CN102788991A (en) * 2012-07-16 2012-11-21 中煤科工集团西安研究院 High-density rapid detection method based on z-component transmitting channel waves
CN103217719A (en) * 2013-04-11 2013-07-24 中国矿业大学 Method of advanced detection of breaking-loss wing coal seam of coal road based on single offset pair observation system
CN103235333A (en) * 2013-04-11 2013-08-07 中国矿业大学 Coal roadway fault advanced-detection method based on reflecting channel wave signals
CN104166164A (en) * 2014-08-08 2014-11-26 山东科技大学 Three-component multiwave reflection three-dimensional earthquake advanced detection method for coal roadway tunneling geologic structure
CN105572745A (en) * 2015-12-10 2016-05-11 北京中矿大地地球探测工程技术有限公司 Seismic prospecting method of three-component slot waves under coalmine well
CN106019371A (en) * 2016-05-13 2016-10-12 中国矿业大学 Outburst coal seam roadway minor fault advanced qualitative forecast method
CN105911588A (en) * 2016-05-16 2016-08-31 中国矿业大学 Stope face small collapse column detection method based on diffracted wave
CN106680868A (en) * 2017-01-05 2017-05-17 南京君诺鑫成网络科技有限公司 In-seam seismic multi-attribute parameter tomographic method
CN107831530A (en) * 2017-10-25 2018-03-23 山东科技大学 High seam reflects slot wave forward probe method along bottom plate or top plate coal driving
CN208506263U (en) * 2018-07-25 2019-02-15 北京思凯维科地球物理信息技术有限公司 A kind of wave detector and advanced forecasting system for tunnel geological prediction

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
基于SH槽波采空区探测技术研究;丁卫;《中州煤炭》;20161031(第10期);第128-132页 *

Also Published As

Publication number Publication date
CN110531415A (en) 2019-12-03

Similar Documents

Publication Publication Date Title
Griffiths et al. Applied geophysics for geologists and engineers: the elements of geophysical prospecting
Urosevic et al. Targeting nickel sulfide deposits from 3D seismicreflection data at Kambalda, Australia
CN104215934B (en) A kind of method utilizing uphole geophone to carry out fracturing micro-seismic monitoring
Götz et al. Vertical seismic profiling using a daisy‐chained deployment of fibre‐optic cables in four wells simultaneously–Case study at the Ketzin carbon dioxide storage site
CN105676279A (en) Earthquake reflection data collection method with concentric-circle equivalent shot-geophone distance
CN110531415B (en) Three-dimensional small fault advanced detection method utilizing influence of surrounding rock loosening ring
Sloan et al. Void detection using near-surface seismic methods
Brodic et al. Three-component seismic land streamer study of an esker architecture through S-and surface-wave imaging
Mari et al. 3D seismic imaging of a near-surface heterogeneous aquifer: a case study
Richter et al. Comparison of pneumatic impact and magnetostrictive vibrator sources for near surface seismic imaging in geotechnical environments
Zhou et al. Seeing coal‐seam top ahead of the drill bit through seismic‐while‐drilling
EP2917857B1 (en) Fracture characterization from refraction travel time data
CN110850472B (en) Variable offset distance advanced fault detection method based on shock wave excitation seismic source
Reiser et al. Imaging the high-temperature geothermal field at Krafla using vertical seismic profiling
CA2018421C (en) Continuity logging using differenced signal detection
CN110531419B (en) Collapse column advanced detection method using Love surface waves
Wawerzinek et al. Performance of high-resolution impact and vibration sources for underground seismic exploration of clay formations at the scale of a rock laboratory
Hammer et al. Seismic reflection imaging of thin, kimberlite dykes and sills: exploration and deposit characterization of the Snap Lake dyke, Canada
WO2017103126A1 (en) Drilling head system with integrated acoustic source and arm which is equipped with electrodynamic detectors
Lüth et al. Seismic tomography and monitoring in underground structures: developments in the Freiberg Reiche Zeche underground lab (Freiberg, Germany) and their application in underground construction (SOUND)
Kahn* et al. Eagle ford microseismic acquisition geometry benchmark
Kepic et al. Enhancing the seismoelectric method via a virtual shot gather
Gritto Subsurface void detection using seismic tomographic imaging
Riedel et al. Seismic imaging of the Kylylahti Cu-Au-Zn ore deposit using conventional and DAS VSP measurements supported by 3D full-waveform seismic modeling
Tomio et al. High-frequency seismic measurements using a piezoelectric type accelerometer array

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant