Method, system and equipment for establishing complex geological surfaces based on precise planar cutting diagrams

Abstract

The present invention belongs to the field of three-dimensional geological modeling and discloses a method, system and equipment for establishing a complex geological surface based on a precise planar cut map. The method for establishing a complex geological surface includes: setting a modeling range and establishing a range clipping box; importing raw data by classification and setting a credibility level; selecting a data reference range and reference data points, projecting the data points onto a planar cut surface, and constructing planar cut points for various types of data on the planar cut surface; performing data correction on each planar cut surface; establishing a precise planar cut map and fitting a complex three-dimensional model surface; setting the credibility of the model surface, and the precise model surface can be corrected for two-dimensional sections later. The present invention corrects data by screening highly credible exploration data, constructs a complex geological model surface based on a precise planar cut map, reduces the number of data crossovers and errors during the fitting process, improves modeling efficiency and model accuracy, solves the problem of difficulty and low accuracy in constructing distorted complex geological bodies in actual engineering, and reduces the difficulty of revising the model surface later.

Description

Technical Field

[0001] The present invention belongs to the field of three-dimensional geological modeling, and in particular relates to a method, system, device and terminal for establishing a complex geological surface based on a precise planar cut map. Background Art

[0002] Currently, 3D geological modeling is a synthetic discipline based on data and information analysis, or rather, a discipline that integrates various disciplines. 3D modeling can construct 3D geological models based on existing 2D exploration data, point cloud data, and field geological analysis data. 3D geological modeling has two main functions: one is to provide a foundational model for numerical simulations, and the other is to facilitate overall reservoir evaluation. The creation of complex 3D geological surfaces not only visually demonstrates the uniquely complex topography and geological structural features of large-scale projects, but also provides a clear understanding of the 3D structural relationships between different complex geological surfaces and key structures. Geological analysis allows for accurate analysis of the impact of complex geological structures on projects. Simultaneously with construction, data control points are continuously added through information updates during the construction phase. Complex geological surfaces are updated by combining new and existing data, further deepening understanding of unique geological structures. Geological analysis enables proactive assessment of potential geological hazards.

[0003] Through the above analysis, the problems and defects of the existing technology are as follows:

[0004] (1) Most existing 3D visualization models are simple, with low fitting data volume, large data gaps, and are inferred based on the equipartition interpolation algorithm over a large range. They have a certain 3D schematic function, but they cannot reflect the complexity and particularity of engineering geology.

[0005] (2) The establishment of existing complex three-dimensional models requires the fusion processing of a large amount of data, which has high requirements for human-computer interaction and computer processors. The workload is huge, the modeling efficiency is low, and the accuracy mostly does not meet the requirements of reverse calibration of two-dimensional data. Summary of the Invention

[0006] In response to the problems existing in the prior art, the present invention provides a method for establishing a complex geological surface based on a precise planar cutting map.

[0007] The present invention is implemented as follows: a method for constructing complex geological surfaces based on precise planar cuts. Different types of control points are projected onto a designated planar cut according to a specific method. Precise deviation corrections are performed on each type of planar cut data by screening highly reliable exploration data, and a highly reliable complex geological model surface is constructed using the precise planar cut. This method primarily includes: selecting highly reliable control point data as the primary basis, discarding obviously erroneous two-dimensional information, and reducing the workload and error rate of later computer fitting calculations; projecting different control points onto the designated planar cut data, and constructing precise planar cuts for each elevation using the projection results and data deviation corrections of various data points from the initial stage; and, when large blank areas in the data are present, extending and supplementing the blank areas with data points.

[0008] Furthermore, the data projection method combines the surrounding geological knowledge to propose a specific method for projecting different control points onto a specified horizontal tangent plane. This makes each elevation horizontal tangent line constructed later have certain specificity and similarity with adjacent elevation horizontal tangent lines, which is the core basis for accurately fitting complex model surfaces in the later stage.

[0009] Data correction: Based on the projection results of each cross-section data point, by retaining the exploration data with higher credibility, the problem of inconsistency between data with different credibility is resolved, avoiding the formation of some erroneous pinch-out points in the later fitting. Constructing accurate cross-section maps: Based on the projection results of various types of data points in the early stage and data correction, accurate cross-section maps for each elevation are constructed. To a certain extent, the fitting direction is specified, avoiding erroneous interpolation that violates geological knowledge, achieving accurate fitting of complex model surfaces, and facilitating later revisions.

[0010] Data point extension to supplement blank areas: When the data blank area is large, an effective method of data extension and supplement is proposed.

[0011] Furthermore, the method for establishing a complex geological surface based on a precise planar cut map includes the following steps:

[0012] Step 1: For closed geological surfaces and some non-closed geological surfaces with sufficient data, set the modeling range (X min ,X max ),(Y min ,Y max ),(Z min ,Z max ), establish the range clipping box;

[0013] For non-closed geological surfaces with blank data, partition modeling is performed;

[0014] Step 2: Classify and import the original data, eliminate outlier data with obvious errors, and set the credibility level, Class I and Class II credibility;

[0015] Step 3: Depending on the sparseness of the existing flat-cut map elevation, a flat-cut map with a specific elevation can be added, and the selected elevation is used as the base surface, and the 1 / 2 elevation difference range is extended to the adjacent elevation surface as the data reference range;

[0016] Step 4: Project the data points within the data reference range onto the tangent plane to construct tangent points for each type of data;

