Entity three-dimensional geologic model creation and curve profile generation method

Abstract

The invention relates to a solid three-dimensional geologic model creation and curve profile generation method, which comprises the following steps: in the process of creating a solid three-dimensional geologic model, defining an overall framework and resolution, importing and converting geological drilling data, defining the sequence of formation in the geological drilling data, visualizing the imported geological drilling data, and generating a curve profile of the solid three-dimensional geologic model. And performing interpolation calculation by using an interpolation method, generating a plane of each stratum, filling the whole geologic body, and generating a three-dimensional geologic model. In the step of generating the curve section, a spline curve is defined, slice setting is carried out, a plane equation is established, planes are calculated at discrete points of the curve step by step, intersection calculation is carried out on the three-dimensional geologic model, the three-dimensional model is cut, a normal vector is calculated for each point on the spline curve, the planes are defined, cutting operation is carried out on the model for multiple times, and then the planes are extracted. And obtaining the section of the spline curve in the three-dimensional geologic model. The method optimizes the creation of dredging engineering and geological engineering and the generation of a geological curve profile on the motion trail of the reamer.

Description

Technical Field

[0001] The invention relates to the technical field of dredging engineering or geological engineering, and in particular to a method for creating a solid three-dimensional geological model and generating a curve section. Background Art

[0002] Dredging projects play a vital role in many fields such as waterway regulation, port construction, and water conservancy projects. When carrying out dredging projects, accurately understanding the geological conditions of the construction area has a key impact on the reasonable planning of construction plans, estimating the amount of work, ensuring construction safety, and improving project efficiency. However, the traditional way of presenting geological data is often two-dimensional drawings or simple text descriptions of geological survey reports, which are difficult to intuitively and comprehensively display complex three-dimensional geological structures, bringing many inconveniences to dredging project-related decisions. Therefore, there is an urgent need for an effective method that can create a solid three-dimensional geological model and generate curved profiles to assist in the better implementation of dredging projects. Summary of the invention

[0003] The technical problem to be solved by the present invention is to overcome the deficiencies in the prior art and provide a method for creating a solid three-dimensional geological model and generating a curved profile.

[0004] The present invention is achieved through the following technical solutions:

[0005] A method for creating a solid three-dimensional geological model and generating a curved profile includes creating a solid three-dimensional geological model and generating a curved profile, wherein the "creating a solid three-dimensional geological model" includes the following steps:

[0006] S11. Define the overall framework and resolution, and determine the regions and voxels of the entire solid 3D geological model;

[0007] S12. Import and convert geological drilling data to obtain two CSV format files as basic data for the three-dimensional geological model;

[0008] S13. defining the order in which the strata in the geological drilling data appear, generating a relation array, and assigning a different color to each stratum;

[0009] S14. Visualizing the imported geological drilling data;

[0010] S15. Generate a two-dimensional plane, and use an interpolation method to perform interpolation calculation to generate a plane for each layer of the stratum;

[0011] S16. Fill the entire geological body, fill volume elements between the geological body strata, and generate a three-dimensional geological model;

[0012] The “generating a curve profile” comprises the following steps:

[0013] S21. Determine the specific shape of the curve through a series of control points and define a spline curve;

[0014] S22. Perform slice setting and establish plane equation;

[0015] S23. gradually calculating planes at discrete points of the curve, performing intersection calculations with the three-dimensional geological model, and cutting the three-dimensional model;

[0016] S24. Calculate the normal vector for each point on the spline curve and define a plane, perform multiple plane cutting operations on the model, and then extract the plane to obtain the section of the spline curve in the three-dimensional geological model.

[0017] According to the above technical solution, preferably, in step S11, the frame is defined using an X, Y, and Z range, wherein X represents the length of the frame, Y represents the width of the frame, and Z represents the depth of the frame; the resolution is defined as a*b*c, wherein a represents the number of voxels in the X direction, b represents the number of voxels in the Y direction, and c represents the number of voxels in the Z direction, and the product of the three represents the total number of voxels in the frame.

[0018] According to the above technical solution, preferably, in step S12, the data structure of the geological drilling data requires two CSVs, one including X, Y, Z and corresponding formation attributes, and the other including X, Y, Z and corresponding azimuth, gradient components, poles, and formation attributes.

