Bauxite concealed ore body three-dimensional detection method and system

Abstract

The invention relates to the technical field of geological exploration, in particular to a bauxite concealed ore body three-dimensional detection method and system, and the method comprises the following steps: obtaining a contact zone graph recognition boundary line, calculating an intersection angle matching landform, extracting a separation line, judging an alignment relation, constructing a sideband, extracting a jump line to recognize a deflection direction, and analyzing a connection relation to mark an extension surface. And generating a three-dimensional distribution structure of the concealed ore body. According to the method, the turning area is identified through the angle change of the construction path and the slope boundary, the spreading sideband is extracted by combining the rock stratum alignment state and the breakpoint, the boundary deviation direction is judged based on the angle difference of the turning point of the jump line, and the structure extension trend of the unclosed area is tracked through the difference value of the connection angle and the distance of the discontinuous area. A stable structure graph is formed through continuous construction of a construction trend and a splicing line, clear description of a hidden ore body space distribution form is achieved, continuity and integrity of sideband recognition are enhanced, and the problems of structure misjudgment and model deformation are avoided.

Description

Technical Field

[0001] This invention relates to the field of geological exploration technology, and in particular to a three-dimensional detection method and system for concealed bauxite ore bodies. Background Technology

[0002] The field of geological exploration technology involves the investigation and research of geological information such as the Earth's internal structure, lithology, and mineral resources. Its core aspects include geophysical exploration, geochemical exploration, remote sensing, drilling and sampling, and ore body modeling. This technical field, through the collection and analysis of geological information, determines the spatial distribution, morphological structure, and reserve scale of mineral resources, serving as a fundamental link in mineral resource development and utilization. Geological exploration not only serves mining development but also plays a crucial role in hydrogeology, engineering geology, and environmental geology. As exploration targets gradually shift from outcrop ore bodies to deeper and concealed ore bodies, higher demands are placed on detection accuracy, depth, and three-dimensional imaging capabilities. Among these, traditional three-dimensional detection methods for concealed bauxite ore bodies refer to a class of technical methods for spatial positioning and three-dimensional morphological identification of bauxite ore bodies beneath overburden or not exposed at the surface. Traditional detection methods of this type typically involve geological surveys combined with geophysical techniques such as gravity, magnetics, electrical resistivity, and seismic surveys. Specific procedures include laying out survey lines to collect ground geophysical signals, using gravity gradient changes to reflect density differences and inferring bauxite distribution, identifying aluminum-rich strata using resistivity anomalies, or inverting orebody interfaces using seismic wave propagation characteristics. Additionally, a small number of boreholes can be used to verify geophysical results, improving the reliability and accuracy of interpretation. These methods often rely on two-dimensional inversion results from single or multiple geophysical data sets, supplemented by empirical inferences to construct a three-dimensional orebody model.

[0003] Traditional methods for detecting concealed bauxite ore bodies primarily focus on the collection and analysis of single-point data when identifying the spatial morphology of the ore body, neglecting the structural linkage between tectonics and geomorphology. This makes it difficult to accurately depict the distribution path and boundary variation trends of rock strata. Geophysical signals rely on single response indicators, making it difficult to accurately identify tectonic turning points in complex geological contexts. In particular, there is a lack of continuous analysis methods for identifying abrupt boundaries, structural faults, and the relationship between boundary extensions. This leads to problems such as blurred ore body boundaries, structural splicing errors, and unclear extension directions. Key structural information is often omitted during model construction, making it difficult to form a stable and closed three-dimensional boundary structure, which in turn affects the prediction of the true distribution of deep or concealed ore bodies. Summary of the Invention

[0004] To achieve the above objectives, the present invention employs the following technical solution: a three-dimensional detection method for concealed bauxite ore bodies, comprising the following steps:

[0005] S1: Obtain the distribution map of the contact zone in the bauxite area, identify the stratigraphic boundary line, calculate the intersection angle between the structural path and the slope line, determine whether the intersection angle variation zone forms a turning zone, divide the overlapping path and match the landform, and project to generate the structural distribution plotting zone.

[0006] S2: Based on the constructed distribution plotting zone, extract the interlayer separation line image within the distribution path, determine the rock layer alignment relationship, mark the breakpoint area, construct the structural edge zone and trace the extension boundary, draw the splicing diagram, and generate the rock layer structure distribution edge zone.

[0007] S3: Based on the distribution zone of the rock strata structure, extract the direction of the jump response line, record the sequence of turning points and calculate the angle change range, identify the spatial range of the boundary deviation from the path, determine the deflection direction, and generate a jump line segment turning layer.

[0008] S4: Based on the jump line segment turning layer, determine whether there is a structural extension gap at the breakpoint, analyze the edge segment connection angle and path difference, screen for unclosed expansion directions, mark the offset edge segment extension surface, and generate a missing segment extension partition map.

[0009] S5: Based on the extended partition map of the missing section, identify the angle difference between the turning point and the jump point, trace the continuous splicing zone, connect the splicing line and the boundary surface, and generate the three-dimensional distribution structure of the concealed ore body.

[0010] As a further aspect of the present invention, the structural distribution mapping zone includes the boundary of the contact zone, the distribution of stratigraphic dip boundary lines, the distribution of the intersection angle between the structural path and the slope line, the turning zone of the intersection angle variation area, the overlapping section of the contact zone and the geomorphic zone, and the spatial location projection layer. The rock strata structure distribution edge zone includes the alignment graphics of interlayer boundaries, the fault zone of rock strata combination, the extension boundary of the structural edge zone, the lithological combination splicing block, and the boundary morphology of the rock strata edge zone. The jump line segment turning layer includes the jump direction of the response line, the angle sequence of the turning point, the distribution of the angle change interval, the spatial zone of the boundary offset, and the structural turning connection path. The missing segment extension zoning map includes the location of the structural extension gap, the angle difference of the edge segment connection, the closed missing segment of the distribution direction, and the extension edge segment offset diagram. The three-dimensional distribution structure of the concealed ore body includes the structural edge zone turnback point, the connection relationship of the jump endpoint, the location section of the continuous splicing zone, the splicing line structural surface connection graphics, and the combination morphology of the closed edge zone.

[0011] As a further aspect of the present invention, the determination of whether the intersection angle variation zone forms a turning zone refers to the change in the angle between the structural path and the slope direction line and the change in the angle with the slope direction line, which shows a continuous spatial trend, and the analysis of the structural orientation deflection area.

[0012] As a further aspect of the present invention, the spatial range of the boundary deviation path refers to the continuous spatial region where the rock stratum boundary line deviates in position and direction from the set distribution path, characterizing the geometric response of structural anomalies and fracture effects.

[0013] As a further aspect of the present invention, the specific steps of S1 are as follows:

[0014] S101: Acquire remote sensing images and geological layers of the bauxite area, compare the differences in reflectance of different bands in the images, identify the location of the stratigraphic contact zone boundary, and combine the known contact line markings in the layers to perform spatial aggregation and generate a contact zone boundary line layer.

