A 3D modeling method and system for geotechnical engineering

Abstract

The present invention relates to the field of geotechnical engineering modeling, and particularly to a three-dimensional modeling method and system for geotechnical engineering; determining verification points and data points according to the verification priorities of each drilling point; determining the setting mode of the search area corresponding to the point to be interpolated according to the geological complexity conditions of the geographical reference area where the point to be interpolated is located; under the first analysis condition, determining the adjustment mode of the search area according to the data point difference state of each search sub-area; under the second analysis condition, determining the data point acquisition condition according to the number of data points and the data point distribution uniformity in the search area, and determining whether to adjust the search area range according to the data point acquisition condition; performing interpolation data processing on the point to be interpolated according to the drilling information of the data points in the search area of the point to be interpolated, and performing three-dimensional modeling when the interpolation data processing of all points to be interpolated is completed; the present invention improves the reliability of the generated three-dimensional geological model.

Description

Technical Field

[0001] The present invention relates to the field of geotechnical engineering modeling, and particularly to a three-dimensional modeling method and system for geotechnical engineering. Background Art

[0002] Establishing a three-dimensional geological model can more accurately represent complex geological phenomena and effectively improve its application value. The establishment of a three-dimensional geological model largely depends on borehole data. Therefore, insufficient borehole depth and borehole density result in low modeling credibility. However, the cost of obtaining borehole data is high. Therefore, it is beneficial to use a spatial interpolation algorithm to construct virtual boreholes based on limited data. However, in current three-dimensional modeling methods and systems, the search range for interpolated data processing cannot be optimized according to the complexity of the actual formation, resulting in the failure to consider the influence of geological spatial anisotropy during interpolated data processing, resulting in unqualified accuracy of the obtained interpolation results, and further resulting in low accuracy of the generated three-dimensional geological model.

[0003] Chinese Patent Publication No. CN112419500B discloses a three-dimensional geological model modeling method, which obtains exploration data of multiple exploration holes, extracts borehole data, performs data interpolation based on the borehole data to obtain interpolated data, integrates the interpolated data and the borehole data into modeling data, generates point cloud data through the modeling data, and then generates a formation surface, generates a formation entity through the formation surface, extracts soil layer data from the exploration data, uses the soil layer data as project parameters of the formation entity, and generates a three-dimensional geological model through the formation entity. The formation surface generated by the above solution is smoother and completely controlled by parameters, and has a good effect on the simulation of the formation. However, the above solution has the following problems: the search range for interpolated data processing cannot be optimized according to the complexity of the actual formation, resulting in the failure to consider the influence of geological spatial anisotropy during interpolated data processing, resulting in unqualified accuracy of the obtained interpolation results, and further resulting in the accuracy of the generated three-dimensional geological model not meeting the requirements of practical applications. Summary of the Invention

[0004] Therefore, the present invention provides a three-dimensional modeling method and system for geotechnical engineering to overcome the problem that the search range for interpolated data processing cannot be optimized according to the complexity of the actual formation, resulting in the failure to consider the influence of geological spatial anisotropy during interpolated data processing, resulting in unqualified accuracy of the obtained interpolation results, and further resulting in the accuracy of the generated three-dimensional geological model not meeting the requirements of practical applications.

[0005] To achieve the above object, the present invention provides a three-dimensional modeling method for geotechnical engineering, including:

[0006] Determine the verification priority of each borehole point according to the reference characteristic value of each borehole point, and determine the verification points and data points according to the verification priority of each borehole point;

[0007] Upload the points to be interpolated, determine the geological complex conditions based on the thickness stability reference value and dip angle stability reference value of the geographical reference area where the points to be interpolated are located, and determine the setting method of the search area corresponding to the points to be interpolated according to the geological complex conditions;

[0008] Under the first analysis condition, determine the data point difference state of each search sub-area according to the data point characteristic value, and determine the adjustment method of the search area according to the data point difference state of each search sub-area;

[0009] The adjustment method of the search area includes adjusting the horizontal width of the search sub-area according to the data feature difference degree, or adjusting the reference length of the edge nodes corresponding to the search sub-area in the search area, or adjusting the length of the search sub-area according to the data feature difference degree of the search sub-area;

[0010] Under the second analysis condition, determine the data point acquisition condition according to the number of data points and the data point distribution uniformity in the search area, and determine whether to adjust the search area range according to the data point acquisition condition;

[0011] Perform interpolation data processing on the points to be interpolated according to the borehole information of the data points in the search area of the points to be interpolated, and perform 3D modeling when the interpolation data processing of all the points to be interpolated is completed;

[0012] Record the borehole points except the verification points among all the borehole points as data points.

[0013] Furthermore, the setting method of the search area is determined according to the geological complex conditions, including:

[0014] If the geological complex condition is a type I geological complex condition, set the search area as a circular area composed of several search sub-areas, the center of the circular area is the point to be interpolated, the initial areas of each search sub-area are the same, and the number of search sub-areas is determined according to the thickness stability reference value or dip angle stability reference value;

[0015] If the geological complex condition is a type II geological complex condition, set the search area as a circular area, the center of the circular area is the point to be interpolated, and the area of the circular area is the preset initial area;

[0016] The number of the search sub-areas is positively correlated with the thickness stability reference value or dip angle stability reference value.

