Geological boundary modeling method and device based on function image, equipment and medium

Abstract

The invention discloses a geological boundary modeling method and device based on a function image, equipment and a medium, and relates to the technical field of geological surveying and mapping. The method comprises the following steps: dividing drilling data into continuous multi-section data along the same horizontal direction of a geological section, including first multi-section data and second multi-section data with different drilling data volumes, and respectively performing segmented fitting on the first multi-section data and the second multi-section data by adopting corresponding fitting functions, and finally, the fitting curve of each segment is drawn into a geological boundary. According to the modeling method, the drawing efficiency and smoothness of the geological boundary are improved.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of geological mapping, and in particular to a geological boundary modeling method and device based on function images, equipment and a medium. BACKGROUND

[0002] In engineering geological work, it is often necessary to draw a geological profile on a specified route. The existing geological profile drawing method is usually to draw a polyline by CAD according to the drilling depth value of the geological interface, the depth value of other exploration points, the surface outcrop point and the boundary line point and the like, and to connect the above-mentioned points into a curve, which is the geological boundary, such as the base-cover boundary, the strong-weak weathering boundary, the weak weathering bottom boundary, the soft-hard rock mass boundary and the like.

[0003] The traditional CAD polyline method for drawing a geological profile has the defects of large drawing workload, large randomness of manual line drawing and non-smoothness of the curve line, which leads to a large deviation of the profile from the actual situation and easy errors in subsequent application. Moreover, the efficiency of this drawing method is not high, and it cannot meet the needs of current actual work. SUMMARY

[0004] The present application provides a geological boundary modeling method and device based on function images to solve the problems of large workload, non-smoothness of the curve and large deviation from the actual situation in the existing CAD polyline method for drawing a geological profile boundary.

[0005] The present application is implemented by the following technical solutions: In a first aspect, the present application provides a geological boundary modeling method based on function images, comprising: obtaining drilling data of a same geological profile; dividing the drilling data into continuous multi-segment data along the same direction of the geological profile, the multi-segment data comprising first multi-segment data and second multi-segment data; wherein, the number of drilling data contained in the first multi-segment data is within a first range, and the number of drilling data contained in the second multi-segment data is within a second range, and the first range and the second range do not have an intersection; performing curve fitting on the first multi-segment data and the second multi-segment data by using different fitting functions respectively to obtain a fitting curve corresponding to each segment of the multi-segment data; sequentially and smoothly connecting the multi-segment fitting curves to obtain a geological boundary of the geological profile.

[0006] The present invention's mapping method first continuously segments the measured borehole data of a geological profile into two different data types, first and second multi-segment data, based on a preset segmentation range. Different fitting functions are then applied to each type of segmented data, resulting in a better fitting effect for each segmented data type. The segmented geological boundaries obtained from the function fitting are smoother than polylines, and the appropriate fitting function is selected based on the amount of segmented data, significantly reducing manual work and improving mapping efficiency.

[0007] Furthermore, different fitting functions are used to perform curve fitting on the first multiple segments of data and the second multiple segments of data, respectively, including: performing curve fitting on the first multiple segments of data using a first fitting function, and performing curve fitting on the second multiple segments of data using a second fitting function.

[0008] Furthermore, the first plurality of segments of data include 4 or 5 drilling data, and the second plurality of segments of data include 2 or 3 drilling data.

[0009] Furthermore, the first fitting function is expressed as:

[0010] in, 、 、 、 、 is a parameter to be determined, " represents the multiplication operation, and each segmented data is fitted in its own local coordinate system. L Indicates the length of the segment to be fitted, representing the horizontal coordinate axis, x is the lateral position of the drill hole in the local coordinate system, x Values ​​range from 0 to L , H is the vertical total height of the segment to be fitted, and takes a negative value when going downward from the origin, representing the vertical axis. z is the vertical position of the drilling depth in the local coordinate system, and z takes a value of 0~ H .

[0011] Furthermore, performing curve fitting on the first plurality of segments of data using a first fitting function includes: If the first plurality of segments of data contains 5 drilling data, then during fitting, take m =0.5, n =1,0.9< a <1, adjust the base a , iterative calculation to find p and r , so that the fitting curve passes through the boundary points of the three middle boreholes excluding the endpoints, and draws the fitting curve of the first plurality of segments of data; If the first plurality of segments of data contain 4 drilling data, then during fitting, take m =0.5, n =1, p =1,0.9< a <1, adjust the base a , iterative calculation to find r The value of is set so that the fitting curve passes through the boundary points of the two middle boreholes except the endpoints, and the fitting curve of the first plurality of segments of data is drawn.

