CN116755154A - Geological structure identification method and device, electronic equipment and storage medium - Google Patents

Geological structure identification method and device, electronic equipment and storage medium Download PDF

Info

Publication number
CN116755154A
CN116755154A CN202311017100.0A CN202311017100A CN116755154A CN 116755154 A CN116755154 A CN 116755154A CN 202311017100 A CN202311017100 A CN 202311017100A CN 116755154 A CN116755154 A CN 116755154A
Authority
CN
China
Prior art keywords
value
variance
ith row
jth column
row
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
Application number
CN202311017100.0A
Other languages
Chinese (zh)
Other versions
CN116755154B (en
Inventor
廉玉广
刘波
李文
邱浩
李江华
马志超
王国库
张�林
崔金亮
华可
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
CCTEG China Coal Research Institute
Original Assignee
CCTEG China Coal Research Institute
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by CCTEG China Coal Research Institute filed Critical CCTEG China Coal Research Institute
Priority to CN202311017100.0A priority Critical patent/CN116755154B/en
Publication of CN116755154A publication Critical patent/CN116755154A/en
Application granted granted Critical
Publication of CN116755154B publication Critical patent/CN116755154B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • 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/307—Analysis for determining seismic attributes, e.g. amplitude, instantaneous phase or frequency, reflection strength or polarity

Abstract

The invention discloses a geological structure identification method, a geological structure identification device, electronic equipment and a storage medium, and belongs to the technical field of coal mine seismic exploration. The method comprises the following steps: acquiring an Ehrlich phase attribute value set of a transmission groove wave of a target area acquired by a seismic channel; acquiring a variance of the Ehrlich-phase attribute value set; and carrying out geological structure identification on the target area based on the variance body. Therefore, the geological structure identification can be carried out on the target area by considering the variance of the Ehrlich-phase attribute value set of the transmission groove wave, the accuracy of geological structure identification based on the transmission groove wave is improved, and the method is suitable for coal mine seismic exploration scenes.

Description

Geological structure identification method and device, electronic equipment and storage medium
Technical Field
The invention relates to the technical field of coal mine seismic exploration, in particular to a geological structure identification method, a device, electronic equipment and a storage medium.
Background
In order to improve the safety and efficiency of coal mining, it is necessary to perform coal mine seismic exploration, such as detecting geological structures of faults, collapse columns, broken zones of coal and rock masses, stress concentration areas, and the like. At present, geological structure identification is mostly carried out based on transmission groove waves, however, the method for geological structure identification based on transmission groove waves in the related art has the problem of low accuracy.
Disclosure of Invention
The present invention aims to solve at least to some extent one of the technical problems in the above-described technology.
Therefore, an object of the present invention is to provide a geological structure identification method, which obtains an ehrling phase attribute value set of a transmission channel wave of a target area acquired by a seismic trace, obtains a variance of the ehrling phase attribute value set, and performs geological structure identification on the target area based on the variance. Therefore, the geological structure identification can be carried out on the target area by considering the variance of the Ehrlich-phase attribute value set of the transmission groove wave, the accuracy of geological structure identification based on the transmission groove wave is improved, and the method is suitable for coal mine seismic exploration scenes.
A second object of the invention is to propose a device for identifying geological formations.
A third object of the present invention is to propose an electronic device.
A fourth object of the present invention is to propose a computer readable storage medium.
An embodiment of a first aspect of the present invention provides a method for identifying a geological structure, including: acquiring an Ehrlich phase attribute value set of a transmission groove wave of a target area acquired by a seismic channel; acquiring a variance of the Ehrlich-phase attribute value set; and carrying out geological structure identification on the target area based on the variance body.
According to the geological structure identification method provided by the embodiment of the invention, the Ehrlich-phase attribute value set of the transmission channel wave of the target area acquired by the seismic channel is obtained, the variance of the Ehrlich-phase attribute value set is obtained, and the geological structure identification is carried out on the target area based on the variance. Therefore, the geological structure identification can be carried out on the target area by considering the variance of the Ehrlich-phase attribute value set of the transmission groove wave, the accuracy of geological structure identification based on the transmission groove wave is improved, and the method is suitable for coal mine seismic exploration scenes.
In addition, the geological structure identification method according to the embodiment of the present invention may further have the following additional technical features:
in one embodiment of the present invention, the identifying the geological structure of the target area based on the variance includes: obtaining the Ehrlich phase attribute characteristics of the seismic channels based on the variance; and carrying out geological structure identification on the target area based on the Ehrlich phase attribute characteristics.
In one embodiment of the invention, the variance includes (N-d) a first variance value of row L and column, and the Ehrlich phase attribute feature includes an Ehrlich phase amplitude intensity coefficient;
the obtaining the Airy phase attribute characteristics of the seismic channel based on the variance comprises the following steps: acquiring absolute values of a first variance value of an ith row and an ith column and a variance value of a first variance value of an ith row and an (j+1) th column as first absolute values of the ith row and the jth column; obtaining the ehrlichia phase amplitude intensity coefficient based on the first absolute value of the (N-d) row (L-1) column; wherein N, L is a positive integer, d and i are positive integers not greater than N, and j is a positive integer not greater than L.
In one embodiment of the present invention, the obtaining the ehrlichia phase amplitude intensity coefficient based on the first absolute value of the (N-d) row (L-1) column includes: determining the minimum value between the absolute value of the first variance value of the ith row and the jth column and the absolute value of the first variance value of the ith row and the jth (j+1) column as a first minimum value of the ith row and the jth column; acquiring the ratio of the first absolute value of the ith row and the jth column to the first minimum value of the ith row and the jth column as a first attribute of the ith row and the jth column; and obtaining the sum value of the first attribute of the row (N-d) and the column (L-1) as the Ehrlich phase amplitude intensity coefficient.
In one embodiment of the invention, the variance includes a first variance value of (N-d) row L column, and the Ehrlich phase attribute feature includes an Ehrlich phase bias;
the obtaining the Airy phase attribute characteristics of the seismic channel based on the variance comprises the following steps: acquiring the absolute value of the difference value of the first difference value of the ith row and the jth column and the absolute value of the difference value of the first difference value of the (i+1) th row and the jth column as a second absolute value of the ith row and the jth column; obtaining the ehrling phase deviation based on the second absolute value of the (N-d-1) row and the L column; wherein N, L is a positive integer, d and i are positive integers not greater than N, and j is a positive integer not greater than L.
In one embodiment of the present invention, the obtaining the ehrling phase deviation based on the second absolute value of the (N-d-1) row and the column includes: determining the minimum value between the absolute value of the first variance value of the ith row and the jth column and the absolute value of the first variance value of the (i+1) th row and the jth column as a second minimum value of the ith row and the jth column; acquiring a ratio of a second absolute value of the ith row and the jth column to a second minimum value of the ith row and the jth column as a second attribute of the ith row and the jth column; and obtaining the sum value of the second attribute of the L columns of the (N-d-1) row as the Egypt deviation degree.
In one embodiment of the present invention, the identifying the geological structure of the target area based on the ehrling phase attribute features includes: identifying a target interval in which the ehrling phase attribute features are located from a plurality of candidate intervals; and acquiring a geological structure identification result corresponding to the target interval as a geological structure identification result of the target area.
In one embodiment of the present invention, the ehrling phase attribute value set includes N rows and L columns of the ehrling phase attribute values, and the acquiring the ehrling phase attribute value set of the transmission channel wave of the target area acquired by the seismic trace includes: performing time-frequency analysis on the transmission channel waves of the target area acquired by the seismic channels to obtain a time frequency spectrum of the transmission channel waves of the target area; taking the amplitude of the transmission groove wave of the target area at the ith sampling time point and the jth sampling frequency of the time spectrum as an Egyptian attribute value of the ith row and the jth column; wherein N is the number of sampling time points, L is the number of sampling frequencies, i is a positive integer not greater than N, and j is a positive integer not greater than L.
In one embodiment of the present invention, the variance includes (N-d) a first variance of L columns, d is a positive integer not greater than N, and the obtaining the variance of the ehrling phase attribute value set includes: acquiring the sum of the Egyptian attribute values of the ith row to the (i+d-1) th row under the jth column as a target value of the ith row and the jth column; acquiring the ratio of the target value of the ith row and the jth column to d, and taking the ratio as the expectation of the ith row and the jth column; obtaining a first square sum of the difference value between the target value of the ith row and the jth column and each target value of the ith row, and taking the ratio of the first square sum to L as a second square difference value of the ith row and the jth column; obtaining a second square sum of the target value of the ith row and the jth column and each expected difference value of the ith row, and taking the ratio of the second square sum to L as a third difference value of the ith row and the jth column; and obtaining a first variance value of the ith row and the jth column based on the second variance value of the ith row and the jth column and the third variance value of the ith row and the jth column.
In one embodiment of the present invention, the obtaining the first variance value of the ith row and the jth column based on the second variance value of the ith row and the jth column and the third variance value of the jth column includes: obtaining a sum of the second variance value of the ith row and the jth column and the third variance value of the ith row and the jth column as a fourth variance value of the ith row and the jth column; and carrying out normalization processing on the fourth difference value of the ith row and the jth column to obtain a first difference value of the ith row and the jth column.
An embodiment of a second aspect of the present invention provides a geological structure identification device, including: the first acquisition module is used for acquiring an Ehrlich phase attribute value set of a transmission groove wave of a target area acquired by the seismic channel; the second acquisition module is used for acquiring the variance of the Ehrlich phase attribute value set; and the identification module is used for carrying out geological structure identification on the target area based on the variance body.
According to the geological structure identification device, an Ehrlich-phase attribute value set of a transmission groove wave of a target area acquired by a seismic channel is obtained, a variance of the Ehrlich-phase attribute value set is obtained, and geological structure identification is carried out on the target area based on the variance. Therefore, the geological structure identification can be carried out on the target area by considering the variance of the Ehrlich-phase attribute value set of the transmission groove wave, the accuracy of geological structure identification based on the transmission groove wave is improved, and the method is suitable for coal mine seismic exploration scenes.
In addition, the geological structure identification device according to the embodiment of the invention may further have the following additional technical features:
in one embodiment of the present invention, the identification module is further configured to: obtaining the Ehrlich phase attribute characteristics of the seismic channels based on the variance; and carrying out geological structure identification on the target area based on the Ehrlich phase attribute characteristics.
In one embodiment of the invention, the variance includes (N-d) a first variance value of row L and column, and the Ehrlich phase attribute feature includes an Ehrlich phase amplitude intensity coefficient;
the identification module is further configured to: acquiring absolute values of a first variance value of an ith row and an ith column and a variance value of a first variance value of an ith row and an (j+1) th column as first absolute values of the ith row and the jth column; obtaining the ehrlichia phase amplitude intensity coefficient based on the first absolute value of the (N-d) row (L-1) column; wherein N, L is a positive integer, d and i are positive integers not greater than N, and j is a positive integer not greater than L.
In one embodiment of the present invention, the identification module is further configured to: determining the minimum value between the absolute value of the first variance value of the ith row and the jth column and the absolute value of the first variance value of the ith row and the jth (j+1) column as a first minimum value of the ith row and the jth column; acquiring the ratio of the first absolute value of the ith row and the jth column to the first minimum value of the ith row and the jth column as a first attribute of the ith row and the jth column; and obtaining the sum value of the first attribute of the row (N-d) and the column (L-1) as the Ehrlich phase amplitude intensity coefficient.
In one embodiment of the invention, the variance includes a first variance value of (N-d) row L column, and the Ehrlich phase attribute feature includes an Ehrlich phase bias;
The identification module is further configured to: acquiring the absolute value of the difference value of the first difference value of the ith row and the jth column and the absolute value of the difference value of the first difference value of the (i+1) th row and the jth column as a second absolute value of the ith row and the jth column; obtaining the ehrling phase deviation based on the second absolute value of the (N-d-1) row and the L column; wherein N, L is a positive integer, d and i are positive integers not greater than N, and j is a positive integer not greater than L.
In one embodiment of the present invention, the identification module is further configured to: determining the minimum value between the absolute value of the first variance value of the ith row and the jth column and the absolute value of the first variance value of the (i+1) th row and the jth column as a second minimum value of the ith row and the jth column; acquiring a ratio of a second absolute value of the ith row and the jth column to a second minimum value of the ith row and the jth column as a second attribute of the ith row and the jth column; and obtaining the sum value of the second attribute of the L columns of the (N-d-1) row as the Egypt deviation degree.
In one embodiment of the present invention, the identification module is further configured to: identifying a target interval in which the ehrling phase attribute features are located from a plurality of candidate intervals; and acquiring a geological structure identification result corresponding to the target interval as a geological structure identification result of the target area.
In one embodiment of the present invention, the ehrling phase attribute value set includes N rows and L columns of the ehrling phase attribute values, and the first obtaining module is further configured to: performing time-frequency analysis on the transmission channel waves of the target area acquired by the seismic channels to obtain a time frequency spectrum of the transmission channel waves of the target area; taking the amplitude of the transmission groove wave of the target area at the ith sampling time point and the jth sampling frequency of the time spectrum as an Egyptian attribute value of the ith row and the jth column; wherein N is the number of sampling time points, L is the number of sampling frequencies, i is a positive integer not greater than N, and j is a positive integer not greater than L.
In one embodiment of the present invention, the variance body includes (N-d) a first variance value of L columns, d is a positive integer not greater than N, and the second obtaining module is further configured to: acquiring the sum of the Egyptian attribute values of the ith row to the (i+d-1) th row under the jth column as a target value of the ith row and the jth column; acquiring the ratio of the target value of the ith row and the jth column to d, and taking the ratio as the expectation of the ith row and the jth column; obtaining a first square sum of the difference value between the target value of the ith row and the jth column and each target value of the ith row, and taking the ratio of the first square sum to L as a second square difference value of the ith row and the jth column; obtaining a second square sum of the target value of the ith row and the jth column and each expected difference value of the ith row, and taking the ratio of the second square sum to L as a third difference value of the ith row and the jth column; and obtaining a first variance value of the ith row and the jth column based on the second variance value of the ith row and the jth column and the third variance value of the ith row and the jth column.
In one embodiment of the present invention, the second obtaining module is further configured to: obtaining a sum of the second variance value of the ith row and the jth column and the third variance value of the ith row and the jth column as a fourth variance value of the ith row and the jth column; and carrying out normalization processing on the fourth difference value of the ith row and the jth column to obtain a first difference value of the ith row and the jth column.
An embodiment of a third aspect of the present invention provides an electronic device, including: a memory, a processor and a computer program stored on the memory and executable on the processor, the processor implementing the method for identifying geological formations according to the embodiments of the first aspect of the invention when executing the program.
According to the electronic equipment, a computer program stored in a memory is executed through a processor, an Ehrlich-phase attribute value set of a transmission groove wave of a target area acquired by a seismic channel is obtained, a variance body of the Ehrlich-phase attribute value set is obtained, and geological structure identification is carried out on the target area based on the variance body. Therefore, the geological structure identification can be carried out on the target area by considering the variance of the Ehrlich-phase attribute value set of the transmission groove wave, the accuracy of geological structure identification based on the transmission groove wave is improved, and the method is suitable for coal mine seismic exploration scenes.
An embodiment of a fourth aspect of the present invention proposes a computer readable storage medium having stored thereon a computer program which, when executed by a processor, implements a method for identifying a geological structure according to an embodiment of the first aspect of the present invention.
The computer readable storage medium of the embodiment of the invention acquires the Ehrlich-phase attribute value set of the transmission slot wave of the target area acquired by the seismic channel by storing a computer program and executing the computer program by a processor, acquires a variance of the Ehrlich-phase attribute value set, and carries out geological structure identification on the target area based on the variance. Therefore, the geological structure identification can be carried out on the target area by considering the variance of the Ehrlich-phase attribute value set of the transmission groove wave, the accuracy of geological structure identification based on the transmission groove wave is improved, and the method is suitable for coal mine seismic exploration scenes.
Additional aspects and advantages of the invention will be set forth in part in the description which follows and, in part, will be obvious from the description, or may be learned by practice of the invention.
Drawings
The foregoing and/or additional aspects and advantages of the invention will become apparent and readily appreciated from the following description of the embodiments, taken in conjunction with the accompanying drawings, in which:
FIG. 1 is a flow chart of a method of identifying geologic structures according to one embodiment of the invention;
FIG. 2 is a flow chart of a method of identifying geologic structures according to another embodiment of the invention;
FIG. 3 is a schematic diagram of a device for identifying geologic structures according to one embodiment of the invention;
fig. 4 is a schematic structural view of an electronic device according to an embodiment of the present invention.
Detailed Description
Embodiments of the present invention are described in detail below, examples of which are illustrated in the accompanying drawings, wherein like or similar reference numerals refer to like or similar elements or elements having like or similar functions throughout. The embodiments described below by referring to the drawings are illustrative and intended to explain the present invention and should not be construed as limiting the invention.
The following describes a geological structure identification method, a geological structure identification device, an electronic device and a storage medium according to an embodiment of the present invention with reference to the accompanying drawings.
FIG. 1 is a flow chart of a method of identifying geologic structures according to one embodiment of the invention.
As shown in fig. 1, a method for identifying a geological structure according to an embodiment of the present invention includes:
s101, acquiring an Ehrlich phase attribute value set of a transmission groove wave of a target area acquired by a seismic channel.
The Airy phase (Airy phase) refers to a wave corresponding to a period of a maximum or minimum value of the group velocity in the dispersion curve. The Ehrlich phase is the vibration phase with the largest amplitude in the wave train of the surface waves, and is also called Ehrlich vibration phase and Ehrlich equal.
The target area for seismic trace acquisition is not excessively limited, and may include a stope face of a coal mine, a coal seam, and the like.
It should be noted that the ehrling phase attribute value set may include a plurality of ehrling phase attribute values, for example, the ehrling phase attribute value may include an amplitude of a time spectrum of the transmitted slot wave of the target region.
In one embodiment, the ehrling phase attribute value set includes N rows and L columns of ehrling phase attribute values, and the obtaining of the ehrling phase attribute value set of the transmission channel wave of the target area acquired by the seismic trace includes performing time-frequency analysis on the transmission channel wave of the target area acquired by the seismic trace to obtain a time spectrum of the transmission channel wave of the target area, and using the i-th sampling time point of the time spectrum and the j-th amplitude of the transmission channel wave of the target area under the j-th sampling frequency as the ehrling phase attribute value of the i-th row and the j-th column. Wherein N is the number of sampling time points, L is the number of sampling frequencies, i is a positive integer not greater than N, and j is a positive integer not greater than L.
For example, the ehrlich phase attribute value set may be represented by a matrix S of N rows and L columns.
Wherein, ehrlich phase attribute values of the ith row and jth column of the matrix SRefers to the amplitude of the transmitted slot wave of the target area at the ith sampling time point and the jth sampling frequency of the time spectrum, and the Egyptian attribute value of the ith row of the matrix S +.>To the point ofRefers to the amplitude of the transmitted slot wave of the target area at the ith sampling time point of the time spectrum, the Egyptian attribute value of the jth column of the matrix S>To->Refers to the amplitude of the transmitted slot wave of the target region at the jth sampling frequency of the time spectrum.
S102, acquiring a variance of the Ehrlich phase attribute value set.
It should be noted that, the variance body includes a plurality of variance values, and the variance body for obtaining the ehrling phase attribute value set may be implemented by any variance body algorithm in the related art, which is not limited herein.
In one embodiment, obtaining a variance of the set of ehrling phase attribute values includes obtaining variance values of a set number of the ehrling phase attribute values in the set of the ehrling phase attribute values as variance values in the variance. It should be noted that the set number is not excessively limited. For example, the variance values of the plurality of ehrling phase attribute values of the ith row and the variance values of the plurality of ehrling phase attribute values of the jth column may be acquired as the variance values in the variance volume.
In one embodiment, the variance includes a first variance value of (N-d) row and L column, d is a positive integer not greater than N, obtaining a variance of the ehrling phase attribute value set, including obtaining a sum of values of the ehrling phase attribute values of the ith row to (i+d-1) th row under the jth column as a target value of the ith row and the jth column, obtaining a ratio of the target value of the ith row and the jth column to d as a desired value of the ith row and the jth column, obtaining a first sum of squares of differences between the target value of the ith row and the jth column and each target value of the ith row, and obtaining a second sum of squares of the target value of the ith row and the jth column and each desired difference of the ith row as a second sum of squares of the target value of the ith row and the jth column, and obtaining a third difference value of the ith row and the jth column based on the second sum of squares of the second difference value of the ith row and the jth column.
D is the variance volume segment length.
In some examples, deriving the first variance value of the ith row and the jth column based on the second variance value of the ith row and the jth column and the third variance value of the ith row and the jth column includes weighted summing the second variance value of the ith row and the jth column and the third variance value of the ith row and the jth column as the first variance value of the ith row and the jth column.
In some examples, based on the second variance value of the ith row and the jth column and the third variance value of the ith row and the jth column, obtaining the first variance value of the ith row and the jth column includes obtaining a sum of the second variance value of the ith row and the jth column and the third variance value of the ith row and the jth column as a fourth variance value of the ith row and the jth column, and normalizing the fourth variance value of the ith row and the jth column to obtain the first variance value of the ith row and the jth column.
For example, the calculation process of the first variance value of the ith row and jth column is as follows:
D=
wherein, the liquid crystal display device comprises a liquid crystal display device,target value for row i, column j,/>For the expectations of row i and column j, < >>Is the first sum of squares->Second variance value for ith row and jth column,/>Is the second sum of squares->Third variance value for ith row and jth column,/>Fourth aspect of ith row and jth columnDifference of->The first variance value of the ith row and jth column.
S103, carrying out geological structure identification on the target area based on the variance body.
The result of identifying the geologic structure is not limited too much. For example, it may include the presence or absence of a geologic formation, a geologic formation of a target area, and the like. The geological structure of the target area may include faults, subsidence posts, coal and rock mass fracture zones, stress concentration zones, and the like, among others.
In one embodiment, performing geologic structure identification of the target region based on the variance body includes inputting the variance body into an identification model and outputting geologic structure identification of the target region from the identification model. It should be noted that the recognition model is not limited too much, and for example, a deep learning model may be included.
In one embodiment, the method includes performing geologic structure identification on a target region based on a variance, including obtaining a target value based on a plurality of variance values in the variance, and performing geologic structure identification on the target region based on the target value.
In some examples, the target value is obtained based on a plurality of variance values in the variance volume, including at least one operation such as addition, subtraction, multiplication, division, and the like, on the plurality of variance values in the variance volume. For example, a target value is obtained by weighted-averaging a plurality of variance values in the variance volume.
In some examples, performing the geologic structure identification on the target region based on the target value includes identifying a second section in which the target value is located from among the plurality of first sections, and obtaining a geologic structure identification result corresponding to the second section as the geologic structure identification result of the target region.
In summary, according to the geological structure identification method provided by the embodiment of the invention, the Ehrlich-phase attribute value set of the transmission channel wave of the target area acquired by the seismic channel is obtained, the variance of the Ehrlich-phase attribute value set is obtained, and the geological structure identification is performed on the target area based on the variance. Therefore, the geological structure identification can be carried out on the target area by considering the variance of the Ehrlich-phase attribute value set of the transmission groove wave, the accuracy of geological structure identification based on the transmission groove wave is improved, and the method is suitable for coal mine seismic exploration scenes.
Fig. 2 is a flow chart of a method of identifying geologic structures according to another embodiment of the invention.
As shown in fig. 2, the method for identifying a geological structure according to an embodiment of the present invention includes:
s201, acquiring an Ehrlich phase attribute value set of a transmission groove wave of a target area acquired by a seismic channel.
S202, acquiring a variance of the Ehrlich phase attribute value set.
The relevant content of steps S201-S202 can be seen in the above embodiments, and will not be described here again.
S203, obtaining the Ehrlich phase attribute characteristics of the seismic channel based on the variance.
The number of the ehrling phase attribute features is at least one.
In one embodiment, the variance includes (N-d) a first variance value of row L and column, and based on the variance, obtaining the hallway attribute feature of the seismic trace may include the following possible embodiments:
mode 1 obtains an absolute value of a difference between a first variance value of an ith row and an ith column and a first variance value of an ith (j+1) th column as a first absolute value of the ith row and the jth column, and obtains an ehrlich phase amplitude intensity coefficient based on the first absolute value of the (N-d) row and the (L-1) th column.
In an embodiment of the present invention, the ehrling phase attribute feature includes an ehrling phase magnitude intensity coefficient, which is used to characterize a lateral variation of the ehrling phase. Wherein N, L is a positive integer, d and i are positive integers not greater than N, and j is a positive integer not greater than L.
In some examples, the Ehrlich phase amplitude intensity coefficient is positively correlated with the sum of the first absolute values of the (N-d) row (L-1) columns.
In some examples, deriving the ehrling phase magnitude intensity coefficient based on the first absolute value of the (N-d) row (L-1) column includes determining a minimum value between the absolute value of the first variance value of the ith row and jth column, and obtaining a ratio of the first absolute value of the ith row and jth column to the first minimum value of the ith row and jth column as the first attribute of the ith row and jth column, and obtaining the sum value of the first attribute of the (N-d) row and (L-1) column as the ehrling phase magnitude intensity coefficient.
For example, the ehrling phase amplitude intensity coefficient is calculated as follows:
=/>
wherein, the liquid crystal display device comprises a liquid crystal display device,first variance value for ith row and jth column,/>For the first variance value of the (j+1) th column of the ith row,/th column>First absolute value of ith row and jth column,>for the first attribute of the ith row and jth column,/->Is the Ehrlich phase amplitude intensity coefficient, +.>For absolute value arithmetic sign ++>The sign is calculated for the minimum value.
Mode 2, obtaining the absolute value of the difference between the first variance value of the ith row and the jth column and the first variance value of the (i+1) th row and the jth column as the second absolute value of the ith row and the jth column, and obtaining the ehrlich phase offset based on the second absolute value of the (N-d-1) th row and the L column.
In an embodiment of the present invention, the ehrling phase attribute feature includes an ehrling phase deviation, and the ehrling phase amplitude intensity coefficient is used to characterize a longitudinal variation feature of the ehrling phase. Wherein N, L is a positive integer, d and i are positive integers not greater than N, and j is a positive integer not greater than L.
In some examples, the Ehrlich phase bias is positively correlated with the sum of the second absolute values of (N-d-1) rows and columns.
In some examples, deriving the ehrling phase offset based on the second absolute value of the (N-d-1) th row and column includes determining a minimum value between the absolute value of the first variance value of the (i+1) th row and column, as the second minimum value of the (i) th row and column, obtaining a ratio of the second absolute value of the (i) th row and column to the second minimum value of the (i) th row and column, as the second attribute of the (i) th row and column, and obtaining a sum of the second attribute of the (N-d-1) th row and column, as the ehrling phase offset.
For example, the ehrling phase deviation is calculated as follows:
=/>
wherein, the liquid crystal display device comprises a liquid crystal display device,first variance value for ith row and jth column,/>For the first variance value of the j-th column of row (i+1), the value of +.>Second absolute value of ith row and jth column,>for the second property of the ith row and jth column,/->For the Ehrlich phase deviation degree- >For absolute value arithmetic sign ++>The sign is calculated for the minimum value.
S204, carrying out geological structure identification on the target area based on the Ehrlich phase attribute characteristics.
In one embodiment, the geological structure identification is performed on the target area based on the ehrling phase attribute features, including identifying a target area where the ehrling phase attribute features are located from a plurality of candidate areas, and acquiring a geological structure identification result corresponding to the target area as a geological structure identification result of the target area.
In some examples, the ehrling phase attribute feature may be divided into M third sections in advance, the third sections are in one-to-one correspondence with the geologic structure identification result, and a difference portion between the q-th third section and the rest of the third sections other than the q-th third section is used as the q-th candidate section. Wherein the difference part is positioned in the q-th third interval, M is a positive integer greater than 1, and q is a positive integer not greater than M. Therefore, in the method, when the overlapped part exists among the plurality of third sections, only the difference part of the third sections is reserved and used as the candidate section, so that the overlapped part does not exist among the plurality of candidate sections, and the geological structure identification result is unique.
In one embodiment, the number of ehrling phase attribute features is H, and the geological structure identification is performed on the target area based on the ehrling phase attribute features, including performing geological structure identification on the target area based on the t-th ehrling phase attribute features to obtain a t-th geological structure identification result of the target area, and determining the geological structure identification result with the largest occurrence number from the H geological structure identification results of the target area as a final geological structure identification result of the target area. Wherein H is a positive integer, and t is a positive integer not greater than H. Therefore, in the method, a plurality of Ehrlich-phase attribute characteristics can be comprehensively considered, the geological structure identification is carried out on the target area, and the accuracy of the geological structure identification is improved.
In one embodiment, the method further comprises determining a boundary of the geologic structure of the target region based on the ehrlichia-phase attribute feature if the geologic structure identification of the target region indicates that the geologic structure is present in the target region.
In some examples, determining the boundary of the geologic structure of the target region based on the ehrling attribute features includes inputting the ehrling attribute features into a recognition model, and outputting the boundary of the geologic structure of the target region from the recognition model. It should be noted that the recognition model is not limited too much, and for example, a deep learning model may be included.
In summary, according to the geological structure identification method provided by the embodiment of the invention, the Ehrlich-phase attribute characteristics of the seismic channel are obtained based on the variance, and the geological structure identification is performed on the target area based on the Ehrlich-phase attribute characteristics. Therefore, the variance body can be subjected to feature extraction to obtain the Ehrlich-phase attribute features of the seismic channel so as to identify the geological structure of the target area, the accuracy of geological structure identification based on the transmission groove wave is improved, and the method is suitable for the seismic exploration scene of the coal mine.
In order to achieve the above embodiment, the invention further provides a device for identifying a geological structure.
FIG. 3 is a schematic diagram of a device for identifying geologic structures according to one embodiment of the invention.
As shown in fig. 3, a geological structure identification device 100 according to an embodiment of the present invention includes: a first acquisition module 110, a second acquisition module 120, and an identification module 130.
A first acquisition module 110, configured to acquire an ehrling phase attribute value set of a transmission channel wave of a target area acquired by a seismic trace;
a second obtaining module 120, configured to obtain a variance of the ehrlichia phase attribute value set;
and the identification module 130 is used for carrying out geological structure identification on the target area based on the variance.
In one embodiment of the present invention, the identification module 130 is further configured to: obtaining the Ehrlich phase attribute characteristics of the seismic channels based on the variance; and carrying out geological structure identification on the target area based on the Ehrlich phase attribute characteristics.
In one embodiment of the invention, the variance includes (N-d) a first variance value of row L and column, and the Ehrlich phase attribute feature includes an Ehrlich phase amplitude intensity coefficient;
the identification module 130 is further configured to: acquiring absolute values of a first variance value of an ith row and an ith column and a variance value of a first variance value of an ith row and an (j+1) th column as first absolute values of the ith row and the jth column; obtaining the ehrlichia phase amplitude intensity coefficient based on the first absolute value of the (N-d) row (L-1) column; wherein N, L is a positive integer, d and i are positive integers not greater than N, and j is a positive integer not greater than L.
In one embodiment of the present invention, the identification module 130 is further configured to: determining the minimum value between the absolute value of the first variance value of the ith row and the jth column and the absolute value of the first variance value of the ith row and the jth (j+1) column as a first minimum value of the ith row and the jth column; acquiring the ratio of the first absolute value of the ith row and the jth column to the first minimum value of the ith row and the jth column as a first attribute of the ith row and the jth column; and obtaining the sum value of the first attribute of the row (N-d) and the column (L-1) as the Ehrlich phase amplitude intensity coefficient.
In one embodiment of the invention, the variance includes a first variance value of (N-d) row L column, and the Ehrlich phase attribute feature includes an Ehrlich phase bias;
the identification module 130 is further configured to: acquiring the absolute value of the difference value of the first difference value of the ith row and the jth column and the absolute value of the difference value of the first difference value of the (i+1) th row and the jth column as a second absolute value of the ith row and the jth column; obtaining the ehrling phase deviation based on the second absolute value of the (N-d-1) row and the L column; wherein N, L is a positive integer, d and i are positive integers not greater than N, and j is a positive integer not greater than L.
In one embodiment of the present invention, the identification module 130 is further configured to: determining the minimum value between the absolute value of the first variance value of the ith row and the jth column and the absolute value of the first variance value of the (i+1) th row and the jth column as a second minimum value of the ith row and the jth column; acquiring a ratio of a second absolute value of the ith row and the jth column to a second minimum value of the ith row and the jth column as a second attribute of the ith row and the jth column; and obtaining the sum value of the second attribute of the L columns of the (N-d-1) row as the Egypt deviation degree.
In one embodiment of the present invention, the identification module 130 is further configured to: identifying a target interval in which the ehrling phase attribute features are located from a plurality of candidate intervals; and acquiring a geological structure identification result corresponding to the target interval as a geological structure identification result of the target area.
In one embodiment of the present invention, the ehrling phase attribute value set includes N rows and L columns of the ehrling phase attribute values, and the first obtaining module 110 is further configured to: performing time-frequency analysis on the transmission channel waves of the target area acquired by the seismic channels to obtain a time frequency spectrum of the transmission channel waves of the target area; taking the amplitude of the transmission groove wave of the target area at the ith sampling time point and the jth sampling frequency of the time spectrum as an Egyptian attribute value of the ith row and the jth column; wherein N is the number of sampling time points, L is the number of sampling frequencies, i is a positive integer not greater than N, and j is a positive integer not greater than L.
In one embodiment of the present invention, the variance includes (N-d) a first variance value of L columns, d is a positive integer not greater than N, and the second obtaining module 120 is further configured to: acquiring the sum of the Egyptian attribute values of the ith row to the (i+d-1) th row under the jth column as a target value of the ith row and the jth column; acquiring the ratio of the target value of the ith row and the jth column to d, and taking the ratio as the expectation of the ith row and the jth column; obtaining a first square sum of the difference value between the target value of the ith row and the jth column and each target value of the ith row, and taking the ratio of the first square sum to L as a second square difference value of the ith row and the jth column; obtaining a second square sum of the target value of the ith row and the jth column and each expected difference value of the ith row, and taking the ratio of the second square sum to L as a third difference value of the ith row and the jth column; and obtaining a first variance value of the ith row and the jth column based on the second variance value of the ith row and the jth column and the third variance value of the ith row and the jth column.
In one embodiment of the present invention, the second obtaining module 120 is further configured to: obtaining a sum of the second variance value of the ith row and the jth column and the third variance value of the ith row and the jth column as a fourth variance value of the ith row and the jth column; and carrying out normalization processing on the fourth difference value of the ith row and the jth column to obtain a first difference value of the ith row and the jth column.
It should be noted that, details not disclosed in the geological structure identification device in the embodiment of the present invention are referred to details disclosed in the geological structure identification method in the embodiment of the present invention, and are not described herein again.
In summary, the geological structure identification device of the embodiment of the invention acquires the Ehrlich-phase attribute value set of the transmission channel wave of the target area acquired by the seismic channel, acquires the variance of the Ehrlich-phase attribute value set, and carries out geological structure identification on the target area based on the variance. Therefore, the geological structure identification can be carried out on the target area by considering the variance of the Ehrlich-phase attribute value set of the transmission groove wave, the accuracy of geological structure identification based on the transmission groove wave is improved, and the method is suitable for coal mine seismic exploration scenes.
In order to implement the above embodiments, as shown in fig. 4, an embodiment of the present invention proposes an electronic device 200, including: the method for identifying the geological structure comprises a memory 210, a processor 220 and a computer program stored in the memory 210 and capable of running on the processor 220, wherein the processor 220 realizes the identification method of the geological structure when executing the program.
According to the electronic equipment, a computer program stored in a memory is executed through a processor, an Ehrlich-phase attribute value set of a transmission groove wave of a target area acquired by a seismic channel is obtained, a variance body of the Ehrlich-phase attribute value set is obtained, and geological structure identification is carried out on the target area based on the variance body. Therefore, the geological structure identification can be carried out on the target area by considering the variance of the Ehrlich-phase attribute value set of the transmission groove wave, the accuracy of geological structure identification based on the transmission groove wave is improved, and the method is suitable for coal mine seismic exploration scenes.
In order to achieve the above embodiments, an embodiment of the present invention proposes a computer-readable storage medium having stored thereon a computer program which, when executed by a processor, implements the above-described method of identifying a geological structure.
The computer readable storage medium of the embodiment of the invention acquires the Ehrlich-phase attribute value set of the transmission slot wave of the target area acquired by the seismic channel by storing a computer program and executing the computer program by a processor, acquires a variance of the Ehrlich-phase attribute value set, and carries out geological structure identification on the target area based on the variance. Therefore, the geological structure identification can be carried out on the target area by considering the variance of the Ehrlich-phase attribute value set of the transmission groove wave, the accuracy of geological structure identification based on the transmission groove wave is improved, and the method is suitable for coal mine seismic exploration scenes.
Furthermore, the terms "first," "second," and the like, are used for descriptive purposes only and are not to be construed as indicating or implying a relative importance or implicitly indicating the number of technical features indicated. Thus, a feature defining "a first" or "a second" may explicitly or implicitly include one or more such feature. In the description of the present invention, the meaning of "a plurality" is two or more, unless explicitly defined otherwise.
In the description of the present specification, a description referring to terms "one embodiment," "some embodiments," "examples," "specific examples," or "some examples," etc., means that a particular feature, structure, material, or characteristic described in connection with the embodiment or example is included in at least one embodiment or example of the present invention. In this specification, schematic representations of the above terms are not necessarily directed to the same embodiment or example. Furthermore, the particular features, structures, materials, or characteristics described may be combined in any suitable manner in any one or more embodiments or examples. Furthermore, the different embodiments or examples described in this specification and the features of the different embodiments or examples may be combined and combined by those skilled in the art without contradiction.
While embodiments of the present invention have been shown and described above, it will be understood that the above embodiments are illustrative and not to be construed as limiting the invention, and that variations, modifications, alternatives and variations may be made to the above embodiments by one of ordinary skill in the art within the scope of the invention.

Claims (13)

1. A method of identifying a geological structure, comprising:
acquiring an Ehrlich phase attribute value set of a transmission groove wave of a target area acquired by a seismic channel;
acquiring a variance of the Ehrlich-phase attribute value set;
and carrying out geological structure identification on the target area based on the variance body.
2. The method of claim 1, wherein the identifying the geologic structure of the target area based on the variance volume comprises:
obtaining the Ehrlich phase attribute characteristics of the seismic channels based on the variance;
and carrying out geological structure identification on the target area based on the Ehrlich phase attribute characteristics.
3. The method of claim 2, wherein the variance volume comprises (N-d) rows L columns of first variance values, the ehrling phase attribute feature comprising an ehrling phase magnitude intensity coefficient;
the obtaining the Airy phase attribute characteristics of the seismic channel based on the variance comprises the following steps:
Acquiring absolute values of a first variance value of an ith row and an ith column and a variance value of a first variance value of an ith row and an (j+1) th column as first absolute values of the ith row and the jth column;
obtaining the ehrlichia phase amplitude intensity coefficient based on the first absolute value of the (N-d) row (L-1) column; wherein, the liquid crystal display device comprises a liquid crystal display device,
n, L is a positive integer, d and i are positive integers not greater than N, and j is a positive integer not greater than L.
4. A method according to claim 3, wherein said deriving said ehrlichia phase amplitude intensity coefficient based on a first absolute value of (N-d) row (L-1) columns comprises:
determining the minimum value between the absolute value of the first variance value of the ith row and the jth column and the absolute value of the first variance value of the ith row and the jth (j+1) column as a first minimum value of the ith row and the jth column;
acquiring the ratio of the first absolute value of the ith row and the jth column to the first minimum value of the ith row and the jth column as a first attribute of the ith row and the jth column;
and obtaining the sum value of the first attribute of the row (N-d) and the column (L-1) as the Ehrlich phase amplitude intensity coefficient.
5. The method of claim 2, wherein the variance volume comprises first variance values of (N-d) rows L columns, the ehrling phase attribute feature comprising an ehrling phase bias;
The obtaining the Airy phase attribute characteristics of the seismic channel based on the variance comprises the following steps:
acquiring the absolute value of the difference value of the first difference value of the ith row and the jth column and the absolute value of the difference value of the first difference value of the (i+1) th row and the jth column as a second absolute value of the ith row and the jth column;
obtaining the ehrling phase deviation based on the second absolute value of the (N-d-1) row and the L column; wherein, the liquid crystal display device comprises a liquid crystal display device,
n, L is a positive integer, d and i are positive integers not greater than N, and j is a positive integer not greater than L.
6. The method of claim 5, wherein the deriving the ehrling phase bias based on the second absolute value of (N-d-1) row L column comprises:
determining the minimum value between the absolute value of the first variance value of the ith row and the jth column and the absolute value of the first variance value of the (i+1) th row and the jth column as a second minimum value of the ith row and the jth column;
acquiring a ratio of a second absolute value of the ith row and the jth column to a second minimum value of the ith row and the jth column as a second attribute of the ith row and the jth column;
and obtaining the sum value of the second attribute of the L columns of the (N-d-1) row as the Egypt deviation degree.
7. The method of claim 2, wherein the identifying the geologic structure of the target region based on the ehrlichia phase attribute features comprises:
Identifying a target interval in which the ehrling phase attribute features are located from a plurality of candidate intervals;
and acquiring a geological structure identification result corresponding to the target interval as a geological structure identification result of the target area.
8. The method of any of claims 1-7, wherein the set of ehrling phase attribute values includes N rows and L columns of ehrling phase attribute values, the acquiring the set of ehrling phase attribute values of the transmission channel wave of the target region of seismic trace acquisition comprising:
performing time-frequency analysis on the transmission channel waves of the target area acquired by the seismic channels to obtain a time frequency spectrum of the transmission channel waves of the target area;
taking the amplitude of the transmission groove wave of the target area at the ith sampling time point and the jth sampling frequency of the time spectrum as an Egyptian attribute value of the ith row and the jth column; wherein, the liquid crystal display device comprises a liquid crystal display device,
n is the number of sampling time points, L is the number of sampling frequencies, i is a positive integer not greater than N, and j is a positive integer not greater than L.
9. The method of claim 8, wherein the variance includes (N-d) rows and columns of first variance values, d being a positive integer no greater than N, the obtaining the variance of the set of ehrlite property values comprising:
Acquiring the sum of the Egyptian attribute values of the ith row to the (i+d-1) th row under the jth column as a target value of the ith row and the jth column;
acquiring the ratio of the target value of the ith row and the jth column to d, and taking the ratio as the expectation of the ith row and the jth column;
obtaining a first square sum of the difference value between the target value of the ith row and the jth column and each target value of the ith row, and taking the ratio of the first square sum to L as a second square difference value of the ith row and the jth column;
obtaining a second square sum of the target value of the ith row and the jth column and each expected difference value of the ith row, and taking the ratio of the second square sum to L as a third difference value of the ith row and the jth column;
and obtaining a first variance value of the ith row and the jth column based on the second variance value of the ith row and the jth column and the third variance value of the ith row and the jth column.
10. The method of claim 9, wherein the obtaining the first variance value of the ith row and the jth column based on the second variance value of the ith row and the jth column and the third variance value of the ith row and the jth column comprises:
obtaining a sum of the second variance value of the ith row and the jth column and the third variance value of the ith row and the jth column as a fourth variance value of the ith row and the jth column;
And carrying out normalization processing on the fourth difference value of the ith row and the jth column to obtain a first difference value of the ith row and the jth column.
11. A geological structure identification device, comprising:
the first acquisition module is used for acquiring an Ehrlich phase attribute value set of a transmission groove wave of a target area acquired by the seismic channel;
the second acquisition module is used for acquiring the variance of the Ehrlich phase attribute value set;
and the identification module is used for carrying out geological structure identification on the target area based on the variance body.
12. An electronic device, comprising: memory, a processor and a computer program stored on the memory and executable on the processor, which processor, when executing the program, implements a method for identifying a geological structure according to any of claims 1-10.
13. A computer-readable storage medium, on which a computer program is stored, characterized in that the program, when being executed by a processor, implements a method of identifying a geological structure according to any one of claims 1-10.
CN202311017100.0A 2023-08-14 2023-08-14 Geological structure identification method and device, electronic equipment and storage medium Active CN116755154B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202311017100.0A CN116755154B (en) 2023-08-14 2023-08-14 Geological structure identification method and device, electronic equipment and storage medium

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202311017100.0A CN116755154B (en) 2023-08-14 2023-08-14 Geological structure identification method and device, electronic equipment and storage medium

Publications (2)

Publication Number Publication Date
CN116755154A true CN116755154A (en) 2023-09-15
CN116755154B CN116755154B (en) 2023-12-12

Family

ID=87949975

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202311017100.0A Active CN116755154B (en) 2023-08-14 2023-08-14 Geological structure identification method and device, electronic equipment and storage medium

Country Status (1)

Country Link
CN (1) CN116755154B (en)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB201321602D0 (en) * 2012-12-06 2014-01-22 Roxar Software Solutions As Method and system for presenting seismic information
CN106680868A (en) * 2017-01-05 2017-05-17 南京君诺鑫成网络科技有限公司 In-seam seismic multi-attribute parameter tomographic method
CN106772637A (en) * 2016-11-29 2017-05-31 北京中矿大地地球探测工程技术有限公司 One kind is based on slot wave forward probe TVSP scan methods
CN108931812A (en) * 2017-05-24 2018-12-04 河南理工大学 A method of deformation coal is known based on mutually sentencing in transmission beam method slot wave angstrom
CN112904424A (en) * 2021-01-21 2021-06-04 中国科学院地理科学与资源研究所 Method and device for determining position of abnormal body and readable storage medium

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB201321602D0 (en) * 2012-12-06 2014-01-22 Roxar Software Solutions As Method and system for presenting seismic information
CN106772637A (en) * 2016-11-29 2017-05-31 北京中矿大地地球探测工程技术有限公司 One kind is based on slot wave forward probe TVSP scan methods
CN106680868A (en) * 2017-01-05 2017-05-17 南京君诺鑫成网络科技有限公司 In-seam seismic multi-attribute parameter tomographic method
CN108931812A (en) * 2017-05-24 2018-12-04 河南理工大学 A method of deformation coal is known based on mutually sentencing in transmission beam method slot wave angstrom
CN112904424A (en) * 2021-01-21 2021-06-04 中国科学院地理科学与资源研究所 Method and device for determining position of abnormal body and readable storage medium

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
乐勇;王伟;申青春;吕东亮;: "槽波地震勘探技术在工作面小构造探测中的应用", 煤田地质与勘探, vol. 41, no. 04, pages 74 - 77 *
王丹: "三维地震多属性解释在东庞矿采空区识别方法", 煤炭与化工, vol. 41, no. 5, pages 51 - 54 *

Also Published As

Publication number Publication date
CN116755154B (en) 2023-12-12

Similar Documents

Publication Publication Date Title
CN111164462B (en) Artificial source surface wave exploration method, surface wave exploration device and terminal equipment
CN105069248B (en) Large landslide polylith sliding scale
CN109146142A (en) A kind of prediction technique of stope drift active workings mine pressure impression data
CN107765302A (en) Inversion method when time-domain single-frequency waveform independent of source wavelet is walked
US10302758B2 (en) Method and device for detecting discontinuous body with ground penetrating radar
CN104155691A (en) Converted wave anisotropy velocity analysis method and device
JP2014119406A (en) Abnormality detection system and abnormality detection method of columnar structure
CN112231974B (en) Deep learning-based method and system for recovering seismic wave field characteristics of rock breaking seismic source of TBM (Tunnel boring machine)
CN104166680A (en) Parallel vulnerability mining method based on open source library and text mining
CN103869362A (en) Method and equipment for obtaining body curvature
CN116755154B (en) Geological structure identification method and device, electronic equipment and storage medium
JPWO2019049263A1 (en) Simple profiling (SPM) method conversion program and simple profiling (SPM) method conversion method
US20160245938A1 (en) Method and device for increasing frequency of seismic digital signal
CN106154328A (en) A kind of micro-seismic event identification pick-up method based on cross-correlation and system thereof
CN105259578A (en) Formation velocity determining method based on vertical seismic profile data
CN109581504B (en) Method and system for calculating polarity of direct longitudinal wave in ground micro-seismic record
CN111273352A (en) Intelligent detection method and device for geological structure and electronic equipment
CN114814951A (en) Method, device, equipment and storage medium for acquiring slowness of fast and slow transverse waves
CN114706128A (en) Microseism data noise suppression method and system, storage medium and seismic information processing equipment
CN104133239A (en) Method for defining earthquake single-shot record target stratum time window
CN114354762A (en) Coal rock instability destruction precursor information identification method
Dindar et al. Vulnerability index assessment using neural networks (Viann): A case study of Nicosia, Cyprus
CN109143330A (en) A kind of optimization controlled source obtains the method and system of single-shot quality
CN111915844B (en) Method and device for evaluating coal rock stability by analyzing vibration signal through cepstrum coefficient
CN112433251B (en) Correction method and device for stacked profile of Vertical Seismic Profile (VSP) corridor of inclined stratum

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