[0017] Furthermore, the geological complex conditions are determined according to the thickness stability reference value and dip angle stability reference value, and the geological complex conditions include:

[0018] A type of complex geological conditions where the thickness stability reference value is greater than the preset thickness stability reference value or the dip angle stability reference value is greater than the preset dip angle stability reference value;

[0019] A type of complex geological conditions where the thickness stability reference value is less than or equal to the preset thickness stability reference value or the dip angle stability reference value is less than or equal to the preset dip angle stability reference value.

[0020] Furthermore, under the first analysis condition, the characteristic values of each data point in each search sub-region are detected, the data point difference state of each search sub-region is determined according to the data point characteristic values, and the adjustment method of the search region is determined according to the data point difference state of each search sub-region. The data point difference state includes that the data characteristic difference degree of the search sub-region is less than the preset characteristic difference degree, or the data characteristic difference degree within the search sub-region is greater than or equal to the preset characteristic difference degree;

[0021] The first analysis condition is that the geological complex condition is a type of complex geological condition and the selection of the setting method of the search region is completed.

[0022] Furthermore, if the data characteristic difference degree within the search sub-region is greater than or equal to the preset characteristic difference degree, it is determined by detecting the distribution condition of the differential data points to increase the horizontal width of the search sub-region according to the data characteristic difference degree, or to increase the reference length of the edge node corresponding to the search sub-region in the search region according to the data characteristic difference degree;

[0023] The increased value of the horizontal width is positively correlated with the data characteristic difference degree;

[0024] The increased value of the reference length of the edge node is positively correlated with the data characteristic difference degree.

[0025] Furthermore, the distribution condition of the differential data points is determined according to the distribution of the differential data points. The distribution condition of the differential data points includes that the differential data points are all distributed on the edge of the sub-region, or the distance between the differential data points and the point to be interpolated is greater than the preset differential distance;

[0026] If the differential data points are all distributed in the edge area of the sub-region, the horizontal width of the search sub-region is increased according to the data characteristic difference degree;

[0027] If the distance between the differential data points and the point to be interpolated is greater than the preset differential distance, the reference length of the edge node corresponding to the search sub-region in the search region is increased according to the data characteristic difference degree.

[0028] Furthermore, if the data characteristic difference degree of the search sub-region is less than the preset characteristic difference degree, the length of the search sub-region is increased according to the data characteristic difference degree of the search sub-region;

[0029] The increased value of the length of the search sub-region is negatively correlated with the degree of data feature difference.

[0030] Further, under the second analysis condition, the data point acquisition condition is determined according to the number of data points and the distribution uniformity of data points in the search region. If the data point acquisition condition is a first type of acquisition condition, the radius length of the search region is increased and adjusted according to the number of data points or the distribution uniformity of data points.

[0031] The increased value of the radius length of the search region is negatively correlated with the number of data points or the distribution uniformity of data points.

[0032] The second analysis condition is that the geological complexity condition is a second type of geological complexity condition and the selection of the setting method of the search region is completed.

[0033] The data point acquisition conditions include a first type of acquisition condition and a second type of acquisition condition.

[0034] The first type of acquisition condition is that the number of data points is less than the preset number of data points, or the distribution uniformity of data points is less than the preset distribution uniformity.

[0035] The second type of acquisition condition is that the number of data points is greater than or equal to the preset number of data points and the distribution uniformity of data points is greater than or equal to the preset distribution uniformity.

[0036] Further, the verification priority of the borehole points is determined according to the reference characteristic values of each borehole point, and the borehole points with the preset verification quantity are selected as verification points according to the verification priority rule, and the other borehole points are recorded as data points.

[0037] The reference characteristic value is determined according to whether the data measured at the borehole point is complete and whether there are anomalies. The reference characteristic value c = ω W ×W + ω G ×G, where W is the data integrity of the borehole point, G is the data difference degree of the borehole point, and ω W and ω G are weight coefficients.

[0038] The increased value of the verification priority is positively correlated with the reference characteristic value of the borehole point.

[0039] The present invention also provides a system for a three-dimensional modeling method for geotechnical engineering, including:

[0040] A data display unit for uploading, displaying and storing the borehole information of each borehole point, where the borehole information includes the formation thickness and formation dip angle of each formation in the borehole point, the compression modulus and the foundation bearing capacity obtained during the borehole process.

[0041] A point selection unit, which is connected to the data display unit, is used to determine the verification priority of borehole points according to the reference characteristic values of each borehole point, and determine verification points and data points according to the verification priorities of each borehole point;

[0042] A data analysis unit, which is connected to the data display unit and the point selection unit, is used to determine the geological complexity conditions according to the thickness stability reference value and dip angle stability reference value of the geographical reference area where the interpolation point is located, and determine the setting method of the search area corresponding to the interpolation point according to the geological complexity conditions;

[0043] A first analysis unit, which is connected to the data analysis unit, is used to determine the data point difference state of each search sub-area according to the data point characteristic values, and determine the adjustment method of the search area according to the data point difference state of each search sub-area;

[0044] A second analysis unit, which is connected to the data analysis unit, is used to determine the data point acquisition conditions according to the number of data points and the data point distribution uniformity in the search area, and determine whether to adjust the search area range according to the data point acquisition conditions;