[0012] Furthermore, the second fitting function is expressed as:

[0013] in, is a parameter to be determined, " represents the multiplication operation, and each segment is fitted in its own local coordinate system. L Indicates the length of the fitting segment and represents the horizontal coordinate axis. x is the lateral position of the drill hole in the local coordinate system, x Values ​​range from 0 to L , H The vertical total height of the fitting segment is negative when it is downward from the origin, representing the vertical axis. z is the vertical position of the drilling depth in the local coordinate system, and z takes a value of 0~ H。

[0014] Furthermore, performing curve fitting on the second plurality of segments of data using a second fitting function includes: If the second multi-segment data contains three drilling data, then when fitting, the local coordinates of the middle drilling holes except the endpoints are brought into the second fitting function to obtain a value, and draw a fitting curve of the second plurality of segments of data; If the second multi-segment data contains two drilling data, then during fitting, the slope of the common point is obtained based on the fitting curve of the multi-segment data adjacent to the second multi-segment data, and the slope of the common point is obtained based on the slope of the common point. a values, and draw a fitting curve for the second plurality of segments of data.

[0015] A second aspect of the present invention provides a geological boundary modeling device based on a function image, comprising: Input module, used to obtain drilling data of the same geological section; The segmentation module is used to divide the drilling data into a plurality of continuous segments along the same horizontal direction of the geological section, wherein the plurality of segments include a first plurality of segments and a second plurality of segments; wherein, The number of drilling data included in the first plurality of segments of data is within a first range, the number of drilling data included in the second plurality of segments of data is within a second range, and the first range and the second range do not intersect; a curve fitting module, performing curve fitting on the first plurality of segments of data and the second plurality of segments of data using different fitting functions, respectively, to obtain a fitting curve corresponding to each segment of the plurality of segments of data; The combination module is used to smoothly connect multiple segments of the fitting curve in sequence to obtain the geological boundary of the geological section.

[0016] The third aspect of the present invention provides an electronic device, comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein when the processor executes the computer program, the method for geological boundary modeling based on function images as described in any one of the first aspects of the present invention is implemented.

[0017] A fourth aspect of the present invention provides a computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the geological boundary modeling method based on function images as described in any one of the first aspects of the present invention.

[0018] Compared with the prior art, the present invention has the following advantages and beneficial effects: the drilling data is reasonably segmented, and suitable fitting functions are used for different data amounts, thereby improving the degree of drawing automation and reducing manual workload; the segmentation method and fitting function of the present invention are used to improve the smoothness and accuracy of the drawn segmented curves, which is beneficial to the subsequent application of geological profiles. BRIEF DESCRIPTION OF THE DRAWINGS

[0019] In order to more clearly illustrate the technical solutions of the exemplary embodiments of the present invention, the following briefly introduces the drawings required for use in the examples. It should be understood that the following drawings only illustrate certain embodiments of the present invention and should not be considered as limiting the scope. A person of ordinary skill in the art can also derive other relevant drawings based on these drawings without inventive effort. In the drawings: Figure 1 It is a schematic diagram of geological boundaries drawn by CAD polylines; Figure 2 This is a flow chart of a geological boundary modeling method based on function images according to an embodiment of the present invention; Figure 3 This is a continuous segmented schematic diagram of an embodiment of the present invention; Figure 4 is a schematic diagram of a second segment fitting curve according to an embodiment of the present invention; Figure 5 is a schematic diagram of a fitting curve of the third segment in an embodiment of the present invention; Figure 6 is a schematic diagram of a fitting curve of a first segment in an embodiment of the present invention; Figure 7 It is a comparison diagram of the geological boundary obtained by fitting the geological interface function of the present invention and the geological boundary drawn by the multi-segment line; Figure 8 This is a schematic diagram of a fitting curve of 5 drilling segments according to an embodiment of the present invention; Figure 9 It is a structural schematic diagram of a geological boundary modeling device based on function images according to an embodiment of the present invention. DETAILED DESCRIPTION