[0015] S102: Call the contact zone boundary line layer, calculate the angle between the structural path and the slope direction based on the stratigraphic dip angle parameter on the boundary line, and filter the boundary segment locations where the angle abruptly exceeds the structural variation angle threshold to obtain a set of abrupt boundary segments;

[0016] S103: Based on the set of abrupt boundary segments, the path location is compared with the contact zone boundary segment, and combined with the zoning boundary information in the geomorphic zone layer, the segment distribution coding relationship is established to generate the structural distribution plotting zone.

[0017] As a further aspect of the present invention, the specific steps of S2 are as follows:

[0018] S201: Based on the plotted area in the structured distribution plotting zone, extract the rock layer image data within the path through gradient threshold, locate the gray-scale abrupt change edge in the image and extract the interlayer separation line, combine the consistency of the orientation of the rock layer boundary on both sides, calculate the offset value, and obtain the rock layer alignment offset relationship group.

[0019] S202: Call the rock strata alignment offset relationship group, retrieve the continuity of rock strata combination in the distribution path, judge whether there is a strike fracture or interruption area at the offset of the rock strata combination, record the position range of the breakpoint segment in the path, and establish a set of structural combination side zone segments.

[0020] S203: Based on the set of structural combination sideband segments, track the boundary extension direction and retrieve the corresponding lithological image identifier, perform a stitching operation on the image area surrounding the combination sideband and draw a closed graphic to generate a rock stratum structure distribution sideband.

[0021] As a further aspect of the present invention, the specific steps of S3 are as follows:

[0022] S301: Based on the layer information within the distribution zone of the rock strata, extract the edge lines of the jump response area in the boundary layer, perform line segment tracking on the direction of continuous pixels, and mark the position nodes where the direction change is greater than the response direction change threshold to obtain the jump response turning point sequence.

[0023] S302: Call the jump response turning point sequence, calculate the angle change value between adjacent turning points, classify the angle change range into segments, identify continuous segments where the boundary direction deviates from the distribution path direction, locate the area position by combining the spatial coordinate range, and obtain the boundary deviation spatial block.

[0024] S303: Based on the boundary deviation spatial block, retrieve the trend of the angle change of the turning point in the block, and connect the continuously changing nodes according to the direction of the direction difference to construct a continuous structural line segment sequence. Then, mark the entire line segment set with a spatial layer to generate a jump line segment turning layer.

[0025] As a further aspect of the present invention, the specific steps of S4 are as follows:

[0026] S401: Based on the discontinuous segments in the transition line segment turning layer, determine whether there is spatial splicing misalignment at the coordinate positions of the breakpoints of adjacent transition line segments in the layer, and calculate the angle between the spatial projection spacing between the breakpoints and the construction extension direction vector, filter the set of breakpoints with corresponding gaps in the structural path, and obtain the structural extension gap index set.

[0027] S402: Call the structure extension gap index set, extract the direction vectors of the two jumping edge segments, calculate the difference in connection angle and the deviation of the response path distance, determine the path segments whose direction difference is higher than the turning tolerance angle threshold and whose response path distance deviation is greater than the connection limit, and establish a list of unclosed direction segments.

[0028] S403: Based on the list of unclosed directional segments, track the extension trend of the corresponding jump line segments in the layer offset segments, extract the overlapping area of ​​the segment extrapolation path and the response range, construct the layer block encoding range corresponding to the spliced ​​fragment, and generate the missing segment extension partition map.

[0029] As a further aspect of the present invention, the specific steps of S5 are as follows:

[0030] S501: Based on the partition range marked in the extended partition map of the missing section, extract the position of the turning point and the coordinates of the endpoint of the jump line segment in the direction of the construction side strip, calculate the connection angle value between the two, and record the connection pairs that are greater than the edge turning angle threshold to obtain the structural turning connection angle sequence.

[0031] S502: Call the structure turning connection angle sequence, track the position index of the corresponding turning line segment in the direction of the construction edge where the connection angle is located, mark the continuous splicing path, and extend the end of the splicing path to the position of the surface structure line in the boundary direction layer to establish an extended splicing construction path set;

[0032] S503: Based on the extended splicing construction path set, determine whether the path combination forms a closed structural loop in the three-dimensional direction, perform spatial volume construction operation on the side strip construction path that meets the closed geometric conditions, extract the spatial boundary of the construction body, and generate the three-dimensional distribution structure of the concealed ore body.

[0033] A three-dimensional detection system for concealed bauxite ore bodies includes:

[0034] The structural distribution plotting and analysis module is used to achieve S1: obtain the distribution range map of the contact zone in the bauxite area, identify the stratigraphic dip boundary line, calculate the intersection angle between the structural path and the slope line, determine whether the intersection angle variation zone forms a turning zone, divide the overlapping path and match the landform, and project to generate the structural distribution plotting zone.

[0035] The rock strata structure edge band construction module is used to realize S2: according to the structural distribution plotting band, extract the interlayer separation line image within the distribution path, determine the rock strata alignment relationship, mark the breakpoint area, construct the structural edge band and trace the extension boundary, draw the splicing diagram, and generate the rock strata structure distribution edge band.

[0036] The jump response and turning analysis module is used to implement S3: Based on the distribution edge of the rock strata structure, extract the direction of the jump response line, record the sequence of turning points and calculate the angle change range, identify the spatial range of the boundary deviation from the path, determine the deflection direction, and generate the jump line segment turning layer.

[0037] The structural missing distribution detection module is used to implement S4: based on the jump line segment turning layer, determine whether there is a structural extension gap at the breakpoint, analyze the edge segment connection angle and path difference, screen for unclosed distribution directions, mark the offset edge segment extension surface, and generate a missing segment extension partition map.

[0038] The concealed ore body three-dimensional structure reconstruction module is used to achieve S5: based on the extended partition map of the missing section, identify the angle difference between the turning point and the jump point, trace the continuous splicing zone, connect the splicing line and the boundary surface, and generate the three-dimensional distribution structure of the concealed ore body.

[0039] Compared with the prior art, the advantages and positive effects of the present invention are as follows:

[0040] In this invention, turning zones are identified by the changes in the angle of the construction path and slope boundary. The distribution side zones are extracted by combining the alignment state of the rock strata and the discontinuity points. The direction of boundary deviation is determined based on the angle difference of the turning point of the jump line. The structural extension trend of the unclosed area is tracked by the connection angle and distance difference of the discontinuous area. A stable structural pattern is formed by the continuous construction of the structural direction and splicing line. This achieves a clear depiction of the spatial distribution morphology of the concealed ore body, enhances the continuity and integrity of the side zone identification, and avoids structural misjudgment and model deformation problems. Attached Figure Description

[0041] To more clearly illustrate the technical solutions in the embodiments of the present invention, the accompanying drawings used in the description of the embodiments will be briefly introduced below. Obviously, the accompanying drawings described below are only some embodiments of the present invention. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0042] Figure 1 This is a schematic diagram of the steps of the present invention;