[0045] A processing unit, which is connected to the data analysis unit, the first analysis unit and the second analysis unit, performs interpolation data processing on the interpolation point according to the borehole information of the data points in the search area of the interpolation point, and performs three-dimensional modeling when the interpolation data processing of all interpolation points is completed.

[0046] Compared with the prior art, the beneficial effect of the present invention is that the technical solution of the present invention determines the setting method of the search area corresponding to the interpolation point according to the actual geological complexity, and after the setting method of the search area is selected, adjusts the search area, so that the interpolation data processing process fully takes into account the influence of geological spatial anisotropy, making the finally obtained interpolation result more accurate, thereby improving the reliability of the three-dimensional geological model generated by the present invention.

[0047] Furthermore, in the present invention, the geological complexity conditions of the area corresponding to each interpolation point are determined according to the thickness stability reference value and dip angle stability reference value, and different setting methods of the search area are selected for each interpolation point, which not only ensures the accuracy of the interpolation processing process for areas with high geological complexity, but also ensures the processing efficiency of the interpolation processing process for areas with low geological complexity.

[0048] Furthermore, in the present invention, the adjustment method for the corresponding search area is determined according to the data feature difference degree in each search sub-area, and the adjustment method for the search area is determined according to the distribution condition of the difference data points, so that the search areas for each interpolation point and the ranges of each search sub-area are more in line with the actual processing process, making the interpolation results of each interpolation point more accurate, and thus improving the reliability of the three-dimensional geological model generated by the present invention.

[0049] Furthermore, in the present invention, the verification priority of the borehole points is determined according to the data integrity and data abnormality degree of the borehole points, and the verification points are selected according to the verification priority, ensuring that the selected verification points are more representative, making the evaluation of the interpolation results for each interpolation point more accurate, and thus making the modeling data more accurate, and improving the reliability of the three-dimensional geological model generated by the present invention. Brief Description of the Drawings

[0050] Figure 1 It is a schematic diagram of the three-dimensional modeling method for geotechnical engineering of the present invention;

[0051] Figure 2 It is a schematic diagram of the search area corresponding to a type of complex geological condition of the present invention;

[0052] Figure 3 It is a flowchart of the method for determining the adjustment method of the search area according to the data point difference state of each search sub-area of the present invention;

[0053] Figure 4 It is a unit connection diagram of the three-dimensional modeling system for geotechnical engineering of the present invention;

[0054] In the figure: search area 1, search sub-area 2, edge node 3, short side 4 of the search sub-area, long side 5 of the search sub-area. Detailed Embodiments

[0055] In order to make the purpose and advantages of the present invention more clear and understandable, the present invention will be further described below in conjunction with embodiments; it should be understood that the specific embodiments described herein are only used to explain the present invention and are not used to limit the present invention.

[0056] The preferred embodiments of the present invention will be described below with reference to the accompanying drawings. Those skilled in the art should understand that these embodiments are only used to explain the technical principles of the present invention and do not limit the protection scope of the present invention.

[0057] It should be noted that in the description of the present invention, the terms indicating the direction or positional relationship such as "upper", "lower", "left", "right", "inner", "outer", etc. are based on the direction or positional relationship shown in the drawings. This is only for the convenience of description and does not indicate or imply that the device or element must have a specific orientation, be constructed and operated in a specific orientation. Therefore, it should not be construed as a limitation to the present invention.

[0058] In addition, it should be noted that in the description of the present invention, unless otherwise clearly specified and defined, the terms "installation", "connection", and "coupling" should be understood in a broad sense. For example, it can be a fixed connection, a detachable connection, or an integral connection; it can be a mechanical connection or an electrical connection; it can be directly connected or indirectly connected through an intermediate medium, and it can be the communication inside two elements. For those skilled in the art, the specific meanings of the above terms in the present invention can be understood according to specific circumstances.

[0059] Please refer to Figure 1 as shown in the figure, which is a schematic diagram of the three-dimensional modeling method for geotechnical engineering of the present invention. The present invention provides a three-dimensional modeling method for geotechnical engineering, including:

[0060] Determine the verification priority for each borehole point according to the reference characteristic values of each borehole point, and determine the verification points and data points according to the verification priorities of each borehole point;

[0061] Upload the points to be interpolated, determine the geological complexity conditions according to the thickness stability reference value and dip angle stability reference value of the geographical reference area where the points to be interpolated are located, and determine the setting method of the search area corresponding to the points to be interpolated according to the geological complexity conditions;

[0062] Under the first analysis condition, determine the data point difference state of each search sub-area according to the data point characteristic values, and determine the adjustment method of the search area according to the data point difference state of each search sub-area;

[0063] The adjustment method of the search area includes adjusting the lateral width of the search sub-area according to the data characteristic difference degree, or adjusting the reference length of the edge nodes corresponding to the search sub-areas in the search area, or adjusting the length of the search sub-area according to the data characteristic difference degree of the search sub-area;

[0064] Under the second analysis condition, determine the data point acquisition conditions according to the number of data points and the data point distribution uniformity in the search area, and determine whether to adjust the search area range according to the data point acquisition conditions;