[0020] In order to make the objectives, technical solutions and advantages of the present invention more clearly understood, the present invention is further described in detail below in conjunction with examples and drawings. The exemplary embodiments of the present invention and their descriptions are only used to explain the present invention and are not intended to limit the present invention.

[0021] It should be noted that the terms "including" and "having" and any variations thereof in the specification and claims of the present invention and the above-mentioned drawings are intended to cover non-exclusive inclusions. For example, a process, method, system, product or device that includes a series of steps or units is not necessarily limited to other steps or units inherent in the device.

[0022] The terms used in various embodiments of the present invention are only used to describe the purpose of specific embodiments and are not intended to limit the various embodiments of the present invention. As used herein, the singular form is intended to also include the plural form, unless the context clearly indicates otherwise. Unless otherwise limited, all terms used here (including technical terms and scientific terms) have the same meaning as those of ordinary skill in the art generally understood by the various embodiments of the present invention. The terms (such as those defined in generally used dictionaries) will be interpreted as having the same meaning as the contextual meaning in the relevant technical field and will not be interpreted as having idealized meaning or too formal meaning, unless clearly defined in various embodiments of the present invention.

[0023] The embodiments of the present invention provide a geological boundary modeling method, device, equipment and medium based on function images, which are suitable for drawing various boundary lines of geological profiles, are conducive to improving the degree of drawing automation, reducing manual workload, and improving the smoothness and accuracy of drawing decomposition curves, so as to facilitate subsequent applications.

[0024] like Figure 1The figure shows a schematic diagram of a geological boundary drawn using CAD polylines. This section has seven known borehole data points, ZK01 to ZK07, and two surface outcrops, O1 and O2. The polyline appears slightly smooth in areas with dense boreholes, but is noticeably uneven in areas with sparse boreholes and large elevation differences. This not only makes the drawn section less aesthetically pleasing, but also brings inconvenience to data processing for modeling geological interfaces and subsequent applications.

[0025] When building a three-dimensional model of a geological body using the skeleton section method, a series of skeleton sections must be created. Points on the geological boundary lines of each skeleton section are extracted, and the geological boundary lines of the sections are formed through all these points. A spatial three-dimensional surface is then constructed from the series of skeleton sections, ultimately yielding a geological interface. To ensure that the geological boundary lines of the geological cross-sections and skeleton sections used in the three-dimensional geological model are smooth and aesthetically pleasing, conform to geological laws, and are free of distortion or sudden changes, the present invention provides an efficient and widely applicable method for modeling geological boundaries based on function graphs.

[0026] like Figure 2 FIG. 1 is a flow chart of a method for modeling geological boundaries based on function images proposed in the present invention, comprising the following steps: S1, obtain the drilling data of the same geological section.

[0027] S2, dividing the drilling data into multiple continuous segments along the same horizontal direction of the geological section, wherein the multiple segments include a first multiple segment of data and a second multiple segment of data having different amounts of drilling data.

[0028] S3, performing curve fitting on the first multi-segment data and the second multi-segment data using different fitting functions, respectively, to obtain segmented geological boundaries corresponding to each segment of the multi-segment data; S4, sequentially and smoothly connecting multiple segmented geological boundaries to obtain the geological boundaries of the geological profile.

[0029] This embodiment models the geological boundary of a geological section to be modeled, and it is necessary to obtain the drilling data on the section in advance. It can be understood that when constructing a geological section, the two endpoints of the section also include two outcrops on the left and right. Therefore, the drilling data to be obtained also include two outcrops on the left and right. Due to the complexity of mountainous terrain and geological conditions, the geological interface often crosses the river valley to form a U-shaped or V-shaped surface. Therefore, it is characterized by a large difference in elevation between the two ends of the drilling holes and a dense distribution of drilling holes in the middle. The general fitting function has a poor fitting effect on the overall data. The effect of numerical fitting is related to the amount of data, the data itself and the fitting function. Different fitting data amounts and data characteristics are suitable for different fitting functions, so it is necessary to reasonably segment the data.