[0043] Figure 2 This is a detailed schematic diagram of S1 of the present invention;

[0044] Figure 3 This is a detailed schematic diagram of S2 of the present invention;

[0045] Figure 4 This is a detailed schematic diagram of S3 of the present invention;

[0046] Figure 5 This is a detailed schematic diagram of S4 of the present invention;

[0047] Figure 6 This is a detailed schematic diagram of S5 of the present invention;

[0048] Figure 7 This is a system module diagram of the present invention. Detailed Implementation

[0049] The technical solution of the present invention will now be described with reference to the accompanying drawings.

[0050] In embodiments of the present invention, words such as "exemplarily," "for example," etc., are used to indicate that something is an example, illustration, or description. Any embodiment or design described as "exemplary" in the present invention should not be construed as being more preferred or advantageous than other embodiments or designs. Specifically, the use of the word "exemplary" is intended to present the concept in a concrete manner. Furthermore, in embodiments of the present invention, the meaning expressed by "and / or" can be both, or either one.

[0051] In the embodiments of this invention, the terms "image" and "picture" may sometimes be used interchangeably. It should be noted that, without emphasizing the distinction between them, they convey the same meaning. Similarly, the terms "of," "corresponding (relevant)," and "corresponding" may sometimes be used interchangeably. It should be noted that, without emphasizing the distinction between them, they convey the same meaning.

[0052] In this embodiment of the invention, sometimes a subscript such as W1 may be written in a non-subscript form such as W1. When the difference is not emphasized, the meaning they express is the same.

[0053] To make the technical problems, technical solutions and advantages of the present invention clearer, a detailed description will be given below in conjunction with the accompanying drawings and specific embodiments.

[0054] Please see Figure 1 This invention provides a three-dimensional detection method for concealed bauxite ore bodies, comprising the following steps:

[0055] S1: Obtain the distribution range map of the contact zone in the bauxite area, identify the boundary line of the stratigraphic dip direction, calculate the intersection angle between the structural distribution path and the surface slope line, determine whether the intersection angle variation area forms a boundary turning zone, divide the contact zone overlapping path and match the distribution pattern of the landform zone, project it to the spatial layer, and generate the structural distribution plotting zone.

[0056] S2: Based on the plotted area in the structural distribution plotting zone, extract the interlayer separation line image within the structural distribution path, determine the alignment relationship of the rock layer boundary in the distribution, mark the rock layer combination fault area within the same distribution trend, establish the structural combination edge zone and trace the position of the extended boundary, draw the lithological combination splicing block graphic, and generate the rock layer structure distribution edge zone.

[0057] S3: Within the distribution zone of rock strata, extract the direction of the jump response lines in the boundary layer within the zone, record the sequence of turning points of the sudden change in response lines, calculate the angle change range of the turning points, identify the spatial range of the boundary deviation path, determine the direction of boundary deflection based on the trend of turning difference, connect the structural turning lines, and generate the turning layer of jump line segments.

[0058] S4: Based on the transition line segment to the discontinuous segment in the layer, determine whether there is a structural extension gap in the spatial splicing of the layer discontinuity point, analyze the connection angle relationship between the two ends of the discontinuous edge segment and the distance difference of the response path, screen the distribution direction of the closed figure that has not been formed, mark the edge extension surface with the offset position, and generate the missing segment extension partition map.

[0059] S5: Based on the marked areas in the extended partition map of the missing section, identify the difference in connection angle between the turning point and the jump endpoint of the side zone direction, track the position of the continuous splicing zone in the construction side zone turning line segment, connect the extended splicing line and the boundary direction layer surface, determine whether a stable closed side zone structure graphic is formed, and generate the three-dimensional distribution structure of the concealed ore body.

[0060] The structural distribution mapping includes the boundary of the contact zone, the distribution of stratigraphic dip boundaries, the distribution of the intersection angle between the structural path and the slope line, the turning zones of the intersection angle variation zone, the overlapping sections of the contact zone and the geomorphic zone, and the spatial location projection layer. The rock strata structure distribution edge zone includes the alignment graphics of interlayer boundaries, the fault zone of rock strata combination, the extension boundary of the structural edge zone, the lithological combination splicing block, and the boundary morphology of the rock strata edge zone. The jump line segment turning layer includes the jump direction of the response line, the angle sequence of the turning point, the distribution of the angle change interval, the spatial zone of the boundary offset, and the structural turning connection path. The missing segment extension zoning map includes the location of the structural extension gap, the angle difference of the edge segment connection, the closed missing segment of the distribution direction, and the extension edge segment offset diagram. The three-dimensional distribution structure of the concealed ore body includes the structural edge zone turnback point, the connection relationship of the jump endpoint, the location section of the continuous splicing zone, the splicing line structural surface connection graphics, and the combination morphology of the closed edge zone.

[0061] Please see Figure 2 The specific steps of S1 are as follows:

[0062] S101: Acquire remote sensing images and geological layers of the bauxite area, compare the differences in reflectance of different bands in the images, identify the location of the stratigraphic contact zone boundary, and combine the known contact line markings in the layers to perform spatial aggregation and generate a contact zone boundary line layer.

[0063] Remote sensing image data files of the bauxite mining area were acquired. Landsat 8 satellite imagery was selected, and reflectance data from four bands—blue (B2), red (B4), near-infrared (B5), and shortwave infrared (B7)—were retrieved. The image data were standardized to a spatial resolution of 30 meters. Radiometric calibration and atmospheric correction were performed to adjust the raw digital values ​​to the range of surface reflectance values. Each band was independently normalized, unifying data acquired from different time periods or sensors to a comparable scale. Different bands were combined differentially, and the reflectance differences between bands were compared pixel-by-pixel. Pixel-by-pixel analysis of the reflectance differences formed by any two bands was used to identify pixel areas with significant reflectance differences. Using known contact line vector data from the geological layer, coordinate point sets were extracted to form contact line sample bands. A sample band was generated around this contact line. Within a 50-meter buffer zone, the average reflectance difference of each pixel within the buffer zone is extracted from the remote sensing image. This average reflectance difference is then compared with the average reflectance difference of the corresponding band in the adjacent area outside the buffer zone. Areas with an average difference greater than 0.1 are marked as contact boundary change areas. Spatial raster connectivity analysis is then used to aggregate the difference areas, removing non-connected patches smaller than 1 hectare. Vectorized boundary extraction is performed on the remaining patches to form preliminary contact zone boundary segments. Finally, the generated boundary segments are spatially overlaid with the original geological layer contact line. The spatial overlap length of each segment is compared. If the overlap length exceeds 80% of the original contact line segment length, it is preferentially identified as an accurate contact boundary segment and its information is retained. Meanwhile, areas that do not reach this proportion are further analyzed to avoid missing potential clues of concealed ore bodies. Adjacent contact segments are then merged to form a continuous boundary line layer.