[0065] Perform interpolation data processing for the points to be interpolated according to the borehole information of the data points in the search area of the points to be interpolated, and perform three-dimensional modeling when the interpolation data processing of all the points to be interpolated is completed.

[0066] Mark the borehole points among all the borehole points except the verification points as data points.

[0067] The present invention is used for three-dimensional geological modeling. An interpolation algorithm is used to perform interpolation processing on the points to be interpolated according to the data of known borehole points. The target modeling area is the entire range of the area where three-dimensional geological modeling needs to be carried out. The user uploads the positions of the points to be interpolated. It can be understood that if there are sub-areas with sparse and scattered borehole points, it is necessary to turn the finite and discrete data in this sub-area into uniform and continuous data. By performing interpolation processing on the points to be interpolated, the risk of low modeling credibility caused by insufficient borehole depth and borehole density can be reduced. When the interpolation data processing of each point to be interpolated in the present invention is completed, three-dimensional modeling is carried out. Three-dimensional modeling is content that is easily understood by those skilled in the art and will not be elaborated here. A method for three-dimensional modeling is provided, including verifying the interpolation results for each point to be interpolated, and recording the verified interpolation results and the known borehole point data as a modeling database. Import the modeling database into three-dimensional modeling software, create a three-dimensional grid and perform numerical settings on the grids at corresponding positions. After the preliminary construction of the three-dimensional geological modeling is completed, adjust the details of the modeling results, adjust the smoothness of the grid, and the adjusted modeling result is the finally obtained three-dimensional geological model.

[0068] Please refer to Figure 2 As shown, it is a schematic diagram of the search area corresponding to a type of geological complex condition of the present invention. The initial search area corresponding to the point to be interpolated is a circular area 1, and the initial shapes of the search sub-areas 2 are all rectangles. The search sub-areas 2 are centered on the point to be interpolated and arranged at equal angles around the base point. The edge nodes 3 are the intersections of the edges of each search sub-area and the search area. Each search sub-area 2 includes two long sides 5 and two short sides 4.

[0069] Please refer to Figures 1 to 3 As shown;

[0070] Determine the setting method of the search area according to the geological complex conditions, including:

[0071] If the geological complex condition is a type of geological complex condition, set the search area as a circular area composed of several search sub-areas. The center of the circular area is the point to be interpolated. The initial areas of each search sub-area are the same, and the number of search sub-areas is determined according to the thickness stability reference value or the dip angle stability reference value;

[0072] If the geological complex condition is a type of geological complex condition, set the search area as a circular area. The center of the circular area is the point to be interpolated, and the area of the circular area is the preset initial area;

[0073] The number of the search sub-areas is positively correlated with the thickness stability reference value or the dip angle stability reference value.

[0074] Among them, the initial area is the area of the search region corresponding to each interpolation point before the search region is adjusted. A method for obtaining the preset initial area is provided. The average value of the areas of the search regions that meet the user's requirements for the accuracy of the interpolation result in the historical records is recorded as the preset initial area.

[0075] Specifically, the geological complexity conditions are determined based on the thickness stability reference value and the dip angle stability reference value. The geological complexity conditions include:

[0076] A type of geological complexity condition where the thickness stability reference value is greater than the preset thickness stability reference value or the dip angle stability reference value is greater than the preset dip angle stability reference value;

[0077] A type of geological complexity condition where the thickness stability reference value is less than or equal to the preset thickness stability reference value or the dip angle stability reference value is less than or equal to the preset dip angle stability reference value.

[0078] Among them, the thickness stability reference value n is the number of data points in the target modeling area, h i is the reference value of the formation thickness measured at the i-th data point. The reference value of the formation thickness is the average value of the formation thicknesses at each borehole point, h0 is the average value of the reference values of the formation thicknesses measured at each data point. The dip angle stability reference value n is the number of data points in the target modeling area, j i is the reference value of the formation dip angle measured at the i-th data point. The reference value of the formation dip angle is the average value of the formation dip angles at each borehole point, j0 is the average value of the reference values of the formation dip angles measured at each data point. The formation thickness is determined according to the logging curve obtained by lowering the resistivity logging tool into the borehole point. According to the change characteristics of the logging curve, the interface of the formation is identified to obtain the formation thickness of each layer. The formation dip angle logging tool is used to measure the dip direction and dip angle of the formation at each borehole point, and the dip angle is recorded as the formation dip angle;

[0079] The user can set the preset thickness stability reference value and the preset dip angle stability reference value according to the actual needs and historical records. The higher the user's requirement for the modeling accuracy of the target modeling area, the smaller the values of the preset thickness stability reference value and the preset dip angle stability reference value. A method for obtaining the preset thickness stability reference value is provided. The average value of the thickness stability reference values that meet the user's requirements for the formation stability of the target modeling area in the historical records is set as the preset formation thickness stability reference value. A method for obtaining the preset dip angle stability reference value is provided. The average value of the dip angle stability reference values that meet the user's requirements for the formation stability of the target modeling area in the historical records is set as the preset dip angle stability reference value.

[0080] Specifically, under the first analysis condition, the eigenvalue of each data point in each search sub-region is detected, the data point difference state of each search sub-region is determined according to the eigenvalue of the data point, and the adjustment method of the search region is determined according to the data point difference state of each search sub-region. The data point difference state includes that the data feature difference degree of the search sub-region is less than the preset feature difference degree, or the data feature difference degree within the search sub-region is greater than or equal to the preset feature difference degree;

