Longitudinal section stratum boundary line identification method based on aeronautical geophysical prospecting magnetic susceptibility result map

Abstract

The invention discloses a vertical section stratum boundary line identification method based on an aeronautical geophysical prospecting susceptibility result map, which can quickly and efficiently delineate a vertical section stratum boundary line (including stratum lithology division and fault / fracture zone boundary division) by utilizing an aeronautical geophysical prospecting susceptibility result in a tunnel engineering geological survey stage in a complex and hard mountainous area, improve the survey efficiency and reduce the construction cost. The method has the advantages of simple process, low cost and low exploration cost investment, is especially suitable for detecting target bodies with large magnetic differences or geologic bodies with components of minerals (especially biotite, hornblende, magnetite and the like) in rocks changed due to different metamorphic degrees, and can be used for detecting the content of the rock according to the magnetic susceptibility value and the abnormal distribution form and by comprehensively considering factors such as stratum lithology and drilling data. A vertical section stratum lithology boundary or a fault fracture zone boundary can be specifically divided; the method is especially applied to tunnel areas with large differences between different stratum lithology and magnetic susceptibility, and has strong identification capability.

Description

Technical Field

[0001] The invention relates to the technical field of stratum boundary identification, and in particular to a longitudinal section stratum boundary identification method based on an aerial geophysical exploration magnetic susceptibility result map. Background Art

[0002] With the continuous expansion of infrastructure projects in my country, more and more engineering construction projects involve difficult and complex mountainous areas. These areas have undulating terrain, complex geological conditions, and inconvenient transportation. Traditional engineering geological survey methods in these areas often face great challenges. Using airborne geophysical methods to carry out survey work in difficult and complex mountainous areas has the advantages of fast speed, high efficiency, and no influence from natural conditions on the ground; therefore, using airborne geophysical magnetic susceptibility results to identify longitudinal stratigraphic boundaries can effectively improve survey efficiency.

[0003] Traditional methods for identifying stratigraphic boundaries mostly rely on ground artificial geological surveys, drilling sampling, indoor tests, ground geophysical exploration and other means, which have problems such as high cost, long cycle, high personnel input, and limited accuracy. In recent years, with the continuous development of airborne geophysical exploration technology, its inversion results maps have been increasingly used in geological exploration. However, the effective application of airborne geophysical exploration magnetic susceptibility results to the identification of stratigraphic boundaries in geological longitudinal sections is still a blank in the current field of geological exploration. For a long time, the results of airborne geophysical exploration magnetic susceptibility at home and abroad have been mainly used in regional geological surveys, mineral resource exploration, marine geological exploration and other fields. At present, the analysis of airborne geophysical data in China is generally qualitative analysis, and relevant research has been carried out on deep underground unfavorable geology, but there is currently no precedent for the interpretation of stratigraphic boundaries and fault / fracture zone locations based on the distribution of airborne geophysical susceptibility. Summary of the invention

[0004] The present application provides a longitudinal section stratigraphic boundary identification method and system based on the aerial geophysical magnetic susceptibility results map, so as to solve the problems of high cost, long cycle, high manpower input, limited accuracy, etc. in the traditional stratigraphic boundary identification method.

[0005] According to the first aspect, an embodiment provides a method for identifying stratum boundaries in a longitudinal section based on an aerial geophysical magnetic susceptibility result map, the method comprising:

[0006] Collect regional geological data, analyze the distribution characteristics of lithology in the region, and obtain the range of aeromagnetic magnetic susceptibility parameters of various strata and lithologies in the region;

[0007] Generate longitudinal section magnetic susceptibility results map based on aeromagnetic data, and preliminarily determine the boundaries of different stratum lithology zones and the approximate locations of fault / fracture zones by analyzing the longitudinal section magnetic susceptibility results map and combining the magnetic susceptibility parameter value ranges of various stratum lithologies;

[0008] Verify and adjust the preliminarily determined boundaries of different stratigraphic lithology zones and the approximate locations of fault / fracture zones through field investigations;

[0009] Drilling work is carried out on the boundaries of different strata and lithology zones and the locations of faults / fracture zones to further identify strata boundaries, fault / fracture zone locations and the overall trend and location of deep underground;

[0010] Through comprehensive analysis, the stratigraphic boundaries and fault / fracture zone boundaries of the regional longitudinal section are accurately divided.

[0011] Furthermore, we collected regional geological data and analyzed the lithology distribution characteristics in the region, including:

[0012] Collect existing regional geological data, including regional geological maps, regional geological survey reports, nearby drilling data, and nearby engineering geological data, to obtain the stratigraphic lithology distribution types and rock characteristics in the region, including the content of the main components of rock minerals.