[0019] According to the above technical solution, preferably, in step S14, it is determined whether the distribution of the boreholes and the distribution of the strata are consistent with the geological survey data, and whether the arrangement order of the strata is consistent with the setting.

[0020] According to the above technical solution, preferably, in step S15, the interpolation method used is spline interpolation:

[0021]

[0022] Where, j = 1, 2, ..., N, N is the number of points involved in the calculation; λ j are the coefficients obtained by solving the linear equations; r j It is the distance between point (x, y) and the jth point; after the points are encrypted, the Delaunay triangulation method is used to create a two-dimensional plane model.

[0023] According to the above technical solution, preferably, step S16 includes:

[0024] The layers are traversed and connected by interpolation points to form a tetrahedron, and the entire model space is interpolated to generate an initial geological structure model;

[0025] Expand discrete data points into continuous geological interfaces and strata distribution, and use Bayesian reasoning to build a continuous three-dimensional geological model;

[0026] Based on geological observation data, the probability distribution of model parameters is inferred and updated to generate multiple geological models;

[0027] Deeply analyze the potential variability of geological systems and the credibility of predictions, and ultimately generate a continuous three-dimensional geological model.

[0028] According to the above technical solution, preferably, in step S21, the cubic spline interpolation method is used to determine the specific shape of the curve, and at the same time, the curve is encrypted by parameter adjustment to form a smooth curve.

[0029] According to the above technical solution, preferably, in step S22, each slice plane is defined by a plane normal vector and a cutting plane equation, and the plane normal vector is usually taken as the tangent direction of the spline curve at that point.

[0030]

[0031] Among them, the tangent direction is calculated by the first-order derivative of the spline curve,

[0032] Cutting plane equation:

[0033]

[0034] Expands to:

[0035]

[0036] in, is a point on a plane. is the plane normal vector.

[0037] According to the above technical solution, preferably, in step S23, the vertex V is inserted into the plane equation:

[0038]

[0039] Among them, d>0 means the vertex is in the direction of the normal vector, d<0 means the vertex is in the opposite direction, and d=0 means the vertex is on the plane.

[0040] Compute the intersection point by linear interpolation:

[0041]

[0042] in,

[0043]

[0044] d1 and d2 are vertices and The distance relative to the plane,

[0045] A cross section is generated at the same time, and the set of intersection points defines the geometry of the cross section, represented as a polygon.

[0046] The beneficial effects of the present invention are:

[0047] The present invention is applicable to the fields of dredging engineering, ports and waterways, rivers, lakes and seas, and can realize the creation of dredging engineering and geological engineering and the generation of curved profiles of geology on the motion trajectory of the auger, thereby optimizing the processing process of dynamic sectioning of three-dimensional geological models. The whole process is fully automated, has high accuracy and small errors, and can facilitate the relevant departments to rationally plan construction plans, thus having high application and promotion value. BRIEF DESCRIPTION OF THE DRAWINGS

[0048] Figure 1 It is a schematic diagram of the process of the method for creating a solid three-dimensional geological model in the present invention.

[0049] Figure 2 It is a schematic flow chart of the curve profile generating method in the present invention. DETAILED DESCRIPTION

[0050] In order to enable those skilled in the art to better understand the technical solution of the present invention, the present invention is further described in detail below in conjunction with the accompanying drawings and the best embodiment. Based on the embodiments in the invention, all other embodiments obtained by ordinary technicians in the field without making creative work are within the scope of protection of the invention.

[0051] Embodiment 1: The present invention includes creating a solid three-dimensional geological model and generating a curved profile, wherein:

[0052] (1) Create a solid 3D geological model:

[0053] S11. Define the frame and resolution. The frame is defined using the range of X, Y, and Z. X represents the length of the frame, Y represents the width of the frame, and Z represents the depth of the frame. The resolution is a*b*c, where a represents the number of voxels in the X direction, b represents the number of voxels in the Y direction, and c represents the number of voxels in the Z direction. The product of the three represents the total number of voxels in the frame. The more voxels, the higher the resolution.

