CN112904424B - Method and device for determining position of abnormal body and readable storage medium - Google Patents
Method and device for determining position of abnormal body and readable storage medium Download PDFInfo
- Publication number
- CN112904424B CN112904424B CN202110081515.9A CN202110081515A CN112904424B CN 112904424 B CN112904424 B CN 112904424B CN 202110081515 A CN202110081515 A CN 202110081515A CN 112904424 B CN112904424 B CN 112904424B
- Authority
- CN
- China
- Prior art keywords
- common
- point
- center point
- coordinates
- seismic data
- 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
Links
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/303—Analysis for determining velocity profiles or travel times
- G01V1/305—Travel times
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/62—Physical property of subsurface
- G01V2210/622—Velocity, density or impedance
- G01V2210/6222—Velocity; travel time
Abstract
The invention discloses a method, a device and a readable storage medium for determining the position of an abnormal body, wherein the method comprises the following steps: acquiring channel wave observation data, wherein the channel wave observation data comprise a plurality of channels of seismic data, and the coordinates of a receiving point and the coordinates of an excitation point corresponding to each channel of seismic data; grouping the multiple channels of seismic data according to the coordinates of the receiving points and the coordinates of the excitation points corresponding to each channel of seismic data to obtain multiple common-center point gather, wherein the common-center point gather comprises the coordinates of a common center point; performing frequency dispersion analysis on each seismic data in each common central point channel set to obtain Ehrlich phase travel time corresponding to each seismic data; for each common-center-point gather, calculating the distance from the common center point of the common-center-point gather to an abnormal body according to the coordinate of the receiving point, the coordinate of the excitation point and the Ehry phase travel time corresponding to each seismic data; and determining the position of the abnormal body according to the coordinate of the common center point corresponding to each common center point gather and the distance from the common center point to the abnormal body.
Description
Technical Field
The present disclosure relates to the field of channel wave detection technologies, and in particular, to a method and an apparatus for determining a position of an anomaly, and a readable storage medium.
Background
The channel wave detection method is one of effective methods for detecting abnormal bodies, collapse columns, gangue clamping and the like in the underground coal mine. The channel wave detection method comprises a transmission method, a reflection method and a transmission-reflection combined detection method. In the reflection method detection, generally, a method of envelope superposition is used to perform multiple times of superposition processing on reflection groove wave data, then a two-pass time of a reflection homophase axis is obtained from the reflection groove wave data after the superposition processing, and then a depth position of an abnormal body is obtained through time-depth conversion. However, when the reflected channel data is overlapped multiple times by using the envelope overlapping method, almost all phase information is lost, so that the resolution of the overlapped reflected channel data is low, and the position of an abnormal body in the obtained coal seam may be inaccurate.
Disclosure of Invention
The embodiment of the invention provides a method and a device for determining the position of an abnormal body and a readable storage medium, which are used for solving the problems that in the detection of a trough wave reflection method, when reflected trough wave data are subjected to multiple overlapping processing by adopting an envelope overlapping method, the resolution of the overlapped reflected trough wave data is low, and the position of the abnormal body in an obtained coal seam is possibly inaccurate.
In order to solve the above problem, in a first aspect, an embodiment of the present invention provides a method for determining a location of an abnormal body, including: acquiring channel wave observation data, wherein the channel wave observation data comprise a plurality of channels of seismic data, and the coordinates of a receiving point and the coordinates of an excitation point corresponding to each channel of seismic data; grouping the multiple channels of seismic data according to the coordinates of the receiving points and the coordinates of the excitation points corresponding to each channel of seismic data to obtain multiple common-center point gather, wherein the common-center point gather comprises the coordinates of a common center point; performing frequency dispersion analysis on each seismic data in each common central point channel set to obtain Ehrlich phase travel time corresponding to each seismic data; for each common-center-point gather, calculating the distance from the common center point of the common-center-point gather to an abnormal body according to the coordinate of the receiving point, the coordinate of the excitation point and the Ehry phase travel time corresponding to each seismic data; and determining the position of the abnormal body according to the coordinate of the common center point corresponding to each common center point gather and the distance from the common center point to the abnormal body.
Optionally, grouping the multiple channels of seismic data according to the coordinates of the receiving point and the coordinates of the excitation point corresponding to each channel of seismic data to obtain multiple common-center gathers, where the grouping includes: calculating the center point coordinates of the connecting line of the receiving point and the excitation point corresponding to each seismic data according to the coordinates of the receiving point and the excitation point corresponding to each seismic data to obtain a plurality of common center point coordinates; and grouping the multiple channels of seismic data according to the common midpoint coordinates to obtain multiple common midpoint gather.
Optionally, performing frequency dispersion analysis on each seismic data in each common midpoint channel set to obtain airy phase travel corresponding to each seismic data, including: performing frequency dispersion analysis on each seismic data in each common central point channel set to form a corresponding frequency dispersion curve; and extracting the airy phase travel time in each frequency dispersion curve to obtain the airy phase travel time corresponding to each seismic data.
Optionally, for each common midpoint gather, a formula for calculating a distance from the common midpoint of the common midpoint gather to the anomaly includes:
where j is the common midpoint gather number, I is the different receiving trace in the same common midpoint gather (I ═ 1,2,3 …, I), and xiIs the distance between the receiving point and the excitation point corresponding to the ith receiving track, TiFor the airy phase travel corresponding to the ith receive channel, hjFor common midpoint gathersDistance from common center point to anomaly.
Optionally, determining the position of the abnormal body according to the coordinate of the common center point corresponding to each common center point gather and the distance from the common center point to the abnormal body, includes: determining the coordinate of the common reflection point corresponding to each common center point gather according to the coordinate of the common center point corresponding to each common center point gather and the distance from the common center point to the abnormal body; and connecting adjacent common reflection points according to the coordinates of the common reflection points to obtain the position of the abnormal body.
Optionally, determining the coordinate of the common reflection point corresponding to each common midpoint gather according to the coordinate of the common midpoint corresponding to each common midpoint gather and the distance from the common midpoint to the abnormal body, including: marking each common central point in a preset coordinate system according to the coordinate of the common central point corresponding to each common central point gather; making a vertical line segment perpendicular to the connection line of the common central points from each common central point, wherein the length of the vertical line segment is the distance from the corresponding common central point to the abnormal body; determining the coordinates of the end points of each vertical line segment; and determining the coordinates of the common reflection point corresponding to each common center point gather according to the coordinates of the end points of each vertical line segment.
Optionally, the method for determining the position of the abnormal body further comprises: when the included angle between the connecting line of two adjacent common reflection points and the connecting line of two corresponding adjacent common central points is larger than a preset angle, drawing a circle by taking the two common central points as circular points and the distance from the corresponding common central point to the abnormal body as a radius; and determining the positions of the abnormal bodies corresponding to the two adjacent concentric points according to the tangent line segments of the two circles.
In a second aspect, an embodiment of the present invention provides an apparatus for determining a location of an abnormal body, including: the acquisition unit is used for acquiring channel wave observation data, wherein the channel wave observation data comprise a plurality of channels of seismic data, and the coordinates of a receiving point and the coordinates of an excitation point which correspond to each channel of seismic data; the grouping unit is used for grouping the multiple channels of seismic data according to the coordinates of the receiving points and the coordinates of the excitation points corresponding to each channel of seismic data to obtain a plurality of common-center point gathers, and each common-center point gather comprises the coordinates of a common center point; the analysis unit is used for performing frequency dispersion analysis on each seismic data in each common central point channel set to obtain the Ehrlich phase travel time corresponding to each seismic data; the computing unit is used for computing the distance from the common center point of each common center point gather to the abnormal body according to the coordinate of the receiving point, the coordinate of the excitation point and the Ehry phase travel time corresponding to each seismic data; and the determining unit is used for determining the position of the abnormal body according to the coordinate of the common center point corresponding to each common center point gather and the distance from the common center point to the abnormal body.
In a third aspect, an embodiment of the present invention provides a computer, including: at least one processor; and a memory communicatively coupled to the at least one processor; wherein the memory stores instructions executable by the at least one processor to cause the at least one processor to perform the method of determining the location of an anomaly as in the first aspect or any of the embodiments of the first aspect.
In a fourth aspect, an embodiment of the present invention provides a computer-readable storage medium, where computer instructions are stored, and the computer instructions are configured to cause a computer to execute the method for determining an abnormal body position according to the first aspect or any implementation manner of the first aspect.
According to the method, the device and the readable storage medium for determining the position of the abnormal body, provided by the embodiment of the invention, multiple channels of seismic data are grouped according to the coordinates of the receiving point and the excitation point corresponding to each channel of seismic data to obtain multiple common-center-point gathers; performing frequency dispersion analysis on each seismic data in each common central point channel set to obtain Ehrlich phase travel time corresponding to each seismic data; for each common-center-point gather, calculating the distance from the common center point of the common-center-point gather to an abnormal body according to the coordinate of the receiving point, the coordinate of the excitation point and the Ehry phase travel time corresponding to each seismic data; the position of the abnormal body is determined according to the coordinate of the common center point corresponding to each common center point gather and the distance from the common center point to the abnormal body, the position of the abnormal body is determined by utilizing the characteristic that the Ehry phase velocity in the channel wave frequency dispersion curve does not change along with the change of the thickness and the propagation distance of the coal bed, and the coordinate of the receiving point, the coordinate of the excitation point and the Ehry phase travel corresponding to each channel of seismic data can be accurately obtained, so that the distance from the common center point corresponding to the common center point gather to the abnormal body can be accurately calculated, and the accurate position of the abnormal body can be obtained. In addition, the method only needs to calculate the distance from the common center point corresponding to the common center point gather to the abnormal body, so the calculation amount is small. And for a plurality of abnormal bodies existing in the coal seam, the method can simultaneously determine the positions of the plurality of abnormal bodies and is suitable for the condition of coal seam thickness variation.
The foregoing description is only an overview of the technical solutions of the present application, and the present application can be implemented according to the content of the description in order to make the technical means of the present application more clearly understood, and the following detailed description of the present application is given in order to make the above and other objects, features, and advantages of the present application more clearly understandable.
Drawings
FIG. 1 is a plot of channel-wave-lux-group velocity in an embodiment of the present invention;
FIG. 2 is a velocity curve of a channel Rayleigh group in an embodiment of the present invention;
FIG. 3 is a schematic diagram of positions of a common center point reflection path, an excitation point and a receiving point according to an embodiment of the present invention;
FIG. 4 is a schematic diagram of the position of an anomaly in an embodiment of the present invention;
FIG. 5 is a flowchart illustrating a method for determining an abnormal body position according to an embodiment of the present invention;
FIG. 6 is a schematic structural diagram of a detection region according to an embodiment of the present invention;
FIG. 7 is a schematic diagram of a dispersion curve according to an embodiment of the present invention;
FIG. 8 is a schematic diagram illustrating the position correction of a non-horizontally oriented anomaly in an embodiment of the present invention;
FIG. 9 is a schematic structural diagram of an apparatus for determining an abnormal body position according to an embodiment of the present invention
Fig. 10 is a schematic diagram of a hardware structure of a computer according to an embodiment of the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are some, but not all, embodiments of the present invention. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
The principle of the invention is as follows: the slot-wave love wave frequency dispersion curve equation is as follows:
in this equation, the inventors have found that when the velocity v of the coal seam isscDensity of coal bed rhoc(c represents the coal seam), the thickness d of the coal seam, and the velocity v of the surrounding rockrDensity of surrounding rock ρr(r represents a surrounding rock) constant, frequency f ═ ω/2 π and group velocity vLOne to one, and there is a group velocity extreme point, i.e., a velocity minimum for the airy phase, as shown in fig. 1. In addition, the group velocity extreme points of the love waves of different coal seam thicknesses are the same, that is, the airy phase velocities of the love waves of different coal seam thicknesses are the same, as shown in fig. 1. This rule also exists for the slot-wave rayleigh wave (as shown in figure 2). That is, under the condition that the physical parameters such as the velocity and density of the surrounding rock and the coal seam are stable, the velocities of the dispersion curve channel wave group velocity extreme points of different coal seam thicknesses are the same (the velocities of the reflection channel wave group velocity extreme points are the same although the reflection paths are different and the reflection time is different).
Therefore, the inventor thought that the airy phase travel time T can be extracted from the slot wave dispersion curveRThen its velocity vR＝lR/TR，lRIndicating the path length (R for reflection). And an equation relationship can be established in each common midpoint gather (CMP) (the relationship between the reflection path of the common midpoint and the positions of the excitation point and the receiving point is shown in fig. 3, and the anomaly is taken as an example in fig. 3, but not limited thereto), and each common midpoint gather establishes an equation:
i is the different receiving channel (I ═ 1,2,3 …, I) in the same common center channel set, TiFor the airy phase travel corresponding to the ith receive channel,/iThe path length from the excitation point to the receiving point corresponding to the ith receiving channel in the gather is shown.
And when the position of the abnormal body is parallel to or not vertical to the trend of the roadway, the distance h between the roadway and the abnormal body and the path l are in the following relation:
x is the offset distance in the common center trace, which can be calculated from the coordinates of the receiving point and the coordinates of the excitation point, and the equation becomes:
wherein j is the common midpoint gather number, xiIs the distance between the receiving point corresponding to the ith receiving track and the excitation point, hjThe distance from the common midpoint of the common midpoint gather to the anomaly.
Solving the equation set of each common-center-point gather to obtain the distance from the common center point of each common-center-point gather to the abnormal body, calculating the coordinate of the common center point of each common-center-point gather, and calculating the h-th distance from the common center point of each common-center-point gather to the abnormal bodyjAnd projecting the CMP co-center point coordinates at corresponding positions, and connecting the end points of the projected line segments to obtain the fault positions, as shown in FIG. 4.
Based on the above principle, an embodiment of the present invention provides a method for determining a location of an abnormal body, as shown in fig. 5, including:
s101, acquiring channel wave observation data, wherein the channel wave observation data comprise a plurality of channels of seismic data, and coordinates of a receiving point and coordinates of an excitation point corresponding to each channel of seismic data; specifically, a tunnel may be established in the detection area, and a receiving point and an excitation point may be arranged on one side of the tunnel to form a reflectometry arrangement, as shown in fig. 6. The excitation points are used for blasting, and the receiving points are used for receiving seismic data. After the seismic data are acquired, the coordinates of the excitation point and the coordinates of the receiving point can be marked on the seismic data, so that each channel of seismic data corresponds to the coordinates of the receiving point and the coordinates of the excitation point. The distance from the excitation point to the receiving point is the offset distance.
S102, grouping the multiple channels of seismic data according to the coordinates of the receiving points and the coordinates of the excitation points corresponding to each channel of seismic data to obtain multiple common-center-point gathers, wherein the common-center-point gathers comprise the coordinates of a common center point; specifically, a plurality of common center points can be calculated according to the coordinates of the receiving points and the coordinates of the excitation points corresponding to the seismic data, and the distance from the common center point to the excitation points is equal to the distance from the common center point to the receiving points. The seismic data with the same common midpoint coordinate position are divided into a group, and a plurality of common midpoint gather sets can be obtained.
S103, performing frequency dispersion analysis on each seismic data in each common central point channel set to obtain Ehrlich phase travel time corresponding to each seismic data; specifically, a dispersion curve can be obtained by performing dispersion analysis on each seismic data in each common midpoint channel set, and then the airy phase travel time corresponding to each seismic data can be obtained.
S104, calculating the distance from the common center point of each common center point gather to an abnormal body according to the coordinate of the receiving point, the coordinate of the excitation point and the Ehry phase travel time corresponding to each seismic data; specifically, for each common midpoint gather, a formula for calculating a distance from a common midpoint of the common midpoint gather to the anomaly includes:
where j is the common midpoint gather number, I is the different receiving trace in the same common midpoint gather (I ═ 1,2,3 …, I), and xiThe distance between the receiving point corresponding to the ith receiving channel and the excitation point can be obtained by calculating the coordinate of the receiving point corresponding to the ith receiving channel and the coordinate of the excitation point, TiFor the ith receiving channelCorresponding airy phase travel time, hjThe distance from the common midpoint of the common midpoint gather to the anomaly. In each common central point channel set, the equation set is an overdetermined equation set, and the solution h of the equation set can be solved by solving the least square solution of the overdetermined equation setj。
And S105, determining the position of the abnormal body according to the coordinate of the common center point corresponding to each common center point gather and the distance from the common center point to the abnormal body. Specifically, each h may bejAnd projecting the CMP common center point coordinates at corresponding positions, and connecting end points of the projected line segments to obtain the fault position.
According to the method for determining the position of the abnormal body, provided by the embodiment of the invention, multiple channels of seismic data are grouped according to the coordinates of the receiving point and the coordinates of the excitation point corresponding to each channel of seismic data, so that multiple common-center point gathers are obtained; performing frequency dispersion analysis on each seismic data in each common central point channel set to obtain Ehrlich phase travel time corresponding to each seismic data; for each common-center-point gather, calculating the distance from the common center point of the common-center-point gather to an abnormal body according to the coordinate of the receiving point, the coordinate of the excitation point and the Ehry phase travel time corresponding to each seismic data; the position of the abnormal body is determined according to the coordinate of the common center point corresponding to each common center point gather and the distance from the common center point to the abnormal body, the position of the abnormal body is determined by utilizing the characteristic that the Ehry phase velocity in the channel wave frequency dispersion curve does not change along with the change of the thickness and the propagation distance of the coal bed, and the coordinate of the receiving point, the coordinate of the excitation point and the Ehry phase travel corresponding to each channel of seismic data can be accurately obtained, so that the distance from the common center point corresponding to the common center point gather to the abnormal body can be accurately calculated, and the accurate position of the abnormal body can be obtained. In addition, the method only needs to calculate the distance from the common center point corresponding to the common center point gather to the abnormal body, so the calculation amount is small. And for a plurality of abnormal bodies existing in the coal seam, the method can simultaneously determine the positions of the plurality of abnormal bodies and is suitable for the condition of coal seam thickness variation.
In an alternative embodiment, in step S102, grouping the multiple channels of seismic data according to the coordinates of the receiving point and the coordinates of the excitation point corresponding to each channel of seismic data to obtain multiple common midpoint gather includes: calculating the center point coordinates of the connecting line of the receiving point and the excitation point corresponding to each seismic data according to the coordinates of the receiving point and the excitation point corresponding to each seismic data to obtain a plurality of common center point coordinates; and grouping the multiple channels of seismic data according to the common midpoint coordinates to obtain multiple common midpoint gather.
Specifically, each seismic data channel corresponds to an excitation point and a receiving point, and the coordinates of the central point of the connecting line between the receiving point and the excitation point corresponding to each seismic data channel can be obtained by calculation. The multi-channel seismic data may have a plurality of different midpoint coordinates, and the seismic data with the same midpoint coordinate is grouped into a group, so as to obtain a common midpoint gather, where the same midpoint coordinate is the common midpoint coordinate of the common midpoint gather.
In the embodiment of the invention, by calculating the central point coordinates of the connecting lines of the receiving points and the excitation points corresponding to each channel of seismic data, the channels of seismic data can be quickly grouped according to the central point coordinates to obtain a plurality of common central point gathers.
In an alternative embodiment, in step S103, performing frequency dispersion analysis on each seismic data in each common midpoint channel set to obtain an airy phase travel time corresponding to each seismic data, including: performing frequency dispersion analysis on each seismic data in each common central point channel set to form a corresponding frequency dispersion curve; and extracting the airy phase travel time in each frequency dispersion curve to obtain the airy phase travel time corresponding to each seismic data.
Specifically, a corresponding dispersion curve may be obtained by performing narrow-band fourier transform or wavelet transform on each seismic data in each common midpoint gather, as shown in fig. 7, where the airy phase is a minimum extreme point of the curve. And then extracting the minimum travel time in the dispersion curve to obtain the corresponding Ehry phase travel time.
In an alternative embodiment, in step S105, determining the position of the abnormal body according to the coordinate of the common center point corresponding to each common center point gather and the distance from the common center point to the abnormal body includes: determining the coordinate of the common reflection point corresponding to each common center point gather according to the coordinate of the common center point corresponding to each common center point gather and the distance from the common center point to the abnormal body; and connecting adjacent common reflection points according to the coordinates of the common reflection points to obtain the position of the abnormal body.
Specifically, as shown in fig. 4, determining the coordinate of the common reflection point corresponding to each common midpoint gather according to the coordinate of the common midpoint corresponding to each common midpoint gather and the distance from the common midpoint to the abnormal body includes: marking the common midpoint (such as M1, M2 and M3) in a preset coordinate system according to the coordinate of the common midpoint corresponding to each common midpoint gather; making a vertical line segment (M1N1, M2N2, M3N3) perpendicular to the connection line of the common center point from each common center point, wherein the length of the vertical line segment is the distance (h1, h2, h3) from the corresponding common center point to the abnormal body; determining coordinates (N1, N2, N3 coordinates) of the end points of each vertical segment; and determining the coordinates of the common reflection point corresponding to each common center point gather according to the coordinates of the end points of each vertical line segment. And then sequentially connecting the adjacent common reflection points to obtain the position of the abnormal body. As shown in fig. 4 by the position of the break layer F.
Because the common reflection points corresponding to each seismic data in the common midpoint gather are the same and the common reflection points are on the abnormal body, the coordinates of the common reflection points corresponding to each common midpoint gather are determined according to the coordinates of the common midpoint corresponding to each common midpoint gather and the distance from the common midpoint to the abnormal body, and the adjacent common reflection points are connected according to the coordinates of the common reflection points, so that the position of the abnormal body can be quickly determined.
In an optional embodiment, the method for determining the location of the anomaly further comprises: when the included angle between the connecting line of two adjacent common reflection points and the connecting line of two corresponding adjacent common central points is larger than a preset angle, drawing a circle by taking the two common central points as circular points and the distance from the corresponding common central point to the abnormal body as a radius; and determining the positions of the abnormal bodies corresponding to the two adjacent concentric points according to the tangent line segments of the two circles.
Specifically, when the fault strike and the roadway included angle are large, dip correction can be adopted. In the specific correction, as shown in fig. 8, two adjacent common center points may be used as dots, a circle may be drawn by using the distance from the corresponding common center point to the abnormal body as a radius, and a tangent line of the two circles may be drawn to obtain positions of the abnormal body corresponding to the two adjacent common center points, where the tangent point is a corrected common reflection point. The corrected angle is θ.
Sinθ＝|hj-hj+1I/Δ x, where Δ x is the spacing between adjacent common center points.
Since the direction of the reflection path is unknown when determining the position of the anomaly, the position of the common reflection point is determined by making a vertical line segment from the common center point of the center point gather. If the determined fault direction is larger than the roadway included angle based on the above, the determined fault direction is indicated to have a certain deviation, and a circle is drawn by taking the two common central points as the circular points and the distance from the corresponding common central point to the abnormal body as the radius; and determining the positions of the abnormal bodies corresponding to the two adjacent common center points according to the tangent line segments of the two circles, and correcting the positions of the abnormal bodies so as to obtain accurate positions of the abnormal bodies.
An embodiment of the present invention further provides a device for determining a location of an abnormal body, as shown in fig. 9, including:
the acquiring unit 201 is configured to acquire channel wave observation data, where the channel wave observation data includes multiple channels of seismic data, and coordinates of a receiving point and coordinates of an excitation point corresponding to each channel of seismic data; the specific implementation manner is described in detail in step S101 of the above embodiment, and is not described again here.
The grouping unit 202 is configured to group the multiple channels of seismic data according to the coordinates of the receiving point and the coordinates of the excitation point corresponding to each channel of seismic data to obtain multiple common-center point gathers, where each common-center point gather includes the coordinates of a common center point; the specific implementation manner is described in detail in step S101 of the above embodiment, and is not described again here.
The analysis unit 203 is configured to perform frequency dispersion analysis on each seismic data in each common midpoint channel set to obtain an airy phase travel time corresponding to each seismic data; the specific implementation manner is described in detail in step S101 of the above embodiment, and is not described again here.
The calculating unit 204 is configured to calculate, for each common midpoint gather, a distance from the common midpoint of the common midpoint gather to the anomaly according to the coordinates of the receiving point and the excitation point corresponding to each seismic data and the airy phase travel time; the specific implementation manner is described in detail in step S101 of the above embodiment, and is not described again here.
The determining unit 205 is configured to determine a position of the abnormal object according to the coordinate of the common center point corresponding to each common center point gather and the distance from the common center point to the abnormal object. The specific implementation manner is described in detail in step S101 of the above embodiment, and is not described again here.
According to the device for determining the position of the abnormal body, provided by the embodiment of the invention, multiple channels of seismic data are grouped according to the coordinates of the receiving point and the coordinates of the excitation point corresponding to each channel of seismic data, so that multiple common-center point gathers are obtained; performing frequency dispersion analysis on each seismic data in each common central point channel set to obtain Ehrlich phase travel time corresponding to each seismic data; for each common-center-point gather, calculating the distance from the common center point of the common-center-point gather to an abnormal body according to the coordinate of the receiving point, the coordinate of the excitation point and the Ehry phase travel time corresponding to each seismic data; the position of the abnormal body is determined according to the coordinate of the common center point corresponding to each common center point gather and the distance from the common center point to the abnormal body, the position of the abnormal body is determined by utilizing the characteristic that the Ehry phase velocity in the channel wave frequency dispersion curve does not change along with the change of the thickness and the propagation distance of the coal bed, and the coordinate of the receiving point, the coordinate of the excitation point and the Ehry phase travel corresponding to each channel of seismic data can be accurately obtained, so that the distance from the common center point corresponding to the common center point gather to the abnormal body can be accurately calculated, and the accurate position of the abnormal body can be obtained. In addition, the method only needs to calculate the distance from the common center point corresponding to the common center point gather to the abnormal body, so the calculation amount is small. And for a plurality of abnormal bodies existing in the coal seam, the method can simultaneously determine the positions of the plurality of abnormal bodies and is suitable for the condition of coal seam thickness variation.
Based on the same inventive concept as the method for determining the position of an abnormal body in the foregoing embodiment, the present invention further provides a computer, as shown in fig. 10, including: a processor 31 and a memory 32, wherein the processor 31 and the memory 32 may be connected by a bus or other means, and the connection by the bus is illustrated in fig. 10 as an example.
The processor 31 may be a Central Processing Unit (CPU). The Processor 31 may also be other general purpose processors, Digital Signal Processors (DSPs), Application Specific Integrated Circuits (ASICs), Field Programmable Gate Arrays (FPGAs) or other Programmable logic devices, discrete Gate or transistor logic devices, discrete hardware components, or combinations thereof.
The memory 32, which is a non-transitory computer readable storage medium, may be used for storing non-transitory software programs, non-transitory computer executable programs, and modules, such as program instructions/modules corresponding to the method for determining the location of an anomaly in an embodiment of the present invention. The processor 31 executes various functional applications and data processing of the processor by executing non-transitory software programs, instructions and modules stored in the memory 32, namely, implements the method for determining the location of the abnormal body in the above method embodiment.
The memory 32 may include a storage program area and a storage data area, wherein the storage program area may store an operating system, an application program required for at least one function; the storage data area may store data created by the processor 31, and the like. Further, the memory 32 may include high speed random access memory, and may also include non-transitory memory, such as at least one magnetic disk storage device, flash memory device, or other non-transitory solid state storage device. In some embodiments, the memory 32 may optionally include memory located remotely from the processor 31, and these remote memories may be connected to the processor 31 via a network. Examples of such networks include, but are not limited to, the internet, intranets, local area networks, mobile communication networks, and combinations thereof.
One or more of the modules described above are stored in the memory 32 and, when executed by the processor 31, perform the method of determining the location of an anomaly as in the embodiment shown in fig. 5.
The details of the computer can be understood by referring to the corresponding related description and effects in the embodiment shown in fig. 5, which are not repeated herein.
It will be understood by those skilled in the art that all or part of the processes of the methods of the embodiments described above can be implemented by a computer program, which can be stored in a computer-readable storage medium, and when executed, can include the processes of the embodiments of the methods described above. The storage medium may be a magnetic Disk, an optical Disk, a Read-Only Memory (ROM), a Random Access Memory (RAM), a Flash Memory (Flash Memory), a Hard Disk (Hard Disk Drive, abbreviated as HDD), a Solid State Drive (SSD), or the like; the storage medium may also comprise a combination of memories of the kind described above.
As will be appreciated by one skilled in the art, embodiments of the present invention may be provided as a method, system, or computer program product. Accordingly, the present invention may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present invention may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.
The present invention is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the invention. It will be understood that each flow and/or block of the flow diagrams and/or block diagrams, and combinations of flows and/or blocks in the flow diagrams and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable information processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable information processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable information processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable information processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
It will be apparent to those skilled in the art that various changes and modifications may be made in the present invention without departing from the spirit and scope of the invention. Thus, if such modifications and variations of the present invention fall within the scope of the claims of the present invention and their equivalents, the present invention is also intended to include such modifications and variations.
Claims (10)
1. A method for determining a location of an anomaly, comprising:
acquiring channel wave observation data, wherein the channel wave observation data comprise a plurality of channels of seismic data, and the coordinates of a receiving point and the coordinates of an excitation point corresponding to each channel of seismic data;
grouping the multiple channels of seismic data according to the coordinates of the receiving points and the coordinates of the excitation points corresponding to each channel of seismic data to obtain multiple common-center-point gather sets, wherein the common-center-point gather sets comprise the coordinates of a common center point;
performing frequency dispersion analysis on each seismic data in each common central point channel set to obtain Ehrlich phase travel time corresponding to each seismic data;
under the condition that physical parameters of surrounding rocks and a coal bed are stable, the speeds of the dispersion curve channel wave group speed extreme points with different coal bed thicknesses are the same, and for each common-center point gather, the distance from the common center point of the common-center point gather to an abnormal body is calculated according to the coordinate of a receiving point, the coordinate of an excitation point and the Ehry phase travel corresponding to each seismic data;
and determining the position of the abnormal body according to the coordinate of the common center point corresponding to each common center point gather and the distance from the common center point to the abnormal body.
2. The method for determining the location of an anomaly according to claim 1, wherein the grouping the multiple channels of seismic data according to the coordinates of the receiving point and the coordinates of the excitation point corresponding to each channel of seismic data to obtain multiple common midpoint gathers comprises:
calculating the center point coordinates of the connecting line of the receiving point and the excitation point corresponding to each seismic data according to the coordinates of the receiving point and the excitation point corresponding to each seismic data to obtain a plurality of common center point coordinates;
and grouping the multiple channels of seismic data according to the common midpoint coordinates to obtain a plurality of common midpoint gather.
3. The method for determining the location of an anomaly according to claim 1, wherein the performing a dispersion analysis on each seismic data in each common midpoint gather to obtain an airy phase travel time corresponding to each seismic data comprises:
performing frequency dispersion analysis on each seismic data in each common central point channel set to form a corresponding frequency dispersion curve;
and extracting the airy phase travel time in each frequency dispersion curve to obtain the airy phase travel time corresponding to each seismic data.
4. The method of claim 1, wherein the formula for calculating the distance from the common center point of the common center point gather to the anomaly for each common center point gather comprises:
where j is the common midpoint gather number, I is the different receiving trace in the same common midpoint gather (I ═ 1,2,3 …, I), and xiIs the distance between the receiving point and the excitation point corresponding to the ith receiving track, TiFor the airy phase travel corresponding to the ith receive channel, hjThe distance from the common midpoint of the common midpoint gather to the anomaly.
5. The method for determining the location of an anomaly according to claim 1, wherein the determining the location of the anomaly according to the coordinates of the common center point corresponding to each common center point gather and the distance from the common center point to the anomaly comprises:
determining the coordinate of the common reflection point corresponding to each common center point gather according to the coordinate of the common center point corresponding to each common center point gather and the distance from the common center point to the abnormal body;
and connecting adjacent common reflection points according to the coordinates of the common reflection points to obtain the position of the abnormal body.
6. The method of claim 5, wherein determining the coordinates of the common reflection point corresponding to each common midpoint gather according to the coordinates of the common midpoint corresponding to each common midpoint gather and the distance from the common midpoint to the anomaly comprises:
marking each common central point in a preset coordinate system according to the coordinate of the common central point corresponding to each common central point gather;
making a vertical line segment perpendicular to a connecting line of the common center points from each common center point, wherein the length of the vertical line segment is the distance from the corresponding common center point to the abnormal body;
determining coordinates of end points of the vertical line segments;
and determining the coordinates of the common reflection point corresponding to each common center point gather according to the coordinates of the end points of each vertical line segment.
7. The method of determining the location of an anomaly according to claim 5, further comprising:
when an included angle between a connecting line of two adjacent common reflection points and a connecting line of two corresponding adjacent common central points is larger than a preset angle, the two common central points are respectively used as round points, and the distance from the corresponding common central points to the abnormal body is used as a radius to draw a circle;
and determining the positions of the abnormal bodies corresponding to the two adjacent common center points according to the tangent line segments of the two circles.
8. An apparatus for determining a location of an anomaly, comprising:
the acquisition unit is used for acquiring channel wave observation data, wherein the channel wave observation data comprises a plurality of channels of seismic data, and the coordinates of a receiving point and the coordinates of an excitation point which correspond to each channel of seismic data;
the grouping unit is used for grouping the multiple channels of seismic data according to the coordinates of the receiving points and the coordinates of the excitation points corresponding to each channel of seismic data to obtain a plurality of common-center point gathers, and each common-center point gather comprises the coordinates of a common center point;
the analysis unit is used for performing frequency dispersion analysis on each seismic data in each common central point channel set to obtain the Ehrlich phase travel time corresponding to each seismic data;
the computing unit is used for computing the distance from the common center point of each common center point gather to the abnormal body according to the coordinate of the receiving point, the coordinate of the excitation point and the Ehry phase travel corresponding to each seismic data for each common center point gather according to the characteristic that the speeds of the frequency dispersion curve channel wave group speed extreme points with different coal seam thicknesses are the same under the condition that the physical property parameters of the surrounding rock and the coal seam are stable;
and the determining unit is used for determining the position of the abnormal body according to the coordinate of the common center point corresponding to each common center point gather and the distance from the common center point to the abnormal body.
9. A computer, comprising:
at least one processor; and a memory communicatively coupled to the at least one processor; wherein the memory stores instructions executable by the at least one processor to cause the at least one processor to perform the method of determining the location of an anomaly of any one of claims 1-7.
10. A computer-readable storage medium storing computer instructions for causing a computer to perform the method for determining the location of an abnormality according to any one of claims 1 to 7.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110081515.9A CN112904424B (en) | 2021-01-21 | 2021-01-21 | Method and device for determining position of abnormal body and readable storage medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110081515.9A CN112904424B (en) | 2021-01-21 | 2021-01-21 | Method and device for determining position of abnormal body and readable storage medium |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112904424A CN112904424A (en) | 2021-06-04 |
| CN112904424B true CN112904424B (en) | 2021-10-15 |
Family
ID=76117920
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110081515.9A Active CN112904424B (en) | 2021-01-21 | 2021-01-21 | Method and device for determining position of abnormal body and readable storage medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112904424B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116755154B (en) * | 2023-08-14 | 2023-12-12 | 煤炭科学技术研究院有限公司 | Geological structure identification method and device, electronic equipment and storage medium |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP2926170A4 (en) * | 2012-11-28 | 2016-07-13 | Exxonmobil Upstream Res Co | Reflection seismic data q tomography |
| CN106772637B (en) * | 2016-11-29 | 2018-08-31 | 北京中矿大地地球探测工程技术有限公司 | A kind of forward probe method of the underground coal mine based on slot wave |
| CN110850472B (en) * | 2019-10-18 | 2021-07-02 | 中国矿业大学 | Variable offset distance advanced fault detection method based on shock wave excitation seismic source |
| CN111551989B (en) * | 2020-05-20 | 2021-02-09 | 中国科学院地理科学与资源研究所 | Transmission channel wave imaging method, equipment and computer readable storage medium |
| CN111929730B (en) * | 2020-08-20 | 2021-04-16 | 中国矿业大学(北京) | Small-scale geological anomaly detection method and device |
-
2021
- 2021-01-21 CN CN202110081515.9A patent/CN112904424B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN112904424A (en) | 2021-06-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US9978161B2 (en) | Supporting a creation of a representation of road geometry | |
| CN100549731C (en) | Demodulator probe secondary localization method | |
| CN107918156B (en) | Detect the method and device of subsea node acquisition earthquake data polarity | |
| CN109828302B (en) | Seismic source positioning method and device based on multiple vibration sensors | |
| CN112904424B (en) | Method and device for determining position of abnormal body and readable storage medium | |
| WO2022052567A1 (en) | Vehicle positioning method and apparatus, vehicle, and storage medium | |
| CN112346120B (en) | Method for calculating variation degree of channel wave frequency dispersion curve and method for identifying abnormal body | |
| CN106597496A (en) | Position correction method and related device | |
| CN103644922A (en) | Method and system for judging deviation from preset action route | |
| CN116892898B (en) | Track error detection method, device and system for agricultural machinery | |
| CN104614554A (en) | Self-amending method of standard error of ship-based wind speed and direction transducer | |
| CN113420004A (en) | Tunnel point cloud data storage method and device, computer equipment and storage medium | |
| CN116907412B (en) | Agricultural machinery row spacing deviation detection method, device and system | |
| CN104965099A (en) | Velocity measurement method and device | |
| CN103925904A (en) | Ultra short base line installation angle offset unbiased estimation method based on symmetric measuring line | |
| CN109188504B (en) | Underground pipe gallery detection method | |
| CN116182862A (en) | Road boundary determination method, device, electronic equipment and storage medium | |
| CN113702928B (en) | Intersection radar calibration method and equipment | |
| CN104570091B (en) | A kind of method and apparatus for obtaining first arrival wave ray | |
| CN104077940A (en) | Vehicle information measuring system and method in driving school examination | |
| CN105487064A (en) | Radar signal processing method and apparatus | |
| CN107333230B (en) | Vehicle angle determination method and device | |
| CN107064942A (en) | A kind of high-precision Doppler Signal Simulation Method | |
| CN104574463A (en) | Computing method and device of approximate parallel curves of curve | |
| CN111624655B (en) | Method and device for determining residual static correction value of first-motion wave |
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 | ||
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20220914 Address after: Room 401, 4th Floor, No. 7 Xueyuan Road, Haidian District, Beijing 100000 Patentee after: BEIJING SIKAIWEIKE PHYSICAL GEOGRAPHY INFORMATION TECHNOLOGY CO.,LTD. Address before: No.11, Datun Road, Chaoyang District, Beijing Patentee before: Institute of Geographic Sciences and Natural Resources Research, CAS |