[0013] Furthermore, the range of magnetic susceptibility parameters of various strata and lithologies in the region were statistically obtained, including:

[0014] According to the lithology of the strata on the regional geological map, the sections are roughly divided, the magnetic survey data of each lithology section collected by airborne geophysical means are analyzed, and the range of magnetic susceptibility parameter values ​​of various types of rocks inverted by airborne magnetic survey data are statistically analyzed.

[0015] Furthermore, a longitudinal section magnetic susceptibility result map is generated based on the aeromagnetic data, specifically including:

[0016] First, the collected aeromagnetic data are corrected, including normal geomagnetic field correction and magnetic diurnal variation correction; then the corrected aeromagnetic data are appropriately leveled; secondly, aeromagnetic potential field conversion processing is carried out, mainly including magnetic field poles, analytical signals, and oblique derivative calculations; then the three-dimensional inversion of the magnetization intensity vector is performed, and the three-dimensional MVI apparent magnetic susceptibility data formed by the inversion is mapped in three dimensions through mapping software, and finally the line is cut along the centerline to extract the contour map of the magnetic susceptibility of the longitudinal section of the line.

[0017] Furthermore, by analyzing the longitudinal section magnetic susceptibility results map and combining the magnetic susceptibility parameter value range of various strata and lithologies, the boundaries of different strata and lithology zones and the approximate locations of faults / fracture zones are preliminarily determined, including:

[0018] The magnitude of the longitudinal section magnetic susceptibility value and the contour morphology characteristics reflect the distribution range and stratigraphic boundaries of the regional strata. The specific identification methods include: based on the statistical range of the aeromagnetic magnetic susceptibility parameter values ​​of various strata lithologies, roughly judging the lithological divisions of different strata according to the difference in longitudinal section magnetic susceptibility values; combining the longitudinal section magnetic susceptibility contour anomalies, including the vertical gradient zone, to reflect the lithological contact boundaries of the area;

[0019] The morphological characteristics of the longitudinal section magnetic susceptibility contour lines reflect the location and boundary of the fault rupture zone. The specific identification methods include: combining the longitudinal section magnetic susceptibility contour line anomalies, including the existence of vertical strip-like characteristics with lower magnetic susceptibility on both sides, which reflects the location and boundary of the fault rupture zone; the axial direction of the magnetic anomaly can reflect the inclination of the fault / rupture zone.

[0020] Furthermore, the preliminary determined lithology division boundaries and approximate locations of fault / fracture zones were verified and adjusted through field investigations, including:

[0021] Tracing the stratigraphic boundaries, determining stratigraphic boundaries through field investigation of rock sequence, lithology changes, and structural evidence;

[0022] Survey the outcrop of the fault surface, conduct a comprehensive analysis through geomorphic signs, lithology changes, rock formation dip changes, and material composition within the fault zone to determine the location and boundaries of the fault, and conduct excavation work for verification when necessary;

[0023] The field geological survey route carries out field verification and supplementary investigation of continuous geological sections along and perpendicular to the geological body boundaries.

[0024] Furthermore, drilling work is carried out on the boundaries of different strata and lithology zones and the locations of faults / fracture zones to further identify the strata boundaries, fault / fracture zone locations and the overall direction and location of the deep underground, including:

[0025] By determining the boundaries of the stratigraphic lithology zoning and the location of the fault / fracture zone and combining the fault dip, drilling holes are arranged on both sides of the boundary line for verification, the stratigraphic lithology and fault zone characteristics on both sides of the boundary line are identified, and the lithology differences on both sides of the boundary line are verified.

[0026] Furthermore, drilling holes were arranged on both sides of the boundary line for verification, including:

[0027] The principles of drilling arrangement are as follows: for the stratigraphic lithology boundary, verification holes are arranged at preset distances on both sides of the intersection of the line centerline and the stratigraphic boundary with the preliminarily determined stratigraphic boundary as the center line; for the fault position and boundary, verification holes are arranged in the vertical strip of lower magnetic susceptibility where the preliminarily determined fault boundary intersects with the line centerline to find out the stratigraphic characteristics within the fault zone; in addition, verification holes are arranged at preset distances outside the two preliminarily determined fault boundaries to find out the difference in lithology inside and outside the fault zone and verify the accuracy of the boundary.

[0028] Furthermore, through comprehensive analysis, the stratigraphic boundaries and fault / fracture zone boundaries of the regional longitudinal section are accurately divided, including:

[0029] By integrating regional geological data, aerial geophysical magnetic susceptibility maps, field survey results, and drilling results, the regional longitudinal section stratigraphic boundaries and fault / fracture zone boundaries are accurately divided.