[0054] S12. Import geological drilling data, including drilling distribution information and stratigraphic location information. The data structure requires two CSVs. One is X, Y, Z and the corresponding stratigraphic attributes, which represent the coordinate points of all layers extracted from the drilling holes, as well as the corresponding stratigraphic attributes. The other is X, Y, Z and the corresponding azimuth, gradient component, pole, and stratigraphic attributes. The azimuth indicates the horizontal inclination direction of the stratigraphic layer, the gradient component indicates the degree of inclination of the stratigraphic layer, and the pole indicates the direction and inclination angle of the stratigraphic vertex.

[0055] S13. Define the stratigraphic order, define the order in which the stratigraphic layers appear in the drilling data, generate a relational array, and assign a different color to each stratigraphic layer.

[0056] S14. Visualize the data, visualize the input data, and observe whether the distribution of the drill holes and the distribution of the strata are consistent with the geological survey data, and whether the arrangement order of the strata is consistent with the settings.

[0057] S15. Generate a two-dimensional plane, and use the interpolation method to perform interpolation calculation, so that the drilling data is interpolated to generate a plane of each layer. The interpolation method used in this example is spline interpolation:

[0058]

[0059] Where, j = 1, 2, ..., N, N is the number of points involved in the calculation; λ j are the coefficients obtained by solving the linear equations; r j It is the distance between point (x, y) and the jth point; after the points are encrypted, the Delaunay triangulation method is used to create a two-dimensional plane model.

[0060] S16. Fill the entire geological body, and use the Delaunay triangulation method to close the sides between the faces to form a closed body. Fill the tetrahedrons between the layers, and use Bayesian reasoning to establish a continuous three-dimensional geological model. Based on the geological observation data, infer and update the probability distribution of the model parameters, generate multiple possible geological models, deeply analyze the potential variability of the geological system and the credibility of the prediction, and finally generate a continuous three-dimensional geological model.

[0061] (2) Generate a curve profile:

[0062] S21. Define a spline curve, which uses a series of control points to determine the specific shape of the curve. In this example, the cubic spline interpolation method is used to determine the shape of the curve, and the curve points are encrypted by parameter adjustment to form a smooth curve.

[0063] S22. Perform slicing settings. Each slice is composed of the intersection of the cutting plane and the model. The normal vector of the cutting plane is usually consistent with the tangent direction of the path. Each slicing plane is defined by the following two parts (a point on the plane is a discrete point on the spline curve).

[0064] The plane normal vector is usually taken as the tangent direction of the spline curve at that point.

[0065]

[0066] Among them, the tangent direction is calculated by the first-order derivative of the spline curve,

[0067] Cutting plane equation:

[0068]

[0069] Expands to:

[0070]

[0071] in, is a point on a plane. is the plane normal vector.

[0072] S23. Calculate the intersection and gradually calculate the plane at the discrete points of the curve. The slices can be used as independent geometric bodies for subsequent operations. Cutting the three-dimensional model is equivalent to calculating the intersection of the model surface and the above plane equation.

[0073] For each unit of the model (such as a triangle), determine whether the vertex is on both sides of the plane:

[0074] Insert vertex V into the plane equation:

[0075]

[0076] Among them, d>0 means the vertex is in the direction of the normal vector, d<0 means the vertex is in the opposite direction, and d=0 means the vertex is on the plane.

[0077] If the two vertices of the triangle are located on both sides of the plane, the intersection point is calculated by linear interpolation:

[0078]

[0079] in,

[0080]

[0081] d1 and d2 are vertices and The distance relative to the plane,

[0082] A cross section is generated at the same time, and the set of intersection points defines the geometry of the cross section, represented as a polygon.

[0083] S24. For each point on the spline curve Calculating the normal vector And define the plane, perform multiple plane cutting operations on the model, and obtain a series of slices along the spline curve. The slicing results contain the cross-section of the model on each plane. They are extracted without changing the original stratigraphic order and properties, which is convenient for analyzing properties such as cross-sectional shape or area.

[0084] Embodiment 2: As attached Figure 1 As shown, the present application provides a method for creating a solid three-dimensional geological model, which specifically includes: S1. defining the overall framework and resolution so as to determine the areas and volume elements of the entire solid three-dimensional geological model; S2. using geological data for conversion to obtain the required format CSV; S3. defining the stratigraphic sequence so that the stratigraphic order can be determined in the correct way; S4. visualizing data, after a series of preliminary data are imported, visualizing the drilling data; S5. generating a two-dimensional plane, interpolating the drilling data to generate a plane for each layer of stratigraphic elements; S6. filling the entire geological body, and filling the volume elements between the stratigraphic layers of the geological body.