[0017] Step 5: Correct the cross-section projection data to create accurate cross-section lines and cross-section diagrams;

[0018] Step 6: Fit the three-dimensional model surface based on the precise horizontal tangent interpolation of each elevation, use the clipping box S1 to perform Boolean operations, and set the model surface credibility;

[0019] Step 7: If the credibility level is lower than Class I, it will be updated in time according to subsequent data, and steps 2 to 7 will be repeated in sequence;

[0020] In step eight, all Class I 3D geological surfaces are gathered in the modeling area, and a vertical plane is established according to the 2D section position to obtain the intersection line with the complex 3D model surface, which is the corrected 2D section line.

[0021] Furthermore, the data within the Class I confidence level are accurate data, and the accurate data include surface control points, adit and borehole control points, and actual excavation section control points;

[0022] The data within the Class II confidence level include sub-precise data and auxiliary points. The sub-precise data are measured occurrence inference profiles on the surface or in the horizontal tunnel, and the auxiliary points are auxiliary points and lines added based on geological analysis and modeling.

[0023] Furthermore, the establishment of the precise horizontal cutting map includes comparing the existing elevation horizontal tangent lines and accurately constraining the Class I and Class II credibility data points, thereby establishing the precise horizontal cutting map.

[0024] Furthermore, the establishment of the precise tangent map also includes, for the newly added tangent lines of the elevation, if the directions of the adjacent tangent lines of the elevation are clear, the projection of some control points can be used as a reference to supplement the data according to the projection method in step 4;

[0025] If there is no actual excavation section near the insufficient data to provide the projection direction, two points Q1 and Q2 with similar X and Y coordinates in the horizontal tangent lines at different elevations are selected to form the projection vector

[0026] If data in some areas are missing, data can be supplemented based on the surface trace control points or the Class I control points in between according to the projection method in step 4;

[0027] Finally, the tangent points of each type of data are reasonably connected into lines, the control points of Class II data that contradict Class I data are deleted, the credibility data points of Class I and Class II are accurately constrained, and an accurate tangent map is established.

[0028] Furthermore, the fitting model surface has different Class I model surface discrimination criteria according to different site types and complexity. Assuming that the model surface n is a geological structure surface in a mountainous area, the maximum spacing between the standard exploration points is R(A n )=50m, then the critical data surface density of model surface n (number / m 2 )for:

[0029]

[0030] Among them, the data area density P(A n ) is the ratio of the number of model surface data points to the model surface area. The Class I model surface is the mapping line density of the surface trace line. Model surface data surface density P(A n )≥P(A ncj ) and the surface density of type I data Model surface.

[0031] Furthermore, the non-closed geological surface includes a lithologic interface, an unloading surface, and a deep and large torsional fault. If the non-closed geological surface has a large data blank area, the process of establishing the non-closed geological surface is as follows:

[0032] S101, partition modeling, divide n modeling areas according to data distribution:

[0033] {(X 1min ,X 1max ),(Y 1min ,Y 1max ),(Z 1min ,Z 1max )},...,{(X nmin ,X nmax ),(Y nmin ,Y nmax ),(Z nmin ,Z nmax )}, and establish the cropping boxes for each area;

[0034] S102, setting the number of blank data areas to i, and adding i blank data area model surfaces one by one;

[0035] S103, after all the model surfaces of the data blank areas are supplemented and fitted, the regional models are merged to obtain the entire model, and the lowest credibility level of the regional model (excluding the data blank areas) is the credibility level of the entire model.

[0036] Furthermore, the process of supplementing the data blank area includes the following steps:

[0037] First, select part of the data in the precise data control area connected to the blank area as the reference basis, select a data reference range equal to 1 / 2 of the data blank area range, and set the data reference range control data point set of each area as M1 and M2. Some special geological control points can be discarded.

[0038] If the blank area of ​​the data is less than 10% of the overall data range, the point set of type I data in M1 and M2 is selected as the priority fitting point, and the blank area is directly fitted;

[0039] If it is greater than 10%, the data extension method is used to extend and supplement the adjacent data blank areas based on the point set data in M1 and M2, and the geological surface of the data blank areas is fitted after the extension and supplement.

[0040] Furthermore, the data extension method process is as follows:

[0041] Select the lowest point nodes J1 and J2 at the boundary between the two areas and the blank, and set the vector In order to avoid contradictions in the extended data, data above 1 / 2 of the reference range is selected for extension. The extension vector of the data in point set M1 is: The extension vector of the data in point set M2 is When the blank area is much larger than 10%, the value of q can be appropriately increased.

[0042] Another object of the present invention is to provide a system for establishing a complex geological surface based on a precise planar section, the system comprising:

[0043] The data classification module sets the credibility level of data according to its source classification, which is used to determine the priority of fitting between data and provide a basis for subsequent data correction.

[0044] The data point projection module projects control points within the data reference range, adjacent geological control points or surface trace control points to the specified horizontal section based on a specific method, providing a basis for the subsequent construction of accurate horizontal section maps.

[0045] The module for establishing precise tangent lines includes data correction and precise tangent line construction. Data correction is used to discard conflicting low-reliability data. Precise tangent line construction is used to precisely and reasonably connect the corrected Class I and Class II credibility data into a line to obtain a precise tangent line, which is the core foundation of the entire system.

[0046] The model surface building module is used to build a complex 3D model surface using the precise planar cut image, including overall modeling and partitioned modeling. Partitioned modeling further includes partitioned modeling and supplementary modeling of data blank areas. The complex 3D model surface is partitioned according to the location of the data blank areas, and reference data points are selected at the boundaries of each area to directly fit or extend the blank area data to fit the model surface.