[0030] According to the second aspect, in one embodiment, a longitudinal section stratigraphic boundary identification system based on an aerial geophysical magnetic susceptibility result map is provided, the system comprising:

[0031] The regional geological data analysis module is used to collect regional geological data, analyze the lithology distribution characteristics in the region, and obtain the range of aeromagnetic magnetic susceptibility parameter values ​​of various strata and lithologies in the region;

[0032] The airborne geophysical prospecting results interpretation module is used to generate longitudinal section magnetic susceptibility results maps based on airborne magnetic survey data. By analyzing the longitudinal section magnetic susceptibility results maps and combining the magnetic susceptibility parameter value ranges of various strata and lithologies, the boundaries of different strata and lithology divisions and the approximate locations of faults / fracture zones can be preliminarily determined;

[0033] Field investigation and verification module, used to verify and adjust the preliminarily determined boundaries of different strata lithology zones and the approximate locations of fault / fracture zones through field investigation;

[0034] The drilling verification module is used to carry out drilling work on the boundaries of different strata lithology zones and the locations of faults / fracture zones, and further determine the strata boundaries, fault / fracture zone locations and the overall direction and location of the deep underground;

[0035] The comprehensive analysis module is used to accurately divide the regional longitudinal section stratigraphic boundaries and fault / fracture zone boundaries through comprehensive analysis.

[0036] The present application provides a method for identifying the longitudinal stratigraphic boundary based on the magnetic susceptibility result map of airborne geophysical prospecting, which is particularly suitable for detecting target bodies with large magnetic differences, or geological bodies whose mineral components (especially biotite, hornblende, magnetite, etc.) contained in rocks with different degrees of metamorphism have changed. According to the magnitude of the magnetic susceptibility value and the abnormal distribution form, and taking into account factors such as stratigraphic lithology and drilling data, the longitudinal stratigraphic lithology boundary or the fault fracture zone boundary can be specifically divided; it is especially used in tunnel areas with large differences in stratigraphic lithology and magnetic susceptibility, and has a strong identification ability. Compared with the prior art, the present invention has the following advantages and effects:

[0037] (1) Highly practical. The results of the aerogeophysical magnetic susceptibility can be used to specifically delineate the lithology boundaries of strata or the boundaries of fault fracture zones in geological surveys in complex and dangerous mountainous areas.

[0038] (2) High survey efficiency and cost savings. The use of airborne geophysical methods can shorten the geological survey period and quickly obtain the lithology and boundary distribution of the tunnel area. Compared with traditional survey methods, it reduces the costs of drilling, geophysical exploration, and field geological surveys.

[0039] (3) It solved the problem of having no survey data to refer to in long sections under complex terrain conditions, ensured reasonable engineering settings, avoided geological risks, and improved the quality of survey and design. BRIEF DESCRIPTION OF THE DRAWINGS

[0040] Figure 1 A flow chart of a method for identifying stratum boundaries in a longitudinal section based on an aerial geophysical magnetic susceptibility result map provided by one embodiment of the present invention;

[0041] Figure 2 A specific implementation flow chart of a method for identifying longitudinal stratigraphic boundaries based on an aerial geophysical magnetic susceptibility result map provided by one embodiment of the present invention;

[0042] Figure 3 A longitudinal section contour map of the magnetic susceptibility of a tunnel obtained by three-dimensional inversion of aeromagnetic MVI in a longitudinal section stratigraphic boundary identification method based on an aerial geophysical magnetic susceptibility result map provided by an embodiment of the present invention;

[0043] Figure 4 The drilling results of ZK1 borehole in a longitudinal section stratigraphic boundary identification method based on an aerogeophysical magnetic susceptibility result map provided by one embodiment of the present invention;

[0044] Figure 5 The drilling results of ZK3 borehole in a longitudinal section stratigraphic boundary identification method based on an aerogeophysical magnetic susceptibility result map provided by one embodiment of the present invention;

[0045] Figure 6 A longitudinal section stratigraphic boundary identification method based on an aerial geophysical magnetic susceptibility result map provided in one embodiment of the present invention is a stratigraphic lithology boundary and fault / fracture zone zoning map of a tunnel longitudinal section;

[0046] Figure 7 A schematic diagram of the logical structure of a longitudinal section stratigraphic boundary identification system based on an aerial geophysical magnetic susceptibility result map is provided as an embodiment of the present invention. DETAILED DESCRIPTION