[0085] Based on the above embodiment, preferably, the overall framework in S1 is [0,1000,0,1000,0,600], which means that the range of the geological body is 0 to 1000 in the X direction, 0 to 1000 in the Y direction, and 0 to 600 in the Z direction, and the resolution is [25,25,20], indicating that the number of voxels within this range is 12500. After the framework is completed, the S2 operation is performed to import the geological drilling data. Its data structure requires two CSVs, one is X, Y, Z and "stratum attribute" fields, which represent the coordinate points of the layers and layers extracted from all the boreholes, as well as the corresponding stratum attributes; the other is X, Y, Z and "azimuth", "gradient component", "pole", "stratum attribute" is used to grasp the overall direction and inclination angle of the stratum, and these two CSV format files are the basic data for forming a three-dimensional geological model.

[0086] Table 1 Drilling distribution information

[0087] Field Name Data Types X float Y float Z float Geological properties Text

[0088] Table 2 Stratigraphic location information

[0089]

[0090]

[0091] Then define the order of the strata, which are sand, sandy soil, clay, and loess. Define the order of the strata, based on the input of various geological location information and drilling information, and visualize them. After confirming that the data is correct, start generating interpolation, using the spline interpolation algorithm based on geological constraints, based on the input of various geological data points,

[0092]

[0093] Where, j = 1, 2, ..., N, N is the number of points involved in the calculation; λ j are the coefficients obtained by solving the linear equations; r j It is the distance between point (x, y) and the jth point; after the points are encrypted, the Delaunay triangulation method is used to create a two-dimensional plane model to generate a two-dimensional plane for each layer of stratum, so that the surface model has been formed.

[0094] The layers are traversed and connected by interpolation points to form a tetrahedron, and the entire model space is interpolated to generate an initial geological structure model. The discrete data points are expanded into continuous geological interfaces and stratigraphic distributions, and Bayesian reasoning is used to establish a continuous three-dimensional geological model. Based on the geological observation data, the probability distribution of the model parameters is inferred and updated to generate multiple possible geological models. The potential variability of the geological system and the credibility of the prediction are deeply analyzed, and finally a continuous three-dimensional geological model is generated.

[0095] As attached Figure 2 As shown, the present application provides a method for generating a curve profile. In one embodiment, the spline interpolation is defined as a semicircle, and its radius is defined as 25. The height of Z, that is, the plane where the arc is located is a plane with Z of 28. The number of points is defined as 100, that is, the semicircle is composed of 100 points. The plane equation in the embodiment is as follows:

[0096]

[0097] Where (i=0,1,2,3...99, N=100),

[0098] These points can also be adjusted by yourself, but more points means more iterations and more resources consumed.

[0099] After defining the basic parameters of the curve, we establish the equation of the plane. The corresponding plane equation is:

[0100]

[0101] in, is a point on the plane (i.e., a curve point ), is the plane normal vector.

[0102] After defining the plane equation, perform intersection calculation with the generated geological body to determine the geological conditions of the point, and then iterate from each point. In this example, the number of points on the curve is 100. The model is cut step by step, and finally the plane is extracted to obtain the profile of the curve in the geological model of this example.

[0103] The present invention is applicable to the fields of dredging engineering, ports and waterways, rivers, lakes and seas, and can realize the creation of dredging engineering and geological engineering and the generation of curved profiles of geology on the motion trajectory of the auger, thereby optimizing the processing process of dynamic sectioning of three-dimensional geological models. The whole process is fully automated, has high accuracy and small errors, and can facilitate the relevant departments to rationally plan construction plans, thus having high application and promotion value.

[0104] The above is only a preferred embodiment of the present invention. It should be pointed out that for ordinary technicians in this technical field, several improvements and modifications can be made without departing from the principle of the present invention. These improvements and modifications should also be regarded as the scope of protection of the present invention.

Claims