[0064] S102: Call the contact zone boundary line layer, calculate the angle between the structural path and the slope direction based on the stratigraphic dip angle parameter on the boundary line, and filter the boundary segment locations where the angle abruptly exceeds the structural variation angle threshold to obtain a set of abrupt boundary segments;

[0065] Node points are extracted at fixed intervals of 100 meters on each boundary segment. The coordinates of the first and last points of each segment are marked sequentially. The azimuth angle of each segment is calculated, and the change in direction between the two endpoints in the two-dimensional coordinate system is used as the boundary azimuth angle parameter. Simultaneously, based on the digital elevation model (DEM) layer data, the slope direction value is extracted within a 10-meter radius of each boundary segment. The mode of the slope direction values ​​within this range is used as the dominant slope direction angle parameter for that segment. The difference between the azimuth angle and the slope direction angle of each boundary segment is calculated to obtain the included angle value, forming an included angle sequence. All included angle values ​​are evaluated, and those with included angles greater than 25 degrees are selected. The segment is considered a structural change abrupt change segment, with 25 degrees as the threshold for structural angle abrupt change. This threshold is set with reference to the stability assessment standards of geological engineering blocks over the years, and is determined by the historical average value of the stratigraphic stability difference area. In the judgment process, each line segment is treated as the object of processing. Line segments with included angles greater than the threshold are marked as abrupt changes. If there are multiple consecutive abrupt change segments, they are merged into a group of abrupt change boundary segments. The start and end coordinates and corresponding attribute field information of each abrupt change boundary segment are recorded, including boundary number, strike angle, slope angle, included angle value, spatial location number, etc. Finally, all abrupt change segments are collected to form an abrupt change boundary segment layer for subsequent path comparison and structural plotting.

[0066] S103: Based on the set of abrupt boundary segments, the path location is compared with the contact zone boundary segment. Combined with the zoning boundary information in the geomorphic zone layer, the segment distribution coding relationship is established to generate the structural distribution plotting zone.

[0067] The coordinates of the center point of the boundary segment are read, and a preset path line layer is invoked accordingly. Following the spatial distance matching principle of a maximum distance not exceeding 20 meters, overlapping path segments exist in the path line layer for each abrupt boundary segment's center point. All boundary segments meeting this distance condition are marked as overlapping path segments. Spatial range identification is then performed on these path segments. The zoning line layer in the geomorphic zone layer is invoked sequentially to perform spatial overlay analysis on the path segments, obtaining the geomorphic zone code information they fall into. A one-to-one correspondence is established between the path segments and the corresponding geomorphic zone zoning, and this relationship is recorded as a path segment-geomorphic zone code lookup table. According to the structural section coding rules, each path segment is assigned a unique number. The number format is set as "structural zone code + section number". The structural zone code is set according to the geological block naming rules. For example, if the current path segment belongs to the area numbered 03 and the structural code is CTB01, then the number is set as CTB01-003. The same naming rules are applied to all path segments. A new coding field is added to the path segment layer for value assignment. A structural distribution plotting zone layer containing information such as unique number, belonging section, abrupt change section location attribute, and spatial geometry is generated. The output layer is used for subsequent analysis or mapping.

[0068] Please see Figure 3 The specific steps of S2 are as follows:

[0069] S201: Based on the plotted area in the structural distribution plotting zone, the rock layer image data within the path is extracted by gradient thresholding, the gray-scale abrupt change edge in the image is located and the interlayer separation line is extracted, and the offset value is calculated by combining the consistency of the rock layer boundary on both sides, and the rock layer alignment offset relationship group is obtained.

[0070] First, all path segments within the plotted area are extracted segment by segment. Each path segment corresponds to a set of rock strata image data. Grayscale image layers are extracted and standardized according to a fixed ratio, normalizing pixel values ​​to a grayscale range of 0 to 255. Then, the row and column coordinate matrix of each image is scanned row by row to check the grayscale value changes. Adjacent pixels with grayscale value differences greater than a set threshold are marked. This threshold is set based on 1.5 times the standard deviation of the grayscale mean of the area image. If the image mean is 128 and the standard deviation is 25, the threshold is set to 37.5. Edge pixels with grayscale jumps exceeding this value are marked as candidate points for rock strata boundaries. A linear connection operation is performed on these candidate points to form grayscale abrupt change edge lines. The rock strata unit outlines on both sides of the abrupt change edge lines are extended 20 meters to both sides of the path direction. The main trend of the extended lines on both sides is read. In the direction of the path, the angle of the two lines is calculated using the coordinate difference method, and a direction consistency judgment is performed. If the difference in the direction of the two lines does not exceed 10 degrees, the direction is considered to be consistent. Then, the horizontal projection distance of the two lines in the path direction is recorded as the offset value. If the value exceeds 5 meters, it is marked as a significant displacement. The offset data of this segment, together with the position coordinates and the direction angle, is recorded as the rock layer alignment offset relationship. For example, in the path segment P15, the two interlayer dividing lines extracted from the rock layer image are located 12 meters to the left and 8 meters to the right of the path centerline, respectively. The directions of the two lines are 85 degrees and 88 degrees, with a difference of 3 degrees and a horizontal projection difference of 4 meters, which meets the conditions for offset. Then, a rock layer alignment offset relationship group is generated for this path segment. The offset direction, offset amount and start and end coordinate range are recorded, and finally the construction of the rock layer alignment offset relationship group is completed.

[0071] S202: Call the rock strata alignment offset relationship group, retrieve the continuity of rock strata combination in the distribution path, judge whether there is a strike fracture or interruption area at the offset of the rock strata combination, record the position range of the breakpoint segment in the path, and establish a set of structural combination side zone segments.

[0072] A continuous retrieval of strata assemblages was performed along the tectonic distribution path. The start and end coordinates, offset direction, and offset value data recorded in each offset relationship were read segment by segment. Based on the path direction, the path segments were divided into segments, and a sequence of strata assemblages within each segment was established. The strata assemblages of adjacent path segments were compared one by one. If the alignment boundary between two consecutive segments was interrupted along the path direction, and the difference in offset direction between adjacent segments exceeded 15 degrees, or the change in offset value exceeded 10 meters, a strata fracture was identified at that location. During the identification process, the difference in offset direction was prioritized. If the offset direction of a segment was 12 degrees east of north, the offset direction of the next segment would be west of south. If the angle is 10 degrees, the directional difference is 22 degrees, which exceeds the set threshold and is identified as a sudden change in direction. Then, the offset value is compared numerically. For example, if the first segment has an offset of 6 meters and the next segment has an offset of 19 meters, the offset change is 13 meters, which also exceeds the 10-meter threshold and is further identified as a faulted section. All segments that meet one of the above two conditions are marked as faults, and their start and end coordinates, path segment numbers, rock layer combination codes and other information are recorded in the distribution path to form a fault identification set. Finally, all fault segments are integrated, and a set of structural combination sideband segments is generated according to the linear order of the path for the next stage of sideband tracking and image region stitching processing.