[0047] The present invention is further described in detail below by specific embodiments in conjunction with the accompanying drawings. Wherein similar elements in different embodiments adopt associated similar element numbers. In the following embodiments, many detailed descriptions are for making the present application better understood. However, those skilled in the art can easily recognize that some features can be omitted in different situations, or can be replaced by other elements, materials, methods. In some cases, some operations related to the present application are not shown or described in the specification, this is to avoid the core part of the present application being overwhelmed by too much description, and for those skilled in the art, it is not necessary to describe these related operations in detail, and they can fully understand the related operations according to the description in the specification and the general technical knowledge in the art.

[0048] In addition, the features, operations or characteristics described in the specification can be combined in any appropriate manner to form various implementations. At the same time, the steps or actions in the method description can also be interchanged or adjusted in a manner that is obvious to those skilled in the art. Therefore, the various sequences in the specification and the drawings are only for the purpose of clearly describing a certain embodiment and are not meant to be a required sequence, unless otherwise specified that a certain sequence must be followed.

[0049] The present invention provides a longitudinal section stratigraphic boundary identification method based on the aerial geophysical magnetic susceptibility result map, which solves the problems of low efficiency, lack of data, and poor accuracy of traditional means of deep underground stratigraphic boundary division. Traditional stratigraphic boundary identification methods mostly rely on ground artificial geological survey and mapping, drilling sampling, ground geophysical prospecting and other means. Especially in difficult and complex mountainous areas, geological surveys often face terrain problems, and personnel surveys, drilling, ground geophysical prospecting and other means are difficult to implement, and it is difficult to effectively obtain geological data. The method of the present invention is suitable for the rapid division of longitudinal section stratigraphic boundaries and fault boundaries in the survey and design stage of deep buried long tunnels in complex and dangerous mountainous areas, providing a basis for the preparation of engineering geological survey reports.

[0050] The first embodiment of the present invention provides a method for identifying the longitudinal stratigraphic boundary based on the aerial geophysical magnetic susceptibility results map, which can be used to quickly and efficiently delineate the longitudinal stratigraphic boundary (including stratigraphic lithology division and fracture boundary division) in the geological survey stage of tunnel engineering in complex and dangerous mountainous areas by using the aerial geophysical magnetic susceptibility results, thereby improving the survey efficiency and reducing the survey cost investment. Figure 1 and Figure 2 Provide detailed explanation.

[0051] like Figure 1 As shown, in step S100, regional geological data are collected, the lithology distribution characteristics in the region are analyzed, and the range of aeromagnetic magnetic susceptibility parameter values ​​of various strata and lithologies in the region are statistically obtained.

[0052] The above steps specifically include:

[0053] S110: Collect existing regional geological data in the tunnel area, including regional geological maps, regional geological survey reports, nearby drilling data, nearby engineering geological data, etc., to obtain the lithology distribution types and rock characteristics of the strata along the tunnel, including the content of the main components of rock minerals (especially biotite, hornblende, magnetite, etc.);

[0054] S120: The magnetic susceptibility is mainly determined by the content of ferromagnetic minerals, especially biotite, hornblende, magnetite, etc. The stratigraphic lithology on the regional geological map is roughly divided into sections. The magnetic survey data of each stratigraphic lithology section collected by airborne geophysical exploration are analyzed, and the range of magnetic susceptibility parameter values ​​inverted from the airborne magnetic survey data of various types of rocks are statistically analyzed.

[0055] The original rock in the study area of ​​a difficult and complex mountain tunnel is mainly intrusive rock. Due to different degrees of metamorphism, the mineral composition changes, especially the content of biotite, amphibole, magnetite and other minerals varies greatly. Therefore, the lithology of the formation can be identified from the distribution characteristics of magnetic susceptibility. The lithology of the tunnel area mainly includes gneissic granodiorite (strong magnetism), granitic gneiss (weak magnetism-non-magnetism), Cretaceous marble intercalated calcareous sandstone (weak magnetism-non-magnetism), fault fracture zone, etc.

[0056] like Figure 1 As shown, in step S200, a longitudinal section magnetic susceptibility result map is generated based on the aeromagnetic data, and the boundaries of different stratigraphic lithology zones and the approximate locations of faults / fracture zones are preliminarily determined by analyzing the longitudinal section magnetic susceptibility result map and combining the magnetic susceptibility parameter value ranges of various stratigraphic lithologies.

[0057] The above steps specifically include:

[0058] S210: First, the collected aeromagnetic data is corrected, including normal geomagnetic field correction, magnetic diurnal variation correction, etc.; then the corrected aeromagnetic data is appropriately leveled; secondly, aeromagnetic potential field conversion processing is carried out, mainly including magnetic field pole, analytical signal, oblique derivative calculation, etc.; then the three-dimensional inversion of the magnetization intensity vector is carried out, and the three-dimensional MVI apparent magnetic susceptibility data formed by the inversion is mapped in three dimensions through mapping software, and finally the tunnel is cut along the centerline to extract the magnetic susceptibility contour map of the tunnel longitudinal section, such as Figure 3 As shown;