[0047] The model surface credibility setting module is used to set the credibility level of the precise and complex three-dimensional model surface based on the data credibility classification results.

[0048] Another object of the present invention is to provide a computer device, which includes a memory and a processor, wherein the memory stores a computer program, and when the computer program is executed by the processor, the processor executes the steps of the method for establishing a complex geological surface based on a precise planar cutting surface.

[0049] Another object of the present invention is to provide an information data processing terminal, which is used to implement the complex geological surface establishment system based on the precise planar section.

[0050] In combination with the above technical solutions and the technical problems solved, the advantages and positive effects of the technical solutions to be protected by the present invention are as follows:

[0051] The present invention sets data credibility and corrects the cross-section data by screening exploration data with high credibility, thereby solving the problem of contradictions between data and avoiding the intersection of different types of control data points when establishing complex model fitting. If new data points are added at a later stage, additional revisions and corrections can be made on the specified cross-section surface, so that the complex geological model surface can be reconstructed quickly and accurately, facilitating secondary revisions.

[0052] The present invention proposes a method for supplementing projection of horizontal tangent data points, which is based on the distribution characteristics of adjacent geological control points or surface trace control points to a certain extent, has certain geological thinking, and is non-traditional interpolation calculation.

[0053] The present invention establishes a precise horizontal tangent line, which, to a certain extent, combines geological knowledge and is equivalent to specifying the fitting direction of the model surface, avoiding the possibility of erroneous interpolation due to direct fitting. It also greatly reduces the workload and number of errors in the later calculations of the computer, and improves the degree of fit between the data control points and the model surface and the efficiency of the later revision of the model.

[0054] The present invention proposes a method for partitioning and modeling models with data blank areas, which reduces the requirements for computer processors and the number of calculation errors to a certain extent; the present invention proposes screening reasonable extended data when establishing data blank areas, reducing the influence of individual special control points in the edge area on the fitting shape of the data blank area. In addition, for models with a data blank area range that is greater than 10% of the overall range, data is extended to supplement the blank area, thereby improving the integration of the data blank area with adjacent model surfaces.

[0055] The expected benefits and commercial value of the technical solution of the present invention after transformation are as follows: exploration data is obtained by screening according to the data credibility level; a three-dimensional geological model is established using the precise flat-cut map to intuitively display the geological information of complex terrain and realize data sharing, providing a basic model for numerical simulation and technical support for the overall evaluation of oil reservoirs; for the geological information of concern in engineering line selection, a three-dimensional geological model is constructed based on data control points, which has powerful visualization functions, dynamically displays geological information in three dimensions and realizes data sharing; and the combination of new and old data to update the three-dimensional geological model can make advance judgments on geological disaster risks.

[0056] The technical solution of the present invention fills a technological gap in the industry at home and abroad: it is an innovation in the method of establishing complex three-dimensional geological surfaces. The present invention projects data points onto specified planar sections based on a specific method to construct various types of data points on each planar section. The present invention corrects existing data by screening highly credible exploration data, establishes a precise planar section map, and thereby fits complex three-dimensional model surfaces. For non-closed geological surfaces with blank data, the present invention selects reference data near the boundary for direct fitting or extends the data to fit a supplementary model surface, and finally integrates it with the precise three-dimensional model surface established by partitioning to obtain an overall model. This method of establishing complex three-dimensional geological surfaces based on precise planar section maps is currently unprecedented at home and abroad.

[0057] The technical solution of the present invention solves a technical problem that people have always wanted to solve but have never been able to solve successfully: based on the precise flat-cut graph, the present invention effectively solves the problems in the prior art such as low amount of fitting data, large data gaps, low modeling efficiency and low accuracy.

[0058] The mathematical model established by this invention has the following advantages: 1. Data source credibility classification prioritizes fitting data, ensuring that measured control points consistently have the primary influence on model morphology; The model surface credibility grading aligns with survey specifications to a certain extent, aligning with actual practice. The establishment of complex geological surfaces requires relatively large amounts of data. By combining the spacing of exploration points during the detailed survey phase, a critical value for the surface density of model surface data control points is defined. Using this critical value as the basis for delineating the model surface credibility grading has certain engineering practical significance. 2. The data range selection and point projection data supplementation method incorporates reasonable assumptions about the surrounding geological environment, incorporates geological thinking, and employs non-traditional interpolation calculations. The establishment of a precise tangent map specifies the fitting direction, to a certain extent avoiding the interpolation errors often found in modeling software. 3. For models with data gaps, the zoning modeling method reduces computer processor requirements and the number of computational errors. By screening data gaps and rationally extending data, the influence of individual control points in marginal areas on the fitted morphology of these gaps is reduced. When the data gap range exceeds 10% of the overall data range, the data extension supplementation method improves the integration of the gap-filled fitting results with adjacent model surfaces. BRIEF DESCRIPTION OF THE DRAWINGS

[0059] Figure 1 This is a flow chart of a method for establishing a complex geological surface based on a precise planar cut diagram provided by an embodiment of the present invention.

[0060] Figure 2 This is a schematic diagram of Class I and Class II credibility control point data provided by an embodiment of the present invention.

[0061] Figure 3 This is a schematic diagram of determining the elevation of a flat-cut map and selecting a data reference range for a non-closed geological surface with sufficient data provided by an embodiment of the present invention.