[0073] S203: Based on the set of structural combination sideband segments, track the boundary extension direction and retrieve the corresponding lithological image identifier, perform splicing operation on the image area around the combination sideband and draw a closed graphic to generate the rock strata structure distribution sideband.

[0074] First, draw the side zone extension lines along the endpoint coordinates of each segment. For each structural combination side zone segment, read its start and end coordinates and offset direction information. Extend the line 50 meters along the offset direction in the endpoint direction as a tracking line. Read the rock layer image data along the extension path direction and identify image regions in the image whose grayscale patterns are similar to the lithological characteristics of the current side zone segment. Image region identification is based on matching the grayscale mean ±20. For example, if the grayscale mean of the current side zone rock layer is 140, then continuous image regions with grayscale values ​​between 120 and 160 are selected as candidate stitching targets in the extension area image. The stitching process is processed according to the overlap of image region boundaries. If the adjacent areas of the image region boundaries overlap, the stitching is performed. If the area exceeds 80 square meters, a splicing operation is performed, connecting the boundary nodes of two image regions. The image regions are connected segment by segment to form a closed sideband structure. When the image regions on both sides of a structural combination sideband segment are completed and the first and last nodes coincide, it is considered to form a closed block. The boundary of the closed shape is smoothed, and the abrupt change in turning angle is reduced by Bézier boundary curve fitting. Finally, the closed shape is assigned an identification code according to the structural combination number, and the newly created structural distribution sideband object attribute record is added to the rock stratum structure layer, including the sideband number, the list of spliced ​​image region numbers, the closed contour area, the set of boundary node coordinates, etc., to complete the generation process of the rock stratum structure distribution sideband.

[0075] Please see Figure 4 The specific steps of S3 are as follows:

[0076] S301: Based on the layer information within the distribution zone of rock strata, extract the edge lines of the jump response area in the boundary layer, perform line segment tracking on the direction of continuous pixels, and mark the position nodes where the direction change is greater than the response direction change threshold to obtain the jump response inflection point sequence.

[0077] First, the spatial range corresponding to each closed shape in the sideband layer is read. Image processing is then performed on the overlapping boundary layers within this spatial range, extracting the grayscale value of each pixel and constructing a grayscale gradient matrix. In this matrix, pixel locations where grayscale abrupt changes exceed a set threshold are identified row by row and marked as transition edge points. This threshold is set to 1.2 times the overall grayscale standard deviation of the layer. For example, if the average grayscale value of a layer within a certain sideband is 135 and the standard deviation is... The threshold for identifying abrupt changes is approximately 5.36. Any pixel whose grayscale difference exceeds this value is identified as an edge-change pixel. Line segment tracing is performed along the continuous arrangement direction between the change pixels. The movement is tracked segment by segment in a sliding window of 5 pixels. For each segment, the direction angle of the line connecting the starting and ending pixels is calculated, and the direction angle sequence is recorded. In the angle sequence, the location where the direction angle difference between adjacent segments is greater than the response direction change threshold is identified as a turning point. This threshold is set to 20 degrees based on the actual rock strata displacement response. During processing, if the direction angle of continuous line segment 1 is 45 degrees and the direction angle of line segment 2 is 70 degrees, the angle difference is 25 degrees, which exceeds the set value. This location is determined to be a turning point of the abrupt response. All abrupt edge lines within the entire layer are processed sequentially, and the locations and spatial coordinates of all turning points that meet the direction change conditions are recorded to construct a sequence of turning points of the abrupt response for subsequent structural deviation analysis.

[0078] S302: Call the transition response inflection point sequence, calculate the angle change value between adjacent inflection points, classify the angle change range into segments, identify continuous segments where the boundary direction deviates from the distribution path direction, locate the region position by combining the spatial coordinate range, and obtain the boundary deviation spatial block;

[0079] For each pair of adjacent turning points, the angle change value is calculated. For each consecutive pair of turning points, the starting and ending angle values ​​of its connection direction are extracted. Based on the starting angle, the relative change of the ending angle is calculated, and the segments are divided according to the absolute value of the change. The division criteria are: 0 to 15 degrees is a smooth segment, 15 to 30 degrees is a transition segment, and above 30 degrees is a deviation segment. If the connection direction of a pair of turning points changes from 60 degrees to 95 degrees, the change angle is 35 degrees, which falls into the category of deviation segment. This segment is marked as a boundary direction change segment. Then, the center coordinate position of each segment of the direction change segment is read and spatially compared with the center point of the line segment of the distribution path. Overlap comparison: If the Euclidean distance between the two is less than 10 meters, it is determined that the change segment has spatial overlap with the distribution path. Then, the start and end positions of the segment on the distribution path are read and recorded as deviation segments in the path. All turning point segments that meet the definition of deviation segment and whose spatial distance does not exceed 10 meters are constructed into a set of boundary direction deviation segments. For each segment in the set, the coordinates of its start and end points, the change angle value, the classification level (smooth, transition, deviation), and the path segment number to which it belongs are recorded. Based on its spatial distribution boundary, the minimum bounding rectangle is generated and defined as the boundary deviation spatial block. A spatial block attribute table is formed and output for subsequent connection processing.

[0080] S303: Based on the boundary deviation spatial block, retrieve the trend of the angle change of the turning point in the block, and connect the continuously changing nodes according to the direction of the direction difference to construct a continuous structural line segment sequence. Then, mark the entire line segment set with spatial layers to generate a jump line segment turning layer.

[0081] The set of turning points in each block is sequentially numbered. The angle change value of every two consecutive numbered turning points is read, and it is determined whether the angle change trend is continuously rising or falling. If three consecutive sets of angle change values ​​are all positive or negative, the segment is marked as a continuous change segment. The spatial coordinates of all turning points in the continuous change segment are extracted, and these coordinate points are connected sequentially according to their numbers to form structural line segments. The distance between any two points in a structural line segment shall not exceed 20 meters. If it exceeds this distance, intermediate fitting points are inserted to supplement the connectivity, forming a continuous structural line segment sequence. After the sequence is formed, it is labeled in the structural layer, and attribute fields such as line segment number, start and end point coordinates, block number, cumulative value of direction change, and number of participating turning points are recorded. All line segment sets are uniformly assigned the layer category identifier "JBL_TURN", completing the generation of the jump structural line segment data layer, which is then superimposed on the original structural distribution layer as an independent structural disturbance reference layer, forming the jump line segment turning layer.

[0082] Please see Figure 5 The specific steps of S4 are as follows:

[0083] S401: Based on the discontinuous segments in the transition line segment turning layer, determine whether there is spatial splicing misalignment at the coordinate positions of the breakpoints of adjacent transition line segments in the layer, and calculate the angle between the spatial projection spacing between the breakpoints and the construction extension direction vector, filter the set of breakpoints with corresponding gaps in the structural path, and obtain the structural extension gap index set.