[0059] S220: The magnitude of the magnetic susceptibility value and the contour morphology of the tunnel longitudinal section can reflect the distribution range and boundary of the strata in the tunnel area. The specific identification methods are as follows: (1) Based on the range of magnetic susceptibility parameter values ​​of various strata and lithologies statistically calculated in S120, the strata and lithologies of different strata can be roughly determined according to the difference in magnetic susceptibility values ​​of the tunnel longitudinal section; (2) Combined with the abnormal magnetic susceptibility contour lines of the tunnel longitudinal section, such as the vertical gradient zone, it reflects the lithology contact boundary of this area;

[0060] S230: The morphological characteristics of the magnetic susceptibility contour lines of the tunnel longitudinal section can reflect the location and boundary of the fault fracture zone. The specific identification method is as follows: Combined with the anomaly of the magnetic susceptibility contour lines of the tunnel longitudinal section, for example, the vertical strip-like characteristics of lower magnetic susceptibility on both sides reflect the location and boundary of the fault fracture zone; the axial direction of the magnetic anomaly can generally reflect the inclination of the fault / fracture zone;

[0061] S240: According to steps S210, S220, and S230, the positions of stratum lithology boundaries and fault / fracture zone boundaries are delineated for the magnetic susceptibility distribution anomaly.

[0062] By processing and calculating the aeromagnetic data, the three-dimensional inversion of the magnetization intensity vector is performed. The inversion results are processed to form a three-dimensional magnetic susceptibility result map, and the magnetic susceptibility contour map of the tunnel longitudinal section is extracted based on this. The boundaries of the stratum lithology divisions and the approximate locations of the fault / fracture zones are delineated according to the differences in magnetic susceptibility distribution, as shown in the figure:

[0063] The magnetic susceptibility of the tunnel varies greatly. In the B-B1 section of the tunnel, a large area of ​​high-strength magnetic susceptibility of 500-3000 (10-5SI) is distributed, mostly in agglomerate, discontinuous, and floating. It is preliminarily judged to be mainly ultrabasic rock or intrusive rock, gneissic granodiorite; the magnetic susceptibility distribution in the A-A1 and C-C1 sections varies less, basically belonging to low magnetic and non-magnetic characteristics, and the magnetic susceptibility is mostly less than 100 (10-5SI). The A-A1 and C-C1 sections are seriously affected by positive metamorphism, and the content of magnetic dark minerals is reduced. From the aeromagnetic susceptibility section, the A-A1 section is divided into a tectonic mélange belt, and the lithology is mainly marble, metasandstone, and schist; the C-C1 section is divided into the granitic gneiss of the Nyainqentanglha Group of the Middle and Lower Proterozoic.

[0064] In addition, the vertical gradient zones of the longitudinal section magnetic susceptibility contour lines reflect the lithological contact boundaries of this area, such as Figure 3 The lithologic boundary between the weakly magnetic sandstone A-A1 and the high magnetic diorite B-B1 is shown in the figure; or the vertical strip-like features with lower magnetic susceptibility on both sides reflect the boundary of the fault fracture zone, based on which the positions and boundaries of the F1, F2, F3 faults and a lithologic contact fracture zone are divided, such as Figure 3 The three faults have high magnetic susceptibility on both sides, and the fault positions show vertical bands with lower magnetic susceptibility.

[0065] like Figure 1 As shown, in step S300, the preliminarily determined boundaries of different stratigraphic lithology zones and the approximate locations of fault / fracture zones are verified and adjusted by conducting field surveys.

[0066] Specifically, field geological survey and verification work is carried out based on the lithology distribution characteristics of the tunnel area determined by S100, the aerial survey magnetic susceptibility parameter values, and the stratigraphic lithology boundaries and fault / fracture zone locations delineated by the S200 aerial geophysical longitudinal section magnetic susceptibility contour lines. The above steps specifically include:

[0067] S310: Tracing the stratigraphic boundaries, determining stratigraphic boundaries through field investigation of rock sequence, lithology changes, and structural evidence;

[0068] S320: Survey the outcrop of the fault surface, conduct a comprehensive analysis through geomorphic signs, lithology changes, rock formation occurrence changes, material composition in the fault zone, etc., to determine the location and boundary of the fault, and conduct excavation work for verification when necessary;

[0069] S330: The field geological survey route should carry out field verification and supplementary investigation of continuous geological sections along and perpendicular to the geological body boundaries.