[0062] Figure 4 It is a schematic diagram of the data tangent point projection method provided by an embodiment of the present invention.

[0063] Figure 5 This is a schematic diagram of a precise horizontal tangent line established by modifying an existing elevation horizontal tangent line provided by an embodiment of the present invention.

[0064] Figure 6 It is a schematic diagram of supplementing data by projecting adjacent elevation tangent control points when data in a partial area is insufficient as provided by an embodiment of the present invention.

[0065] Figure 7 It is a schematic diagram of a method for supplementing data by projecting tangent elevation lines close to the surface based on surface trace control points when data in a partial area is insufficient, as provided by an embodiment of the present invention.

[0066] Figure 8This is a schematic diagram of data correction provided by an embodiment of the present invention.

[0067] Figure 9 The embodiment of the present invention provides a precise planar cut map established by combining data correction and data projection supplementation.

[0068] Figure 10 It is a schematic diagram of the three-dimensional model surface provided by an embodiment of the present invention. The left one is a schematic diagram of the effect of directly fitting the model surface, and the right one is a schematic diagram of the effect of establishing the model surface by the method provided by the present invention.

[0069] Figure 11 It is a schematic diagram of the overall three-dimensional model surface provided by an embodiment of the present invention.

[0070] Figure 12 This is a schematic diagram of selecting a data reference range for a non-closed geological surface with blank data provided by an embodiment of the present invention.

[0071] Figure 13 Schematic diagram of a delineated data extension method provided in an embodiment of the present invention.

[0072] Figure 14 This is a schematic diagram of a complex three-dimensional geological model surface provided by another more complex example provided by an embodiment of the present invention. DETAILED DESCRIPTION

[0073] In order to make the purpose, technical solutions and advantages of the present invention more clearly understood, the present invention is further described in detail below in conjunction with the embodiments. It should be understood that the specific embodiments described herein are only used to explain the present invention and are not intended to limit the present invention.

[0074] In order to enable those skilled in the art to fully understand how to implement the present invention, this section provides an explanatory embodiment that expands on the technical solutions of the claims.

[0075] like Figure 1 As shown, the method for establishing a complex geological surface based on a precise planar cut map provided by an embodiment of the present invention includes the following steps:

[0076] S1: For closed geological surfaces and some non-closed geological surfaces with sufficient data, set the model node coordinates to (X, Y, Z), determine the maximum values ​​of the model node X, Y, and Z coordinates, and determine the modeling range (X min ,X max ),(Y min ,Y max ),(Z min ,Z max ), establish the range clipping box;

[0077] For non-closed geological surfaces with blank data, partition modeling is performed;

[0078] S2: Classify and import the original data, eliminate outlier data, and set the credibility level, Class I and Class II credibility;

[0079] S3: According to the sparseness of the existing flat-cut map elevation, a flat-cut map with a specific elevation can be added, and the selected elevation is used as the base surface to extend 1 / 2 of the elevation difference range to the adjacent elevation surface as the data reference range, such as Figure 3 As shown, the left figure shows the selection of the horizontal tangent elevation: based on the distribution of data points and the elevation surface of the existing horizontal tangent line, the elevation surface for establishing a precise horizontal tangent line is finally delineated; the right figure shows the selection of the data reference range of the horizontal tangent plane 2210m: the elevation difference between the horizontal tangent plane 2210m and the adjacent horizontal tangent elevation is 25m and 15m respectively, so the data reference range is extended upward and downward by 12.5m and 7.5m respectively based on the horizontal tangent plane 2210m;

[0080] S4: Projecting the data points within the data reference range onto the horizontal tangent plane to construct horizontal tangent points for each elevation data;

[0081] S5: Correct the horizontal section projection data and establish accurate horizontal tangent lines at each elevation;

[0082] S6: Fitting the three-dimensional model surface based on the precise horizontal tangent interpolation of each elevation, performing Boolean operations using the clipping box S1, and setting the model surface credibility;

[0083] S7: Determine whether reconstruction is required based on the model surface credibility. If the credibility level is lower than Class I, timely update is required based on subsequent data, and steps S2-S7 are repeated in sequence.

[0084] S8: All Class I 3D geological surfaces are gathered in the modeling area. A vertical plane is established according to the 2D section position to obtain the intersection line of the complex model surface. After the intersection line is derived, the 2D section correction is performed.

[0085] Class I credibility data is accurate data, including surface trace control points, exploration drilling control points, adit and underground cavern exploration control points, and actual excavation profile control points;

[0086] The Class II credibility data includes sub-precise data and auxiliary points. The sub-precise data are measured occurrence inference profiles on the surface or in the horizontal tunnel, and the auxiliary points are auxiliary points and lines added based on geological analysis and modeling.

[0087] The principle of projecting data points within the reference range of S4 onto the tangent plane and constructing the tangent points of various types of data is as follows. The main cases are as follows:

[0088] like Figure 4As shown in the figure, if there is no actual excavation section near the data point range, the projection points within the data range are directly projected vertically onto the plane, and the plane equation is set to Z = C, where C is the plane elevation. The required projection point is S1 (X0, Y0, Z0), and the projection vector of point S1 onto the plane Z = C is Projection point coordinates S1'(X0, Y0, C-Z0); if there is a real excavation section near the data point range, select the direction vector of the adjacent real excavation section as the projection direction, set the surface outcrop point of the selected real excavation section as D1(X1, Y1, Z1), the intersection point with the horizontal section as D2(X2, Y2, Z2), the lowest control point as D3(X3, Y3, Z3), set the vector

