Geological profile map automatic generation method and device based on drilling data and storage medium
Through the automatic generation method based on drilling data, the problems of low efficiency and low accuracy of traditional geological profile maps are solved, and the automatic generation of high-precision geological profile maps is realized, which is suitable for linear engineering geological analysis.
Patent Information
- Application Number
- CN202411969539.8
- Authority / Receiving Office
- CN · China
- Patent Type
- Applications(China)
- Current Assignee / Owner
- Filing Date
- 2024-12-30
- Publication Date
- 2025-05-30
AI Technical Summary
The traditional geological profile map generation method relies on manual operations, is inefficient and prone to errors. Especially in linear engineering, due to the positional offset between the drilling hole and the profile line, the geological profile map is inaccurately drawn.
Automatic generation method based on drilling data is adopted, including collecting and preprocessing drilling data, spatial interpolation is performed through Kriging interpolation, virtual drilling points are generated, and the formations between adjacent drilling holes are connected using the geotechnical layer connection relationship inference algorithm. Smooth curve fitting is used to optimize the curve shape of the geological layer, and finally a high-precision geological profile map is automatically generated.
It realizes automatic generation of accurate geological profiles, reduces manual intervention, improves work efficiency, can more accurately reflect the actual engineering geological profile, and is suitable for large-scale linear engineering geological analysis.
Smart Images