[0084] First, extract the coordinate information of the termination and start nodes of all jump line segments in the layer, and establish spatial coordinate pairs of adjacent line segment breakpoints in sequence. Then, extract the differences in the horizontal and vertical directions between the coordinates of each pair of breakpoints, and use the sum of the squares of the coordinate differences to calculate the relative position distance to determine whether there is a spatial misalignment between adjacent line segments. Then, based on the direction vector information of each jump line segment, decompose the direction vector into X and Y components, and calculate the angle difference between the angles of the direction vectors of the line segments to which each pair of breakpoints belongs. Determine whether there is a combination of breakpoints with a vector angle greater than the construction extension tolerance. This tolerance angle is set to 20 degrees, based on the statistical parameters of the regional construction line distribution tolerance. If a breakpoint pair has a greater than 20 degrees, the tolerance angle is set to 20 degrees. If the angle between the two points is 23 degrees, it meets the directional deviation standard. Then, combined with the previously calculated spatial distance, a screening is performed. If the projected spacing exceeds 12 meters (this value is the spatial tolerance threshold for judging misconnection of the on-site path), the breakpoint is judged to have a splicing misalignment phenomenon. For the breakpoint pair that meets the above two conditions, its starting node ID, ending node ID, spatial spacing value, vector angle value, the number of the jump line segment to which it belongs, and the directional consistency mark value are recorded in the layer. The breakpoint is then marked as the extended fracture gap point in the constructed path. All the fracture points that meet the standard are combined in spatial order to form a fracture point set, and a structural extended gap index set containing the above attribute fields is constructed to guide the subsequent path extension and closure logic analysis.

[0085] S402: Call the structural extension gap index set, extract the direction vectors of the two-end transition segments, calculate the difference in connection angle and the deviation of the response path distance, determine the path segments whose direction difference is higher than the turning tolerance angle threshold and whose response path distance deviation is greater than the connection limit, and establish a list of unclosed direction segments.

[0086] The direction vectors of each segment are extracted sequentially according to the jump edge segment corresponding to the breakpoint. The principal direction angle value is then extracted for each pair of direction vectors. The angle between these pairs of vectors is calculated. If the angle difference exceeds the turning tolerance angle threshold (set to 18 degrees), a directional abrupt change is considered established. Combined with the statistical frequency setting of geological structural abrupt change segments, while judging the angle difference, the connection distance of the actual response path segment between each pair of breakpoints is measured and compared with the standard path segment length benchmark (set to 15 meters). If the deviation is greater than 6 meters, the response path distance is considered to exceed the limit. For combinations that meet the criteria of a directional difference greater than 18 degrees and a path connection segment length error greater than 6 meters, an abnormal connection segment is recorded. At the same time, information such as the path segment number, path node start and end identifiers, specific directional difference value, path distance difference value, and layer number of the associated jump line segment is extracted to form an unclosed path structure marker entry. All path segments that meet the conditions are combined to form an unclosed directional segment list, and a unique identification ID is assigned to each record to ensure that it can be located and tracked in the path construction closure logic verification and segment backfilling operation. This list serves as the basic input source for the subsequent construction of extended blocks.

[0087] S403: Based on the list of unclosed directional segments, track the extension trend of the corresponding jump line segments in the offset edge segments in the layer, extract the overlapping area of ​​the edge segment extrapolation path and the response range, construct the layer block encoding range corresponding to the spliced ​​fragment, and generate the missing segment extension partition map.

[0088] The process reads the termination position of the corresponding transition segment in each record, and extends the direction based on the direction vector of the transition segment at that point. The default extension length is set to 40 meters. The extension process extrapolates the coordinates segment by segment according to the unit direction component, forming a continuous directional extrapolation path. Then, a response space buffer zone with a radius of 30 meters centered on this extrapolation path is constructed in the layer. It is then determined whether this buffer zone overlaps with other transition segments, structural boundaries, or image region borders in the layer. Spatial encoding extraction is performed on all overlapping areas to extract the partition layer number, region encoding identifier, coordinate set of the overlapping boundary, and overlap length value. Records shorter than 15 meters or with an overlap area of ​​less than 150 square meters are discarded. Only candidate segments that meet the spatial overlap criteria are retained. For each retained segment, a layer block encoding operation is performed, uniformly labeling the segment ID, overlap area block number, extension direction code, associated jump segment number, and matching confidence level (divided into high, medium, and low levels according to the overlap length ratio). Finally, a spatial block data set corresponding to the missing segment is generated, and a missing segment extension partition map is established. This map has the characteristics of structural splicability, clear encoding, and complete associated jump structure logic, and can be used as a direct input layer in subsequent path structure restoration and sideband completion scenarios.

[0089] Please see Figure 6 The specific steps of S5 are as follows:

[0090] S501: Based on the partition range marked in the extended partition map of the missing section, extract the position of the turning point and the coordinates of the endpoint of the jump line segment in the direction of the construction side strip, calculate the connection angle value between the two, and record the connection pairs that are greater than the edge turning angle threshold to obtain the structural turning connection angle sequence.

[0091] First, the boundary coordinate set in the corresponding layer of each partition is read, the set of structural line segments of the sideband direction is extracted, and the main direction angle is calculated according to the start and end coordinates of each line segment. The nodes with turning changes in the sideband direction are marked as the turning point positions. Then, the coordinates of the end points of all jumping line segments in the jumping line segment turning layer are extracted. A one-to-one pairing relationship is established between the endpoint of each jumping line segment and the coordinates of the adjacent turning point. The direction angle of the connecting line constructed between the turning point and the endpoint of the jumping line segment is calculated and compared with the direction angle of the corresponding sideband direction segment to obtain the angle difference. It is determined whether the difference exceeds the edge turning angle judgment threshold. The threshold is set to 22 degrees. The reasonable criterion is determined by referring to the statistical angle range of multiple structural turning points in the reference area. When the difference of a pair of connecting angles is 28 degrees, it meets the screening condition and is recorded. The start number, end number, direction angle difference, spatial coordinate pair, the code of the sideband to which it belongs, and the jumping structure number of the pair of connecting angles are extracted. All records that meet the connection angle greater than the set threshold are combined into a structural turning connection angle sequence for subsequent extension processing in the scenario of constructing path direction continuity judgment.

[0092] S502: Call the structural turning connection angle sequence, track the position index of the corresponding turning line segment in the direction of the construction edge where the connection angle is located, mark the continuous splicing path, and extend the end of the splicing path to the position of the surface structure line in the boundary direction layer to establish an extended splicing construction path set;

