CN111812731B - Three-dimensional imaging method for fusion of boulder detection resistivity data in subway shield region - Google Patents
Three-dimensional imaging method for fusion of boulder detection resistivity data in subway shield region Download PDFInfo
- Publication number
- CN111812731B CN111812731B CN202010547826.5A CN202010547826A CN111812731B CN 111812731 B CN111812731 B CN 111812731B CN 202010547826 A CN202010547826 A CN 202010547826A CN 111812731 B CN111812731 B CN 111812731B
- Authority
- CN
- China
- Prior art keywords
- data
- dimensional
- resistivity
- subway shield
- hole
- 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
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/38—Processing data, e.g. for analysis, for interpretation, for correction
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention discloses a three-dimensional imaging method for fusing boulder detection resistivity data in a subway shield zone, which comprises the steps of obtaining cross-hole CT and high-density electrical method geophysical prospecting original data in the subway shield zone, and obtaining a cross-hole CT and high-density electrical method two-dimensional resistivity profile after inverting the two original data; extracting cross-hole CT and high-density electrical method data at the same position, and then fusing the two data by using a principal component analysis method; and performing two-dimensional imaging on the resistivity data to obtain a plurality of two-dimensional profile maps of the resistivity data, and converting the two-dimensional coordinates into three-dimensional coordinates to form a three-dimensional model map. The method better depicts the distribution condition of the boulder group in the three-dimensional space, has more obvious relevance among all boulders, and provides guarantee for safe construction and smooth development of the shield.
Description
Technical Field
The invention belongs to the technical field of tunnel advanced prediction, and relates to a three-dimensional imaging method for fusion of resistivity data of boulder detection in a subway shield zone.
Background
The statements in this section merely provide background information related to the present disclosure and may not necessarily constitute prior art.
With the rapid development of the urban rail transit industry, subways meet wide opportunity and development climax in the air, and the shield tunneling technology gradually becomes an important means for subway construction due to the advantages of high tunneling speed, high automation degree and the like, but has poor adaptability to the stratum, so that a cutter head is slightly clamped, deformed or rapidly worn when encountering unexplored boulder distribution in the long-distance tunneling process, a tunnel face is collapsed when the cutter head is seriously damaged, and finally unnecessary casualties and economic losses are caused. Therefore, before shield construction, the occurrence state and distribution range of the boulder group in front of the TBM need to be found out, and guarantee is provided for safe construction and smooth development of the shield.
According to the knowledge of the inventor, the conventional geophysical prospecting method for the boulder in front of the shield machine comprises a ground penetrating radar method, a high-density electrical method, a cross-hole resistivity CT method and the like, but a single geophysical prospecting method has many limitations, false abnormity can be generated during inversion imaging, or good reflection can not be achieved on the distribution range of the boulder, so that combined imaging of multiple geophysical prospecting methods is necessary, and the conventional imaging is usually two-dimensional and cannot reflect the distribution condition of the boulder in space.
Disclosure of Invention
The invention aims to solve the problems and provides a three-dimensional imaging method for detecting resistivity data fusion of boulders in a subway shield zone.
According to some embodiments, the following technical scheme is adopted in the disclosure:
a three-dimensional imaging method for fusion of detection resistivity data of boulders in a subway shield zone comprises the following steps:
acquiring cross-hole CT and high-density electrical method geophysical prospecting original data of a subway shield interval, and after inverting the two original data, acquiring a cross-hole CT and high-density electrical method two-dimensional resistivity profile;
extracting cross-hole CT and high-density electrical method data at the same position, and then fusing the two data by using a principal component analysis method;
and performing two-dimensional imaging on the resistivity data to obtain a plurality of two-dimensional profile maps of the resistivity data, and converting the two-dimensional coordinates into three-dimensional coordinates to form a three-dimensional model map.
As an alternative embodiment, at least two measuring lines in the horizontal direction and the vertical direction are arranged in the subway shield zone, and the two measuring lines are crossed.
As an alternative embodiment, a plurality of high-density measuring lines are symmetrically distributed on the axial line of the tunnel respectively, a plurality of rows of drill holes are distributed at the edge of the detection area and the center of the detection area, and cross-hole resistivity CT detection is implemented by utilizing the drill holes.
As an alternative embodiment, the cross-hole CT method resistivity data points are subjected to coordinate transformation by taking the resistivity data points obtained by high-density electrical detection as a reference.
As an alternative embodiment, in the process of utilizing the principal component analysis method, the resistivity sample data obtained by the cross-hole method and the high-density electrical method are standardized; then solving a covariance matrix between the centralized variables, and judging whether the deviation change trends of the two variables are consistent; then, obtaining eigenvalues and eigenvectors of covariance, arranging the eigenvalues in a sequence from small to large, selecting the largest one of the eigenvalues, and obtaining the corresponding eigenvector; and finally, projecting the sample points of the centralized data to the eigenvector base with the largest eigenvalue to obtain a fused resistivity data result.
As an alternative embodiment, after the two data are fused, two-dimensional imaging of the resistivity data is performed to obtain a plurality of two-dimensional profile maps of the resistivity data, and the obtained two-dimensional data is converted into three-dimensional data by a coordinate conversion method.
As an alternative embodiment, a kriging interpolation method is used to comprehensively image data points in a plurality of three-dimensional coordinates to form a three-dimensional model map.
A three-dimensional imaging system for detecting resistivity data fusion of boulders in a subway shield zone comprises:
the system comprises a plurality of measuring lines, a plurality of sensors and a controller, wherein the measuring lines are arranged in a subway shield interval and used for carrying out cross-hole CT and high-density electrical geophysical prospecting work to obtain two detection data;
the data inversion module is configured to invert the two kinds of original data to obtain a cross-hole CT and high-density electrical method two-dimensional resistivity profile;
the data fusion module is configured to extract cross-hole CT and high-density electrical method data at the same position and then fuse the two data by using a principal component analysis method;
the imaging module is configured to perform resistivity data two-dimensional imaging, obtain a plurality of resistivity data two-dimensional section maps, convert two-dimensional coordinates into three-dimensional coordinates, and form a three-dimensional model map.
A computer readable storage medium stores a plurality of instructions, and the instructions are suitable for being loaded by a processor of terminal equipment and executing the method for three-dimensional imaging of the fusion of the resistivity data of the orphan rock detection in the subway shield zone.
A terminal device comprising a processor and a computer readable storage medium, the processor being configured to implement instructions; the computer-readable storage medium is used for storing a plurality of instructions, and the instructions are suitable for being loaded by a processor and executing the subway shield interval solitary stone exploration resistivity data fusion three-dimensional imaging method.
Compared with the prior art, the beneficial effect of this disclosure is:
the invention provides a three-dimensional imaging method for fusing boulder exploration resistivity data in a metro shield zone, compared with the prior art, the method fuses cross-hole CT data and high-density data at the same position, solves the problem that cross-hole resistivity CT is poor in imaging near a drill hole, reduces false abnormal body distribution near a high-density electrical method electrode, improves vertical and horizontal resolutions of the high-density electrical method and the cross-hole resistivity CT method, performs three-dimensional imaging on two-dimensional exploration data fusion by using a three-dimensional coordinate conversion formula, better describes the distribution condition of boulder groups in a three-dimensional space, has more obvious relevance among boulders, and provides guarantee for safe construction and smooth development of a shield.
Drawings
The accompanying drawings, which are included to provide a further understanding of the disclosure, illustrate embodiments of the disclosure and together with the description serve to explain the disclosure without limiting the disclosure in any way.
FIG. 1 is a schematic diagram of the position of a geophysical prospecting boulder;
FIG. 2 is a flow chart;
wherein: 1. high-density survey line, 2 cross-hole CT survey line, 3 drilling, 4 boulder, 5 TBM, 6 survey line coincidence.
The specific implementation mode is as follows:
the present disclosure is further described with reference to the following drawings and examples.
It should be noted that the following detailed description is exemplary and is intended to provide further explanation of the disclosure. Unless defined otherwise, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this disclosure belongs.
It is noted that the terminology used herein is for the purpose of describing particular embodiments only and is not intended to be limiting of example embodiments according to the present disclosure. As used herein, the singular forms "a", "an", and/or "the" are intended to include the plural forms as well, unless the context clearly indicates otherwise, and it should be understood that when the terms "comprises" and/or "comprising" are used in this specification, they specify the presence of the features, steps, operations, devices, components, and/or combinations thereof.
As shown in fig. 1, a schematic diagram of a position of a geophysical prospecting boulder according to the present disclosure is used for better positioning a distribution situation of the boulder of each tunnel line, and performing high-density electrical method and cross-hole CT exploration work in front of a TBM 5. 5 high-density measuring lines 1 are symmetrically distributed about the axial line of the tunnel, 2 rows of drill holes 3 are distributed at the edge of the detection area and the center of the detection area, cross-hole resistivity CT detection is implemented by utilizing the drill holes 3, and the arrangement mode is shown as cross-hole CT measuring lines 2. A measuring line superposition part 6 exists between the high-density measuring line 1 and the cross-hole CT measuring line 2, so that data fusion can be carried out, and the distribution range of the boulder can be better detected.
As shown in fig. 2, a specific process of the present disclosure includes:
1. geophysical prospecting work of a high-density electrical method and a cross-hole CT method is carried out in front of tunnel construction, as shown in figure 1, when measuring lines of the two methods are arranged, a superposed part is needed, and thus data fusion work can be carried out. After the original data is subjected to necessary processing, a two-dimensional resistivity map is obtained through inversion, and the two-dimensional resistivity map can be subjected to preliminary analysis to obtain the approximate distribution range of the boulder.
2. And (3) performing coordinate transformation on the cross-hole CT method resistivity data points by taking the resistivity data points obtained by high-density electrical method detection as a reference. Because the positions of the starting points of the measuring lines of the cross-hole CT detection method and the high-density electrical detection method can be different, the initial coordinates of the two methods are different. And extracting resistivity data at the same position at the part where the high-density measuring line is overlapped with the cross-hole CT measuring line after conversion.
3. And determining proper weights of the two components by using a principal component analysis method for data fusion. Principle of principal component analysis: performing data centralization, namely standardizing resistivity sample data obtained by a cross-hole method and a high-density electrical method; then solving a covariance matrix between the centralized variables, and judging whether the deviation change trends of the two variables are consistent; then, characteristic values and characteristic vectors of the covariance are obtained, the characteristic values are arranged in a sequence from small to large, the largest one is selected, and the corresponding characteristic vector is obtained; and finally, projecting the sample points of the centralized data to the eigenvector base with the largest eigenvalue to obtain a fused resistivity data result, which can be regarded as the comprehensive property of the cross-hole method detection data and the high-density method detection data.
4. The extracted two-dimensional data is converted into three-dimensional data by the following formula
Wherein X, Y, Z is the final three-dimensional coordinate with O as the origin, X1 is the horizontal distance from the origin of the survey line to the origin of the coordinate, Y1 is the longitudinal distance from the origin of the survey line to the origin of the coordinate, L is the span hole distance, X ' is the horizontal coordinate of the original data point, Y ' is the depth coordinate of the original data point, the initial value of Z ' is 0, R1 is the position matrix, and R2 is the data point matrix.
5. After the resistivity data points under the three-dimensional coordinate system are obtained, comprehensive imaging can be carried out on the data points under the three-dimensional coordinate system, and the schematic diagram of the solitary stone distribution three-dimensional model is obtained by utilizing a kriging interpolation method.
The following product examples are also provided:
a three-dimensional imaging system for detecting resistivity data fusion of boulders in a subway shield zone comprises:
the system comprises a plurality of measuring lines, a plurality of sensors and a controller, wherein the measuring lines are arranged in a subway shield interval and used for carrying out cross-hole CT and high-density electrical geophysical prospecting work to obtain two detection data;
the data inversion module is configured to invert the two kinds of original data to obtain a cross-hole CT and high-density electrical method two-dimensional resistivity profile;
the data fusion module is configured to extract cross-hole CT and high-density electrical method data at the same position and then fuse the two data by using a principal component analysis method;
the imaging module is configured to perform resistivity data two-dimensional imaging, obtain a plurality of resistivity data two-dimensional section maps, convert two-dimensional coordinates into three-dimensional coordinates, and form a three-dimensional model map.
A computer readable storage medium stores a plurality of instructions, and the instructions are suitable for being loaded by a processor of terminal equipment and executing the method for three-dimensional imaging of the fusion of the resistivity data of the orphan rock detection in the subway shield zone.
A terminal device comprising a processor and a computer readable storage medium, the processor being configured to implement instructions; the computer-readable storage medium is used for storing a plurality of instructions, and the instructions are suitable for being loaded by a processor and executing the subway shield interval solitary stone exploration resistivity data fusion three-dimensional imaging method.
As will be appreciated by one skilled in the art, embodiments of the present disclosure may be provided as a method, system, or computer program product. Accordingly, the present disclosure may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present disclosure 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 so forth) having computer-usable program code embodied therein.
The present disclosure is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the disclosure. 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 programs may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data 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 data 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 data 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.
The above description is only a preferred embodiment of the present disclosure and is not intended to limit the present disclosure, and various modifications and changes may be made to the present disclosure by those skilled in the art. Any modification, equivalent replacement, improvement and the like made within the spirit and principle of the present disclosure should be included in the protection scope of the present disclosure.
Although the present disclosure has been described with reference to specific embodiments, it should be understood that the scope of the present disclosure is not limited thereto, and those skilled in the art will appreciate that various modifications and changes can be made without departing from the spirit and scope of the present disclosure.
Claims (10)
1. A three-dimensional imaging method for fusion of resistivity data of boulder detection in a subway shield zone is characterized by comprising the following steps of: the method comprises the following steps:
acquiring cross-hole CT and high-density electrical method geophysical prospecting original data of a subway shield interval, and inverting the two original data to obtain a cross-hole CT and high-density electrical method two-dimensional resistivity profile;
extracting cross-hole CT and high-density electrical method data at the same position, and then fusing the two data by using a principal component analysis method;
and performing two-dimensional imaging on the resistivity data to obtain a plurality of two-dimensional profile maps of the resistivity data, and converting the two-dimensional coordinates into three-dimensional coordinates to form a three-dimensional model map.
2. The subway shield zone boulder detection resistivity data fusion three-dimensional imaging method as claimed in claim 1, characterized in that: at least two measuring lines in the horizontal direction and the vertical direction are arranged in the subway shield region, and the two measuring lines are crossed.
3. The subway shield zone boulder detection resistivity data fusion three-dimensional imaging method as claimed in claim 1, characterized in that: and respectively and symmetrically arranging a plurality of high-density measuring lines on the central axis of the tunnel, arranging a plurality of rows of drilling holes on the edge of the detection area and the center of the detection area, and performing cross-hole CT resistivity detection by using the drilling holes.
4. The subway shield zone boulder detection resistivity data fusion three-dimensional imaging method as claimed in claim 1, characterized in that: and (3) performing coordinate transformation on the cross-hole CT method resistivity data points by taking the resistivity data points obtained by high-density electrical method detection as a reference.
5. The subway shield zone boulder detection resistivity data fusion three-dimensional imaging method as claimed in claim 1, characterized in that: in the process of utilizing the principal component analysis method, standardizing resistivity sample data obtained by a cross-hole CT method and a high-density electrical method; solving a covariance matrix between the centralized variables, and judging whether the deviation change trends of the two variables are consistent; then, the eigenvalue and the eigenvector of the covariance are solved, the eigenvalues are arranged according to the sequence from small to large, and the largest one is selected to solve the corresponding eigenvector; and finally, projecting the sample points of the centralized data to the eigenvector base with the largest eigenvalue to obtain a fused resistivity data result.
6. The subway shield zone boulder detection resistivity data fusion three-dimensional imaging method as claimed in claim 1, characterized in that: after the two kinds of data are fused, two-dimensional imaging of the resistivity data is carried out firstly, a plurality of two-dimensional profile maps of the resistivity data are obtained, and the obtained two-dimensional data are converted into three-dimensional data through a coordinate conversion method.
7. The subway shield zone boulder detection resistivity data fusion three-dimensional imaging method as claimed in claim 1, characterized in that: and comprehensively imaging a plurality of data points under the three-dimensional coordinates by using a kriging interpolation method to form a three-dimensional model map.
8. A three-dimensional imaging system for detecting resistivity data fusion of boulders in a subway shield region is characterized in that: the method comprises the following steps:
the system comprises a plurality of measuring lines, a plurality of sensors and a controller, wherein the measuring lines are distributed in a subway shield interval and used for carrying out cross-hole CT and high-density electrical geophysical prospecting work to obtain two kinds of detection data;
the data inversion module is configured to invert the two kinds of original data to obtain a cross-hole CT and high-density electrical method two-dimensional resistivity profile;
the data fusion module is configured to extract cross-hole CT and high-density electrical method data at the same position and then fuse the two data by using a principal component analysis method;
the imaging module is configured to perform resistivity data two-dimensional imaging, obtain a plurality of resistivity data two-dimensional section maps, convert two-dimensional coordinates into three-dimensional coordinates, and form a three-dimensional model map.
9. A computer-readable storage medium characterized by: a plurality of instructions are stored, and the instructions are suitable for being loaded by a processor of terminal equipment and executing the method for three-dimensional imaging of the fusion of the detection resistivity data of the orphan rock in the shield zone of the subway shield according to any one of claims 1-7.
10. A terminal device is characterized in that: the system comprises a processor and a computer readable storage medium, wherein the processor is used for realizing instructions; the computer readable storage medium is used for storing a plurality of instructions, and the instructions are suitable for being loaded by a processor and executing the subway shield interval solitary stone exploration resistivity data fusion three-dimensional imaging method as set forth in any one of claims 1-7.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010547826.5A CN111812731B (en) | 2020-06-16 | 2020-06-16 | Three-dimensional imaging method for fusion of boulder detection resistivity data in subway shield region |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010547826.5A CN111812731B (en) | 2020-06-16 | 2020-06-16 | Three-dimensional imaging method for fusion of boulder detection resistivity data in subway shield region |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111812731A CN111812731A (en) | 2020-10-23 |
CN111812731B true CN111812731B (en) | 2021-07-06 |
Family
ID=72845196
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010547826.5A Active CN111812731B (en) | 2020-06-16 | 2020-06-16 | Three-dimensional imaging method for fusion of boulder detection resistivity data in subway shield region |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111812731B (en) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US11549359B2 (en) * | 2020-05-11 | 2023-01-10 | Saudi Arabian Oil Company | Systems and methods to identify and quantify field development opportunities through integration of surface and sub-surface data |
CN112305624B (en) * | 2020-10-30 | 2024-02-06 | 中国地质调查局西安地质调查中心(西北地质科技创新中心) | Non-contact obstacle-crossing high-density measurement method |
CN112799139A (en) * | 2020-12-29 | 2021-05-14 | 河北煤炭科学研究院有限公司 | Three-dimensional detection method and device for water-applying abnormal area |
CN113552635B (en) * | 2021-06-25 | 2023-03-03 | 山东大学 | Three-dimensional advanced forecasting method and system for fixed-source transient electromagnetic data fusion tunnel |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20160124108A1 (en) * | 2014-10-30 | 2016-05-05 | Schlumberger Technology Corporation | Inversion Technique For Fracture Characterization In Highly Inclined Wells Using Multiaxial Induction Measurements |
CN104407394B (en) * | 2014-11-26 | 2015-07-22 | 山东大学 | Subway shield section boulder refining detecting method based on resistivity trans-hole CT (computed tomography) |
US20180106708A1 (en) * | 2015-05-20 | 2018-04-19 | Schlumberger Technology Corporation | Hydraulic fracturability index using high resolution core measurements |
CN105807336A (en) * | 2016-03-11 | 2016-07-27 | 山东大学 | Comprehensive boulder detection model test device suitable for plurality of earth physical detection methods, and test method thereof |
CN105804763B (en) * | 2016-03-11 | 2017-03-29 | 山东大学 | Across the hole CT boulders detection system of advanced 3 D resistivity and detection method that shield excavation machine is carried |
CN110361792B (en) * | 2019-06-24 | 2020-08-25 | 山东大学 | Method, medium and equipment for fusing and imaging geophysical data |
CN110823962B (en) * | 2019-11-14 | 2020-11-06 | 山东大学 | Three-dimensional imaging method and system for landslide mass |
CN110988999A (en) * | 2019-12-02 | 2020-04-10 | 山东大学 | Detection method and system for analyzing pile foundation based on cross-hole resistivity method CT inversion imaging |
-
2020
- 2020-06-16 CN CN202010547826.5A patent/CN111812731B/en active Active
Also Published As
Publication number | Publication date |
---|---|
CN111812731A (en) | 2020-10-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111812731B (en) | Three-dimensional imaging method for fusion of boulder detection resistivity data in subway shield region | |
CN110109895B (en) | Surrounding rock grading combined prediction method suitable for TBM tunneling tunnel and application | |
CN108345049B (en) | Underground engineering unfavorable geological detection multi-method constraint inversion and joint interpretation method | |
WO2021243967A1 (en) | Three-dimensional resistivity tomography method and system | |
CN110119994B (en) | Metal ore blasting pile displacement extraction and prediction method supported by GIS | |
CN112114359B (en) | Dangerous area detection method, system and terminal based on active and passive seismic source signals and readable storage medium | |
US20170146690A1 (en) | Method for generating secondary data in geostatistics using observed data | |
CN111812730B (en) | Resistivity data fusion three-dimensional imaging method and system for landslide detection | |
CN112580165B (en) | Method and system for predicting card jamming of open TBM (tunnel boring machine) through unfavorable geological cutter | |
CN109212592A (en) | A kind of prestack orientation P wave method for detecting anisotropic fracture and device | |
AU2019404613A1 (en) | Automated updating of geological model boundaries for improved ore extraction | |
CN111986321B (en) | Three-dimensional geological modeling method, device, equipment and storage medium | |
CN116524017B (en) | Underground detection, identification and positioning system for mine | |
CN110632667B (en) | Hidden collapse column advanced detection method based on shock wave shock condition | |
CN109630089B (en) | Horizontal well geological structure recognition method and device | |
CN104777513A (en) | Non-continuity boundary detection method for seismic data gradient information | |
CN103670382A (en) | Drill core direction recovery method | |
CN109444960A (en) | Determination method for fault tectonic | |
CN117331124A (en) | Tunnel seismic wave detection data preprocessing method, device and storage device | |
CN115239108B (en) | Weak broken surrounding rock sensing method based on TBM real-time broken rock data | |
CN111273352A (en) | Intelligent detection method and device for geological structure and electronic equipment | |
CN113077546B (en) | Automatic search ellipsoid setting method for mineral resource reserve estimation | |
Brzovic et al. | Integrated photogrammetry and discrete fracture network modeling to determine rock structure around excavations at the El Teniente mine | |
CN105840175A (en) | Three-dimensional fracture modeling method | |
CN206960672U (en) | A kind of wireless controlled passes the Tunnel prediction system of acceleration transducer |
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 |