[0081] The first analysis condition is that the geological complex condition is a type of geological complex condition and the selection of the setting method of the search region is completed.

[0082] Among them, the eigenvalue m is the type of parameters considered for calculating the eigenvalue, A z is the value of the z-th parameter, B z is the weight coefficient corresponding to the z-th parameter. The number of types of parameters considered for calculating the eigenvalue can be set by the user according to actual needs. The user can set the weight coefficient corresponding to each parameter according to historical experience. The user can judge the influence degree of each parameter on the interpolation process according to historical experience. The greater the influence degree of the parameter on the interpolation process, the greater the weight coefficient corresponding to the parameter. In the present invention, the parameters considered for calculating the eigenvalue include the reference value of the change in compression modulus and the difference value of the foundation bearing capacity. Therefore, m takes the value of 2. A method for obtaining the weight coefficient corresponding to each parameter in the present invention is provided. The weight coefficients of the change value of the compression modulus and the difference value of the foundation bearing capacity are 65% and 35% respectively;

[0083] The compression modulus is the ratio of the vertical compressive stress to the vertical total strain when the soil body is completely unable to deform laterally. The change value of the compression modulus is detected at the same downward depth interval during the drilling process, and the average value of the change values of the compression modulus each time is recorded as the reference value of the change in compression modulus. The foundation bearing capacity is the maximum pressure that the foundation soil can bear per unit area. The difference value of the foundation bearing capacity is the absolute value of the difference between the average value of the foundation bearing capacity in the horizontal direction and the average value of the foundation bearing capacity in the vertical direction. The anisotropy difference between the drilling points is evaluated through the change situation of the compression modulus and the difference value of the foundation bearing capacity between different drilling points, so that the interpolation process fully takes into account the geological spatial anisotropy and improves the accuracy of the modeling result.

[0084] The data feature difference degree b is the number of data points in the search sub-region corresponding to each interpolation point to be interpolated, t ais the eigenvalue of the a-th data point, and t0 is the average value of the eigenvalues of the data points in the search sub-region corresponding to each interpolation point to be interpolated. The user can set the preset feature difference degree according to actual needs and historical records. The higher the user's requirement for the accuracy of the interpolation result, the smaller the value of the preset feature difference degree. A method for obtaining the value of the preset feature difference degree is provided. The historical records that meet the user's requirement for the accuracy of the interpolation result under a certain type of complex geological conditions are set as the first type of reference records, and the maximum value of the data feature difference degrees in the first type of reference records is recorded as the preset feature difference degree. It should be understood that it is easy for those skilled in the art to understand that the user determines the accuracy of the interpolation result through the interpolation result verification method, which will not be elaborated here.

[0085] Specifically, if the data feature difference degree in the search sub-region is greater than or equal to the preset feature difference degree, it is determined by detecting the distribution condition of the difference data points to increase and adjust the horizontal width of the search sub-region according to the data feature difference degree, or increase and adjust the reference length of the edge node corresponding to the search sub-region in the search region according to the data feature difference degree;

[0086] The increased value of the horizontal width is positively correlated with the data feature difference degree;

[0087] The increased value of the reference length of the edge node is positively correlated with the data feature difference degree.

[0088] Among them, when in a certain type of complex geological conditions, the search sub-regions are all rectangles, the center points of the short sides of each search sub-region are the interpolation points to be interpolated, the horizontal width of the search sub-region is the shortest distance between the two long sides, the edge nodes are the intersection points of each search sub-region and the edge of the search region, and the reference length of the edge node is the distance length between the edge node and the interpolation point to be interpolated.

[0089] Specifically, the distribution condition of the difference data points includes that the difference data points are all distributed in the edge region of the corresponding search sub-region, or the distances between the difference data points and the interpolation point to be interpolated are all greater than the preset difference distance;

[0090] If the difference data points are all distributed in the edge region of the sub-region, the horizontal width of the search sub-region is increased and adjusted according to the data feature difference degree;

[0091] If the distances between the difference data points and the interpolation point to be interpolated are all greater than the preset difference distance, the reference length of the edge node corresponding to the search sub-region in the search region is increased and adjusted according to the data feature difference degree.

[0092] Among them, the differential data points are data points whose eigenvalue is greater than the average value of the eigenvalues of the data points within the corresponding search sub-region. The value of the preset differential distance can be set by the user according to actual needs and historical records. The higher the user's requirement for the accuracy of the interpolation result, the smaller the value of the preset differential distance. A method for setting the value of the preset differential distance is provided, where the average value of the distances between the differential data points and the point to be interpolated in a class of reference records is set as the preset differential distance;

[0093] The shortest distance from any position within the edge region of the search sub-region to the center line of the search sub-region is greater than the preset center distance. The center line of the search sub-region is a line segment connecting the center points of the two short sides of the search sub-region. The value of the preset center distance can be set by the user according to actual needs and historical records. The higher the user's requirement for the accuracy of the interpolation result, the smaller the value of the preset center distance. A method for setting the value of the preset center distance is provided, and the value of the preset center distance is 1 / 4 of the short side of the search sub-region.