[0093] Read the sideband direction segment number and the jump line segment endpoint number from each connection angle record in the sequence. Based on the number, quickly index the corresponding turn-off point in the main direction structure line segment in the sideband structure layer. Further obtain the line segment's arrangement order in the layer, and determine if it is located at the connection point between consecutive segments in the sideband path structure. If the continuous difference between the starting direction of the connection angle and the direction angle of the subsequent line segment is less than 10 degrees, it is considered to meet the conditions for a continuous splicing path. Record the continuous segment set number, turn-off angle value, and path direction trend of this path, and construct the line segment combination relationship of this splicing path. The end of the path continues to extend along the current structural direction. Read the coordinates of the end point of the direction segment in the side strip direction layer and extend it a certain distance (set to 60 meters) according to the same direction angle to generate an extension line segment. Perform spatial intersection determination with the construction surface line segment in the surface structure line layer. If the shortest distance between the extension segment and the surface structure line is less than 10 meters, it is marked as connected. Record the number, direction angle, and coordinate position of the connection point in the surface structure line. Take the last node in the splicing path and the node of the surface structure line as the beginning and end connection line of the closed path segment to generate a new splicing path segment set. Finally, form an extended splicing construction path set and uniformly identify, number and archive the path set according to the construction number and segment range.

[0094] S503: Based on the extended splicing construction path set, determine whether the path combination forms a closed structural loop in the three-dimensional direction, perform spatial volume construction operation on the side construction path that meets the closed geometric conditions, extract the spatial boundary of the construction body, and generate the three-dimensional distribution structure of the concealed ore body.

[0095] For each spliced ​​path in the path set, the starting point, turning point, and ending point are reconstructed using spatial coordinates. Z-axis elevation values ​​are introduced into 3D space, and elevation values ​​are supplemented for path segment nodes along the vertical direction. If any values ​​are missing, geological profile data from the corresponding location is used for estimation. A 3D polyline structure is constructed. Then, the spatial projection range of neighboring path segments is retrieved for each spliced ​​path. If two path segments intersect in the X and Y planes, and their Z-axis elevation start and end ranges overlap by more than 70%, they are marked as potential 3D closed segments. Adjacent segments are then grouped according to their structural path numbers. Envelope calculations are performed on the path segments of the section to detect whether a circular structural trajectory is formed. If a closed loop structure is formed, the vertex coordinate set of the enclosed space is further extracted to construct the surface structure of the three-dimensional structure. The surfaces are stitched together according to the node docking rules to form a closed structural space entity model. Each closed model is assigned a unique number and its boundary coordinates, volume estimate, surface area, structural path start and end segment numbers, and internal structural zone numbers are recorded. Finally, a three-dimensional distribution structure of the concealed ore body is generated and written into the three-dimensional layer for exploration reference.

[0096] Please see Figure 7A three-dimensional detection system for concealed bauxite ore bodies, including:

[0097] The structural distribution plotting and analysis module is used to achieve S1: obtain the distribution range map of the contact zone in the bauxite area, identify the stratigraphic dip boundary line, calculate the intersection angle between the structural path and the slope line, determine whether the intersection angle variation zone forms a turning zone, divide the overlapping path and match the landform, and project to generate the structural distribution plotting zone.

[0098] The rock strata structure edge band construction module is used to realize S2: based on the structural distribution plotting zone, extract the interlayer separation line image within the distribution path, determine the rock strata alignment relationship, mark the breakpoint area, construct the structural edge band and trace the extension boundary, draw the splicing diagram, and generate the rock strata structure distribution edge band.

[0099] The jump response and turning analysis module is used to implement S3: based on the distribution of rock strata, extract the direction of jump response lines, record the sequence of turning points and calculate the angle change range, identify the spatial range of boundary deviation from the path, determine the deflection direction, and generate a jump line segment turning layer;

[0100] The structural missing distribution detection module is used to implement S4: based on the jump line segment turning layer, it determines whether there is a structural extension gap at the breakpoint, analyzes the edge segment connection angle and path difference, screens unclosed distribution directions, marks the offset edge segment extension surface, and generates a missing segment extension partition map.

[0101] The concealed ore body 3D structure reconstruction module is used to achieve S5: based on the extended partition map of the missing section, identify the angle difference between the turning point and the jump point, trace the continuous splicing zone, connect the splicing line and the boundary surface, and generate the 3D distribution structure of the concealed ore body.

[0102] The above description is merely a specific embodiment of the present invention, but the scope of protection of the present invention is not limited thereto. Any variations or substitutions that can be easily conceived by those skilled in the art within the technical scope disclosed in the present invention should be included within the scope of protection of the present invention. Therefore, the scope of protection of the present invention should be determined by the scope of the claims.

Claims

1. A three-dimensional detection method for concealed bauxite ore bodies, characterized in that, Includes the following steps: S1: Obtain the distribution map of the contact zone in the bauxite area, identify the stratigraphic boundary line, calculate the intersection angle between the structural path and the slope line, determine whether the intersection angle variation zone forms a turning zone, divide the overlapping path and match the landform, and project to generate the structural distribution plotting zone. S2: Based on the constructed distribution plotting zone, extract the interlayer separation line image within the distribution path, determine the rock layer alignment relationship, mark the breakpoint area, construct the structural edge zone and trace the extension boundary, draw the splicing diagram, and generate the rock layer structure distribution edge zone. S3: Based on the distribution zone of the rock strata structure, extract the direction of the jump response line, record the sequence of turning points and calculate the angle change range, identify the spatial range of the boundary deviation from the path, determine the deflection direction, and generate a jump line segment turning layer. S4: Based on the jump line segment turning layer, determine whether there is a structural extension gap at the breakpoint, analyze the edge segment connection angle and path difference, screen for unclosed expansion directions, mark the offset edge segment extension surface, and generate a missing segment extension partition map. S5: Based on the extended partition map of the missing section, identify the angle difference between the turning point and the jump point, trace the continuous splicing zone, connect the splicing line and the boundary surface, and generate the three-dimensional distribution structure of the concealed ore body.

2. The three-dimensional detection method for concealed bauxite ore bodies according to claim 1, characterized in that, The structural distribution mapping zone includes the boundary of the contact zone, the distribution of stratigraphic dip boundary lines, the distribution of the intersection angle between the structural path and the slope line, the turning zones of the intersection angle variation zone, the overlapping sections of the contact zone and the geomorphic zone, and the spatial location projection layer. The rock strata structure distribution zone includes the alignment graphics of interlayer boundaries, the fault zone of rock strata combination, the extension boundary of the structural zone, the lithological combination splicing block, and the boundary morphology of the rock strata zone. The jump line segment turning layer includes the jump direction of the response line, the angle sequence of the turning point, the distribution of the angle change interval, the spatial zone of the boundary offset, and the structural turning connection path. The missing segment extension zoning map includes the location of the structural extension gap, the angle difference of the segment connection, the closed missing segment of the distribution direction, and the extension segment offset diagram. The three-dimensional distribution structure of the concealed ore body includes the structural zone turnback point, the connection relationship of the jump endpoint, the location section of the continuous splicing zone, the splicing line structural surface connection graphics, and the combination morphology of the closed zone.