[0089] If X0>C, then the projection vector of point T1 onto the tangent plane Z=C is Projection point coordinates

[0090] If X0<C, then the projection vector of point T1 onto the tangent plane Z=C is Projection point coordinates

[0091] The construction of the tangent points of each elevation data described in S4 also includes the projection supplement of adjacent elevation tangent lines when data is insufficient. If there is a nearby actual excavation section to provide the projection direction, the adjacent tangent control points are supplemented with data according to the projection method in S4.

[0092] like Figure 6 As shown, if there is no actual excavation section nearby to provide the projection direction, two points Q1 and Q2 with similar X and Y coordinates in the horizontal tangent lines of different elevations are selected to form the projection vector Project the control points on the adjacent elevation horizontal tangent line to the horizontal tangent plane;

[0093] like Figure 7 As shown in Figure 5, if data in some areas are missing and there are no Class I control points to prove that the development trend is unrelated to the surface trace trend, data can also be supplemented based on the surface trace control points or the Class I control points in between according to the projection method in S4.

[0094] The correction of the projection data of the tangent map described in S5 includes deleting the control points of the type II data that are inconsistent with the type I data, such as Figure 8 As shown in the figure, because the adjacent Class I data must be strictly constrained, if this Class II data point is not deleted, a sharp corner point will appear, and there is no similar corner point in the adjacent elevation tangent line. It does not belong to a special geological structure and is therefore discarded.

[0095] like Figure 5 As shown, the establishment of the precise horizontal cutting map described in S5 includes comparing the existing elevation horizontal tangent lines, accurately constraining the Class I and Class II credibility data points, and thus establishing the precise horizontal cutting map.

[0096] The establishment of an accurate tangent map also includes adding tangent lines of new elevations. If the directions of adjacent tangent lines of elevations are clear, some control points can be projected as references to supplement the data according to the projection method in step S4.

[0097] like Figure 9 As shown in the figure, the credibility data points of types I and II are precisely constrained, the tangent points of each type of data are reasonably connected into lines, and an accurate tangent diagram is established.

[0098] The source of Class I data is related to the angle between the constructed geological surface and the ground surface. For steeply inclined geological model surfaces, such as the lithologic interface of intrusive rocks and mixed rocks, the sources of Class I data mainly include surface trace control points, horizontal shafts and underground cavern exploration control points. The actual excavation profile depends on the engineering conditions and it is difficult to have a standard.

[0099] According to the "Geotechnical Engineering Investigation Specifications", the detailed investigation stage of engineering geological mapping can be selected at a scale of 1:500 to 1:2000. The drawing accuracy of geological boundaries and geological observation points should not be less than 3mm on the map. L(A) is the mapping line density of the surface trace line. During detailed surveys, underground caverns in mountainous areas shall be arranged according to geological structures, and the distance between exploration points shall be less than 50m; for urban underground caverns, the distance between exploration points shall be less than 25m for sites with complex rock and soil changes, 25-40m for sites with moderately complex rock and soil changes, and 40-80m for sites with simple rock and soil changes.

[0100] The fitted three-dimensional model surface is as follows Figure 10 As shown in the figure, the model surface has different Class I model surface identification criteria according to different site types and complexity. For example, if the model surface n is a geological structure surface in a mountainous area, the maximum spacing between exploration points is R(A n )=50m, then the critical data surface density of model surface n is:

[0101]

[0102] Among them, the data area density P(A n ) is the ratio of the number of data points on the model surface n to the area of ​​the model surface, R(A n ) is the distance between exploration points on the model surface.

[0103] Nearly horizontal or partially gently dipping geological model surfaces, such as the lithologic interface of sedimentary rocks and the basement-overburden interface between the Quaternary system and bedrock, are mainly derived from Type I data, which are surface drilling exploration control points and surface trace control points. Usually, if such model surfaces are distorted and complex, vertical tangent lines can be established (similar to constructing precise section lines) to achieve the purpose of establishing complex model surfaces.

[0104] The layout of exploration holes is related to the site type and complexity. Different sites and complexities will give a layout spacing range of exploration lines (exploration profiles) and a layout spacing range of exploration points. Assume that the maximum layout spacing of exploration lines is C (A n ), the maximum spacing of exploration points is R(A n ), some special sites will arrange exploration holes according to corner points, and only set the layout spacing range, that is, C(A n ) does not exist. Then the critical data surface density of such a model surface n is:

[0105]

[0106] Among them, if C(A n ) does not exist, then let C(A n )=R(A n ).

[0107] The type I model surface is the mapping line density of the surface trace Model data surface density P(A n )≥P(A ncj ) and the density of type I data Model surface.

[0108] The non-closed geological surface includes lithologic interfaces, unloading surfaces, and deep and large torsional faults. If the non-closed geological surface has a large range of data blank areas, the process of establishing the non-closed geological surface is as follows:

[0109] S101, partition modeling, divide n modeling areas according to data distribution, and determine the range of each modeling area:

[0110] {(X 1min ,X 1max ),(Y 1min ,Y 1max ),(Z 1min ,Z 1max )},...,{(X nmin ,X nmax ),(Y nmin ,Y nmax ),(Z nmin ,Z nmax )}, and establish the cropping boxes for each area;

[0111] S102, setting the number of blank data areas to i, and creating i data blank area model surfaces one by one;