1. A method for creating a solid three-dimensional geological model and generating a curved profile, characterized in that: The method includes creating a solid three-dimensional geological model and generating a curve profile, wherein the "creating a solid three-dimensional geological model" includes the following steps: S11. Define the overall framework and resolution, and determine the regions and voxels of the entire solid 3D geological model; S12. Import and convert geological drilling data to obtain two CSV format files as basic data for the three-dimensional geological model; S13. defining the order in which the strata in the geological drilling data appear, generating a relation array, and assigning a different color to each stratum; S14. Visualizing the imported geological drilling data; S15. Generate a two-dimensional plane, and use an interpolation method to perform interpolation calculation to generate a plane for each layer of the stratum; S16. Fill the entire geological body, fill volume elements between the geological body strata, and generate a three-dimensional geological model; The "generating a curve profile" comprises the following steps: S21. Determine the specific shape of the curve through a series of control points and define a spline curve; S22. Perform slice setting and establish plane equation; S23. gradually calculating planes at discrete points of the curve, performing intersection calculations with the three-dimensional geological model, and cutting the three-dimensional model; S24. Calculate the normal vector for each point on the spline curve and define a plane, perform multiple plane cutting operations on the model, and then extract the plane to obtain the section of the spline curve in the three-dimensional geological model.

2. A method for creating a solid three-dimensional geological model and generating a curved profile according to claim 1, characterized in that: In step S11, the frame is defined using an X, Y, and Z range, where X represents the length of the frame, Y represents the width of the frame, and Z represents the depth of the frame; The resolution is defined as a*b*c, where a represents the number of voxels in the X direction, b represents the number of voxels in the Y direction, and c represents the number of voxels in the Z direction. The product of the three represents the total number of voxels in the frame.

3. A method for creating a solid three-dimensional geological model and generating a curved profile according to claim 2, characterized in that: In step S12, the data structure of the geological drilling data requires two CSVs, one including X, Y, Z and corresponding formation attributes, and the other including X, Y, Z and corresponding azimuth, gradient component, pole, and formation attributes.

4. A method for creating a solid three-dimensional geological model and generating a curved profile according to claim 3, characterized in that: In step S14, it is determined whether the distribution of the boreholes and the distribution of the strata are consistent with the geological survey data, and whether the arrangement order of the strata is consistent with the setting.

5. A method for creating a solid three-dimensional geological model and generating a curved profile according to claim 3, characterized in that: In step S15, the interpolation method used is spline interpolation: Where, j = 1, 2, ..., N, N is the number of points involved in the calculation; λ j are the coefficients obtained by solving the linear equations; r j It is the distance between point (x, y) and the jth point; after the points are encrypted, the Delaunay triangulation method is used to create a two-dimensional plane model.

6. A method for creating a solid three-dimensional geological model and generating a curved profile according to claim 5, characterized in that: Step S16 includes: The layers are traversed and connected by interpolation points to form a tetrahedron, and the entire model space is interpolated to generate an initial geological structure model; Expand discrete data points into continuous geological interfaces and strata distribution, and use Bayesian reasoning to build a continuous three-dimensional geological model; Based on geological observation data, the probability distribution of model parameters is inferred and updated to generate multiple geological models; Deeply analyze the potential variability of geological systems and the credibility of predictions, and ultimately generate a continuous three-dimensional geological model.

7. The method for creating a solid three-dimensional geological model and generating a curved profile according to any one of claims 1 to 6, characterized in that: In step S21, the cubic spline interpolation method is used to determine the specific shape of the curve, and the points of the curve are encrypted by parameter adjustment to form a smooth curve.

8. A method for creating a solid three-dimensional geological model and generating a curved profile according to claim 7, characterized in that: In step S22, each slice plane is defined by a plane normal vector and a cutting plane equation. The plane normal vector is usually taken as the tangent direction of the spline curve at that point. Among them, the tangent direction is calculated by the first-order derivative of the spline curve, Cutting plane equation: Expands to: in, is a point on a plane. is the plane normal vector.

9. A method for creating a solid three-dimensional geological model and generating a curved profile according to claim 8, characterized in that: In step S23, the vertex V is inserted into the plane equation: Among them, d>0 means the vertex is in the direction of the normal vector, d<0 means the vertex is in the opposite direction, and d=0 means the vertex is on the plane. Compute the intersection point by linear interpolation: in, d1 and d2 are vertices and The distance relative to the plane, A cross section is generated at the same time, and the set of intersection points defines the geometry of the cross section, represented as a polygon.

Images