[0070] like Figure 6 The observation point G-1 in the field survey is the boundary between marble, metasandstone and granodiorite; one side of the observation point G-2 is complete granodiorite, and the other side has fault breccia, crushed rock and other fault zone materials exposed on the surface; therefore, based on the actual survey results, the stratigraphic boundaries circled on the tunnel longitudinal section magnetic susceptibility results map are fine-tuned.

[0071] like Figure 1 As shown, in step S400, drilling work is carried out on the boundaries of different stratum lithology zones and the locations of faults / fracture zones to further determine the stratum boundaries, fault / fracture zone locations and the overall direction and location of the deep underground.

[0072] Specifically, based on the stratigraphic boundaries (strata and fault boundaries) initially determined by S100, S200 and S300 and combined with the fault dips, drilling holes are arranged on both sides of the boundary lines for verification, the stratigraphic lithology and fault zone characteristics on both sides of the boundary lines are identified, and the lithology differences on both sides of the boundary lines are verified. Specifically, it includes:

[0073] S410: The principle of drilling arrangement is as follows: for the stratum lithology boundary, with the initially determined stratum boundary as the center line, a verification hole is arranged about 100m on both sides of the intersection of the line centerline and the stratum boundary; for the fault position and boundary, a verification hole is arranged in the vertical strip with lower magnetic susceptibility where the initially determined fault boundary intersects with the line centerline to find out the stratum characteristics in the fault zone; in addition, a verification hole is arranged 100m outside each of the two initially determined fault boundaries to find out the difference in lithology inside and outside the fault zone and verify the accuracy of the boundary;

[0074] S420: Carry out drilling verification work on different stratigraphic boundaries / fault fracture zones; e.g. Figure 4As shown, the ZK1 borehole reveals marble and metasandstone in the upper part and gneissic granodiorite in the lower part, and the ZK2 full borehole reveals relatively complete to complete gneissic granodiorite; Figure 5 As shown, ZK3 borehole revealed fault breccia, crushed rock, etc. According to the drilling revelations in each lithology segment, the lithology boundary interpretation and actual drilling revelations have a high degree of matching, and the effect is good.

[0075] like Figure 1 As shown, in step S100, the regional longitudinal section stratigraphic boundaries and fault / fracture zone boundaries are accurately divided through comprehensive analysis.

[0076] Specifically, regional geological data, aerial geophysical magnetic susceptibility results, field survey results, and drilling results are integrated to accurately delineate the stratum boundaries and fault / fracture zone boundaries of the tunnel longitudinal section. Figure 6 Through comprehensive analysis, the A-A1 marble and metasandstone formation boundary range, the B-B1 schist gneissic granodiorite formation boundary range, the C-C1 granite gneiss formation boundary range, and the F1, F2, F3 and broken zone formation boundary ranges were delineated.

[0077] The present invention collects regional geological data, analyzes regional lithology distribution characteristics and rock geophysical parameter background values; analyzes the distribution boundaries of different strata lithology and the approximate location of fault / fracture zones by interpreting the aerial geophysical magnetic longitudinal section map; conducts field investigation and verification based on the strata lithology boundaries and fault / fracture zone boundaries initially interpreted by aerial geophysical exploration; conducts drilling work on the strata lithology boundaries and fault / fracture zone locations of different strata, and further identifies the strata boundaries, fault / fracture zone locations and the overall trend and location of the deep underground; and accurately divides the strata boundaries and fault / fracture zone boundaries of the tunnel longitudinal section by integrating regional geological data, aerial geophysical magnetic susceptibility results maps, survey data, and drilling results. The present invention can shorten the geological survey period, quickly obtain the distribution of strata lithology and geological body boundaries in the tunnel area, and reduce the cost of drilling, geophysical exploration, and field geological surveys compared to traditional survey methods. At the same time, it solves the problem of no survey data available for reference in long sections under complex terrain conditions, ensures reasonable engineering settings, avoids geological risks, and improves the quality of survey and design.

[0078] Corresponding to the above-disclosed method for identifying the boundary of a longitudinal section of a stratum based on an aerial geophysical magnetic susceptibility result map, the embodiment of the present invention further discloses a system for identifying the boundary of a longitudinal section of a stratum based on an aerial geophysical magnetic susceptibility result map, such as Figure 7 As shown, it specifically includes:

[0079] The regional geological data analysis module is used to collect regional geological data, analyze the lithology distribution characteristics in the region, and obtain the range of aeromagnetic magnetic susceptibility parameter values ​​of various strata and lithologies in the region;

[0080] The airborne geophysical prospecting results interpretation module is used to generate longitudinal section magnetic susceptibility results maps based on airborne magnetic survey data. By analyzing the longitudinal section magnetic susceptibility results maps and combining the magnetic susceptibility parameter value ranges of various strata and lithologies, the boundaries of different strata and lithology divisions and the approximate locations of faults / fracture zones can be preliminarily determined;