[0094] Specifically, if the data feature difference degree within the search sub-region is less than the preset feature difference degree, the length of the search sub-region is increased and adjusted according to the data feature difference degree within the search sub-region;

[0095] The increased value of the length of the search sub-region has a negative correlation with the data feature difference degree.

[0096] Among them, when in a class of complex geological conditions, the length of the search sub-region is the distance between the center points of the two short sides of the search sub-region respectively.

[0097] Specifically, under the second analysis condition, the data point acquisition condition is determined according to the number of data points and the data point distribution uniformity within the search region. If the data point acquisition condition is a class of acquisition conditions, the radius length of the search region is increased and adjusted according to the number of data points or the data point distribution uniformity;

[0098] The increased value of the radius length of the search region has a negative correlation with the number of data points or the data point distribution uniformity;

[0099] The second analysis condition is that the geological complex condition is a class II geological complex condition and the selection of the setting method of the search region is completed;

[0100] The data point acquisition conditions include a class of acquisition conditions and a class II acquisition conditions;

[0101] A class of acquisition conditions is that the number of data points is less than the preset number of data points, or the data point distribution uniformity is less than the preset distribution uniformity;

[0102] A class II acquisition condition is that the number of data points is greater than or equal to the preset number of data points and the data point distribution uniformity is greater than or equal to the preset distribution uniformity.

[0103] Among them, the number of data points is the data point distribution uniformity when in the second-class geological complex conditions and the search area area is the initial area f1 is the data point distance distribution uniformity, f 10 is the standard value of the data point distance distribution uniformity s is the number of data points in the search area corresponding to each interpolation point to be interpolated, l c is the distance length between the c-th data point and the interpolation point to be interpolated in the corresponding search area, l0 is the average value of the distance lengths between each data point and the interpolation point to be interpolated in the corresponding search area, f2 is the data point angle distribution uniformity. The search area is evenly divided into several fan-shaped areas with the same area, and the corresponding fan angles of each fan-shaped area are the same, f 20 is the standard value of the data point angle distribution uniformity k is the number of fan-shaped areas, p e is the number of data points in the e-th fan-shaped area, p0 is the average value of the number of data points in each fan-shaped area. The user can set the number of fan-shaped areas according to historical records and actual needs. The higher the accuracy requirement of the user for the modeling result, the more the number of fan-shaped areas. A value of the number of fan-shaped areas is provided. The value of the number of fan-shaped areas is 8. The values of the standard value of the data point distance distribution uniformity and the standard value of the data point angle distribution uniformity can be set by the user according to actual needs and historical records. The higher the accuracy requirement of the user for the modeling result, the smaller the values of the standard value of the data point distance distribution uniformity and the standard value of the data point angle distribution uniformity. A method for obtaining the value of the standard value of the data point distance distribution uniformity is provided. The average value of the data point distance distribution uniformity in the historical records that meet the accuracy requirement of the user for the modeling result is set as the standard value of the data point distance distribution uniformity. A method for obtaining the value of the standard value of the data point angle distribution uniformity is provided. The average value of the data point angle distribution uniformity in the historical records that meet the accuracy requirement of the user for the modeling result is set as the standard value of the data point angle distribution uniformity;

[0104] The values of the preset number of data points and the preset distribution uniformity can be set by the user according to actual needs and historical records. The higher the accuracy requirement of the user for the interpolation result, the larger the value of the preset number of data points and the larger the value of the preset distribution uniformity. A method for obtaining the value of the preset number of data points is provided. The historical records that meet the accuracy requirement of the user for the interpolation result when in the second-class geological complex conditions are recorded as the second-class reference records. The average value of the number of data points in the second-class reference records is set as the preset number of data points. A method for obtaining the value of the preset distribution uniformity is provided. The historical records that meet the accuracy requirement of the user for the interpolation result when in the second-class geological complex conditions are recorded as the second-class reference records. The minimum value of the data point distribution uniformity in the second-class reference records is set as the preset distribution uniformity. A value of the preset distribution uniformity is provided. The value of the preset distribution uniformity is 0.7;

[0105] Increase the adjustment of the radius length of the search area to increase the range of the search area, so as to obtain sufficient data points for the interpolation process and improve the accuracy of the interpolation result.

[0106] Among them, after the search area adjustment is completed, the value of the point to be interpolated and the estimated value of the point to be interpolated are obtained according to the following formula D v is the distance between the v-th data point and the point to be interpolated, N is the number of data points within the adjusted search range of the point to be interpolated, P' v is the v-th data point value, σ is the power parameter, and the user can set the power parameter according to historical experience. A value of the power parameter is provided, and the value of the power parameter is 2.

[0107] Specifically, according to the reference characteristic values of each drilling point, the verification priority of the drilling point is determined, and the drilling points with a preset verification quantity are selected as verification points according to the verification priority rule, and the other drilling points are recorded as data points;

[0108] The reference characteristic value is determined according to whether the data measured at the drilling point is complete and whether there are abnormalities. The reference characteristic value c = ω W ×W + ω G ×G, W is the data integrity of the drilling point, G is the data difference degree of the drilling point, ω W and ω G are weight coefficients;

[0109] The increase value of the verification priority is positively correlated with the reference characteristic value of the drilling point.