[0030] In step S2, this embodiment continuously segments the drillhole data on the cross-section along the same direction based on the amount of drillhole data, resulting in two types of multi-segment data. Adjacent segments are continuous, with no drillholes missing. The drillhole data is divided into two types of segments based on pre-set segmentation ranges: the number of drillhole data included in the first multi-segment data falls within a first range, and the number of drillhole data included in the second multi-segment data falls within a second range. The first and second ranges are set to have no intersection, meaning that the amount of drillhole data in the first and second multi-segment data differs, allowing for different fitting functions.

[0031] Continuous segmentation is performed along the same direction, for example, from the outcrop point on the left bank to the outcrop point on the right bank, or from the outcrop point on the right bank to the outcrop point on the left bank. When considering which type of segmentation to divide the next set of unsegmented borehole data into, the elevation difference between adjacent boreholes can be comprehensively considered. When the elevation difference is large, a smaller amount of data is divided into one segment. Conversely, a denser area is divided into a larger amount of data. This allows the fitting curve to be more consistent with the actual data. The segmentation step of this embodiment can also be set to more than two division ranges, divided into two or more segmentation types, and one type is applicable to one fitting function. Dividing into two segmentation types can simplify the segmentation process and the function fitting process, and optimize the connection between multiple fitting curves to avoid mutations at the connection points.

[0032] In step S3, after selecting an appropriate fitting function, the program can be set to automatically fit the function, i.e., automatically solve the function fitting parameters based on the known coordinates. Finally, the fitting curve is plotted using a curve drawing tool to obtain the segmented geological boundaries of each segment. Each segment can be fitted sequentially according to the order of division, or multiple threads can be used to fit the segments simultaneously, and then the geological boundaries of each segment are connected together.

[0033] When fitting each segment, a unified reference coordinate system is used, such as the left bank outcrop or the right bank outcrop as the origin, the horizontal direction along the divided borehole data as the X-axis, the height direction or the reverse direction as the Y-axis, the horizontal position of the borehole as the horizontal coordinate reference value, and the elevation of the borehole as the vertical coordinate reference value. Of course, other coordinate systems can also be used, and any borehole data can be used as the origin of the unified reference coordinate system, which will not affect the final presentation of the geological boundary of the profile.

[0034] When fitting each segment, the local coordinate system of each segment is applied. For example, the horizontal length of the segment to be fitted is used as the horizontal axis, and the relative height of the segment (which can be positive or negative) is used as the vertical axis. After fitting in the local coordinate system, it is converted to a unified reference coordinate system to splice the geological boundaries of each segment.

[0035] Furthermore, the modeling method of the present invention can be used to draw geological boundaries of multiple geological sections and establish a three-dimensional model of the geological body.

[0036] In one embodiment, to optimize function fitting, the first range for segmentation is set to 3-5, and the second range is set to 2. That is, the first multi-segment data includes data from 3 to 5 and borehole data, while the second multi-segment data includes data from two boreholes. During the specific segmentation, the elevation difference and distribution concentration of adjacent boreholes are considered to determine whether the multiple consecutive points are divided into the first multi-segment data or the second multi-segment data.

[0037] The following is a specific example to illustrate the data segmentation process. Figure 3 The continuous segmentation diagram shown in the figure shows the profile to be modeled, consisting of the left bank outcrop O1, the right bank outcrop O2, and the seven borehole data points ZK01 to ZK07 in between. Continuous segmentation is performed along the direction from left bank outcrop O1 to right bank outcrop O2, resulting in three continuous segmented data points: the first segment, the second segment, and the third segment. The first segment contains two data points, namely the left bank outcrop O1 and borehole ZK01, representing the second multi-segment data. The second segment contains five data points, namely boreholes ZK01 to ZK05. The third segment contains four data points, namely boreholes ZK05 to ZK07 and the right bank outcrop O2, representing the first multi-segment data. Curve fitting was performed using the first fitting function for the first multi-segment data, and the second fitting function for the second multi-segment data.

[0038] Furthermore, to achieve smooth connections between multiple fitting curves, the present invention proposes an efficient and widely applicable geological boundary function (such as the base-cover boundary function) for the above-mentioned segmented fitting. Specifically, the general form of the geological boundary function is:

[0039] in, 、 、 、 、 All parameters are pending. " represents the multiplication operation, and each segment is fitted in its own local coordinate system. L Indicates the length of the fitting segment and represents the horizontal coordinate axis. x is the lateral position of the drill hole in the local coordinate system, x Values ​​range from 0 to L , H is the vertical total height of the fitting segment (positive value when upward from the origin, negative value when downward from the origin), represents the vertical axis, z is the longitudinal position of the drilling depth in the local coordinate system, and z takes values ​​from 0 to H In actual calculation, the fitting curve may be larger thanH .

[0040] Below Figure 3 The fitting process of each segment data shown in FIG. 1 illustrates the application of the geological boundary function of the present invention. First, take the second segment as an example. Figure 4 The figure shows the fitting curve of the second segment. The second segment contains 5 data points from boreholes ZK01 to ZK05. Its highest point is the geological boundary point of borehole ZK01, and its lowest point is the boundary point of borehole ZK05. Figure 3 The local coordinate system shown in HL There are three boreholes ZK02, ZK03 and ZK04 in the middle of the second curve, which are suitable for the general form of geological interface function. m =0.5, n =1,0.9< a <1. Adjust the base a , use Excel spreadsheet or calculation program to iteratively calculate p and r Two undetermined coefficients.

[0041] For example, during iterative calculation, first take a The initial value of is 0.96, and the coordinates of holes ZK02, ZK03, and ZK04 are substituted and iterative calculation is performed. p and r , and then continue to adjust within the range of 0.96~1 a until the fitting curve of the geological interface function passes through the geological boundary points of the three boreholes ZK2, ZK3, and ZK4. a Converges to a value of 0.9939, corresponding to p =0.4085, r =-1.5651, that is, the functional form of the fitting curve of the second segment is expressed as: ,in H is the total vertical height of the second segment. The curve corresponding to this function can be drawn using the drawing plug-in tool.

[0042] like Figure 5 The figure shows the fitting curve of the third segment. First, the third segment is mirror-symmetrically processed to establish the local coordinate system of the third segment. HL , at this time, in the local coordinate system of the third segment L Axis direction and the second segment local coordinate system L In addition, a fourth quadrant HL The coordinate system is the local coordinate system of the third segment pair, that is, the intersection point O of the vertical line through the borehole ZK05 and the horizontal line through the right bank outcrop point O2 is taken as the origin of the local coordinate system, which is located in the upper left corner. LThe axis direction is along O-O2, and the local coordinate system of the second segment L The axis direction is the same, H Take negative value, z≤0.

[0043] The second segment contains two boreholes ZK06 and ZK07 in addition to the two endpoints. Since there are only two data, the calculation process takes p =1, the simplified form of the geological interface function is as follows:

[0044] Similarly, take m =0.5, n =1, 0.9<a<1, adjust the base a , you only need to use Excel spreadsheet or calculation program to iteratively calculate r When the function graph passes through the two drilling points ZK06 and ZK07, a Converges to 0.972, at which point the parameter r =1.968. The functional form of the fitting curve of the third segment is expressed as: . Then use the drawing plug-in tool to draw the fitting curve of the second segment, and then use the drawn curve to H The axis is used as the axis of symmetry for mirroring and moved to the correct position in the figure to obtain the geological boundary of the second segment.

[0045] like Figure 6 The figure shows the fitting curve of the first segment. The local coordinate system of the first segment is established as shown in the figure. HL There is no borehole in the middle of this section, and it only contains two endpoints. However, it intersects with the adjacent first segment at borehole ZK01. At this point, the two curves transition smoothly, that is, the derivatives of the two curves at the common point should be roughly equal. Therefore, the geological boundary of the first segment can be simply made based on the derivative of the intersection point on the basis of the second segment curve.

[0046] The first segment is applied with the second fitting function, which is also a simplified form of the geological interface function and is expressed as:

[0047] in, is a parameter to be determined, " represents the multiplication operation, and each segment is fitted in its own local coordinate system. L Indicates the length of the fitting segment and represents the horizontal coordinate axis. x is the lateral position of the drill hole in the local coordinate system, x Values ​​range from 0 to L , H is the vertical total height of the fitting segment, representing the ordinate axis, z is the longitudinal position of the drilling depth in the local coordinate system, and z takes values ​​of 0~H In actual calculation, the local fitting curve may be larger than H .

[0048] Similarly, first measure the horizontal inclination of the geological boundary line of the adjacent second segment at the intersection point ZK01, and then calculate the derivative of the intersection point, that is, the horizontal inclination of the geological boundary line of the first segment in the local coordinate system. x = L The derivative at can be found a Finally, use the drawing plug-in tool to obtain the first geological boundary.