[0081] Field investigation and verification module, used to verify and adjust the preliminarily determined boundaries of different strata lithology zones and the approximate locations of fault / fracture zones through field investigation;

[0082] The drilling verification module is used to carry out drilling work on the boundaries of different strata lithology zones and the locations of faults / fracture zones, and further determine the strata boundaries, fault / fracture zone locations and the overall direction and location of the deep underground;

[0083] The comprehensive analysis module is used to accurately divide the regional longitudinal section stratigraphic boundaries and fault / fracture zone boundaries through comprehensive analysis.

[0084] It should be noted that, for the detailed description of a longitudinal section stratigraphic boundary identification system based on an aerial geophysical magnetic susceptibility results map provided in an embodiment of the present invention, reference can be made to the relevant description of a longitudinal section stratigraphic boundary identification method based on an aerial geophysical magnetic susceptibility results map provided in an embodiment of the present application, which will not be repeated here.

[0085] Those skilled in the art will appreciate that all or part of the functions of the various methods in the above-mentioned embodiments can be implemented by hardware or by computer programs. When all or part of the functions in the above-mentioned embodiments are implemented by computer programs, the program can be stored in a computer-readable storage medium, and the storage medium can include: read-only memory, random access memory, disk, optical disk, hard disk, etc., and the program is executed by a computer to implement the above-mentioned functions. For example, the program is stored in the memory of the device, and when the program in the memory is executed by the processor, all or part of the above-mentioned functions can be implemented. In addition, when all or part of the functions in the above-mentioned embodiments are implemented by computer programs, the program can also be stored in a storage medium such as a server, another computer, disk, optical disk, flash disk or mobile hard disk, and can be downloaded or copied and saved in the memory of the local device, or the system of the local device is updated, and when the program in the memory is executed by the processor, all or part of the functions in the above-mentioned embodiments can be implemented.

[0086] The above specific examples are used to illustrate the present invention, which is only used to help understand the present invention and is not intended to limit the present invention. For those skilled in the art, according to the concept of the present invention, some simple deductions, modifications or substitutions can be made.

Claims

1. A method for identifying the boundary of a longitudinal section of a stratigraphic layer based on an aerial geophysical exploration magnetic susceptibility result map, characterized in that: The method comprises: Collect regional geological data, analyze the distribution characteristics of lithology in the region, and obtain the range of aeromagnetic magnetic susceptibility parameters of various strata and lithologies in the region; Generate longitudinal section magnetic susceptibility results map based on aeromagnetic data, and preliminarily determine the boundaries of different stratum lithology zones and the approximate locations of fault / fracture zones by analyzing the longitudinal section magnetic susceptibility results map and combining the magnetic susceptibility parameter value ranges of various stratum lithologies; Verify and adjust the preliminarily determined boundaries of different stratigraphic lithology zones and the approximate locations of fault / fracture zones through field investigations; Drilling work is carried out on the boundaries of different strata and lithology zones and the locations of faults / fracture zones to further identify strata boundaries, fault / fracture zone locations and the overall trend and location of deep underground; Through comprehensive analysis, the stratigraphic boundaries and fault / fracture zone boundaries of the regional longitudinal section are accurately divided.

2. A method for identifying stratum boundaries in longitudinal sections based on an aerial geophysical magnetic susceptibility result map as claimed in claim 1, characterized in that: Collect regional geological data and analyze the lithology distribution characteristics in the region, including: Collect existing geological data in the region to obtain the stratum lithology distribution types and rock characteristics in the region.

3. A method for identifying stratum boundaries in longitudinal sections based on an aerial geophysical magnetic susceptibility result map as claimed in claim 2, characterized in that: The range of magnetic susceptibility parameters of various strata and lithologies in the region are statistically obtained, including: According to the lithology of the strata on the regional geological map, the sections are roughly divided, the magnetic survey data of each lithology section collected by airborne geophysical means are analyzed, and the range of magnetic susceptibility parameter values ​​of various types of rocks inverted by airborne magnetic survey data are statistically analyzed.

4. A method for identifying stratum boundaries in longitudinal sections based on an aerial geophysical magnetic susceptibility result map as claimed in claim 1, characterized in that: Generate longitudinal magnetic susceptibility results map based on aeromagnetic data, including: First, the collected aeromagnetic data are corrected and processed; then the corrected aeromagnetic data are appropriately leveled; secondly, aeromagnetic potential field conversion is carried out; then the three-dimensional inversion of the magnetic intensity vector is performed, and the three-dimensional MVI apparent magnetic susceptibility data formed by the inversion is mapped in three dimensions, and finally the line is cut along the centerline to extract the contour map of the magnetic susceptibility of the longitudinal section of the line.

