CN112114370B - Multi-scale analysis method for applying aviation geophysical prospecting three-dimensional inversion result to traffic corridor - Google Patents
Multi-scale analysis method for applying aviation geophysical prospecting three-dimensional inversion result to traffic corridor Download PDFInfo
- Publication number
- CN112114370B CN112114370B CN202010595709.6A CN202010595709A CN112114370B CN 112114370 B CN112114370 B CN 112114370B CN 202010595709 A CN202010595709 A CN 202010595709A CN 112114370 B CN112114370 B CN 112114370B
- Authority
- CN
- China
- Prior art keywords
- resistance
- area
- rock
- low
- interpretation
- 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
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
-
- 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/08—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with magnetic or electric fields produced or modified by objects or geological structures or by detecting devices
- G01V3/081—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with magnetic or electric fields produced or modified by objects or geological structures or by detecting devices the magnetic field is produced by the objects or geological structures
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Electromagnetism (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention belongs to the technical field of geophysical exploration, and relates to a multi-scale analysis method for applying an aviation geophysical prospecting three-dimensional inversion result to a traffic corridor. To address the deficiencies of the prior interpretation analysis. The method provided by the invention is to utilize the three-dimensional inversion result of the aerospace geophysical prospecting, perform large-scale analysis and interpretation for geological structure units, perform mesoscale analysis and interpretation for specific structural bands and stratum lithology, perform small-scale analysis and interpretation for engineering geological features such as special rock stratum, joint intensive bands, broken water-rich bands, karst development areas and the like, and obtain the basic geological conditions of the tunnel engineering area through large-medium-small multi-scale comprehensive interpretation and analysis.
Description
Technical field:
the invention belongs to the technical field of geophysical exploration, and relates to a multi-scale analysis method for applying an aviation geophysical prospecting three-dimensional inversion result to a traffic corridor.
The background technology is as follows:
the linear engineering in the mountain area of the plateau is mainly a long and large deeply buried tunnel engineering, the basic geological conditions in the traffic corridor need to be ascertained in the line selection stage, the overall stability and the safety and reliability of the traffic corridor are guaranteed, a basis is provided for scientific planning and timely optimization of a line scheme, data is provided for geological condition evaluation of the tunnel engineering, and the purpose of rapidly and efficiently analyzing geological problems in the tunnel engineering area is achieved through assistance.
For a long time, the aviation geophysical prospecting method at home and abroad is mainly applied to the field of regional geological investigation or prospecting, the obtained aviation geophysical prospecting result data mainly comprise low-precision and small scale, and the gap of the domestic aviation electromagnetic result with high precision and large-medium scale is particularly huge. The analysis of the aviation geophysical prospecting data in China is generally qualitative analysis; the aeromagnetic data inversion is limited to the inversion of a local two-dimensional platy body model, the application is very limited, and the three-dimensional inversion is only in the theoretical research and numerical simulation stage at present; the aviation electromagnetic method is in experimental research stage for a long time, and no production example is seen.
The area range related to the aerospace geophysical prospecting in the field of regional geological investigation and prospecting is limited, the structural units are single, large-scale analysis of multiple structural units is not needed to be studied, stratum, lithology and structural analysis and ore forming condition analysis are generally only carried out on the aerospace geophysical prospecting data, and the method belongs to mesoscale analysis and interpretation.
The traffic corridor belongs to linear engineering, has long length and narrow width, spans more construction units, but tracks stratum, lithology and construction in the direction perpendicular to the direction of the line, has large limitation, and needs large-scale analysis of the construction units as a basis to support mesoscale analysis of the stratum, lithology and construction; meanwhile, the engineering problem finally relates to engineering geological condition judgment, belongs to small-scale analysis, and needs to reasonably judge on the basis of stratum, lithology and structure analysis by taking a structural unit as a background; the multi-scale analysis flow of the aero-geophysical prospecting data from large to small and from whole to detail has obvious geological analysis characteristics of the traffic corridor, and no precedent exists at present.
In the railway route selection stage of a mountain area in the southeast of the Tibetan, an aviation transient electromagnetic method (VTEM) and an aviation magnetotelluric method (ZTEM) are introduced aiming at complex terrain and geological features, and aviation magnetic method (TMI) detection is carried out, so that the high-precision measurement of the full coverage of the aviation electromagnetic method of the 370 km traffic corridor is realized. When aviation geophysical prospecting data are processed, a ZTEM and TMI three-dimensional inversion technology is successfully applied, a high-precision stereoscopic result chart of the resistivity and the magnetic susceptibility of the railway traffic corridor is formed, and brand-new basic data and realization possibility are provided for multi-scale comprehensive analysis of geological conditions of the traffic corridor.
The invention comprises the following steps:
in view of the above, the invention provides a multi-scale analysis method for applying the three-dimensional inversion result of the aerospace geophysical prospecting to the traffic corridor in order to solve the defects existing in the conventional interpretation analysis.
In order to solve the problems in the prior art, the technical scheme of the invention is that the aviation geophysical prospecting three-dimensional inversion result is applied to a multi-scale analysis method of a traffic corridor, and is characterized in that: the method comprises the following steps: and performing large-scale analysis and interpretation on geological structure units, performing mesoscale analysis and interpretation on specific structural bands and formation lithology, performing small-scale analysis and interpretation on engineering geological features such as special rock stratum, joint intensive bands, broken water-rich bands, karst development areas and the like, and obtaining basic geological conditions of the tunnel engineering area through multi-scale comprehensive interpretation and analysis of large, medium and small.
Further, the method comprises the following specific steps:
1) Extracting a typical section three-dimensional inversion magnetic susceptibility result graph of the traffic corridor, and performing macroscopic interpretation and large-scale analysis according to the integral size of the magnetic susceptibility and the strong and weak magnetic body distribution characteristic interpretation to divide a geological structure unit of the corridor and a stratum lithology large area;
2) Extracting a tunnel longitudinal section three-dimensional inversion resistivity result graph of a corresponding region, comparing the three-dimensional susceptibility result graph, calibrating a fracture structure according to a vertical connected resistivity or susceptibility gradient zone and a through type low stop band, dividing lithology rock groups in the tunnel section according to distribution characteristics of high and low resistance regions, and carrying out refinement interpretation and mesoscale analysis;
3) On the basis of large-scale and medium-scale analysis and interpretation, the non-vertical connected low-resistance region of the three-dimensional resistivity result pattern is further refined and analyzed, a special stratum, a joint dense zone, a water-rich zone and a karst region are divided by combining the low-resistance form and the lithology rock group, and fine interpretation and small-scale analysis are performed to obtain the basic geological condition of the tunnel engineering region.
Further, the method for performing macroscopic interpretation and large scale analysis in step 1) is as follows: the magnetic susceptibility result map is inverted three-dimensionally by utilizing a typical section of a traffic corridor, and land block construction units and normal sedimentary rock and negative metamorphic rock large areas are defined according to the characteristic of 'overall nonmagnetic-weak magnetism'; according to the fact that the strong magnetic bodies and the weak magnetic bodies are distributed alternately in the whole, the strong magnetic bodies in the central area are represented by root lumps and comprise extremely strong magnetic bodies, and the strong magnetic bodies in the edge area are represented by root-free lumps or floating bodies, so that the seam belt construction units and the mixed rock large area are defined; according to the characteristic that the whole strong magnetic bodies and the weak magnetic bodies are alternately distributed, the strong magnetic bodies are presented in a rooted block, and the characteristic that the strong magnetic bodies are not contained is used for defining a magma arc construction unit and invading a rock mass area.
Further, the method for performing the refined interpretation and the mesoscale analysis in the step 2) comprises the following steps: comparing the tunnel longitudinal section three-dimensional inversion resistivity result graph with the longitudinal section three-dimensional inversion susceptibility result graph which is analyzed and interpreted in a large scale by the step 1), and calibrating interlayer fracture according to the characteristic of vertical connected resistivity or susceptibility gradient zone (abrupt zone); calibrating the internal faults according to the characteristics of 'through type low stop band from the earth surface to the deep'; dividing lithology rock groups of the land block construction units according to the characteristic that the high-resistance area and the low-resistance area are separated by an obvious vertical communicated gradient zone; dividing lithology rock groups of the suture belt construction units according to the characteristic that the high-resistance area and the low-resistance area are mutually cut by irregular boundaries, and the resistivity lump and the magnetic susceptibility lump approximately correspond to each other; and dividing lithology rock groups of the magma arc construction units according to the characteristic that the high-resistance area and the low-resistance area are mutually wrapped by smooth boundary lines and the resistivity block and the magnetic susceptibility block have no corresponding rule.
Further, the specific method for performing fine interpretation and small scale analysis in step 3): the tunnel longitudinal section three-dimensional inversion resistivity result graph which is analyzed through the mesoscale analysis is completed through the steps 1) and 2) is utilized to analyze the 'non-vertical connected' low-resistance morphology in each lithology rock group in detail, the sedimentary rock area and the negative metamorphic rock area are continuously and directionally spread in a large range, and the low-resistance area which is irregular in shape but can see the spreading trend and trend is translated into special stratum such as carbonaceous shale, carbonaceous slate, mudstone and the like; the lens-shaped or beaded distribution of the intrusion rock area, the metamorphic hard rock area and the deposition harder rock area in a smaller range appears as a low-resistance ring in the whole high-resistance area or an ultra-low-resistance ring in the whole low-resistance area, and the low-resistance area of which the spreading trend and the trend are difficult to distinguish is interpreted as a joint dense belt and a broken water rich belt; the low-resistance region with the dissoluble rock region being distributed in a bead shape in the whole high-resistance block, and having no spreading trend but a certain vertical zonal characteristic is interpreted as a karst development region.
Compared with the prior art, the invention has the following advantages:
1) Based on aviation electromagnetic and aviation magnetic combined detection data, the invention forms a set of 'big-middle-small' multi-scale comprehensive analysis interpretation method by utilizing three-dimensional inversion magnetic susceptibility and resistivity achievements, and achieves the purpose of rapidly and efficiently analyzing tunnel engineering geology, thereby providing basis for line selection and scheme optimization of long and large deep buried tunnels in complicated and difficult mountain areas, controlling tunnel engineering geology risks and reducing engineering construction cost.
Description of the drawings:
FIG. 1 is a contour plot of the three-dimensional inversion susceptibility of a traffic corridor (closure strip segment) representative section TMI;
FIG. 2 is a contour plot of TMI three-dimensionally inverted susceptibility of a typical section of a traffic corridor (block);
FIG. 3 is a contour plot of the three-dimensional inversion susceptibility of a traffic corridor (magma arc segment) representative section TMI;
FIG. 4 is a three-dimensional inverted resistivity contour plot of a tunnel profile ZTEM (joint strip segment);
FIG. 5 is a three-dimensional inverted resistivity contour plot of a tunnel profile ZTEM;
fig. 6 is a three-dimensional inversion resistivity contour plot of a tunnel profile ZTEM (magma arc segment).
In the figure: HT 1-HT 9 medium-high susceptibility regions (blocks); HE 1-HE 5 high resistivity regions (blocks); LE1 to LE15 low resistivity regions (blocks); 10/f11/f12 general fault numbering; F17/F18/F19 regional cleavage number.
The specific embodiment is as follows:
the present invention will be described in further detail with reference to the drawings and examples, in order to make the objects, technical solutions and advantages of the present invention more apparent. It should be understood that the specific embodiments described herein are for purposes of illustration only and are not intended to limit the scope of the invention.
The multi-scale analysis method of the aviation geophysical prospecting three-dimensional inversion result applied to the traffic corridor comprises the following steps: and performing large-scale analysis and interpretation on geological structure units, performing mesoscale analysis and interpretation on specific structural bands and formation lithology, performing small-scale analysis and interpretation on engineering geological features such as special rock stratum, joint intensive bands, broken water-rich bands, karst development areas and the like, and obtaining basic geological conditions of the tunnel engineering area through multi-scale comprehensive interpretation and analysis of large, medium and small.
The method comprises the following specific steps:
1) Extracting a three-dimensional magnetic susceptibility result map of a typical section of a traffic corridor, dividing a gallery geological structure unit and a stratum lithology large area according to the integral size of the magnetic susceptibility and the strong and weak magnetic body distribution characteristic interpretation, performing macroscopic interpretation and large-scale analysis, and providing a geological structure background and lithology analysis basis for the mesoscale analysis;
the method for carrying out macroscopic interpretation and large-scale analysis comprises the following steps: the magnetic susceptibility result map is inverted three-dimensionally by utilizing a typical section of a traffic corridor, and land block construction units and normal sedimentary rock and negative metamorphic rock large areas are defined according to the characteristic of 'overall nonmagnetic-weak magnetism'; according to the fact that the strong magnetic bodies and the weak magnetic bodies are distributed alternately in the whole, the strong magnetic bodies in the central area are represented by root lumps and comprise extremely strong magnetic bodies, and the strong magnetic bodies in the edge area are represented by root-free lumps or floating bodies, so that the seam belt construction units and the mixed rock large area are defined; according to the characteristic that the whole strong magnetic bodies and the weak magnetic bodies are alternately distributed, the strong magnetic bodies are presented in a rooted block, and the characteristic that the strong magnetic bodies are not contained is used for defining a magma arc construction unit and invading a rock mass area.
2) Extracting a tunnel longitudinal section three-dimensional inversion resistivity result graph of a corresponding region, comparing the three-dimensional susceptibility result graph, calibrating a fracture structure according to a vertical connected resistivity or susceptibility gradient zone (abrupt zone) and a through low-stop zone, dividing lithology rock groups in a tunnel section according to distribution characteristics of high and low resistance zones (blocks), performing refined interpretation and mesoscale analysis, and providing specific physical characteristic reference for small-scale analysis;
the method for carrying out refinement interpretation and mesoscale analysis comprises the following steps: comparing the tunnel longitudinal section three-dimensional inversion resistivity result graph with the longitudinal section three-dimensional inversion susceptibility result graph which is analyzed and interpreted in a large scale by the step 1), and calibrating interlayer fracture according to the characteristic of vertical connected resistivity or susceptibility gradient zone (abrupt zone); calibrating the internal faults according to the characteristics of 'through type low stop band from the earth surface to the deep'; dividing lithology rock groups of the land block construction units according to the characteristic that the high-resistance area and the low-resistance area are separated by an obvious vertical communicated gradient zone; dividing lithology rock groups of the suture belt construction units according to the characteristic that the high-resistance area and the low-resistance area are mutually cut by irregular boundaries, and the resistivity lump and the magnetic susceptibility lump approximately correspond to each other; and dividing lithology rock groups of the magma arc construction units according to the characteristic that the high-resistance area and the low-resistance area are mutually wrapped by smooth boundary lines and the resistivity block and the magnetic susceptibility block have no corresponding rule.
3) On the basis of large-scale and medium-scale analysis and interpretation, the non-vertical connected low-resistance area in the three-dimensional resistivity result graph is further refined and analyzed, a special stratum, a joint dense zone, a water-rich zone, a karst area and the like are divided by combining the low-resistance form and the lithology rock group, fine interpretation and small-scale analysis are performed, basic geological conditions of a tunnel engineering area are obtained, and reliable data support is provided for tunnel surrounding rock division and deep engineering geological problem assessment.
The method for fine interpretation and small-scale analysis comprises the following steps: the tunnel longitudinal section three-dimensional inversion resistivity result graph which is analyzed through the mesoscale analysis is completed through the steps 1) and 2) is utilized to analyze the 'non-vertical connected' low-resistance morphology in each lithology rock group in detail, the sedimentary rock area and the negative metamorphic rock area are continuously and directionally spread in a large range, and the low-resistance area which is irregular in shape but can see the spreading trend and trend is translated into special stratum such as carbonaceous shale, carbonaceous slate, mudstone and the like; the lens-shaped or beaded distribution of the intrusion rock area, the metamorphic hard rock area and the deposition harder rock area in a smaller range appears as a low-resistance ring in the whole high-resistance area or an ultra-low-resistance ring in the whole low-resistance area, and the low-resistance area of which the spreading trend and the trend are difficult to distinguish is interpreted as a joint dense belt and a broken water rich belt; the low-resistance region with the dissoluble rock region being distributed in a bead shape in the whole high-resistance block, and having no spreading trend but a certain vertical zonal characteristic is interpreted as a karst development region.
Examples: taking a complex mountain railway engineering of a certain altitude as an example:
a multi-scale analysis method for applying the three-dimensional inversion result of the aviation geophysical prospecting to the traffic corridor comprises the following steps:
1. large scale macroscopic interpretation geological structure unit and lithology large area
The geological structure units through which the traffic corridor passes mainly comprise plots (relatively complete sedimentary or metamorphic areas which are not cut or cracked since the geological history evolves), seam strips (mixed rock strips which are mainly formed by marine land transition and comprise marine ridge expansion strips and are subjected to multiple times of sea-land transition in the geological history evolvement process), magma arcs (large mountain bodies or mountain systems which are mainly formed by invaded rocks and are formed by strong magma activities in the crust lifting and mountain making processes). The invention mainly analyzes the three geological structure units and lithology large areas contained in the geological structure units according to the macroscopic interpretation of the aeromagnetic MVI three-dimensional inversion result so as to realize large-scale analysis.
1.1 land elements
The ground is in a relatively stable environment in the geological structure process, and the stratum is mainly formed by normal sedimentary rocks and negative metamorphic rocks (metamorphic rocks are formed by the sedimentary rocks) and is characterized by no magnetism and weak magnetism. Taking the f 10-f 11 sections in FIG. 1 and FIG. 2 as an example, the three-dimensional inversion susceptibility contour map of the typical section of the gallery is wholly represented as non-weakly magnetic characteristics of 0-50, reflecting the relatively stable diagenetic environment and lithology large area characteristics mainly comprising sedimentary rocks and negative metamorphic rocks; the magnetic susceptibility contour map is in a gentle form and has no response to the size fracture structure (F27-F30).
1.2 seam tape Unit
In the joint belt unit, strong magnetic bodies and weak magnetic bodies are distributed alternately in many ways, the strong magnetic bodies in the edge area of the joint belt are usually represented by ungrooved lumps or floating bodies, and the strong magnetic bodies in the central area are represented by rooted lumps and contain extremely strong magnetic body lumps. Taking segments f 11-f 19 in fig. 1 as an example: from F11 to F18, the formation magnetic susceptibility gradually increases, and the overall low magnetic susceptibility comprises a small amount of medium-high magnetic susceptibility (HT 1, HT 2) and gradually transits to the overall medium-high magnetic susceptibility (HT 3) which comprises a large range of high magnetic susceptibility (HT 4), so that the magnetic change characteristic of transition from the edge of the suture belt to the center is reflected; the extremely high magnetic susceptibility (HT 6) between F18 and F19 shows root ferromagnetic bodies on the section, reflects the existence of superbasic rock blocks such as olive rock, snake green rock and the like, and proves the development history of the middle ridge of the suture belt; the high magnetic susceptibility and the low magnetic susceptibility are distributed alternately, and the medium-high magnetic susceptibility is presented by a root-free block or a floating body (HT 1, HT2 and HT 4), so that the characteristics of a mixed rock zone formed by repeatedly kneading seam-zone stratum rock bodies through a multi-stage structure are reflected; the response of the susceptibility contours to fracture is relatively pronounced in the central region of the seam tape, predominantly as a through-type gradient tape (F12, F18, F19) with high and low susceptibility transitions, while there is no pronounced response in the edge regions of the seam tape (F11).
1.3 magma arc unit
In terms of space, the magma arc unit is mainly composed of weak magnetic bodies and strong magnetic bodies, the electrodeless strong magnetic bodies are arranged, the strong magnetic bodies and the weak magnetic bodies are distributed alternately, and the strong magnetic bodies are usually represented by rooted blocks and rooted floating bodies. Taking the susceptibility contour line of fig. 3 as an example: the medium-high magnetic susceptibility (HT 8, HT 9) is in a shape of a root cluster, the (HT 7) is in a shape of a cap with a root, and the low magnetic susceptibility is cut and surrounded by the medium-high magnetic susceptibility block, so that the spatial distribution characteristics of an invaded body, an invaded stratum, different phases of invaded bodies and different basic rock bodies are reflected; the medium-high magnetic susceptibility represents neutral flash rock formed by late invasion or contemporaneous invasion, the low magnetic susceptibility represents meta-acid granite formed by early invasion, and the non-magnetism of the tunnel inlet and outlet sections represents early deposition and deterioration stratum; the magnetic susceptibility contour map is uniform and flat as a whole, and does not react exactly to the position of the fracture structure (F39, F40), but reacts well to the boundary fault (F41) between the invaded body and the sedimentary rock.
2. Mesoscale analysis and interpretation structural band and formation lithology rock group
The aeromagnetic has no obvious reaction or inaccurate reaction position on most of the breaks of the land block, the seam tape edge area and the magma arc area, and the aeromagnetic has obvious reaction on the area break area and the water-rich fault, so that the structural units and the lithology large area which are analyzed in a large scale can be combined, and the structural tape and the stratum lithology rock group can be analyzed by utilizing the avionic three-dimensional inversion resistivity result to realize the mesoscale analysis.
2.1 fracture Structure
The rule of the fracture structure on the aviation electromagnetic three-dimensional inversion resistivity section chart is as follows: inter-layer faults (i.e., faults that develop between two sets of formations or lithologies) are primarily represented by resistivity gradient bands (abrupt bands) that communicate from the earth's surface to deep, such as F11, F17, F18, F27, F28, F29, F30, F40; part of interlayer faults are simultaneously represented as connected resistivity and magnetic susceptibility gradient bands (abrupt bands), such as F17, F18 and F41; the in-situ faults (i.e. faults developing in the same set of strata) are mainly represented by low stop bands from the earth surface to deep through, such as F10, F19 and F39.
2.2 lithology rock group
The tunnel aviation electromagnetic three-dimensional inversion resistivity profile of the plot section (figure 5) shows a form that high resistance and low resistance are alternately arranged in an obvious vertical rate gradient zone, the position of the gradient zone is the development position of a fracture structure, and according to 2.1 analysis, the fracture zone is a lithology rock group boundary. Avionics high resistance zone lithology is limestone (HE 4), quartz sandstone (HE 5), etc. Avionics low resistance zone lithology is carbonaceous shale (LE 9), slate, mudstone (LE 7), etc.
The tunnel electromagnetic three-dimensional inversion resistivity profile of the seam tape section (fig. 4) shows that the high resistance and the low resistance are distributed in a form of being alternately cut by regular dividing lines and mutually wrapped, and the gradient tape is coarser relative to the land block section. The bulk high resistance zone is characterized by granite and flash rock (HE 1) invading along the edge zone of the joint strip. The high-resistance and low-resistance interphase region distribution lithology is a mixed rock zone which mainly contains various schists and is mixed with superbasic rock blocks and quartz rocks in the joint zone, the lithology of the high-resistance block is generally quartz rocks, quartz schists, gneiss (HE 2 and HE 3) and the like, and the lithology of the low-resistance block is generally mica schists, green mud schists, superbasic rocks and the like (LE 4 and LE 5).
The tunnel aviation electromagnetic three-dimensional inversion resistivity profile of the magma arc section (figure 6) shows a distribution form that high resistance and low resistance are mutually wrapped by smooth fuzzy boundary lines, and through-type gradient bands are rarely arranged. Comparing and analyzing the figures 3 and 6, the resistivity and lithology subareas are not necessarily related, wherein the high magnetic resistance, the low magnetic resistance, the high magnetic resistance, the low magnetic resistance and the low magnetic resistance are all corresponding blocks. The low resistance zone is typically an altered zone (LE 11), a fault zone (LE 12) or a joint dense zone (LE 13) that invades the rock.
3. Small scale detailed interpretation of specific strata, joint dense zones or broken water-rich zones, karst development zones
Based on the large-scale and medium-scale analysis, various low-resistance forms are analyzed in detail by utilizing an aviation electromagnetic three-dimensional inversion resistivity section chart, and rock physical characteristics related to engineering geology are interpreted to realize small-scale analysis.
3.1 Special formations
The low-resistance area caused by the special stratum is continuously and directionally spread in a large range, the shape is irregular, but the spreading trend and trend of the low-resistance area can be seen, as shown by LE7, LE9 and LE16 in figure 5, the physical characteristics of dwarfism Gong Tang groups of carbonaceous slate, chalk-system doney carbon shale and dwarfism marl mudstone are respectively reflected.
3.2 dense joint zones or broken water-rich zones
The low-resistance areas caused by the joint dense bands and the water-rich bands are distributed in a small range in a lens shape or a bead shape, and appear in a low-resistance ring in the whole high-resistance area or an ultra-low-resistance ring in the whole low-resistance area, so that the spreading trend and the trend of the low-resistance rings are difficult to distinguish, such as LE6 and LE8 in FIG. 5, which respectively reflect the broken water-rich areas in the dwarf series La Gong Tang group sandstone and the near series Zong Bai Qun sandstone mudjacking.
3.3 karst development zone
According to the analysis of 2.2, the karst is characterized by overall high resistance, karst cave, solution gap and other karst phenomena develop in the karst region, and the resulting low resistance is generally distributed in the overall high resistance block in the form of a bead, and has certain vertical banding characteristics, and has positive correlation with the development of plane of the patterns in the region, such as LE17 in FIG. 5.
Taking a mountain area traffic corridor of Qinghai-Tibet plateau as an example, the corridor spans the seam area of the Yangjiang, the Athletic and the like, passes through the magma arc area of the Bersla and the like, and passes through the land areas of the Changdu, the Gangzusi and the like. And extracting a 370 km traffic corridor aviation electromagnetic and magnetic three-dimensional inversion result diagram, performing multi-scale interpretation analysis according to the method, rapidly and efficiently analyzing tunnel engineering geological problems, and providing a basis for line scheme optimization and tunnel geological condition evaluation. The steps for solving the problems are as follows:
1. through large-scale analysis and interpretation, the full-line 16-seat long deep buried tunnel engineering is divided into different geological structure units, stratum large areas where each tunnel passes through and main geological problems possibly encountered are analyzed, suggestions are provided for optimizing a line scheme, and the tunnel engineering scheme is guided to be optimized.
2. Aiming at different physical characteristics of different geological structure units, three-dimensional inversion magnetic susceptibility and resistivity achievements of each tunnel section are comprehensively analyzed, interlayer fracture is calibrated by vertical communication type resistivity or magnetic susceptibility gradient bands (abrupt bands), fracture is calibrated in a through type low-stop band calibration layer, and stratum lithology names, alteration bands in invaded rocks and the like are calibrated in sequence according to the combination characteristics of high-low-resistance areas (blocks) in each geological structure unit. And a reference is provided for tunnel surrounding rock division and engineering geological problem analysis.
3. The low-resistance areas (blocks) in the tunnel section resistivity result chart are analyzed in detail, the low-resistance areas which are continuously and directionally spread in a larger range are marked as special strata such as carbonaceous shale, mudstone and the like, the low-resistance areas which are distributed in a lens shape or a bead shape in a smaller range are marked as joint dense bands or broken water-rich bands, and the low-resistance areas (points) which are distributed in the high-resistance areas in the bead shape in the soluble rock areas are marked as karst development areas. On the basis, surrounding rock grades are divided in detail, and concrete engineering geological problems such as rock burst, soft rock large deformation, gushing water, harmful gas and the like possibly encountered by a tunnel are analyzed by combining the tunnel burial depth, so that basis is provided for engineering scheme design and engineering construction measures.
The foregoing description is only of the preferred embodiments of the present invention, and is not intended to limit the scope of the present invention.
Claims (1)
1. The multi-scale analysis method for the aviation geophysical prospecting three-dimensional inversion result applied to the traffic corridor is characterized by comprising the following steps of: the method utilizes the three-dimensional inversion result of the aerospace geophysical prospecting, carries out large-scale analysis and interpretation aiming at a geological structure unit, carries out mesoscale analysis and interpretation aiming at the lithology of a specific structural zone and a stratum, carries out small-scale analysis and interpretation aiming at engineering geological features of a special stratum, a joint intensive zone, a broken water-rich zone and a karst development zone, and obtains basic geological conditions of a tunnel engineering zone through multi-scale comprehensive interpretation analysis of large, medium and small;
the method comprises the following specific steps:
1) Extracting a typical section three-dimensional inversion magnetic susceptibility result graph of the traffic corridor, and performing macroscopic interpretation and large-scale analysis according to the integral size of the magnetic susceptibility and the strong and weak magnetic body distribution characteristic interpretation to divide a geological structure unit of the corridor and a stratum lithology large area;
the specific method of the step 1) is as follows: the magnetic susceptibility result map is inverted three-dimensionally by utilizing a typical section of a traffic corridor, and land block construction units and normal sedimentary rock and negative metamorphic rock large areas are defined according to the characteristic of 'overall nonmagnetic-weak magnetism'; according to the fact that the strong magnetic bodies and the weak magnetic bodies are distributed alternately in the whole, the strong magnetic bodies in the central area are represented by root lumps and comprise extremely strong magnetic bodies, and the strong magnetic bodies in the edge area are represented by root-free lumps or floating bodies, so that the seam belt construction units and the mixed rock large area are defined; according to the fact that the strong magnetic bodies and the weak magnetic bodies are distributed alternately, the strong magnetic bodies are presented in a rooted block, and the characteristic of the strong magnetic bodies is not contained, so that a magma arc construction unit and an intrusion rock large area are defined;
2) Extracting a tunnel longitudinal section three-dimensional inversion resistivity result graph of a corresponding region, comparing the three-dimensional susceptibility result graph, calibrating a fracture structure according to a vertical connected resistivity or susceptibility gradient zone and a through type low stop band, dividing lithology rock groups in the tunnel section according to distribution characteristics of high and low resistance regions, and carrying out refinement interpretation and mesoscale analysis;
the specific method of the step 2) is as follows: comparing the tunnel longitudinal section three-dimensional inversion resistivity result graph with the longitudinal section three-dimensional inversion susceptibility result graph which is analyzed and interpreted in a large scale by the step 1), and calibrating interlayer fracture according to the characteristics of vertical connected resistivity or susceptibility gradient bands; calibrating the internal faults according to the characteristics of 'through type low stop band from the earth surface to the deep'; dividing lithology rock groups of the land block construction units according to the characteristic that the high-resistance area and the low-resistance area are separated by an obvious vertical communicated gradient zone; dividing lithology rock groups of the suture belt construction units according to the characteristic that the high-resistance area and the low-resistance area are mutually cut by irregular boundaries, and the resistivity lump and the magnetic susceptibility lump approximately correspond to each other; dividing lithology rock groups of the magma arc construction units according to the characteristic that the high-resistance area and the low-resistance area are mutually wrapped by smooth boundary lines and the resistivity block and the magnetic susceptibility block have no corresponding rule;
3) On the basis of large-scale and medium-scale analysis and interpretation, further refining and analyzing a non-vertical connected low-resistance region in a three-dimensional resistivity result graph, dividing a special stratum, a joint dense zone, a broken water-rich zone and a karst development zone by combining a low-resistance form and a lithology rock group, and carrying out fine interpretation and small-scale analysis to obtain basic geological conditions of a tunnel engineering region;
the specific method of the step 3) is as follows: the tunnel longitudinal section three-dimensional inversion resistivity result graph which is analyzed through the mesoscale analysis is completed through the steps 1) and 2) is utilized to analyze the 'non-vertical connected' low-resistance form in each lithology rock group in detail, the sedimentary rock area and the negative metamorphic rock area are continuously and directionally spread in a large range, and the low-resistance area which is irregular in shape but can see the spreading trend and trend is translated into special stratum of carbonaceous shale, carbonaceous slate and mudstone; the lens-shaped or beaded distribution of the intrusion rock area, the metamorphic hard rock area and the deposition harder rock area in a smaller range appears as a low-resistance ring in the whole high-resistance area or an ultra-low-resistance ring in the whole low-resistance area, and the low-resistance area of which the spreading trend and the trend are difficult to distinguish is interpreted as a joint dense belt and a broken water rich belt; the low-resistance region with the dissoluble rock region being distributed in a bead shape in the whole high-resistance block, and having no spreading trend but a certain vertical zonal characteristic is interpreted as a karst development region.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010595709.6A CN112114370B (en) | 2020-06-28 | 2020-06-28 | Multi-scale analysis method for applying aviation geophysical prospecting three-dimensional inversion result to traffic corridor |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010595709.6A CN112114370B (en) | 2020-06-28 | 2020-06-28 | Multi-scale analysis method for applying aviation geophysical prospecting three-dimensional inversion result to traffic corridor |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112114370A CN112114370A (en) | 2020-12-22 |
| CN112114370B true CN112114370B (en) | 2023-10-10 |
Family
ID=73799050
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010595709.6A Active CN112114370B (en) | 2020-06-28 | 2020-06-28 | Multi-scale analysis method for applying aviation geophysical prospecting three-dimensional inversion result to traffic corridor |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112114370B (en) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113609645B (en) * | 2021-07-07 | 2024-01-26 | 杨建云 | Gold mine exploration target area prediction method |
| CN115618702B (en) * | 2022-10-21 | 2023-05-09 | 河海大学 | Method for generating two-dimensional ultrahigh volume fraction mixed rock numerical simulation model by adopting overlapping contact cutting algorithm |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR20100018811A (en) * | 2008-08-07 | 2010-02-18 | 에스케이건설 주식회사 | Construction simulation system and method for tunnel and underground storage cavern considering geological condition |
| CN108345049A (en) * | 2018-02-12 | 2018-07-31 | 山东大学 | Underground engineering unfavorable geology detects multi-method constraint inverting and joint interpretation method |
| CN110231659A (en) * | 2019-06-03 | 2019-09-13 | 中铁二院工程集团有限责任公司 | Buried depth changes railway tunnel aviation electromagnetic detection method greatly |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CA2802722C (en) * | 2010-07-27 | 2023-04-04 | Exxonmobil Upstream Research Company | Inverting geophysical data for geological parameters or lithology |
-
2020
- 2020-06-28 CN CN202010595709.6A patent/CN112114370B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR20100018811A (en) * | 2008-08-07 | 2010-02-18 | 에스케이건설 주식회사 | Construction simulation system and method for tunnel and underground storage cavern considering geological condition |
| CN108345049A (en) * | 2018-02-12 | 2018-07-31 | 山东大学 | Underground engineering unfavorable geology detects multi-method constraint inverting and joint interpretation method |
| CN110231659A (en) * | 2019-06-03 | 2019-09-13 | 中铁二院工程集团有限责任公司 | Buried depth changes railway tunnel aviation electromagnetic detection method greatly |
Non-Patent Citations (3)
| Title |
|---|
| 基于CSAMT和3D重磁地球物理模型综合解译相山石洞地区深部地质结构;王峰 等;地球物理学进展;第31卷(第6期);2657-2663 * |
| 大井坡航空地球物理试验场典型剖面综合地球物理解释;张婉 等;吉林大学学报(地球科学版);第46卷(第2期);569-580 * |
| 航空和地面综合地球物理方法开展干热岩探测研究;赵丛 等;矿产勘查;第9卷(第5期);938-946 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112114370A (en) | 2020-12-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Deng et al. | Two distinct strike-slip fault networks in the Shunbei area and its surroundings, Tarim Basin: Hydrocarbon accumulation, distribution, and controlling factors | |
| CN103454678B (en) | Method and system for determining seismic slice isochronism | |
| CN109324355B (en) | Pegmatite type rare metal prospecting method | |
| CN112114370B (en) | Multi-scale analysis method for applying aviation geophysical prospecting three-dimensional inversion result to traffic corridor | |
| CN105717551A (en) | Blind ore space positioning method for volcanic rock type uranium mine | |
| CN113608278B (en) | Sandstone-type uranium ore positioning method in red heterolayer of sedimentary basin | |
| Behlau et al. | Geological and tectonic investigations in the former Morsleben salt mine (Germany) as a basis for the safety assessment of a radioactive waste repository | |
| Basson et al. | Structural analysis and 3D modelling of major mineralizing structures at the Phalaborwa copper deposit | |
| CN112711078A (en) | Method for identifying favorable sandstone type uranium mineralization sand body in deep sedimentary basin | |
| CN115128698B (en) | Method for determining depth-of-extension pattern of rock-magma hot-liquid type multi-metal ore field or ore deposit rock-control and ore-control structure | |
| Carlos et al. | Imaging iron ore from the Quadrilátero Ferrífero (Brazil) using geophysical inversion and drill hole data | |
| Benham et al. | The function of baffles within heavy oil reservoir and their impact on field development: A case study from kuwait | |
| CN111679340A (en) | Ore seeking method for black rock series vanadium-molybdenum ore | |
| CN113156531B (en) | Hidden manganese ore bed exploration method | |
| Harris et al. | Volcanic setting of the Bajo de la Alumbrera porphyry Cu-Au deposit, Farallón Negro volcanics, northwest Argentina | |
| Srivastava et al. | Detailed gravity studies over Jabera—Damoh region of the Vindhyan basin (Central India) and crustal evolution | |
| de Macêdo Filho et al. | Unraveling the tectonic setting and crystallization history of the Equatorial Atlantic Magmatic Province | |
| CN115577213B (en) | Water invasion direction prediction method | |
| CN114594518B (en) | Fine stratum contrast method for complex fault blocks in later development period based on well-seismic alternation | |
| Wilson | Geological appraisal of the Mootwingee area, New South Wales | |
| CN108732639B (en) | Shale stratum isochronous comparison method and system | |
| Kalicki | Krzysztof Żurek | |
| Ahmed et al. | Seismic facies modelling of potwar basin using seismic and well log data | |
| CN117310826A (en) | Geophysical method for delineating favorable target area of ore formation and prospecting | |
| Warrington | Inside Mid-Triassic fluvial deposits: a legacy of mining sediment-hosted ore bodies in Cheshire, UK |
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 |