[0110] Among them, the data integrity of the drilling point = the number of rock layers collected at the drilling point / the total number of rock layers in the target modeling area, and the data difference degree of the drilling point μ0 is the average value of the characteristic values of each drilling point in the target modeling area, r is the number of drilling points in the target modeling area, μ r is the characteristic value of the r-th drilling point in the target modeling area, ω W and ω G are weight coefficients. The user can set the weight coefficients corresponding to the data integrity of the drilling point and the data difference degree of the drilling point according to historical experience. The user can judge the influence degree of the data integrity of the drilling point and the data difference degree on the verification result according to historical experience. The greater the influence degree of the parameter on the verification result, the greater the weight coefficient corresponding to the parameter.

[0111] Please refer to Figure 4 shown, which is the unit connection diagram of the three-dimensional modeling system for geotechnical engineering of the present invention. The present invention also provides a system for a three-dimensional modeling method for geotechnical engineering, including:

[0112] A data display unit is used to upload, display, and store the drilling information of each drilling point. The drilling information includes the rock layer thickness and dip angle of each rock layer within the drilling point, the compression modulus and foundation bearing capacity obtained during the drilling process.

[0113] A point selection unit is connected to the data display unit and is used to determine the verification priority of each drilling point according to the reference characteristic values of each drilling point, and determine the verification points and data points according to the verification priority of each drilling point.

[0114] A data analysis unit is connected to the data display unit and the point selection unit and is used to determine the geological complex conditions according to the thickness stability reference value and dip angle stability reference value of the geographical reference area where the point to be interpolated is located, and determine the setting method of the search area corresponding to the point to be interpolated according to the geological complex conditions.

[0115] A first analysis unit is connected to the data analysis unit and is used to determine the data point difference state of each search sub-area according to the data point characteristic values, and determine the adjustment method of the search area according to the data point difference state of each search sub-area.

[0116] A second analysis unit is connected to the data analysis unit and is used to determine the data point acquisition conditions according to the number of data points and the distribution uniformity of data points within the search area, and determine whether to adjust the search area range according to the data point acquisition conditions.

[0117] A processing unit is connected to the data analysis unit, the first analysis unit, and the second analysis unit, and performs interpolation data processing on the point to be interpolated according to the drilling information of the data points within the search area of the point to be interpolated, and performs three-dimensional modeling when the interpolation data processing of all points to be interpolated is completed.

[0118] So far, the technical solution of the present invention has been described in conjunction with the preferred embodiments shown in the drawings. However, it is easy for those skilled in the art to understand that the protection scope of the present invention is obviously not limited to these specific embodiments. Without departing from the principle of the present invention, those skilled in the art can make equivalent changes or substitutions to the relevant technical features, and the technical solutions after these changes or substitutions will fall within the protection scope of the present invention.

[0119] The above are only the preferred embodiments of the present invention and are not used to limit the present invention; for those skilled in the art, the present invention can have various changes and modifications. Any modification, equivalent substitution, improvement, etc. made within the spirit and principle of the present invention shall be included within the protection scope of the present invention.

Claims

1. A three-dimensional modeling method for geotechnical engineering, characterized in that, Including: Determine the verification priority for each drilling point according to the reference eigenvalue of each drilling point, and determine the verification points and data points according to the verification priority of each drilling point; Upload the points to be interpolated, determine the geological complexity conditions according to the thickness stability reference value and dip angle stability reference value of the geographical reference area where the points to be interpolated are located, and determine the setting method of the search area corresponding to the points to be interpolated according to the geological complexity conditions; Under the first analysis condition, determine the data point difference status of each search sub-area according to the data point eigenvalue, and determine the adjustment method of the search area according to the data point difference status of each search sub-area; The adjustment method of the search area includes adjusting the lateral width of the search sub-area according to the data feature difference degree or adjusting the reference length of the edge node corresponding to the search sub-area in the search area, or adjusting the length of the search sub-area according to the data feature difference degree of the search sub-area; Under the second analysis condition, determine the data point acquisition condition according to the number of data points and the data point distribution uniformity in the search area, and determine whether to adjust the search area range according to the data point acquisition condition; Perform interpolation data processing on the points to be interpolated according to the drilling information of the data points in the search area of the points to be interpolated, and perform 3D modeling when the interpolation data processing of all points to be interpolated is completed; Mark the drilling points other than the verification points among all the drilling points as data points; The differential data point distribution conditions include that the differential data points are all distributed on the edge of the sub-area, or the distance between the differential data points and the points to be interpolated is greater than the preset differential distance; If the differential data points are all distributed in the edge area of the sub-area, increase and adjust the lateral width of the search sub-area according to the data feature difference degree; If the distance between the differential data points and the points to be interpolated is greater than the preset differential distance, increase and adjust the reference length of the edge node corresponding to the search sub-area in the search area according to the data feature difference degree.

2. The three-dimensional modeling method for geotechnical engineering according to claim 1, characterized in that Determine the setting method of the search area according to the geological complexity conditions, including: If the geological complexity condition is a type I geological complexity condition, set the search area as a circular area composed of several search sub-areas, the center of the circular area is the point to be interpolated, the initial areas of each search sub-area are the same, and the number of search sub-areas is determined according to the thickness stability reference value or dip angle stability reference value; If the geological complexity condition is a type II geological complexity condition, set the search area as a circular area, the center of the circular area is the point to be interpolated, and the area of the circular area is the preset initial area; The number of the search sub-areas is positively correlated with the thickness stability reference value or dip angle stability reference value.