[0112] S103, after all the model surfaces of the data blank area are supplemented and fitted, the regional models are merged to obtain the entire model. The lowest credibility level of the regional model (excluding the data blank area) is the credibility level of the entire model, such as Figure 11As shown, the confidence level of this model surface is the lowest confidence level of the model surfaces controlled by precise data on both sides.

[0113] The embodiment of the present invention provides a method for supplementing the blank area of ​​data, including the following steps:

[0114] like Figure 12 As shown, first select part of the data in the precise data control area connected to the blank area as a reference, and the selected data reference range is equal to 1 / 2 of the data blank area range;

[0115] If the data blank area is less than 10% of the overall data range, select the Class I data within the data reference range as the priority fitting point and directly fit to fill the blank area;

[0116] If it is greater than 10%, the data extension method is used to extend and supplement the adjacent data blank areas based on the data within the data reference range, and the geological surface of the data blank areas is fitted after extension and supplementation.

[0117] The data extension method process is as follows:

[0118] like Figure 13 As shown, select the lowest point node of the boundary between the two areas and the blank, and set the vector In order to avoid contradictions in the extended data, data above 1 / 2 of the reference range is selected for extension. The M1 extension direction vector is: The M2 extension direction vector is When the blank area is much larger than 10%, the value of q can be appropriately increased.

[0119] In order to prove the creativity and technical value of the technical solution of the implementation regulations of this invention, this section provides application examples of the claimed technical solution on specific products or related technologies.

[0120] The following is the lithologic boundary surface between granite and migmatite in a dam area. Granite has a distinct regional zonal distribution. Due to different tectonic movements within the same rock belt at different geological periods, the boundary morphology is mostly irregular, such as strips, veins, dendrites, and ellipses. Compared with the examples in the specific steps, this example is more complex and cannot be accurately achieved by direct fitting. Figure 13 is a schematic diagram of the example fitting results.

[0121] The lithologic boundaries in this example have vein-like and dendritic development characteristics, and develop to a specified elevation with complex morphology, but are not exposed on the surface. Modeling requires combining geological knowledge with modeling, otherwise it is easy to form an erroneous connection with the ground surface.

[0122] The method for establishing a precise plane-cut surface to fit a complex geological model surface proposed in the present invention has several outstanding advantages in this example: ① Direct fitting cannot usually be generated directly in one step, mainly because it cannot achieve the effect of basic consistency between the data control points and the model surface. In particular, the geological surface in this case has complex data at the dendritic and vein-shaped locations, and the direct fitting has poor consistency and accuracy. Local adjustments must be made later, but local adjustments usually form irregular protrusions, which may not achieve the ideal effect and cause model errors and distortions. In comparison, the complex three-dimensional geological model surface established based on the precise plane-cut map can be generated by fitting the plane tangent line in one step. The consistency between the data control points and the model surface is high, the model accuracy is high, and the ideal effect can be achieved. ② The data points in this case are complex and intertwined. If contradictory data control points appear, direct fitting may result in erroneous pinch-out corners. In comparison, the establishment of a precise plane-cut map basically eliminates erroneous pinch-out corners because the early data correction basically eliminates similar situations. ③ The geological model surface in this case is distorted and complex, and it is necessary to combine geological knowledge to construct the model. Direct fitting will result in erroneous interpolation that violates geological knowledge. In comparison, modeling based on precise cross-sections specifies the fitting direction to a certain extent, avoiding most of the possibility of erroneous interpolation.

[0123] The modeling software used in this paper is BM_GeoModelerS2019, a geological three-dimensional modeling software with point and line drawing and surface fitting functions. It can generate simple geological model surfaces with one click based on geological control points and surface occurrence. However, the fitting effect of overly complex geological surfaces is similar to that of other software, and direct fitting will result in unreasonable irregular protrusions. The method proposed in this paper mainly uses functions such as "point projection" and "surface construction based on multiple non-closed curves". By establishing precise horizontal tangents for each elevation and controlling the fitting direction, it can eliminate unfounded irregular protrusions and erroneous supplements formed by automatic interpolation of fitting.

[0124] The method for establishing a complex geological surface based on a precise flat-cut map provided in an application embodiment of the present invention is applied to a computer device, wherein the computer device includes a memory and a processor, wherein the memory stores a computer program. When the computer program is executed by the processor, the processor executes the steps of the method for establishing a complex geological surface based on a precise flat-cut map.

[0125] The complex geological surface establishment method based on the precise planar cut map provided in the application embodiment of the present invention is applied to an information data processing terminal, and the information data processing terminal is used to implement the complex geological surface establishment system based on the precise planar cut map.

[0126] It should be noted that the embodiments of the present invention can be implemented by hardware, software, or a combination of software and hardware. The hardware portion can be implemented using dedicated logic; the software portion can be stored in a memory and executed by an appropriate instruction execution system, such as a microprocessor or dedicated design hardware. Those skilled in the art will understand that the above-mentioned devices and methods can be implemented using computer-executable instructions and / or contained in processor control code, for example, such as a carrier medium such as a disk, CD or DVD-ROM, a programmable memory such as a read-only memory (firmware), or a data carrier such as an optical or electronic signal carrier. The device and its modules of the present invention can be implemented by hardware circuits such as very large-scale integrated circuits or gate arrays, semiconductors such as logic chips, transistors, etc., or programmable hardware devices such as field programmable gate arrays, programmable logic devices, etc., can also be implemented by software executed by various types of processors, or can be implemented by a combination of the above-mentioned hardware circuits and software, such as firmware.