[0049] Thus, the geological boundaries of the first, second and third segments are smoothly connected together to obtain the geological boundaries of the section to be modeled, such as Figure 7 The figure shows a comparison between the geological boundary obtained by fitting the geological interface function of the present invention and the geological boundary drawn by multi-segment lines. It can be seen that the geological boundary drawn by the interface function is smoother and improves the deformed mutation situation.

[0050] The following is another fitting example to illustrate the fitting effect of the geological interface function of this embodiment. Figure 8 The figure shows the fitting curve of 5 drilling segments, the local coordinate system HL As shown in the figure, the two end points of the segment are drilled in the local coordinate system of the segment. HL The coordinates in the are (0,200) and (240,0). The three holes in the middle are in the segmented local coordinate system. HL The coordinates in are (50,60), (127,70), and (210,20), respectively. The first fitting function is applicable, which is the general form of the geological interface function: , when calculating m =0.5, n =1,0.9< a <1. The geological interface curve must pass through the two end holes (0,200), (240,0) and the three middle holes. According to this requirement, iterative calculation is obtained. a =0.98062, p =93.01, r =-95.81, the function expression of the fitting curve is: .like Figure 8 It can be seen that the geological interface function fitting method can be used to make geological boundaries with complex concavity and convexity changes, and it has a wide range of applications.

[0051] When using the geological interface function fitting method implemented in this embodiment to make a geological profile or a skeleton profile for three-dimensional modeling, the drilling data is first reasonably segmented. Each segment only needs to contain 4 to 5 boreholes or 2 to 3 boreholes. The amount of fitting data is small and the calculation is efficient. Different fitting function forms are used for segments with different data amounts, so that the fitting effect of each segment is good, and adjacent segments of different types are smoothly connected using derivatives to optimize the smoothness at the connection points.

[0052] For the case where a segment contains 4 to 5 boreholes, the fitting function (1) is applied. When there are 3 boreholes in the middle of a segment, the three middle boreholes are used to find the unknown parameters in the local coordinate system through the iteration method. a 、 p 、 r , and then draw the fitting curve; when there are two drill holes in the middle of a section, take p = 1, and iteratively calculated in the local coordinate system using the two middle boreholes. a and r , make a fitting curve. For the case where the segment contains 2~3 drill holes, the fitting function (2) is applied. When there is a drill hole in the middle of a segment, the local coordinates of the middle drill hole are substituted into the function to obtain a When there is no drill hole in a certain section, the derivative of the section and the known adjacent section at the common point is roughly equal, and the fitting curve is obtained. a value, and then make a fitting curve.

[0053] The process of adjusting the parameters of the function fitting method implemented in this paper is simple and feasible. a A value within the small range of 0.9 to 1 is sufficient, and it is easy to get a suitable a The value only needs to be adjusted in rare cases if the calculation does not converge. a The value is slightly greater than 1. This calculation process can be performed in Excel. A standard calculation table is pre-compiled to easily and quickly calculate the value. a 、 p 、 r Three parameters to improve fitting efficiency.

[0054] The second aspect of the present invention provides a geological boundary modeling device based on function images, such as Figure 9 The schematic diagram of the device structure shown includes: Input module 100, used to obtain drilling data of the same geological section; The segmentation module 200 is used to divide the drilling data into multiple continuous segments along the same horizontal direction of the geological section, wherein the multiple segments include first multiple segments of data and second multiple segments of data; wherein, The number of drilling data included in the first plurality of segments of data is within a first range, the number of drilling data included in the second plurality of segments of data is within a second range, and the first range and the second range do not intersect; The curve fitting module 300 performs curve fitting on the first plurality of segments of data and the second plurality of segments of data using different fitting functions to obtain a fitting curve corresponding to each plurality of segments of data; The combination module 400 smoothly connects multiple fitting curves in sequence to obtain the geological boundary of the geological section.

[0055] In one embodiment, the segmentation module divides the drilling data into the first multiple segments of data or the second multiple segments of data according to a preset division range, wherein the first range is set to 4-5 and the second range is set to 2-3.