5. A method for identifying stratum boundaries in longitudinal sections based on an aerial geophysical magnetic susceptibility result map as claimed in claim 1, characterized in that: By analyzing the longitudinal section magnetic susceptibility results map and combining the magnetic susceptibility parameter value range of various stratum lithologies, the boundaries of different stratum lithology divisions and the approximate locations of fault / fracture zones are preliminarily determined, including: Determine the distribution range of regional strata lithology and strata boundaries based on the magnitude of the longitudinal section magnetic susceptibility value and the contour line morphology characteristics; The location and boundary of the fault fracture zone are determined based on the morphological characteristics of the longitudinal section magnetic susceptibility contour lines.

6. A method for identifying stratum boundaries in longitudinal sections based on an aerial geophysical magnetic susceptibility result map as claimed in claim 1, characterized in that: Through field investigation, the preliminarily determined boundaries of different strata and lithology zones and the approximate locations of fault / fracture zones are verified and adjusted, including: Tracing the stratigraphic boundaries and determining stratigraphic boundaries through field surveys; Survey the surface outcrops of the fault, determine the location and boundaries of the fault through comprehensive analysis, and conduct excavation work to verify if necessary; The field geological survey route carries out field verification and supplementary investigation of continuous geological sections along and perpendicular to the geological body boundaries.

7. A method for identifying stratum boundaries in longitudinal sections based on an aerial geophysical magnetic susceptibility result map as claimed in claim 1, characterized in that: Drilling work is carried out on the boundaries of different strata and lithology zones and the locations of faults / fracture zones to further identify the strata boundaries, fault / fracture zone locations and the overall direction and location of the deep underground, including: By determining the boundaries of the stratigraphic lithology zoning and the location of the fault / fracture zone and combining the fault dip, drilling holes are arranged on both sides of the boundary line for verification, the stratigraphic lithology and fault zone characteristics on both sides of the boundary line are identified, and the lithology differences on both sides of the boundary line are verified.

8. A method for identifying stratum boundaries in longitudinal sections based on an aerial geophysical magnetic susceptibility result map as claimed in claim 7, characterized in that: Drill holes are arranged on both sides of the boundary line for verification, including: The principles of drilling arrangement are: For the stratum lithology boundary, the initially determined stratum boundary is taken as the center line, and verification holes are arranged at preset distances on both sides of the intersection of the line center line and the stratum boundary; For the location and boundary of the fault, verification holes are arranged in the vertical strip with lower magnetic susceptibility where the fault boundary determined initially intersects with the center line of the line to identify the stratigraphic characteristics within the fault zone; In addition, verification holes are arranged at a preset distance outside the two initially determined fault boundaries to find out the differences in lithology inside and outside the fault zone and to verify the accuracy of the boundaries.

9. A method for identifying stratum boundaries in longitudinal sections based on an aerial geophysical magnetic susceptibility result map as claimed in claim 1, characterized in that: Through comprehensive analysis, the stratigraphic boundaries and fault / fracture zone boundaries of the regional longitudinal section are accurately divided, including: By integrating regional geological data, aerial geophysical magnetic susceptibility maps, field survey results, and drilling results, the regional longitudinal section stratigraphic boundaries and fault / fracture zone boundaries are accurately divided.

10. A longitudinal section stratigraphic boundary identification system based on the aerial geophysical magnetic susceptibility results map, characterized in that: The system comprises: The regional geological data analysis module is used to collect regional geological data, analyze the lithology distribution characteristics in the region, and obtain the range of aeromagnetic magnetic susceptibility parameter values ​​of various strata and lithologies in the region; The airborne geophysical prospecting results interpretation module is used to generate longitudinal section magnetic susceptibility results maps based on airborne magnetic survey data. By analyzing the longitudinal section magnetic susceptibility results maps and combining the magnetic susceptibility parameter value ranges of various strata and lithologies, the boundaries of different strata and lithology divisions and the approximate locations of faults / fracture zones can be preliminarily determined; Field investigation and verification module, used to verify and adjust the preliminarily determined boundaries of different strata lithology zones and the approximate locations of fault / fracture zones through field investigation; The drilling verification module is used to carry out drilling work on the boundaries of different strata lithology zones and the locations of faults / fracture zones, and further determine the strata boundaries, fault / fracture zone locations and the overall direction and location of the deep underground; The comprehensive analysis module is used to accurately divide the regional longitudinal section stratigraphic boundaries and fault / fracture zone boundaries through comprehensive analysis.

Images