[0127] The embodiments of the present invention have achieved some positive results during the development or use process, and indeed have great advantages over the existing technology. The following content describes them in conjunction with data, charts, etc. from the experimental process.

[0128] like Figure 10 As shown, the left picture is a model surface established based on a precise horizontal cutting map. Combined with geological knowledge, the horizontal tangent line correction and fitting are used to make the control points and the model surface accurately match. The right picture is a model surface directly fitted based on the control points. There are several obvious problems: ① The surface control points are difficult to completely match in the distortion. ② In areas where there are few and scattered control points, the fitting has irregular protrusions, which are incorrect supplements formed by automatic interpolation. ③ There are large differences between adjacent elevation control points, and the fitting consistency is poor. Most control points are not on the surface. In summary, the method proposed in the present invention for establishing a complex three-dimensional surface model based on a precise horizontal cutting map has obvious advantages.

[0129] The above description is only a specific embodiment of the present invention, but the scope of protection of the present invention is not limited thereto. Any modifications, equivalent substitutions and improvements made by any technician familiar with this technical field within the technical scope disclosed by the present invention and within the spirit and principles of the present invention should be covered by the scope of protection of the present invention.

Claims

1. A method for establishing a complex geological surface based on a precise planar cut map, characterized in that: Project different types of control points onto the data projection of the specified horizontal section, and build accurate horizontal section maps of each elevation through the projection results of various types of data points in the early stage and data correction; Correct the deviation of various types of horizontal cutting map data by screening highly reliable exploration data, and construct highly reliable complex geological model surfaces through accurate horizontal cutting maps; By filtering the data points at the boundary and extending them to the data blank area, a model for a larger data blank area is constructed. The method for establishing a complex geological surface based on a precise planar cut map comprises the following steps: Step 1: For closed geological surfaces and some non-closed geological surfaces with sufficient data, set the model node coordinates to (X, Y, Z), set the modeling range (X min ,X max ),(Y min ,Y max ),(Z min ,Z max ), establish the range clipping box; For non-closed geological surfaces with blank data, partition modeling is performed; Step 2: Classify and import the original data, eliminate the pinch-out data, and set the credibility level, Class I credibility and Class II credibility; Step 3: Based on the sparseness of the existing horizontal cutting map elevation, or by adding a horizontal cutting map with a specific elevation, the selected elevation is used as the base surface and the 1 / 2 elevation difference range is extended to the adjacent elevation surface as the data reference range; Step 4: Project the data points within the data reference range onto the tangent plane to construct tangent points for each type of data; Step 5: Correct the cross-section projection data to create accurate cross-section lines and cross-section diagrams; Step 6: Fit the three-dimensional model surface based on the precise horizontal tangent interpolation of each elevation, use the clipping box S1 to perform Boolean operations, and set the model surface credibility; Step 7: If the credibility level is lower than Class I, it will be updated in time according to subsequent data, and steps 2 to 7 will be repeated in sequence; Step 8: All Class I 3D geological surfaces are gathered in the modeling area, and a vertical plane is established according to the 2D section position to obtain the intersection line with the complex 3D model surface, which is the corrected 2D section line. The Class I credibility data is accurate data, including surface control points, adit and borehole control points, and actual excavation profile control points; The Class II credibility data includes sub-precise data and auxiliary points. The sub-precise data are measured occurrence inference profiles on the surface or in the adit. The auxiliary points are auxiliary points and lines added based on geological analysis and modeling. The method of projecting data points within the data reference range to the plane is as follows: if there is no actual excavation section near the data point range, the projection points within the data range are directly projected vertically to the plane, and the plane equation is set to Z = C, C is the plane elevation, and the required projection point is S1 (X0, Y0, Z0). The projection vector of point S1 projected to the plane Z = C is The projection point coordinates are S1'(X0, Y0, C-Z0); if there is a real excavation section near the data point range, the direction vector of the adjacent real excavation section is selected as the projection direction. The surface outcrop point of the selected real excavation section is set as D1(X1, Y1, Z1), the intersection point with the horizontal tangent plane is set as D2(X2, Y2, Z2), and the lowest control point is set as D3(X3, Y3, Z3); Set Vector If X0>C, then the projection vector of point T1 onto the tangent plane Z=C is Projection point coordinates If X0<C, then the projection vector of point T1 onto the tangent plane Z=C is Projection point coordinates The fitting three-dimensional model surface has different Class I model surface discrimination criteria according to different site types and complexity; Since the source of Class I data is related to the angle between the geological surface and the ground, the source of Class I data for steeply inclined geological model surfaces mainly includes surface trace control points, horizontal tunnel and underground cave exploration control points. Assuming that the model surface n is a geological structural surface in a mountainous area, according to the standard, the maximum distance between exploration points is R(A n )=50m, then the critical data surface density of model surface n is: Among them, the data area density P(A n ) is the ratio of the number of data points on the model surface n to the area of ​​the model surface, S(A n ) is the area of ​​the model surface n, R(A n ) is the maximum spacing of n exploration points on the model surface, P(A n ) and the model surface area S(A n ) is inversely proportional to; For nearly horizontal or partially gently dipping geological model surfaces, the main source of Class I data is surface drilling exploration control points and surface trace control points. If such model surfaces are distorted and complex, vertical tangent lines can be established to achieve the purpose of building complex model surfaces. The layout of exploration holes is related to the site type and complexity. Different sites and complexities will give a layout spacing range of exploration lines and a layout spacing range of exploration points. Assume that the maximum layout spacing of exploration lines is C(A n ), the maximum spacing of exploration points is R(A n ), some special sites will arrange exploration holes according to corner points, and only set the layout spacing range, that is, C(A n ) does not exist; then the critical data surface density of such model surface n is: Among them, if C(A n ) does not exist, then let C(A n )=R(A n ); The type I model surface is the mapping line density of the surface trace / m, model surface data density P(A n )≥P(A ncj ) and the density of type I data Model surface.