[0056] Furthermore, the curve fitting module calls a preset fitting function to fit each segmented data, calls a first fitting function to perform curve fitting on the first plurality of segments of data, and calls a second fitting function to perform curve fitting on the second plurality of segments of data.

[0057] Furthermore, the curve fitting module includes a function fitting submodule and a drawing submodule. The function fitting submodule is used to solve the function's undetermined parameters based on the fitting function and the segmented data. The specific solution process is described in the previous embodiment. The drawing submodule is used to draw the curve based on the fitting function, i.e., the solved undetermined parameters.

[0058] Furthermore, the above device also includes a setting module for setting data such as the first range, the second range, and the fitting function.

[0059] Furthermore, the above device also includes a display module for displaying visual data such as fitting curves and geological boundaries.

[0060] In a third aspect, embodiments of the present invention further provide an electronic device comprising a processor and a memory, where the number of processors may be one or more. The memory, as a computer-readable storage medium, can be used to store software programs, computer-executable programs, and modules. The processor executes the software programs, instructions, and modules stored in the memory to perform various functional applications and data processing of the electronic device, thereby implementing the function graph-based geological boundary modeling method of any of the aforementioned embodiments of the present invention.

[0061] The memory may primarily include a program storage area and a data storage area. The program storage area may store an operating system and at least one application required for a function; the data storage area may store data generated based on the use of the terminal, etc. Furthermore, the memory may include high-speed random access memory and may also include non-volatile memory, such as at least one disk storage device, flash memory device, or other non-volatile solid-state memory device. In some instances, the memory may further include memory remotely located relative to the processor, and these remote memories may be connected to the electronic device via a network. Examples of the aforementioned networks include, but are not limited to, the Internet, an intranet, a local area network, a mobile communication network, and combinations thereof.

[0062] In a fourth aspect, an embodiment of the present invention further provides a computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the geological boundary modeling method based on function images according to any embodiment of the present invention.

[0063] The computer storage media of the embodiments of the present invention may employ any combination of one or more computer-readable media. The computer-readable medium may be a computer-readable signal medium or a computer-readable storage medium. The computer-readable storage medium may be, for example, but not limited to, an electrical, magnetic, optical, electromagnetic, infrared, or semiconductor system, device, or component, or any combination thereof. More specific examples (a non-exhaustive list) of computer-readable storage media include: an electrical connection having one or more wires, a portable computer disk, a hard disk, a random access memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or flash memory), an optical fiber, a portable compact disk read-only memory (CD-ROM), an optical storage device, a magnetic storage device, or any suitable combination thereof. In the present invention, a computer-readable storage medium may be any tangible medium containing or storing a program that can be used by or in conjunction with an instruction execution system, device, or component.

[0064] A computer-readable signal medium may include a data signal propagated in baseband or as part of a carrier wave, which carries computer-readable program code. Such propagated data signals may take various forms, including but not limited to electromagnetic signals, optical signals, or any suitable combination thereof. A computer-readable signal medium may also be any computer-readable medium other than a computer-readable storage medium that can transmit, propagate, or transport a program for use by or in conjunction with an instruction execution system, apparatus, or device.

[0065] In a fifth aspect, an embodiment of the present invention further provides a computer program product, which, when run on a computer, enables the computer to execute the geological boundary (such as base-cover boundary) modeling method based on function images of any of the above embodiments of the present invention.

[0066] The specific implementation methods described above further illustrate the objectives, technical solutions and beneficial effects of the present invention in detail. It should be understood that the above description is only a specific implementation method of the present invention and is not intended to limit the scope of protection of the present invention. Any modifications, equivalent substitutions, improvements, etc. made within the spirit and principles of the present invention should be included in the scope of protection of the present invention.

Claims

1. A geological boundary modeling method based on function images, characterized in that: include: Obtaining drilling data of the same geological section; Dividing the drilling data into a continuous plurality of segments of data along the same direction of the geological section, wherein the plurality of segments of data includes a first plurality of segments of data and a second plurality of segments of data; wherein, The number of drilling data included in the first plurality of segments of data is within a first range, the number of drilling data included in the second plurality of segments of data is within a second range, and the first range and the second range do not intersect; Performing curve fitting on the first plurality of segments of data and the second plurality of segments of data using different fitting functions, respectively, to obtain a fitting curve corresponding to each segment of the plurality of segments of data; The multiple segments of the fitting curves are connected smoothly in sequence to obtain the geological boundary of the geological section.

