CN114814946A - Tomography method based on transmission channel wave adjacent channel centroid frequency - Google Patents
Tomography method based on transmission channel wave adjacent channel centroid frequency Download PDFInfo
- Publication number
- CN114814946A CN114814946A CN202210432816.6A CN202210432816A CN114814946A CN 114814946 A CN114814946 A CN 114814946A CN 202210432816 A CN202210432816 A CN 202210432816A CN 114814946 A CN114814946 A CN 114814946A
- Authority
- CN
- China
- Prior art keywords
- frequency
- channel
- wave
- centroid frequency
- seismic source
- 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.)
- Granted
Links
- 230000005540 biological transmission Effects 0.000 title claims abstract description 70
- 238000000034 method Methods 0.000 title claims abstract description 34
- 238000003325 tomography Methods 0.000 title claims abstract description 21
- 239000003245 coal Substances 0.000 claims abstract description 30
- 238000001514 detection method Methods 0.000 claims abstract description 13
- 238000003384 imaging method Methods 0.000 claims abstract description 12
- 238000004364 calculation method Methods 0.000 claims abstract description 10
- 238000001228 spectrum Methods 0.000 claims abstract description 9
- 230000008569 process Effects 0.000 claims abstract description 6
- 239000011159 matrix material Substances 0.000 claims description 9
- 238000005070 sampling Methods 0.000 claims description 9
- 238000005516 engineering process Methods 0.000 claims description 3
- 238000000605 extraction Methods 0.000 claims description 3
- 230000008859 change Effects 0.000 abstract description 8
- 238000005065 mining Methods 0.000 description 4
- 230000008878 coupling Effects 0.000 description 2
- 238000010168 coupling process Methods 0.000 description 2
- 238000005859 coupling reaction Methods 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 1
- 230000006978 adaptation Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000014759 maintenance of location Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000009394 selective breeding Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. analysis, for interpretation, for correction
- G01V1/30—Analysis
- G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
Abstract
The invention discloses a tomography method based on transmission channel wave adjacent channel centroid frequency, which comprises the steps of firstly obtaining the centroid frequency of a channel wave signal according to the statistical rule of seismic signal amplitude-frequency spectrum, and then obtaining the centroid frequency change characteristic through an algorithm of the adjacent channel centroid frequency, wherein the algorithm for calculating the centroid frequency is combined with the centroid frequency of the original channel wave signal, the attenuation change of the channel wave centroid frequency can be accurately estimated, and f does not need to be determined in the calculation process S Value, thereby overcoming seismic source variability and f s And (3) artificially selecting conditions such as influence on an imaging result caused by improper selection, then carrying out linear tomography according to the attenuation change rule of the trough wave centroid frequency to obtain the distribution of the frequency shift quantity M of the trough wave centroid frequency in the coal seam working surface, determining whether a geological structure exists in the detection area or not according to the distribution, and determining the position of the geological structure with high precision if the geological structure exists.
Description
Technical Field
The invention relates to a seismic tomography method, in particular to a tomography method based on transmission groove wave adjacent channel centroid frequency.
Background
In recent years, the intelligent unmanned coal mining performance of the coal industry is remarkable, however, the coal mining conditions in China are complex, and the transparent geological conditions for coal mining need to be found and reconstructed urgently. There are many small geological structures (small faults, collapse columns and the like) which affect safe production in the coal mining process, and the channel wave exploration is one of effective means for detecting the small structures (as shown in figure 1). At present, the influence of geological abnormal bodies on the working surface of a coal bed on the energy attenuation of transmitted trough waves is mainly utilized in trough wave exploration, but in the trough wave signal acquisition process, due to factors such as wave detector coupling and the like, serious noise interference generally exists, and the amplitude retention of actually measured signals is low. This phenomenon may cause an excessive energy difference of the channel wave signal, and reduce the applicability of directly using the channel wave energy for imaging, so a data processing method with higher stability and stronger applicability needs to be developed. At present, the method of utilizing the characteristic imaging of the frequency domain of the channel wave is a better channel wave data processing idea.
The channel wave frequency domain characteristic imaging adopts a centroid frequency migration method, and the main principle is that the statistical relationship of high-frequency and low-frequency energy of a seismic signal frequency spectrum is utilized, so that the interference of noise and detector coupling factors is small, and the stability is better. In the existing calculation process of the centroid frequency shift method, f is given S An empirical value, but since the source signal is not acquired in transmitted channel exploration, the frequency expectation f of the channel source S The method cannot be accurately calculated and can only be assigned by experience. At the same time, the method faces f S The estimation is difficult and the difference between different seismic sources cannot be considered, so that the frequency shift characteristic of the coal tank which is obtained finally can only be a relative value. Thus how the source frequency f can be avoided S Because the influence on the imaging result caused by improper manual selection, the tomography of the coal seam working face with high precision can be carried out, and the method is one of the research directions of the industry.
Disclosure of Invention
Aiming at the problems in the prior art, the invention provides a tomography method based on the transmission channel wave adjacent channel centroid frequency without determining f S Value, avoiding source frequency f S Because of the influence of improper artificial selection on the imaging result, the coal seam can be accurately and precisely controlledThe working face is subjected to tomography.
In order to achieve the purpose, the invention adopts the technical scheme that: a tomography method based on transmission channel wave adjacent channel centroid frequency comprises the following specific steps:
A. collecting transmission channel wave signals and forming a plurality of shot-sharing gather: a plurality of seismic source points are arranged on the side wall of the roadway on one side of the coal seam working surface at equal intervals in a row, a plurality of detectors are arranged on the side wall of the roadway on the other side at equal intervals in a row, the plurality of seismic source points and the plurality of detectors are all positioned on the same horizontal plane, and the plurality of detectors are all connected with a seismometer to form a coal seam transmission channel wave exploration and observation system; transmitting transmission channel wave signals to the coal seam from each seismic source point in sequence, wherein the transmission channel wave signals of each seismic source point are received by each detector after being transmitted by the coal seam, and the seismometer forms a shot point-sharing gather by receiving the transmission channel wave signals of the same seismic source point by the plurality of detectors respectively, so that each seismic source point forms a shot point-sharing gather;
B. and (3) performing centroid frequency extraction on each transmission channel slot wave signal received by each detector: according to the statistical relationship of the amplitude-frequency spectrum of the seismic signals, extracting the centroid frequency of each channel of transmission channel wave signals excited by each seismic source acquired by each detector to obtain the centroid frequency of each channel of transmission channel wave signals;
C. calculating the attenuation parameters of the centroid frequency of the adjacent channel according to the centroid frequency of each channel of transmission channel wave signals: the step is to calculate the centroid frequency of the adjacent channel to obtain the attenuation parameter on the basis of solving the centroid frequency of each channel of channel wave signal, the algorithm of the step can accurately estimate the attenuation change of the centroid frequency of the channel wave, and the differences of seismic sources and f are overcome s The method comprises the following steps of artificially selecting conditions such as influence on an imaging result caused by improper selection, and calculating the attenuation parameter of the centroid frequency of an adjacent channel: firstly selecting a seismic source point, then determining the detector which is the shortest from the seismic source point as an initial detector, setting the center of mass frequency of a transmission trough wave excited by the seismic source point received by the detector as an initial value, and acquiring the transmission trough excited by the same seismic source point by the detectors which are the nearest from the detector on two sides of the initial detectorThe wave centroid frequency and the initial value are respectively subjected to adjacent channel calculation, specifically:
a) the attenuation propagation distance of the frequency of the transmitted slot wave is known to be linear, and there are:
in the formula: f. of S 、Respectively exciting the centroid frequency and the variance of the transmitted channel wave signals for the seismic source points; alpha is alpha o Is the slot wave frequency attenuation coefficient; f. of R Receiving the centroid frequency of the single channel transmission channel wave signal of the seismic source point for the detector;
Kα 0 L R =f S -f R (3)
wherein L is R The propagation distance from each detector to the seismic source point;
c) according to equation (3), there is the result of the ith trace in the common shot trace set for each detector:
Kα 0 L i =f s -f i ,i=1,2,.....,N
in a transmission slot wave exploration observation system, adjacent propagation paths are close, so that the centroid frequency attenuation coefficients of slot waves in the two propagation paths are set to be the same, and the following steps are provided:
Kα i (L i+1 -L i )=(f s -f i+1 )-(f s -f i )=f i -f i+1
then, the centroid frequency f of the transmission channel wave signal received by each detector at each seismic source point is obtained according to the formula i And a propagation distance L i And calculating the obtained frequency attenuation coefficient K alpha i Defined as the result on the ith propagation path, i.e.
In the formula: m i The attenuation parameter of the ith channel wave centroid frequency in each common shot point channel set is in a linear relation with the propagation path and the frequency attenuation coefficient, and the parameter is obtained by calculating adjacent channel wave signals of the common shot point channel sets;
according to the formula, the attenuation parameters of the centroid frequency of the channel wave of each channel of the transmission channel wave signal adjacent to each seismic source point can be obtained;
D. performing tomography according to the acquired data: and D, performing inversion on the attenuation parameters of the centroid frequency of each adjacent channel obtained in the step C by adopting a tomography technology, and estimating the distribution condition of the M value in the working surface of the coal bed, thereby realizing geological imaging of the detection area.
Further, the method for solving the centroid frequency of each transmission channel wave signal in the step B comprises the following steps:
let the time signal of the transmitted channel wave be A (t) i ) N is the number of sampling points, the sampling interval is Δ T, and the sampling time T is Δ T (N-1);
secondly, the transmission channel wave signal is processed by FFT to obtain a frequency-amplitude spectrum R (f) j ),
In the formula (f) R The center of mass frequency of the single channel transmission channel wave signal;
through the above formula, each transmission channel wave signal in each detector can be obtained.
Furthermore, in the step D, a discrete image reconstruction technique is adopted during inversion calculation, and the expression is
In the formula,. DELTA.f i The total frequency shift quantity of the transmission slot wave; m is the frequency shift amount in the corresponding grid; d ij The length of the ray in the jth grid on the ith ray is taken as the length of the ray in the jth grid on the ith ray; n is the total number of rays; k is the number of grids; at this time, a grid matrix a (N × K) can be established in the detection area, and the following matrix equation is formed:
M=AΔf -1 (6)
and (4) solving the coefficient matrix of the formula (6), namely obtaining the distribution information of the coal bed M value in the detection area.
Compared with the prior art, the method firstly obtains the centroid frequency of the tank wave signal according to the statistical rule of the amplitude-frequency spectrum of the seismic signal, and then obtains the centroid frequency change characteristic through the algorithm of the centroid frequency of the adjacent channel, wherein the algorithm of the centroid frequency is combined with the centroid frequency of the original tank wave signal, the attenuation change of the centroid frequency of the tank wave can be accurately estimated, and f of each seismic source point does not need to be determined in the calculation process S Value, thereby overcoming seismic source variability and f s And (3) artificially selecting conditions such as influence on an imaging result caused by improper selection, then carrying out linear tomography according to the attenuation change rule of the trough wave centroid frequency to obtain the distribution of the frequency shift quantity M of the trough wave centroid frequency in the coal seam working surface, determining whether a geological structure exists in the detection area or not according to the distribution, and determining the position of the geological structure with high precision if the geological structure exists.
Drawings
FIG. 1 is a schematic diagram of the detection of a transmitted channel wave exploration observation system according to the present invention;
FIG. 2 is a graph showing the attenuation law of the centroid frequency of the transmitted channel wave in the present invention;
FIG. 3 is a graph of the results of tomography of the present invention.
Detailed Description
The present invention will be further explained below.
Example (b): the method of the invention is adopted to carry out actual detection on a certain coal mine working face, and the specific steps are as follows:
A. collecting transmission channel wave signals and forming a plurality of shot-sharing gather: arranging a plurality of seismic source points at a row of equal intervals of 10m on the side wall of the roadway on one side of the coal seam working surface, arranging a plurality of detectors at a row of equal intervals of 10m on the side wall of the roadway on the other side, wherein the plurality of seismic source points and the plurality of detectors are all positioned on the same horizontal plane, and the plurality of detectors are all connected with a seismometer to form a coal seam transmission channel wave exploration and observation system; transmitting transmission channel wave signals to the coal seam from each seismic source point in sequence, wherein the transmission channel wave signals of each seismic source point are received by each detector after being transmitted by the coal seam, and a plurality of transmission channel wave signals received by each detector form a shot-point-sharing gather according to the fact that a plurality of detectors receive the transmission channel wave signals of the same seismic source point respectively, so that each seismic source point forms a shot-point-sharing gather as shown in figure 1;
B. and (3) performing centroid frequency extraction on each transmission channel slot wave signal received by each detector: according to the statistical relationship of the amplitude-frequency spectrum of the seismic signals, extracting the centroid frequency of each channel of transmission channel wave signals excited by each seismic source acquired by each detector to obtain the centroid frequency of each channel of transmission channel wave signals; the method for solving the centroid frequency of each transmission channel wave signal comprises the following steps:
let the time signal of the transmitted channel wave be A (t) i ) N is the number of sampling points, the sampling interval is Δ T, and the sampling time T is Δ T (N-1);
secondly, the transmission channel wave signal is subjected to FFT to obtain a frequency-amplitude spectrum R (f) j ),
In the formula (f) R The center of mass frequency of the single channel transmission channel wave signal;
through the above formula, each transmission channel wave signal in each detector can be obtained.
C、Calculating the attenuation parameters of the centroid frequency of the adjacent channel according to the centroid frequency of each channel of transmission channel wave signals: the step is to calculate the centroid frequency of the adjacent channel to obtain the attenuation parameter on the basis of solving the centroid frequency of each channel of channel wave signals, the algorithm of the step can accurately estimate the attenuation change of the centroid frequency of the channel waves, and the differences of seismic sources and f are overcome s The method comprises the following steps of artificially selecting conditions such as influence on an imaging result caused by improper selection, and calculating attenuation parameters of the centroid frequency of adjacent channels: firstly, selecting a seismic source point, then determining a detector which is the shortest from the seismic source point as an initial detector, setting the center of mass frequency of a transmission trough wave excited by the seismic source point received by the detector as an initial value, and respectively carrying out adjacent channel calculation on the center of mass frequency of the transmission trough wave excited by the same seismic source point, acquired by detectors at two sides of the initial detector, which are nearest to the detector, and the initial value, wherein the method specifically comprises the following steps:
a) the attenuation propagation distance of the frequency of the transmitted slot wave is known to be linear, and there are:
in the formula: f. of S 、Respectively is the centroid frequency and variance of the transmitted channel wave signal when the seismic source point is excited; alpha is alpha o Is the slot wave frequency attenuation coefficient; f. of R Receiving the centroid frequency of the single channel transmission channel wave signal of the seismic source point for the detector;
Kα 0 L R =f S -f R (3)
wherein L is R The propagation distance from each detector to the seismic source point;
c) according to equation (3), there is the result of the ith trace in the common shot trace set for each detector:
Kα 0 L i =f s -f i ,i=1,2,.....,N
in a transmission slot wave exploration observation system, adjacent propagation paths are close, so that the centroid frequency attenuation coefficients of slot waves in the two propagation paths are set to be the same, and the following steps are provided:
Kα i (L i+1 -L i )=(f s -f i+1 )-(f s -f i )=f i -f i+1
then, the centroid frequency f of the transmission channel wave signal received by each detector at each seismic source point is obtained according to the formula i And a propagation distance L i And calculating the obtained frequency attenuation coefficient K alpha i Defined as the result on the ith propagation path, i.e.
In the formula: m i The attenuation parameter of the ith channel wave centroid frequency in each common shot point channel set is in a linear relation with the propagation path and the frequency attenuation coefficient, and the parameter is obtained by calculating adjacent channel wave signals of the common shot point channel sets;
according to the formula, the attenuation parameters of the centroid frequency of the channel wave of each channel of the transmission channel wave signal excited by each seismic source point can be obtained, as shown in FIG. 2;
D. performing tomography according to the acquired data: c, inverting the attenuation parameters of the centroid frequencies of the adjacent channels acquired in the step C by adopting a tomography technology, and estimating the distribution condition of M values in the working surface of the coal bed, thereby realizing geological imaging of the detection area; wherein the reconstruction technique of discrete image is adopted in the inversion calculation, and the expression is
In the formula,. DELTA.f i The total frequency shift quantity of the transmission slot wave; m is the frequency shift in the corresponding gridAn amount; d ij The length of the ray in the jth grid on the ith ray is taken as the length of the ray in the jth grid on the ith ray; n is the total number of rays; k is the number of grids; at this time, a grid matrix a (N × K) can be established in the detection area, and the following matrix equation is formed:
M=AΔf -1 (6)
solving the coefficient matrix of the formula (6) can obtain the distribution information of the coal seam M values in the detection area, as shown in fig. 3.
The above description is only of the preferred embodiments of the present invention, and it should be noted that: it will be apparent to those skilled in the art that various modifications and adaptations can be made without departing from the principles of the invention and these are intended to be within the scope of the invention.
Claims (3)
1. A tomography method based on transmission channel wave adjacent channel centroid frequency is characterized by comprising the following specific steps:
A. collecting transmission channel wave signals and forming a plurality of shot-sharing gather: a plurality of seismic source points are arranged on the side wall of the roadway on one side of the coal seam working surface at equal intervals in a row, a plurality of detectors are arranged on the side wall of the roadway on the other side at equal intervals in a row, the plurality of seismic source points and the plurality of detectors are all positioned on the same horizontal plane, and the plurality of detectors are all connected with a seismometer to form a coal seam transmission channel wave exploration and observation system; transmitting transmission channel wave signals to the coal seam from each seismic source point in sequence, wherein the transmission channel wave signals of each seismic source point are received by each detector after being transmitted by the coal seam, and the seismometer forms a shot point-sharing gather by receiving the transmission channel wave signals of the same seismic source point by the plurality of detectors respectively, so that each seismic source point forms a shot point-sharing gather;
B. and (3) performing centroid frequency extraction on each transmission channel slot wave signal received by each detector: according to the statistical relationship of the amplitude-frequency spectrum of the seismic signals, extracting the centroid frequency of each channel of transmission channel wave signals excited by each seismic source acquired by each detector to obtain the centroid frequency of each channel of transmission channel wave signals;
C. calculating the attenuation parameters of the centroid frequency of the adjacent channel according to the centroid frequency of each channel of transmission channel wave signals: the attenuation parameter calculation process of the centroid frequency of the adjacent track is as follows: firstly, selecting a seismic source point, then determining a detector which is the shortest from the seismic source point as an initial detector, setting the center of mass frequency of a transmission trough wave excited by the seismic source point received by the detector as an initial value, and respectively carrying out adjacent channel calculation on the center of mass frequency of the transmission trough wave excited by the same seismic source point, acquired by detectors at two sides of the initial detector, which are nearest to the detector, and the initial value, wherein the method specifically comprises the following steps:
a) the attenuation propagation distance of the frequency of the transmitted slot wave is known to be linear, and there are:
in the formula: f. of S 、Respectively is the centroid frequency and variance of the transmitted channel wave signal when the seismic source point is excited; alpha is alpha o Is the slot wave frequency attenuation coefficient; f. of R Receiving the centroid frequency of the single channel transmission channel wave signal of the seismic source point for the detector;
Kα 0 L R =f S -f R (3)
wherein L is R The propagation distance from each detector to the seismic source point;
c) according to equation (3), there is the result of the ith trace in the common shot trace set for each detector:
Kα 0 L i =f s -f i ,i=1,2,.....,N
in a transmission slot wave exploration observation system, adjacent propagation paths are close, so that the centroid frequency attenuation coefficients of slot waves in the two propagation paths are set to be the same, and the following steps are provided:
Kα i (L i+1 -L i )=(f s -f i+1 )-(f s -f i )=f i -f i+1
then, the centroid frequency f of the transmission channel wave signal received by each detector at each seismic source point is obtained according to the formula i And a propagation distance L i And calculating the obtained frequency attenuation coefficient K alpha i Defined as the result on the ith propagation path, i.e.
In the formula: m i The attenuation parameter of the ith channel wave centroid frequency in each common shot point channel set is in a linear relation with the propagation path and the frequency attenuation coefficient, and the parameter is obtained by calculating adjacent channel wave signals of the common shot point channel sets;
according to the formula, the attenuation parameters of the centroid frequency of the channel wave of each channel of the transmission channel wave signal adjacent to each seismic source point can be obtained;
D. performing tomography according to the acquired data: and C, inverting the attenuation parameters of the centroid frequencies of the adjacent channels acquired in the step C by adopting a tomography technology, and estimating the distribution condition of the M values in the working surface of the coal bed, thereby realizing geological imaging of the detection area.
2. The tomography method based on the transmission channel adjacent channel centroid frequency in the step B as claimed in claim 1, wherein the method for calculating the centroid frequency of each transmission channel signal in the step B is as follows:
let the time signal of the transmitted channel wave be A (t) i ) N is the number of sampling points, the sampling interval is Δ T, and the sampling time T is Δ T (N-1);
secondly, the transmission channel wave signal is subjected to FFT to obtain a frequency-amplitude spectrum R (f) j ),
In the formula (f) R The center of mass frequency of the single channel transmission channel wave signal;
through the formula, each transmission channel wave signal in each detector can be obtained.
3. The tomography method based on the centroid frequency of the adjacent channel of the transmission slot wave as claimed in claim 1, wherein the inversion calculation in the step D adopts a discrete image reconstruction technique, and the expression is
In the formula,. DELTA.f i The total frequency shift quantity of the transmission slot wave; m is the frequency shift amount in the corresponding grid; d ij The length of the ray in the jth grid on the ith ray is taken as the length of the ray in the jth grid on the ith ray; n is the total number of rays; k is the number of grids; at this time, a grid matrix a (N × K) can be established in the detection area, and the following matrix equation is formed:
M=AΔf -1 (6)
and (4) solving the coefficient matrix of the formula (6), namely obtaining the distribution information of the coal bed M value in the detection area.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210432816.6A CN114814946B (en) | 2022-04-24 | 2022-04-24 | Tomography method based on centroid frequency of adjacent channel of transmission groove wave |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210432816.6A CN114814946B (en) | 2022-04-24 | 2022-04-24 | Tomography method based on centroid frequency of adjacent channel of transmission groove wave |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114814946A true CN114814946A (en) | 2022-07-29 |
CN114814946B CN114814946B (en) | 2024-04-05 |
Family
ID=82507137
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210432816.6A Active CN114814946B (en) | 2022-04-24 | 2022-04-24 | Tomography method based on centroid frequency of adjacent channel of transmission groove wave |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114814946B (en) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5144591A (en) * | 1991-01-02 | 1992-09-01 | Western Atlas International, Inc. | Method for determining geometry of subsurface features while drilling |
US20110273961A1 (en) * | 2010-05-05 | 2011-11-10 | Wenyi Hu | Q Tomography Method |
CN111077572A (en) * | 2019-12-17 | 2020-04-28 | 安徽理工大学 | Quantitative coal thickness prediction method based on inversion of transmission groove wave frequency dispersion curve |
CN111551989A (en) * | 2020-05-20 | 2020-08-18 | 中国科学院地理科学与资源研究所 | Transmission channel wave imaging method, equipment and computer readable storage medium |
CN113107599A (en) * | 2021-04-14 | 2021-07-13 | 山东科技大学 | Amplitude ratio imaging method for adjacent channels of transmission channel waves of hidden structure in coal face |
-
2022
- 2022-04-24 CN CN202210432816.6A patent/CN114814946B/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5144591A (en) * | 1991-01-02 | 1992-09-01 | Western Atlas International, Inc. | Method for determining geometry of subsurface features while drilling |
US20110273961A1 (en) * | 2010-05-05 | 2011-11-10 | Wenyi Hu | Q Tomography Method |
CN111077572A (en) * | 2019-12-17 | 2020-04-28 | 安徽理工大学 | Quantitative coal thickness prediction method based on inversion of transmission groove wave frequency dispersion curve |
CN111551989A (en) * | 2020-05-20 | 2020-08-18 | 中国科学院地理科学与资源研究所 | Transmission channel wave imaging method, equipment and computer readable storage medium |
CN113107599A (en) * | 2021-04-14 | 2021-07-13 | 山东科技大学 | Amplitude ratio imaging method for adjacent channels of transmission channel waves of hidden structure in coal face |
Non-Patent Citations (2)
Title |
---|
何文欣;: "槽波地震勘探在煤层构造探测中的应用", 煤炭技术, no. 02, 10 February 2017 (2017-02-10) * |
田瀚等: "透射槽波相邻道质心频率的层析成像方法", 《煤田地质与勘探》, 13 July 2022 (2022-07-13) * |
Also Published As
Publication number | Publication date |
---|---|
CN114814946B (en) | 2024-04-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108064348B (en) | Seismic travel time tomography inversion method based on two-point ray tracing | |
CN111158049B (en) | Seismic reverse time migration imaging method based on scattering integration method | |
CN109669212B (en) | Seismic data processing method, stratum quality factor estimation method and device | |
CN102590862B (en) | Prestack time migration method for compensating absorptive attenuation | |
CN107219554B (en) | The automatic obtaining method of the Value of residual static correction of land seismic data | |
CN103926622A (en) | Method for suppressing multiple waves based on L1 norm multichannel matched filtering | |
CN105607124A (en) | Seismic-wave near-surface stratum quality factor compensation method and device | |
CN103424777A (en) | Method for increasing seismic imaging resolution ratio | |
CN104330826A (en) | A method for removing various noises under the condition of complex surface | |
CN104614769A (en) | Beam-forming filtering method for suppressing seismic surface waves | |
CN113156427A (en) | Inversion method and device of ground penetrating radar data | |
CN104570116A (en) | Geological marker bed-based time difference analyzing and correcting method | |
CN112487613A (en) | Method and device for determining travel time of stratum wave | |
CN113640881B (en) | Multi-offset-distance two-dimensional transverse high-resolution transient surface wave detection method | |
CN114236612A (en) | Signal processing method for earthquake background noise of water area in land | |
CN105298482B (en) | The separation method of construction back wave and borehole wall direct wave by a kind of well | |
CN102323618A (en) | Coherent noise suppression method based on fractional order Fourier transformation | |
CN114814946A (en) | Tomography method based on transmission channel wave adjacent channel centroid frequency | |
CN114415234B (en) | Method for determining shallow surface transverse wave speed based on active source surface wave dispersion and H/V | |
CN112882101B (en) | Random noise attenuation method and device for pre-stack seismic data | |
CN110888158B (en) | Full waveform inversion method based on RTM constraint | |
CN104142513B (en) | Method for improving chromatography inversion accuracy in seismic exploration | |
CN112526611A (en) | Method and device for extracting surface seismic wave quality factor | |
Han et al. | Automatic time picking and velocity determination on full waveform sonic well logs | |
CN113138419B (en) | Method and device for extracting downlink wavelet and attenuation parameters |
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 |