Figure CN120070613A_ABST
Abstract
Description
Technical Field
[0001] The present invention relates to the fields of geological exploration and engineering design, and particularly to an automatic generation method, device, and storage medium for geological profiles based on borehole data. Background Art
[0002] During the construction of linear projects such as highways and railways, geological exploration is a very important link. As basic geological data, borehole data usually needs to be used to generate geological profiles to describe the distribution of underground rock and soil layers. However, most traditional methods for generating geological profiles rely on manual operations, which are not only inefficient but also prone to errors due to human factors. With the increase in geological exploration data, the automatic generation of accurate geological profiles has become a requirement for improving work efficiency and ensuring data accuracy.
[0003] A geological profile is obtained by taking a profile line on the plane of the project area, based on the engineering geological exploration data on this profile line, dividing the profile strata, drawing stratification lines, and marking relevant geological elements and symbols. The drawing process of traditional geological profiles is as follows: multiple boreholes are sequentially selected according to the plane profile line, the connections of the same strata between each borehole are inferred based on the borehole stratum distribution information, and then the divided areas are filled with different filling patterns to obtain the engineering geological profile division on the two-dimensional map.
[0004] Geological profiles of linear projects such as highways and railways are divided into two categories. One is the longitudinal profile arranged along the center line of the line, and the other is the transverse profile perpendicular to the line at the specified mileage stake. Linear projects generally arrange boreholes staggered along the center line of the line. Therefore, the boreholes do not coincide with the profile line, that is, there is a positional offset on the plane between the boreholes and the profile line. This leads to situations such as inconsistent elevation of the borehole opening and the ground line, and inconsistent stratum information in the borehole data associated with the profile line, and cannot truly reflect the actual engineering geological profile. In addition, due to the wide range of linear projects, large amounts of borehole data, and updated exploration information, geological profiles need to be continuously corrected and redrawn, which brings a huge workload to engineering personnel. Therefore, how to automatically and intelligently generate accurate geological profiles has become an urgent problem to be solved in the current field of geological exploration. Summary of the Invention
[0005] The purpose of the present invention is to address the problem that when drilling currently, the boreholes do not coincide with the profile line, resulting in inaccurate drawing of geological profiles. In view of this deficiency, an automatic generation method, device, and storage medium for geological profiles based on borehole data are proposed.
[0006] To achieve the above objectives, the present invention adopts the following technical solutions:
[0007] A method, device, and storage medium for automatically generating a geological section diagram based on drilling data, comprising the following specific steps:
[0008] S1. Collect drilling data and preprocess the drilling data;
[0009] S2. Perform spatial interpolation on the depth data of rock and soil layers based on the mathematical model of Kriging interpolation;
[0010] S3. Predict and generate virtual drilling points along the section line based on the spatial interpolation;
[0011] S4. Infer the connection of rock and soil layers between adjacent drill holes based on the rock and soil layer connection relationship inference algorithm, and optimize the curve shape of the geological layer on the section diagram by using smooth curve fitting for the rock and soil layer lines;
[0012] S5. Automatically generate a high-precision geological section diagram based on the drilling data and the rock and soil layer lines.
[0013] As a further preference of the present invention, the drilling data includes: drilling coordinates, drilling depth, rock and soil layer types, and depths of rock and soil layers with different attributes.
[0014] As a further preference of the present invention, the S2 includes:
[0015] S21. Statistically analyze the depth value data of the same type of rock and soil layer based on the drilling data to form a depth data set of the same type of rock and soil layer {Z(x 1 ), Z(x 2 ), …, Z(x n )}, where x i is the spatial coordinate of the rock and soil layer under the drill hole, and Z(x i ) is the variable value corresponding to the depth value of the rock and soil layer;
[0016] S22. Calculate the variogram value γ(x i , x j ) based on the depth data set of the same type of rock and soil layer; the variogram describes the spatial autocorrelation between data points, that is, as the spatial distance increases, the difference between variable values.
[0017] Among them, Z(x i ) and Z(x j ) are any two points in the depth data set of the same type of rock and soil layer, and Z(x i ) and Z(x j ) form a data point pair, k is the record of the data point pair, and N is the number of data point pairs;
[0018] S23. According to
[0019]
[0020] Among them, γ(x i , x j ) is a known data point x i and x j The variance function value between γ(x 0 , x i ) is the point to be estimated x 0 With the known point x i The variance function value between them is used to calculate the weight λ of the data point 1 ,λ 2 ,λ 3 …λ n and the Lagrange multiplier μ.
[0021] S24, according to Among them, λ i is the weight of the i-th data point, and the specified spatial position x is calculated. 0 The depth of the rock layer.
[0022] As a further preferred embodiment of the present invention, S3 includes:
[0023] S31, based on the method of projecting the actual drilling point onto an arbitrary curve, calculate the minimum distance d(t) between the actual drilling point and the cutting curve and the curve parameter t, that is, the plane coordinate position of the virtual drilling point;
[0024] S32. According to r(t)=(x(t), y(t)), obtain the virtual drilling plane coordinate point P'(x', y').
[0025] As a further preferred embodiment of the present invention, S4 includes:
[0026] S41, adding a virtual continuous rock and soil layer with a thickness of 0 at the top and bottom of each borehole respectively;
[0027] S42, selecting two adjacent boreholes, connecting the upper and lower boundary points of the rock and soil layers with the same attribute in the two adjacent boreholes, and determining the distribution mode of each rock and soil layer according to geological rules;
[0028] The distribution modes of rock and soil layers include: continuous distribution and discontinuous distribution;
[0029] S43, determining whether the distribution mode of the rock and soil layers is continuous distribution;
[0030] S44, if the distribution is continuous, connecting the depth position points of the rock and soil layers with the same attribute between two adjacent boreholes;
[0031] S45. If there are multiple rock-soil layer connection lines with the same attribute at the depth position point of the rock-soil layer with the same attribute, the rock-soil layer connection line with the smallest slope is selected as the connection line of the continuous rock-soil layer;
[0032] S46. If there are no multiple connecting lines of the same property rock and soil layer at the depth position point of the same property rock and soil layer, it is determined that there is a discontinuous distribution of the rock and soil layer;
[0033] S47. Determine the connection relationship of the rock and soil layer according to the sorting relationship of the same property rock and soil layer between two adjacent boreholes,
[0034] The connection relationships of the discontinuous distribution rock and soil layer include pinch-out, lens body, missing, and discontinuous missing;
[0035] S48. Connect the lines according to the connection relationship of the discontinuous distribution rock and soil layer, and jump to S45;
[0036] S49. Connect all adjacent two boreholes in sequence to form the rock and soil layer line on the cross-section diagram.
[0037] As a further preference of the present invention, the S47 includes:
[0038] Select any one of two adjacent boreholes, and determine whether there is a continuously distributed rock and soil layer between this borehole and other adjacent boreholes;
[0039] If it exists, connect the lines;
[0040] When the m-th layer of one borehole and the (i + 1)-th layer and the (i - 1)-th layer of another borehole are of the same property rock and soil layer, and the i-th layer of another borehole is of another property rock and soil layer, it is a pinch-out connection. Connect the lower critical line of the (i + 1)-th layer with the lower critical line of the m-th layer, connect the upper critical line of the (i - 1)-th layer with the upper critical line of the m-th layer, and connect the upper critical line and the lower critical line of the i-th layer with the pinch-out point respectively;
[0041] When the m-th layer of one borehole and the (i - 1)-th layer of another borehole are of the same property rock and soil layer, the (m + 1)-th layer of one borehole and the (i + 1)-th layer of another borehole are of the same property rock and soil layer, and the property of the i-th layer of another borehole is different from both the (i - 1)-th layer and the (i + 1)-th layer, it is a missing connection. Then connect the upper critical line and the lower critical line of the same property rock and soil layer respectively;
[0042] When the (m - 1)-th layer of one borehole and the (i - 1)-th layer of another borehole are of the same property rock and soil layer, the (m + 1)-th layer of one borehole and the (i + 1)-th layer of another borehole are of the same property rock and soil layer, the property of the i-th layer of another borehole, the property of the m-th layer of one borehole are different from both the (i - 1)-th layer and the (i + 1)-th layer, and the property of the i-th layer is also different from the property of the m-th layer, it is an indirect missing connection. Then connect the lower critical line of the i-th layer with the upper critical line of the m-th layer, and connect the upper critical line and the lower critical line of the same property rock and soil layer respectively.
[0043] As a further preference of the present invention, the pinch-out point is selected according to the pinch-out rule, and the pinch-out rule includes:
[0044] If the thickness of the pinched-out rock and soil layer is less than 2m, this layer can be ignored, and the pinching-out position is determined at 1 / 10d away from the borehole; the pinched-out rock and soil layer is a missing rock and soil layer existing in adjacent boreholes at the current depth;
[0045] If the thickness of the pinched-out rock and soil layer is 2 - 5m, the pinching-out position is 1 / 3d away from the borehole;
[0046] If the thickness of the pinched-out rock and soil layer is 5 - 8m, the pinching-out position is 1 / 2d away from the borehole;
[0047] If the thickness of the pinched-out rock and soil layer is greater than 8m, the pinching-out position is 2 / 3d away from the borehole;
[0048] Wherein, d is the distance between adjacent boreholes.
[0049] As a further preference of the present invention, the S5 includes:
[0050] S51. Using ObjectARX technology to generate borehole custom entities, rock and soil layer line entities, and geological section custom entities;
[0051] S52. Optimizing the color, layer thickness, and line style of the profile diagram;
[0052] S53. Automatically generating a DWG file, which contains the profile diagram, borehole information, and rock and soil layer description.
[0053] A device for automatically generating a geological profile diagram based on borehole data, the device includes one or more processors and a memory, and the memory is used to store one or more programs; when the one or more programs are executed by the processor, the processor can implement a method for automatically generating a geological profile diagram based on borehole data.
[0054] A computer storage medium, in which at least one program instruction is stored, and the at least one program instruction is used to be loaded and executed by a processor to implement a method for automatically generating a geological profile diagram based on borehole data.
[0055] A method, device, and storage medium for automatically generating a geological profile diagram based on borehole data proposed by the present invention, compared with the prior art, have the following beneficial effects:
[0056] 1. The present invention can effectively convert borehole data into accurate exploration line profile diagrams, and the whole process has a high degree of automation, reducing manual intervention and improving work efficiency;
[0057] 2. The present invention uses spatial interpolation to predict virtual boreholes, addresses the positional offset between actual boreholes and profile lines on a plane, and is more in line with the actual profile situation;
[0058] 3. The self-identification algorithm for formation connection relationships is adopted to connect formations between adjacent boreholes, and a smooth curve fitting is used to optimize the curve shape of the formation on the geological layer profile diagram, realizing the automatic connection of layers between adjacent boreholes;
[0059] 4. The present invention allows users to flexibly adjust parameters such as drawing styles, scales, and colors according to specific project requirements, meeting the needs of different linear engineering projects such as highways and railways;
[0060] 5. The present invention can process a large amount of borehole data and is applicable to geological analysis of large-scale linear projects (such as highways and railways). BRIEF DESCRIPTION OF THE DRAWINGS
[0061] Figure 1 is a schematic flow diagram of the present invention;
[0062] Figure 2 is a schematic diagram of the connection relationships of 3 types of discontinuous rock and soil layers. DETAILED DESCRIPTION OF THE EMBODIMENTS
[0063] The present invention will be specifically introduced below in conjunction with the accompanying drawings and specific embodiments.
[0064] The present invention combines a mathematical model and computer-aided design (CAD) technology to realize the automatic generation of geological profile diagrams based on borehole data, mainly for geological analysis and design of linear projects such as highways and railways.
[0065] Embodiment 1: In combination with Figure 1 - Figure 2 , a method, device, and storage medium for automatically generating a geological profile diagram based on borehole data include the following specific steps:
[0066] S1. Collect borehole data and preprocess the borehole data.
[0067] The borehole data includes: borehole coordinates, borehole depth, rock and soil layer types, and depths of different types of rock and soil layers with different attributes.
[0068] S2. Perform spatial interpolation on the rock and soil layer depth data based on the mathematical model of Kriging interpolation.
[0069] S21. Assume that there are n known data points of the same type of rock and soil layer under all known boreholes, and statistically analyze the depth value data of the same type of rock and soil layer based on the borehole data to form a depth data set of the same type of rock and soil layer {Z(x 1 ), Z(x 2 ), …, Z(x n )}, where x iis the spatial coordinate of the rock and soil layer under the borehole, Z(x i ) is the variable value corresponding to the depth value of the rock and soil layer.
[0070] S22, based on the same type of rock and soil layer depth data set, calculate the variance function value γ(x i ,x j ); the variogram describes the spatial autocorrelation between data points, that is, the difference between variable values as the spatial distance increases.
[0071] Among them, Z(x i ) and Z(x j ) are any two points of the same type of soil layer depth data set, and Z(x i ) and Z(x j ) form data point pairs, k is the data point pair record, and N is the number of data point pairs.
[0072] S23. According to
[0073]
[0074] Among them, γ(x i , x j ) is a known data point x i and x j The variance function value between γ(x 0 , x i ) is the point to be estimated x 0 With the known point x i The variance function value between them is used to calculate the weight λ of the data point 1 ,λ 2 ,λ 3 …λ n and the Lagrange multiplier μ.
[0075] S24, according to Among them, λ i is the weight of the i-th data point, and the specified spatial position x is calculated. 0 The depth of the rock layer.
[0076] S3. Predict and generate virtual drilling points along the section line based on spatial interpolation.
[0077] The cutting line is any continuous plane curve that satisfies the center line of the linear engineering line. The virtual drilling referred to in this method is the projection of the actual drilling position to the corresponding position of the cutting line.
[0078] S31. Based on the method of projecting the actual drilling point onto an arbitrary curve, the minimum distance d(t) between the actual drilling point and the cutting curve and the curve parameter t, that is, the plane coordinate position of the virtual drilling point, are calculated.
[0079] S32. According to r(t)=(x(t), y(t)), obtain the virtual drilling plane coordinate point P'(x', y').
[0080] The plane coordinate position of the virtual drilling point is solved by projecting the point onto an arbitrary curve. The minimum distance between the actual drilling point and the cutting curve is solved. The cutting curve is composed of three types of curve units: several polylines, arcs, and straight lines. During the solution process, the cutting curve is split into the smallest curve units for solution.
[0081] Let the curve unit be C, defined by the parametric equation r(t) = (x(t), y(t)), where t is the curve parameter. The actual drilling plane coordinate point P(x 0 ,y 0 ) to any point (x(t), y(t)) on the curve unit:
[0082]
[0083] Minimize d(t), that is, solve the following optimization problem:
[0084] min t ((x(t)-x 0 ) 2 +(y(t)-y 0 ) 2 ),
[0085] By taking the derivative of d(t) with respect to t and setting the derivative to zero, we can get the parameter t* corresponding to the minimum distance. Substituting the parameter t* into the curve equation r(t), we can get the coordinate point P′(x′,y′) of the virtual drilling plane.
[0086] S4. The rock and soil layers between adjacent boreholes are connected based on the rock and soil layer connection relationship inference algorithm, and the curve shape of the geological layer is optimized by smooth curve fitting to form the rock and soil layer line on the profile.
[0087] S41. Add a virtual continuous rock and soil layer with a thickness of 0 at the top and bottom of each borehole.
[0088] S42. Select two adjacent boreholes, connect the upper and lower boundary points of the rock and soil layers with the same properties in the two adjacent boreholes, and determine the distribution of each rock and soil layer according to geological rules.
[0089] Starting from the starting hole A of the section line, find the next adjacent hole B of the hole.
[0090] The method of judging the distribution pattern of rock and soil layers according to geological rules is used to determine the distribution pattern of each rock and soil layer in boreholes A and B.
[0091] The distribution patterns of rock and soil layers are divided into continuous distribution and discontinuous distribution. The discontinuous distribution is further divided into four cases: pinch-out, lens, missing, and intermittent missing. The phenomenon where the rock and soil layer is continuous in the middle borehole and missing at both ends is regarded as a lens.
[0092] Continuous distribution: The rock and soil layer is continuously distributed in adjacent boreholes. Connect the upper and lower demarcation points of the rock and soil layer respectively, and the connecting line is the upper and lower demarcation line of the rock and soil layer.
[0093] Missing connection: A certain rock and soil layer in a borehole is missing at the corresponding layer in its adjacent borehole. The rock and soil layer pinches out at the adjacent borehole. Connect the upper and lower demarcation points of the rock and soil layer and the pinch-out point respectively.
[0094] Intermittent missing connection: The rock and soil layers in two adjacent boreholes are both discontinuous, and the missing rock and soil layers in the two boreholes are not the same rock and soil layer, then it is regarded as intermittent missing.
[0095] Pinch-out connection: The lithologies of the upper and lower rock and soil layers in a borehole are the same and continuous in adjacent boreholes, while the rock and soil layer is missing in the adjacent borehole. This situation is regarded as the pinch-out connection of the rock and soil layer. It is necessary to reasonably determine its pinch-out point, and connect the upper and lower demarcation points of the rock and soil layer and the pinch-out point respectively.
[0096] S43. Judge whether the distribution pattern of the rock and soil layer is continuous distribution.
[0097] The judgment principle of continuous rock and soil layer is that during the deposition process of the rock and soil layer, it basically remains horizontal and the rock and soil layers do not intersect. The gentle connection line of the rock and soil layer can be used to determine the continuously distributed rock and soil layer.
[0098] S44. If it is continuous distribution, connect the depth position points of the rock and soil layers with the same attribute between two adjacent boreholes. The depth position points of the rock and soil layers with the same attribute are the bottom ends of the rock and soil layers with the same attribute.
[0099] S45. If there are multiple connection lines of the rock and soil layers with the same attribute at this depth position point of the rock and soil layer with the same attribute, select the connection line of the rock and soil layer with the smallest slope as the connection line of the continuous rock and soil layer.
[0100] S46. If there are no multiple connection lines of the rock and soil layers with the same attribute at this depth position point of the rock and soil layer with the same attribute, then judge that there are discontinuous distributed rock and soil layers.
[0101] S47. Judge the connection relationship of the rock and soil layer according to the sorting relationship of the rock and soil layers with the same attribute in two adjacent boreholes.
[0102] Select any one of the two adjacent boreholes and judge whether there are continuously distributed rock and soil layers in this borehole and other adjacent boreholes;
[0103] If there are, then make a connection;
[0104] When the m-th layer of one borehole and the (i + 1)-th and (i - 1)-th layers of another borehole are of the same property rock and soil layer, and the i-th layer of another borehole is of another property rock and soil layer, it is a pinch-out connection. Connect the lower critical line of the (i + 1)-th layer with the lower critical line of the m-th layer, connect the upper critical line of the (i - 1)-th layer with the upper critical line of the m-th layer, and connect the upper critical line and the lower critical line of the i-th layer with the pinch-out point respectively.
[0105] The pinch-out point is selected according to the pinch-out rule, and the pinch-out rule includes:
[0106] If the thickness of the pinched-out rock and soil layer is less than 2m, then this layer can be ignored, and the pinch-out position is determined at 1 / 10d away from this borehole; the pinched-out rock and soil layer is a missing rock and soil layer existing in adjacent boreholes at the current depth;
[0107] If the thickness of the pinched-out rock and soil layer is 2 - 5m, then the pinch-out position is 1 / 3d away from this borehole;
[0108] If the thickness of the pinched-out rock and soil layer is 5 - 8m, then the pinch-out position is 1 / 2d away from this borehole;
[0109] If the thickness of the pinched-out rock and soil layer is greater than 8m, then the pinch-out position is 2 / 3d away from this borehole;
[0110] Where d is the distance between adjacent boreholes.
[0111] When the m-th layer of one borehole and the (i - 1)-th layer of another borehole are of the same property rock and soil layer, the (m + 1)-th layer of one borehole and the (i + 1)-th layer of another borehole are of the same property rock and soil layer, and the property of the i-th layer of another borehole is different from both the (i - 1)-th layer and the (i + 1)-th layer, it is a missing connection, then connect the upper critical line and the lower critical line of the rock and soil layers with the same property respectively.
[0112] When the (m - 1)-th layer of one borehole and the (i - 1)-th layer of another borehole are of the same property rock and soil layer, the (m + 1)-th layer of one borehole and the (i + 1)-th layer of another borehole are of the same property rock and soil layer, the property of the i-th layer of another borehole, the property of the m-th layer of one borehole are different from both the (i - 1)-th layer and the (i + 1)-th layer, and the property of the i-th layer is also different from the property of the m-th layer, it is an indirect missing connection, then connect the lower critical line of the i-th layer with the upper critical line of the m-th layer, and connect the upper critical line and the lower critical line of the rock and soil layers with the same property respectively.
[0113] S48. Make connections according to the connection relationship of the discontinuous distribution rock and soil layers, and skip to S45.
[0114] S49. Make connections for all adjacent two boreholes in sequence to form the rock and soil layer line on the cross-section diagram.
[0115] S5. Automatically generate a high-precision geological cross-section diagram based on the borehole data and the rock and soil layer line.
[0116] S51. Generate custom drilling entities, rock and soil layer line entities, and custom geological section entities using ObjectARX technology.
[0117] S52. Optimize the color, layer thickness, and line style of the cross-section diagram.
[0118] S53. Automatically generate a DWG file that contains the cross-section diagram, drilling information, and rock and soil layer descriptions.
[0119] Embodiment 2: A device for automatically generating a geological cross-section diagram based on drilling data, the device includes one or more processors and a memory, and the memory is used to store one or more programs; when the one or more programs are executed by the processor, the processor can implement a method for automatically generating a geological cross-section diagram based on drilling data.
[0120] Embodiment 3: A computer storage medium, in which at least one program instruction is stored, and the at least one program instruction is used to be loaded and executed by a processor to implement a method for automatically generating a geological cross-section diagram based on drilling data.
[0121] The foregoing has shown and described the basic principles, main features, and advantages of the present invention. Those skilled in the art should understand that the above embodiments do not limit the present invention in any form. Any technical solutions obtained by using equivalent replacements or equivalent transformations fall within the protection scope of the present invention.
Claims
1. A method for automatically generating a geological profile based on drilling data, characterized in that: The specific steps include: S1, collecting drilling data and preprocessing the drilling data; S2, spatial interpolation of rock and soil layer depth data based on the mathematical model of Kriging interpolation; S3, predicting and generating virtual drilling points along the section line based on spatial interpolation; S4, connecting the rock and soil layers between adjacent boreholes based on the rock and soil layer connection relationship inference algorithm, and optimizing the rock and soil layer line on the curve shape profile of the geological layer by smooth curve fitting; S5. Automatically generate high-precision geological profiles based on drilling data and rock and soil layer lines.
2. The method for automatically generating a geological profile based on drilling data according to claim 1, characterized in that: The drilling data include: drilling coordinates, drilling depth, rock and soil layer type, and depth of rock and soil layers of different attributes.
3. The method for automatically generating a geological profile based on drilling data according to claim 2, characterized in that: The S2 includes: S21, based on the drilling data, the depth value data of the same type of rock and soil layer are counted to form a depth data set of the same type of rock and soil layer {Z(x1), Z(x2),…, Z(x n )}, where x i is the spatial coordinate of the rock and soil layer under the borehole, Z(x i ) is the variable value corresponding to the depth value of the rock and soil layer; S22, based on the same type of rock and soil layer depth data set, calculate the variance function value γ(x i ,x j ); Among them, Z(x i ) and Z(x j ) are any two points of the same type of soil layer depth data set, and Z(x i ) and Z(x j ) form data point pairs, k is the data point pair record, and N is the number of data point pairs; S23. According to Among them, γ(x i , x j ) is a known data point x i and x j The variance function value between γ(x0, x i ) is the distance between the estimated point x0 and the known point x i The variance function value between them is used to calculate the weights of the data points λ1,λ2,λ3…λ n and the Lagrange multiplier μ. S24, according to Among them, λ i is the weight of the i-th data point, and the depth of the rock and soil layer at the specified spatial position x0 is calculated.
4. The method, device and storage medium for automatically generating a geological profile based on drilling data according to claim 2, characterized in that: The S3 includes: S31, based on the method of projecting the actual drilling point onto an arbitrary curve, calculate the minimum distance d(t) between the actual drilling point and the cutting curve and the curve parameter t, that is, the plane coordinate position of the virtual drilling point; S32. According to r(t)=(x(t), y(t)), obtain the virtual drilling plane coordinate point P'(x', y').
5. The method for automatically generating a geological profile based on drilling data according to claim 2, characterized in that: The S4 includes: S41, adding a virtual continuous rock and soil layer with a thickness of 0 at the top and bottom of each borehole respectively; S42, selecting two adjacent boreholes, connecting the upper and lower boundary points of the rock and soil layers with the same attribute in the two adjacent boreholes, and determining the distribution mode of each rock and soil layer according to geological rules; The distribution modes of rock and soil layers include: continuous distribution and discontinuous distribution; S43, determining whether the distribution mode of the rock and soil layers is continuous distribution; S44, if the distribution is continuous, connecting the depth position points of the rock and soil layers with the same attribute between two adjacent boreholes; S45. If there are multiple rock-soil layer connection lines with the same attribute at the depth position point of the rock-soil layer with the same attribute, the rock-soil layer connection line with the smallest slope is selected as the connection line of the continuous rock-soil layer; S46. If there are no multiple connecting lines of rock and soil layers with the same attribute at the depth position point of the rock and soil layer with the same attribute, it is determined that a non-continuously distributed rock and soil layer exists; S47, judging the connection relationship of rock and soil layers according to the order relationship of rock and soil layers with the same attributes in two adjacent boreholes, The connection relationship of non-continuously distributed rock and soil layers includes pinch-out, lens, missing and discontinuous missing; S48, connect the lines according to the connection relationship of the non-continuously distributed rock and soil layers, and skip S45; S49. Connect all adjacent boreholes in sequence to form rock and soil layer lines on the cross-section diagram.
6. The method for automatically generating a geological profile based on drilling data according to claim 5, characterized in that: The S47 includes: Select one of the two adjacent boreholes to determine whether there is a continuously distributed rock and soil layer between the borehole and the other adjacent boreholes; If it exists, connect it; When the mth layer of a borehole and the i+1th layer and i-1th layer of another borehole are rock and soil layers of the same attribute, and the ith layer of the other borehole is a rock and soil layer of another attribute, it is a pinch-out connection, connecting the lower critical line of the i+1th layer with the lower critical line of the mth layer, connecting the upper critical line of the i-1th layer with the upper critical line of the mth layer, and connecting the upper critical line and lower critical line of the ith layer to the pinch-out point respectively; When the mth layer of a borehole and the i-1th layer of another borehole are the same attribute rock and soil layers, the m+1th layer of a borehole and the i+1th layer of another borehole are the same attribute rock and soil layers, and the attribute of the i-th layer of the other borehole is different from that of the i-1th layer and the i+1th layer, it is a missing connection, then the upper critical line and the lower critical line of the rock and soil layers with the same attribute are connected respectively; When the m-1th layer of a borehole and the i-1th layer of another borehole are rock and soil layers with the same attributes, the m+1th layer of a borehole and the i+1th layer of another borehole are rock and soil layers with the same attributes, the attributes of the i-th layer of the other borehole and the m-th layer of a borehole are different from the i-1th layer and the i+1th layer, and the attributes of the i-th layer are also different from the attributes of the m-th layer, it is an indirect missing connection, then the lower critical line of the i-th layer is connected with the upper critical line of the m-th layer, and the upper critical line and lower critical line of the rock and soil layers with the same attributes are connected respectively.
7. The method for automatically generating a geological profile based on drilling data according to claim 6, characterized in that: The pinch-out point is selected according to the pinch-out rule, and the pinch-out rule includes: If the thickness of the pinched-out rock and soil layer is less than 2m, the layer can be ignored, and the pinch-out position is determined by the distance of 1 / 10d from the borehole; the pinched-out rock and soil layer is a missing rock and soil layer in the adjacent borehole at the current depth; If the thickness of the rock and soil layer where the pinch-out occurs is 2 to 5 m, the pinch-out position is 1 / 3 d away from the borehole; If the thickness of the rock and soil layer where the pinch-out occurs is 5 to 8 m, the pinch-out position is 1 / 2 d away from the borehole; If the thickness of the rock and soil layer where the pinch-out occurs is greater than 8m, the pinch-out position is 2 / 3d away from the borehole; Where d is the distance between adjacent boreholes.
8. The method for automatically generating a geological profile based on drilling data according to claim 2, characterized in that: The S5 includes: S51. Use ObjectARX technology to generate custom entities for drilling holes, rock and soil layer lines, and geological profiles; S52, optimizing the color, layer thickness and line style of the cross-section; S53, the DWG file generated automatically, the DWG file contains the profile, drilling information, and rock and soil layer description.
9. A device for automatically generating geological profiles based on drilling data, characterized in that: The device includes one or more processors and a memory, wherein the memory is used to store one or more programs; when the one or more programs are executed by the processor, the processor is enabled to implement a method for automatically generating a geological profile based on drilling data as described in any one of claims 1-8.
10. A computer storage medium, characterized in that: The computer storage medium stores at least one program instruction, and the at least one program instruction is used to be loaded and executed by the processor to implement the method for automatically generating a geological profile based on drilling data as described in any one of claims 1-8.