2. The method for geological boundary modeling based on function graph according to claim 1, characterized in that: The curve fitting is performed on the first plurality of segments of data and the second plurality of segments of data using different fitting functions respectively, including: curve fitting is performed on the first plurality of segments of data using a first fitting function, and curve fitting is performed on the second plurality of segments of data using a second fitting function.

3. The method for geological boundary modeling based on function graph according to claim 2, characterized in that: The first plurality of segments of data include 4 or 5 drilling data, and the second plurality of segments of data include 2 or 3 drilling data.

4. The method for geological boundary modeling based on function graph according to claim 3, characterized in that: The first fitting function is expressed as: in, 、 、 、 、 is a parameter to be determined, " represents the multiplication operation, and each segmented data is fitted in its own local coordinate system. L Indicates the length of the segment to be fitted, representing the horizontal coordinate axis, x is the lateral position of the drill hole in the local coordinate system, x Values ​​range from 0 to L , H is the vertical total height of the segment to be fitted, and takes a negative value when going downward from the origin, representing the longitudinal coordinate axis. z is the longitudinal position of the drilling depth in the local coordinate system, and z takes a value of 0~ H .

5. The method for geological boundary modeling based on function graph according to claim 4, characterized in that: The performing curve fitting on the first plurality of segments of data using a first fitting function includes: If the first plurality of segments of data contains 5 drilling data, then during fitting, take m =0.5, n =1,0.9< a <1, adjust the base a , iterative calculation to find p and r , so that the fitting curve passes through the geological boundary points of the three middle boreholes except the endpoints, and draws the fitting curve of the first plurality of segments of data; If the first plurality of segments of data contains 4 drilling data, then during fitting, take m =0.5, n =1, p =1,0.9< a <1, adjust the base a , iterative calculation to find r The value of is set so that the fitting curve passes through the boundary points of the two middle boreholes except the endpoints, and the fitting curve of the first plurality of segments of data is drawn.

6. The method for geological boundary modeling based on function graph according to claim 3, characterized in that: The second fitting function is expressed as: in, is a parameter to be determined, " represents the multiplication operation, and each segment is fitted in its own local coordinate system. L Indicates the length of the fitting segment and represents the horizontal coordinate axis. x is the lateral position of the drill hole in the local coordinate system, x Values ​​range from 0 to L , H The vertical total height of the fitting segment is negative when it is downward from the origin, representing the vertical axis. z is the vertical position of the drilling depth in the local coordinate system, and z takes a value of 0~ H。 7. The method for geological boundary modeling based on function graph according to claim 6, characterized in that: The performing curve fitting on the second plurality of segments of data using a second fitting function includes: If the second multi-segment data contains three drilling data, then when fitting, the local coordinates of the middle drilling holes except the endpoints are brought into the second fitting function to obtain a value, and draw a fitting curve of the second plurality of segments of data; If the second multi-segment data contains 2 drilling data, then during fitting, the slope of the common point is calculated based on the fitting curve of the multi-segment data adjacent to the second multi-segment data, and the slope of the common point is calculated based on the slope of the common point. a values, and draw a fitting curve for the second plurality of segments of data.

8. A geological boundary modeling device based on function images, characterized in that: include: Input module, used to obtain drilling data of the same geological section; The segmentation module is used to divide the drilling data into a plurality of continuous segments along the same horizontal direction of the geological section, wherein the plurality of segments include a first plurality of segments and a second plurality of segments; wherein, The number of drilling data included in the first plurality of segments of data is within a first range, the number of drilling data included in the second plurality of segments of data is within a second range, and the first range and the second range do not intersect; a curve fitting module, performing curve fitting on the first plurality of segments of data and the second plurality of segments of data using different fitting functions, respectively, to obtain a fitting curve corresponding to each segment of the plurality of segments of data; The combination module is used to smoothly connect multiple segments of the fitting curve in sequence to obtain the geological boundary of the geological section.

9. An electronic device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein: When the processor executes the computer program, the function image-based geological boundary modeling method described in any one of claims 1 to 7 is implemented.

10. A computer-readable storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the function image-based geological boundary modeling method described in any one of claims 1 to 7 is implemented.

Images