3. The three-dimensional detection method for concealed bauxite ore bodies according to claim 1, characterized in that, The determination of whether the angle variation zone forms a turning zone refers to the change in the angle between the structural path and the slope direction line, which shows a continuous spatial trend, and the analysis of the structural orientation deflection area.

4. The three-dimensional detection method for concealed bauxite ore bodies according to claim 1, characterized in that, The spatial range of the boundary deviation path refers to the continuous spatial region where the rock stratum boundary line deviates in position and direction from the set distribution path, characterizing the geometric response of structural anomalies and fracture effects.

5. The three-dimensional detection method for concealed bauxite ore bodies according to claim 1, characterized in that, The specific steps of S1 are as follows: S101: Acquire remote sensing images and geological layers of the bauxite area, compare the differences in reflectance of different bands in the images, identify the location of the stratigraphic contact zone boundary, and combine the known contact line markings in the layers to perform spatial aggregation and generate a contact zone boundary line layer. S102: Call the contact zone boundary line layer, calculate the angle between the structural path and the slope direction based on the stratigraphic dip angle parameter on the boundary line, and filter the boundary segment locations where the angle abruptly exceeds the structural variation angle threshold to obtain a set of abrupt boundary segments; S103: Based on the set of abrupt boundary segments, the path location is compared with the contact zone boundary segment, and combined with the zoning boundary information in the geomorphic zone layer, the segment distribution coding relationship is established to generate the structural distribution plotting zone.

6. The three-dimensional detection method for concealed bauxite ore bodies according to claim 1, characterized in that, The specific steps of S2 are as follows: S201: Based on the plotted area in the structured distribution plotting zone, extract the rock layer image data within the path through gradient threshold, locate the gray-scale abrupt change edge in the image and extract the interlayer separation line, combine the consistency of the orientation of the rock layer boundary on both sides, calculate the offset value, and obtain the rock layer alignment offset relationship group. S202: Call the rock strata alignment offset relationship group, retrieve the continuity of rock strata combination in the distribution path, judge whether there is a strike fracture or interruption area at the offset of the rock strata combination, record the position range of the breakpoint segment in the path, and establish a set of structural combination side zone segments. S203: Based on the set of structural combination sideband segments, track the boundary extension direction and retrieve the corresponding lithological image identifier, perform a stitching operation on the image area surrounding the combination sideband and draw a closed graphic to generate a rock stratum structure distribution sideband.

7. The three-dimensional detection method for concealed bauxite ore bodies according to claim 1, characterized in that, The specific steps for S3 are as follows: S301: Based on the layer information within the distribution zone of the rock strata, extract the edge lines of the jump response area in the boundary layer, perform line segment tracking on the direction of continuous pixels, and mark the position nodes where the direction change is greater than the response direction change threshold to obtain the jump response turning point sequence. S302: Call the jump response turning point sequence, calculate the angle change value between adjacent turning points, classify the angle change range into segments, identify continuous segments where the boundary direction deviates from the distribution path direction, locate the area position by combining the spatial coordinate range, and obtain the boundary deviation spatial block. S303: Based on the boundary deviation spatial block, retrieve the trend of the angle change of the turning point in the block, and connect the continuously changing nodes according to the direction of the direction difference to construct a continuous structural line segment sequence. Then, mark the entire line segment set with a spatial layer to generate a jump line segment turning layer.

8. The three-dimensional detection method for concealed bauxite ore bodies according to claim 1, characterized in that, The specific steps of S4 are as follows: S401: Based on the discontinuous segments in the transition line segment turning layer, determine whether there is spatial splicing misalignment at the coordinate positions of the breakpoints of adjacent transition line segments in the layer, and calculate the angle between the spatial projection spacing between the breakpoints and the construction extension direction vector, filter the set of breakpoints with corresponding gaps in the structural path, and obtain the structural extension gap index set. S402: Call the structure extension gap index set, extract the direction vectors of the two jumping edge segments, calculate the difference in connection angle and the deviation of the response path distance, determine the path segments whose direction difference is higher than the turning tolerance angle threshold and whose response path distance deviation is greater than the connection limit, and establish a list of unclosed direction segments. S403: Based on the list of unclosed directional segments, track the extension trend of the corresponding jump line segments in the layer offset segments, extract the overlapping area of ​​the segment extrapolation path and the response range, construct the layer block encoding range corresponding to the spliced ​​fragment, and generate the missing segment extension partition map.

9. The three-dimensional detection method for concealed bauxite ore bodies according to claim 1, characterized in that, The specific steps of S5 are as follows: S501: Based on the partition range marked in the extended partition map of the missing section, extract the position of the turning point and the coordinates of the endpoint of the jump line segment in the direction of the construction side strip, calculate the connection angle value between the two, and record the connection pairs that are greater than the edge turning angle threshold to obtain the structural turning connection angle sequence. S502: Call the structure turning connection angle sequence, track the position index of the corresponding turning line segment in the direction of the construction edge where the connection angle is located, mark the continuous splicing path, and extend the end of the splicing path to the position of the surface structure line in the boundary direction layer to establish an extended splicing construction path set; S503: Based on the extended splicing construction path set, determine whether the path combination forms a closed structural loop in the three-dimensional direction, perform spatial volume construction operation on the side strip construction path that meets the closed geometric conditions, extract the spatial boundary of the construction body, and generate the three-dimensional distribution structure of the concealed ore body.

10. A three-dimensional detection system for concealed bauxite ore bodies, characterized in that, The system is used to implement the three-dimensional detection method for concealed bauxite ore bodies according to any one of claims 1-9, and the system comprises: The structural distribution plotting and analysis module is used to achieve S1: obtain the distribution range map of the contact zone in the bauxite area, identify the stratigraphic dip boundary line, calculate the intersection angle between the structural path and the slope line, determine whether the intersection angle variation zone forms a turning zone, divide the overlapping path and match the landform, and project to generate the structural distribution plotting zone. The rock strata structure edge band construction module is used to realize S2: according to the structural distribution plotting band, extract the interlayer separation line image within the distribution path, determine the rock strata alignment relationship, mark the breakpoint area, construct the structural edge band and trace the extension boundary, draw the splicing diagram, and generate the rock strata structure distribution edge band. The jump response and turning analysis module is used to implement S3: Based on the distribution edge of the rock strata structure, extract the direction of the jump response line, record the sequence of turning points and calculate the angle change range, identify the spatial range of the boundary deviation from the path, determine the deflection direction, and generate the jump line segment turning layer. The structural missing distribution detection module is used to implement S4: based on the jump line segment turning layer, determine whether there is a structural extension gap at the breakpoint, analyze the edge segment connection angle and path difference, screen for unclosed distribution directions, mark the offset edge segment extension surface, and generate a missing segment extension partition map. The concealed ore body three-dimensional structure reconstruction module is used to achieve S5: based on the extended partition map of the missing section, identify the angle difference between the turning point and the jump point, trace the continuous splicing zone, connect the splicing line and the boundary surface, and generate the three-dimensional distribution structure of the concealed ore body.

Images