3. The three-dimensional modeling method for geotechnical engineering according to claim 2, characterized in that, The geological complexity conditions are determined according to the thickness stability reference value and dip angle stability reference value, and the geological complexity conditions include: A type I geological complexity condition where the thickness stability reference value is greater than the preset thickness stability reference value or the dip angle stability reference value is greater than the preset dip angle stability reference value; A type II geological complexity condition where the thickness stability reference value is less than or equal to the preset thickness stability reference value or the dip angle stability reference value is less than or equal to the preset dip angle stability reference value.

4. The 3D modeling method for geotechnical engineering according to claim 3, wherein Under the first analysis condition, the eigenvalue of each data point in each search sub-region is detected, the data point difference state of each search sub-region is determined according to the eigenvalue of the data point, and the adjustment method of the search region is determined according to the data point difference state of each search sub-region. The data point difference state includes that the data feature difference degree of the search sub-region is less than the preset feature difference degree, or the data feature difference degree within the search sub-region is greater than or equal to the preset feature difference degree; The first analysis condition is that the geological complexity condition is a type of geological complexity condition and the selection of the setting method of the search region is completed.

5. The three-dimensional modeling method for geotechnical engineering according to claim 4, wherein If the data feature difference degree within the search sub-region is greater than or equal to the preset feature difference degree, the horizontal width of the search sub-region is increased according to the data feature difference degree by detecting the distribution condition of the differential data points, or the reference length of the edge node corresponding to the search sub-region in the search region is increased according to the data feature difference degree; The increase value of the horizontal width is positively correlated with the data feature difference degree; The increase value of the reference length of the edge node is positively correlated with the data feature difference degree.

6. The three-dimensional modeling method for geotechnical engineering according to claim 5, characterized in that, If the data feature difference degree of the search sub-region is less than the preset feature difference degree, the length of the search sub-region is increased according to the data feature difference degree of the search sub-region; The increase value of the length of the search sub-region is negatively correlated with the data feature difference degree.

7. The three-dimensional modeling method for geotechnical engineering according to claim 6, wherein Under the second analysis condition, the data point acquisition condition is determined according to the number of data points and the data point distribution uniformity in the search region. If the data point acquisition condition is a type of acquisition condition, the radius length of the search region is increased according to the number of data points or the data point distribution uniformity; The increase value of the radius length of the search region is negatively correlated with the number of data points or the data point distribution uniformity; The second analysis condition is that the geological complexity condition is a type of geological complexity condition and the selection of the setting method of the search region is completed; The data point acquisition condition includes a type of acquisition condition and a type II acquisition condition; A type of acquisition condition is that the number of data points is less than the preset number of data points, or the data point distribution uniformity is less than the preset distribution uniformity; The type II acquisition condition is that the number of data points is greater than or equal to the preset number of data points and the data point distribution uniformity is greater than or equal to the preset distribution uniformity.

8. The three-dimensional modeling method for geotechnical engineering according to claim 7, characterized in that, The verification priority of the borehole points is determined according to the reference eigenvalue of each borehole point, and the borehole points with a preset verification number are selected as verification points according to the verification priority rule, and the other borehole points are recorded as data points; The reference characteristic value is determined according to whether the data measured at the drilling point is complete and whether there are abnormalities. The reference characteristic value c = ω W ×W + ω G ×G, where W is the integrity of the data at the drilling point, G is the difference degree of the data at the drilling point, and ω W and ω G are weight coefficients; The increase value of the verification priority is positively correlated with the reference eigenvalue of the borehole point.

9. A system applying the three-dimensional modeling method for geotechnical engineering according to any one of claims 1 to 8, characterized in that, Including: A data display unit for uploading, displaying and storing the borehole information of each borehole point. The borehole information includes the rock layer thickness and dip angle of each rock layer in the borehole point, the compression modulus and foundation bearing capacity obtained during the borehole process; A point position screening unit connected to the data display unit for determining the verification priority of the borehole points according to the reference eigenvalue of each borehole point, and determining the verification points and data points according to the verification priority of each borehole point; A data analysis unit, which is connected to the data display unit and the point selection unit, is used to determine geological complex conditions according to the thickness stability reference value and dip angle stability reference value of the geographical reference area where the point to be interpolated is located, and determine the setting method of the search area corresponding to the point to be interpolated according to the geological complex conditions; A first analysis unit, which is connected to the data analysis unit, is used to determine the data point difference state of each search sub-area according to the data point characteristic value, and determine the adjustment method of the search area according to the data point difference state of each search sub-area; A second analysis unit, which is connected to the data analysis unit, is used to determine the data point acquisition conditions according to the number of data points and the data point distribution uniformity in the search area, and determine whether to adjust the search area range according to the data point acquisition conditions; A processing unit, which is connected to the data analysis unit, the first analysis unit and the second analysis unit, performs interpolation data processing on the point to be interpolated according to the borehole information of the data points in the search area of the point to be interpolated, and performs three-dimensional modeling when the interpolation data processing of all points to be interpolated is completed.

Images