2. The method for establishing a complex geological surface based on a precise planar cutting diagram according to claim 1, wherein: Combined with the knowledge of the surrounding geology, a specific method of projecting different data control points onto the designated horizontal tangent plane is proposed, so that each elevation horizontal tangent line constructed later has certain characteristics and similarities with the adjacent elevation horizontal tangent lines. According to the projection results of the data points on each plane, data correction is performed by retaining the exploration data with higher credibility and discarding the low credibility data with obvious errors; The fitting direction is specified to a certain extent by fitting the complex model surface through the precise horizontal cutting diagram of each elevation.

3. The method for establishing a complex geological surface based on a precise planar cutting diagram according to claim 1, wherein: The establishment of the precise horizontal tangent map includes comparing the existing elevation horizontal tangent lines, accurately constraining the Class I and Class II credibility data points, and thus establishing the precise horizontal tangent map; The establishment of the precise tangent map also includes, for the newly added tangent lines of the elevation, if the directions of the adjacent tangent lines of the elevation are clear, using the projections of some control points as references to supplement the data; If there is no actual excavation section near the data gap to provide the projection direction, two points Q1 and Q2 with similar X and Y coordinates in the horizontal tangent lines at different elevations are selected to form the projection vector If data in some areas are missing, data supplementation shall be carried out based on the surface trace control points or the Class I control points therebetween; Finally, the tangent points of each type of data are reasonably connected into lines, the control points of type II data that contradict type I data are deleted, the credibility data points of types I and II are accurately constrained, and an accurate tangent map is established.

4. The method for establishing a complex geological surface based on a precise planar cutting diagram according to claim 1, wherein: The non-closed geological surface includes lithologic interfaces, unloading surfaces, and deep and large torsional faults. If the non-closed geological surface has a large range of data blank areas, the process of establishing the non-closed geological surface is as follows: S101, partition modeling, set the model node coordinates to (X, Y, Z), and divide n modeling areas according to data distribution: {(X 1min ,X 1max ),(Y 1min ,Y 1max ),(Z 1min ,Z 1max )},...,{(X nmin ,X nmax ),(Y nmin ,Y nmax ),(Z nmin ,Z nmax )}, and establish the cropping boxes for each area; S102, setting the number of blank data areas to i, and creating i data blank area model surfaces one by one; S103, after all the model surfaces of the data blank areas are supplemented and fitted, the regional models are merged to obtain the entire model, and the lowest credibility level of the regional models except the data blank areas is the credibility level of the entire model.

5. The method for establishing a complex geological surface based on a precise planar cutting diagram according to claim 4, characterized in that: The process of supplementing the data blank area includes the following steps: First, select part of the data in the precise data control area connected to the blank area as the reference basis, select a data reference range equal to 1 / 2 of the data blank area range, and set the data reference range control data point set of each area as M1 and M2. Some special geological control points can be discarded. If the blank area of ​​the data is less than 10% of the overall data range, the point set of type I data in M1 and M2 is selected as the priority fitting point, and the blank area is directly fitted; If it is greater than 10%, then the data extension method is used to extend and supplement the adjacent data blank areas based on the point set data in M1 and M2, and the geological surface of the data blank areas is fitted after the extension and supplement; The data extension method process is as follows: Select the lowest point node of the boundary between the two areas and the blank, and set the vector In order to avoid contradictions in the extended data, data above 1 / 2 of the reference range is selected for extension. The extension direction vector of the data in point set M1 is: The extension direction vector of the data in point set M2 is When the blank area is much larger than 10%, the value of q can be appropriately increased.

6. A system for establishing a complex geological surface based on a precise planar cut map, which implements the method according to any one of claims 1 to 5, characterized in that: The complex geological surface establishment system based on the precise planar cutting map includes: The data classification module sets the credibility level of data according to its source classification, which is used to determine the priority of fitting between data and provide a basis for subsequent data correction; The data point projection module projects the control points within the data reference range, adjacent geological control points or surface trace control points to the specified horizontal section based on a specific method, providing a basis for the subsequent construction of accurate horizontal section maps; The module for establishing precise tangent maps includes data correction and precise tangent line construction. Data correction is used to discard conflicting low-reliability data. Precise tangent line construction is used to precisely and reasonably connect corrected Class I and Class II credibility data into lines to establish a precise tangent map, which is the core foundation of the entire system. A model surface establishment module is used to establish a complex three-dimensional model surface using the precise planar cut map, including overall modeling and partition modeling; the partition modeling further includes partition modeling and supplementary establishment of data blank area models, and the complex three-dimensional model surface is established based on the data blank area partitions, and reference data points at the boundaries are selected to directly fit or extend the blank area data to fit the model surface; The model surface credibility setting module is used to set the credibility level of the precise and complex three-dimensional model surface based on the data credibility classification results.

7. An information data processing terminal, characterized in that: The information data processing terminal is used to implement the complex geological surface establishment system based on the precise cutting map as described in claim 6.

Images