Method for constructing lenticular body three-dimensional model

Abstract

The invention discloses a method for constructing a lenticular body three-dimensional model, and relates to the technical field of three-dimensional geological modeling. The method comprises the following steps: reading spatial data and geological data; marking the exposure times and sequence of each lithologic section of the lenticular body, and calculating a horizontal maximum boundary; generating an initial three-dimensional model of each lens lithologic section; regenerating a three-dimensional model to adjust the size of the lenticular body; performing model inspection, and regenerating a three-dimensional model; providing a list or screening each lithologic section of the lenticules with the same lithology, selecting lithologic sections to be combined, and generating an inner boundary and an outer boundary located on the top surface and the bottom surface of the lenticules; generating top and bottom curved surfaces of each lens body; determining an interpolation function and a function parameter in a space, generating a side curved surface, and finally forming an integrated three-dimensional model; and carrying out model inspection, and regenerating a three-dimensional entity model. According to the method, automatic modeling and manual intervention are combined, and the accuracy and geological authenticity of the lenticular body three-dimensional model are remarkably improved.

Description

Technical Field

[0001] The invention relates to the technical field of three-dimensional geological modeling, and in particular to a method for constructing a three-dimensional model of a lens. Background Art

[0002] Lens modeling is an important area of ​​3D modeling under complex geological conditions. Traditional 3D modeling techniques often have difficulty accurately representing complex stratigraphic sequence relationships and the spatial distribution of lenses. Manual modeling using borehole geological data of lenses is not only inefficient, but also difficult to maintain consistency and accuracy. More importantly, in the process of automated modeling, the problem of spatial multi-solutions is often encountered, that is, the same geological data set may produce multiple reasonable 3D models.

[0003] Traditional modeling methods often have problems when dealing with complex situations such as one-to-many or many-to-many where a single lens appears in multiple boreholes. These methods usually ignore the construction of the inner boundary of the lens, resulting in significant defects in the scientificity and accuracy of the model. In addition, when multiple solutions appear in the lens modeling process, traditional methods lack effective automated rules and manual intervention mechanisms to optimize and verify the model, which ultimately limits the practicality and reliability of the model results and affects the efficiency of actual geological work. Summary of the invention

[0004] The present invention aims at solving the problems existing in the prior art and provides a method for constructing a three-dimensional model of a lens body.

[0005] To achieve the above purpose, the technical solution adopted by the present invention is as follows:

[0006] A method for constructing a three-dimensional model of a lens body, comprising the following steps: step S1, reading the spatial data and geological data of the lens body lithology segments in all boreholes stored in a database; step S2, marking the number and order of exposure of each lens body lithology segment at the spatial position of all boreholes according to the geological data; calculating the maximum horizontal boundary allowed for each lens body lithology segment according to the spatial data; step S3, generating the initial three-dimensional model of each lens body lithology segment according to the maximum horizontal boundary, the top elevation and bottom elevation data of the lens body lithology segment in the borehole, and automatically naming it according to the lithology and order; step S4, selecting the initial three-dimensional model of each lens body lithology segment, setting a value less than the maximum horizontal boundary allowed, and regenerating the initial three-dimensional model to adjust the size of the lens body; step S5, checking the initial three-dimensional model, and recording the initial three-dimensional model. Nested relationship or intersection relationship, perform Boolean operation, and regenerate the initial three-dimensional model; when the scale of the lens is larger than the initial three-dimensional model, operate steps S6-S9; step S6, provide a list or screen the various lithology segments of the lens with the same lithology in the cross-section view or three-dimensional view; select the lithology segments to be merged, and generate the inner boundary and outer boundary located on the top and bottom surfaces of the lens; step S7, generate the top and bottom surfaces of each lens by interpolation calculation; determine the interpolation function and function parameters in space according to the lens morphology; step S8, form the side surface of the lens based on the top surface boundary, bottom surface boundary and interpolation function; form a three-dimensional model after the collection based on the top surface, bottom surface and side surface; step S9, perform a three-dimensional model check, record the nested relationship or intersection relationship of the model, perform Boolean operation, and regenerate the three-dimensional solid model.

[0007] Based on the above technical solution, further, in step S1, the reading process is: first determine the storage format and interface of the required data in the database; then write the corresponding query statement according to the selected database to extract the spatial coordinate data and lithology segment data of the required borehole; wherein the spatial coordinate data is the key data of the spatial data, and the lithology segment data is the key data of the geological data.

[0008] Based on the above technical solution, further, in step S2, all boreholes are traversed, and the lithology data in the geological data of each borehole is analyzed; the number of exposures of different lens lithology sections in each borehole is counted, and the order of exposure is recorded.

[0009] Based on the above technical solution, further, for the same lens lithology segment, if it appears multiple times in the same borehole, the number of times it appears is recorded; for different lens lithology segments that appear in the same borehole, the order of their appearance is recorded.

[0010] Based on the above technical solution, further, in step S2, a convex hull algorithm or an Alpha Shapes algorithm is used to calculate the maximum horizontal boundary of each lens; wherein the elevation of the maximum horizontal boundary is within the range of the top and bottom burial depths of the lithologic section.

[0011] Based on the above technical solution, further, in step S3, a parabolic function is generated on each cross section on the horizontal maximum boundary using the top elevation point and the bottom elevation point, and an initial three-dimensional model of the lens is generated along the cross section.

[0012] Based on the above technical solution, further, in step S4, the process of adjusting the size of the lens body is: providing an interface or script to adjust the horizontal boundary of the lens body, and recalculating and generating an adjusted initial three-dimensional model.

[0013] Based on the above technical solution, further, in step S6, the process of screening and merging lithology segments is: developing an interactive interface, using visualization tools to display lists, profiles or three-dimensional views, and merging lithology segments by manually screening and selecting them.

[0014] Based on the above technical solution, further, in step S7, the process of determining the interpolation function and function parameters is: according to the geological characteristics and lens morphology, the most suitable interpolation function is selected, and the parameters are optimized.

[0015] Based on the above technical solution, further, in step S8, several dividing lines are extracted along the top and bottom surface boundaries as sweeping paths; the interpolation function is calculated along each sweeping path to obtain the side surface on the path; all the side surfaces are connected in order to form a complete side surface of the lens body.

[0016] Compared with the prior art, the present invention has the following beneficial effects:

[0017] (1) The present invention proposes a method for constructing a three-dimensional model of a lens. By comprehensively utilizing drilling data and geological data, combined with automated modeling and manual intervention, the present invention solves several problems in the traditional method. Specifically, the present method calculates the maximum horizontal boundary of each lens and generates an initial three-dimensional model based on the number of times the lens lithology segment is exposed. By adjusting the boundaries and Boolean operations, the nesting and intersection relationships of the model are accurately processed. In addition, an interactive interface is provided to allow geological experts to screen and merge lithology segments, regenerate the top and bottom surfaces of the lens, and finally form a high-precision three-dimensional solid model.

[0018] (2) The present invention significantly improves the accuracy and geological authenticity of the lens 3D model by comprehensively utilizing drilling data and geological data, combining automated modeling with manual intervention. Specifically, by marking the number and sequence of exposures of the lens lithology segments and calculating the maximum horizontal boundary, the spatial distribution of the model is ensured to be accurate.

[0019] (3) The present invention can effectively handle complex situations such as one-to-many and many-to-many situations of a single lens in multiple boreholes, and consider the construction of the inner boundary of the lens to generate a more scientific and accurate model. By screening and merging lithological segments, the inner and outer boundaries of the top and bottom surfaces of the lens are generated to ensure the integrity and rationality of the lens model.

[0020] (4) The present invention can effectively handle the spatial multi-solution problem in the modeling process by setting up automated rules and manual intervention mechanisms to ensure the consistency and reliability of the model. When multiple solutions occur, model checking and Boolean operations are performed to record and correct the nested or intersecting relationships of the model.

[0021] (5) The present invention optimizes the model verification and adjustment mechanism, provides an interactive interface and visualization tools, enables geological experts to intuitively screen, view and adjust the model, and enhances the operability and practicality of the model. Visualization tools are used to display cross-sections and three-dimensional views, and an interactive interface is provided for model screening and adjustment, ensuring the practicality and reliability of the model results. BRIEF DESCRIPTION OF THE DRAWINGS

[0022] Figure 1 This is an overall flow chart of constructing a three-dimensional model of a lens body in an embodiment of the present invention;

[0023] Figure 2 The spatial distribution of the boreholes and the exposed positions of the top and bottom surfaces of each lithology section in the borehole in the embodiment of the present invention;

[0024] Figure 3 is the initial three-dimensional model in the embodiment of the present invention;

[0025] Figure 4 An interface for selecting and merging initial three-dimensional models in an embodiment of the present invention;

[0026] Figure 5 The top surface and bottom surface boundaries and interpolation functions recalculated for the merged lens body in the embodiment of the present invention;

[0027] Figure 6 The top curved surface and the bottom curved surface are formed by interpolation calculation in the embodiment of the present invention;

[0028] Figure 7 is a side curved surface formed by sweeping an interpolation function in an embodiment of the present invention;

[0029] Figure 8 It is a three-dimensional solid model finally formed by manual intervention and merging in the embodiment of the present invention. DETAILED DESCRIPTION

[0030] The present invention is further described and illustrated below in conjunction with the accompanying drawings and specific embodiments. The technical features of each embodiment of the present invention can be combined accordingly without conflicting with each other.

[0031] In order to make the above-mentioned purpose, features and advantages of the present invention more obvious and easy to understand, the specific implementation mode of the present invention is described in detail below in conjunction with the accompanying drawings. In the following description, many specific details are set forth to facilitate a full understanding of the present invention. However, the present invention can be implemented in many other ways different from those described herein, and those skilled in the art can make similar improvements without violating the connotation of the present invention. Therefore, the present invention is not limited to the specific embodiments disclosed below. The technical features in each embodiment of the present invention can be combined accordingly without conflicting with each other.

[0032] In this method, unless the context clearly dictates otherwise, the following terms have the following meanings: Lens: A lithologic body with unique geological characteristics that can be identified in one or more boreholes and has a certain spatial distribution and morphology. Lithologic segment: A continuous portion of a specific rock type recorded in a borehole, which can be a single lithology or a complex composed of multiple lithologies. Each lithologic segment may represent a lens independently or a part of a larger lens. Borehole data: includes the location coordinates (x, y, z) of the borehole, as well as the top elevation (stratum_top) and bottom elevation (stratum_bottom) of the lithologic segment recorded in the borehole. Top surface: The upper surface of the lens, usually determined by the top elevation data of the lithologic segment recorded in the borehole. Bottom surface: The lower surface of the lens, usually determined by the bottom elevation data of the lithologic segment recorded in the borehole. Side surface: A vertical or inclined surface connecting the top and bottom boundaries, forming the three-dimensional structure of the lens. Maximum horizontal boundary: the maximum extended boundary of the lens on the horizontal plane determined by spatial analysis, usually calculated using a geometric algorithm (such as a convex hull algorithm or an Alphashapes algorithm).

[0033] Example 1

[0034] Combination Figure 1 As shown, this embodiment provides a method for constructing a three-dimensional model of a lens body. It should be noted that: Figure 2-Figure 8In the figure, A1 is the borehole number; A2 is the borehole space; A3 is the top elevation of the lithology segment; A4 is the bottom elevation of the lithology segment; A5 is the maximum horizontal radius; A6 is the maximum horizontal boundary; A7 is a parabola; A8 is the initial three-dimensional model; A9 is the merge interface; A10 is the lithology list; A11 is the lithology segment lens list box; A12 is the visualization window frame; A13 is the top surface boundary, A14 is the bottom surface boundary; A15 is the interpolation function curve; A16 is the top surface; A17 is the bottom surface; A18 is the side surface; and A19 is the assembled three-dimensional model. The method includes the following steps.

[0035] Step S1, reading the spatial data and geological data of the lens lithology sections in all boreholes stored in the database.

[0036] In this embodiment, when reading the database, SQL Server or other database query languages ​​are used to extract spatial data and geological data, including the spatial coordinates X, Y, and Z of each borehole; the top and bottom burial depths of each lens lithology section. Specifically, the storage format and interface of the relevant data in the database must be determined first, such as using a relational database SQL Server. Then, the corresponding query statements are written according to the selected database to extract the required borehole spatial coordinate data and lithology section data, such as hole_id, x_coordinate, y_coordinate, elevation, lithology, stratum_top, stratum_bottom, etc. The extracted raw data is preprocessed to remove outliers, and the data is converted into a unified format for subsequent use, such as using Python's Pandas library.

[0037] like Figure 2 As shown, Figure 2 The spatial distribution of the boreholes in this embodiment and the exposed positions of the top and bottom surfaces of each lithology segment in the borehole are shown. The information included mainly includes the borehole number A1, namely hole_id (for example, ZK1-ZK4), borehole space information A2 (including coordinates, namely x_coordinate, y_coordinate, elevation, etc.), the top elevation A3 (stratum_top) and the bottom elevation A4 (stratum_bottom) of a lithology segment in the borehole.

[0038] Step S2: Mark the number and order of exposure of each lens lithology segment at the spatial position of all boreholes according to the geological data; calculate the maximum horizontal boundary allowed for each lens lithology segment according to the spatial data.

[0039] In this embodiment, all boreholes are traversed, and the lithology data in the geological data of each borehole are analyzed. The number of exposures of different lens lithology segments in each borehole is counted, and the order of exposure is recorded, for example, using Pandas's data.size() and data.cumcount() functions, and finally a data set containing spatial coordinates, lithology exposure times and order is formed.

[0040] It should be noted that if the same lens lithology segment appears multiple times in the same borehole, the number of times it appears should be recorded; if different lens lithology segments appear in the same borehole, the order in which they appear should be recorded.

[0041] Calculate the maximum horizontal boundary of each lens using the Convex Hull algorithm or the Alpha Shapes algorithm. For example, use the convex_hull class / function in the Shapely library. The calculated maximum horizontal radius A5 or maximum horizontal boundary A6 can be found in Figure 2 It should be further explained that the elevation of the maximum horizontal boundary is within the range of the top and bottom burial depths of the lithologic section.

[0042] Step S3: Generate an initial three-dimensional model A8 of each lens lithology segment according to the maximum horizontal boundary, the top elevation and the bottom elevation of the lens lithology segment in the borehole, and automatically name them according to the lithology and order.

[0043] In this embodiment, a parabola function is generated using the top elevation point and the bottom elevation point on each cross section on the horizontal maximum boundary, and an initial three-dimensional model A8 of the lens is generated along the cross section. The initial three-dimensional model A8 of each lens is automatically named according to the lithology and sequence. For example, the NumPy library and the Open3d library are used. The generated parabola A7 function is shown in Figure 2 .like Figure 3 As shown, Figure 3 The figure shows an initial three-dimensional model A8 of a lens body initially generated from four drilling data in a typical case of an embodiment of the method.

[0044] It should be noted that the process of generating the parabolic function is: generate a lens based on the top and bottom elevation data and the maximum horizontal boundary, extract the top and bottom elevation points of the borehole and the points on the maximum horizontal boundary, and select Gaussian basis functions, polynomial basis functions, etc. as the basis function family of the interpolation function according to the lens shape.

[0045] Step S4: Select the initial three-dimensional model A8 of each lens lithology segment, manually intervene to set a value smaller than the maximum allowed horizontal boundary, and regenerate the initial three-dimensional model A8 to change the lens size.

[0046] In this embodiment, the horizontal boundary value of the initial three-dimensional model A8 of the lens body is adjusted, and the initial three-dimensional model A8 is regenerated to change the size of the lens body. An interface or script is provided to allow the user to adjust the horizontal boundary value of the lens body. The adjusted initial three-dimensional model A8 is recalculated and generated. Specifically, the process of adjusting the size of the lens body is: providing an interface or script to adjust the horizontal boundary of the lens body, and recalculating and generating the adjusted initial three-dimensional model A8.

[0047] It should be further explained that the initial three-dimensional model A8 here is a model formed by the envelope of the top elevation (point), the bottom elevation (point), and the maximum boundary (line). These models are independent or discrete three-dimensional models.

[0048] Step S5, check the initial three-dimensional model A8, record the nesting relationship or intersection relationship of the initial three-dimensional model A8, perform Boolean operation, and regenerate the initial three-dimensional model A8. In this embodiment, a Boolean operation library is used to check the nesting and intersection relationship of the model for correcting the three-dimensional model. For example, the Polygon class and the Unary_Union class / function in the Shapely library are used. It should be noted that the Boolean operation is an existing method and will not be recorded in detail here.

[0049] Furthermore, when the scale of the lens is larger than the initial three-dimensional model A8, that is, a single lens contains multiple lithology segments, the lithology segments need to be merged, and steps S6 to S9 are performed.

[0050] Step S6: Provide a list or select the lithology segments of the lens with the same lithology in the cross-section or three-dimensional view; select the lithology segments to be merged according to the actual situation, and generate the inner and outer boundaries on the top and bottom surfaces of the lens. It should be noted that the generated isolated and merged lens models are displayed here in the form of one or more of a text list, cross-section, or three-dimensional view, such as Figure 4 shown.

[0051] In this embodiment, an interactive interface is developed to allow users to filter and view lens lithology segments of the same lithology in a list, a cross section, or a three-dimensional view. A visualization tool is used to display the cross section or three-dimensional view, for example, using the Plotly library.

[0052] like Figure 4 As shown, Figure 4 The interface A9 for selecting and merging the three-dimensional lens model initially formed in this embodiment, the lithology list box A10 to be viewed and selected, the lithology segment lens list box A11 to be viewed and selected, and the visualization window frame A12 for displaying the section or three-dimensional view are displayed.

[0053] Specifically, the steps of screening and merging lithologic segments include developing an interactive interface that uses visualization tools to display lists, sections, or 3D views, allowing geologists to screen and select lithologic segments for merging. Figure 4 In the interactive interface shown in the figure, select the lithology segment to be merged. Use geometric algorithms (such as Voronoi diagram or polygon Boolean operation) to calculate the inner and outer boundaries after merging. Figure 5 As shown, Figure 5 The top surface boundary A13, the bottom surface boundary A14 and the interpolation function curve A15 of the lens body recalculated after merging in this embodiment are shown.

[0054] Furthermore, the step of generating top and bottom surfaces includes using an interpolation method such as Kriging interpolation to generate the top and bottom surfaces based on the selected lithology segments and boundary data.

[0055] Step S7, generate the top and bottom surfaces of each lens by interpolation calculation; determine the interpolation function and function parameters in space according to the lens shape. Specifically, according to the selected lithology segment and boundary data, use interpolation methods such as Kriging interpolation to generate the top and bottom surfaces. For example, use the SciPy library or dedicated geological modeling software to implement interpolation calculation. Figure 6 As shown, Figure 6 The top curved surface A16 and the bottom curved surface A17 formed by interpolation calculation in this embodiment are shown.

[0056] In this embodiment, the process of determining the interpolation function and parameters is: according to the geological characteristics and lens morphology, combined with experience or other data, the most suitable interpolation function is selected, and the parameters are optimized. Further, the interpolation function and function parameters can select Gaussian basis functions, polynomial basis functions, etc. as the basis function family of the interpolation function according to the lens morphology. According to the basis function, a specific interpolation function form is fitted, and according to the provided interface, the basis function form can be set and the function parameters can be modified.

[0057] Step S8: Based on the top surface boundary, the bottom surface boundary and the interpolation function, the side surface of the lens body is formed; and a three-dimensional model after the top surface, the bottom surface and the side surface are assembled. It should be noted that the boundary here includes the inner boundary and the outer boundary, and only the outer boundary is involved in the example.

[0058] In this embodiment, a sweeping method is used to extract several segmentation lines along the top and bottom surface boundaries as sweeping paths; the interpolation function is calculated along each sweeping path to obtain the side curved surface on the path. All side curved surfaces are connected in order to form a complete side curved surface of the lens body. Figure 7 As shown, Figure 7 The side curved surface A18 formed by sweeping the interpolation function in this embodiment is shown.

[0059] Furthermore, the data of the top surface, bottom surface and side surface are integrated into the same 3D model representation; the surface data is discretized into an approximate entity composed of small planes through triangulation or other methods; operations such as coloring the entity and adding lighting are performed to obtain a 3D solid model of the lens with good rendering effect. Figure 8 As shown, Figure 8 The three-dimensional lens model A19 is shown as a result of the manual merging and finally forming a collection in this embodiment.

[0060] Step S9, check the assembled three-dimensional model. If it forms a nested or intersecting relationship with the three-dimensional models of other lens bodies, the three-dimensional model status should be recorded, Boolean operations should be performed, and the assembled three-dimensional model should be regenerated. Specifically, the Boolean operation library is used to check and record the nesting and intersection relationships of the models; the Polygon class and Unary_Union class / function in the Shapely library are used. It should be noted that the early three-dimensional models are relatively independent or discrete, and ultimately multiple discrete three-dimensional models are assembled to generate a new three-dimensional solid model.

[0061] Finally, it should be noted that the above content is only used to illustrate the technical solution of the present invention, rather than to limit the scope of protection of the present invention. Simple modifications or equivalent substitutions of the technical solution of the present invention by ordinary technicians in this field do not deviate from the essence and scope of the technical solution of the present invention.

Claims

1. A method for constructing a three-dimensional model of a lens, characterized in that: The following steps are involved: Step S1, reading the spatial data and geological data of the lens lithology sections in all boreholes stored in the database; Step S2: Mark the number and order of exposure of each lens lithology segment at the spatial position of all boreholes according to the geological data; calculate the maximum horizontal boundary allowed for each lens lithology segment according to the spatial data; Step S3, generating an initial three-dimensional model of each lens lithology segment according to the horizontal maximum boundary, the top elevation and the bottom elevation data of the lens lithology segment in the borehole, and automatically naming them according to the lithology and order; Step S4, selecting the initial three-dimensional model of each lens lithology segment, setting a value smaller than the maximum allowed horizontal boundary, and regenerating the initial three-dimensional model to adjust the lens size; Step S5, checking the initial three-dimensional model, recording the nesting relationship or intersection relationship of the initial three-dimensional model, performing Boolean operations, and regenerating the initial three-dimensional model; When the scale of the lens is larger than the initial three-dimensional model, the steps S6 to S9 are performed; Step S6, providing a list or screening various lithology segments of lenses of the same lithology in a cross-section or three-dimensional view; selecting the lithology segments to be merged, and generating inner boundaries and outer boundaries located on the top and bottom surfaces of the lenses; Step S7, generating the top and bottom surfaces of each lens body by interpolation calculation; determining the interpolation function and function parameters in space according to the lens body shape; Step S8, forming a side curved surface of the lens body based on the top surface boundary, the bottom surface boundary and the interpolation function; forming a three-dimensional model after being assembled according to the top curved surface, the bottom curved surface and the side curved surface; Step S9: Check the three-dimensional model, record the nesting relationship or intersection relationship of the model, perform Boolean operations, and regenerate the three-dimensional solid model.

2. A method for constructing a three-dimensional model of a lens according to claim 1, characterized in that: In step S1, the reading process is: first determine the storage format and interface of the required data in the database; then write the corresponding query statement according to the selected database to extract the spatial coordinate data and lithology segment data of the required borehole; among which, the spatial coordinate data is the key data of the spatial data, and the lithology segment data is the key data of the geological data.

3. The method for constructing a three-dimensional model of a lens according to claim 1, characterized in that: In step S2, all boreholes are traversed, and the lithology data in the geological data of each borehole is analyzed; the number of times different lens lithology segments are exposed in each borehole is counted, and the order of exposure is recorded.

4. A method for constructing a three-dimensional model of a lens according to claim 3, characterized in that: If the same lens lithology segment appears multiple times in the same borehole, the number of times it appears is recorded; if different lens lithology segments appear in the same borehole, the order of their appearance is recorded.

5. The method for constructing a three-dimensional model of a lens according to claim 1, characterized in that: In step S2, the maximum horizontal boundary of each lens is calculated using a convex hull algorithm or an Alpha Shapes algorithm; wherein the elevation of the maximum horizontal boundary is within the range of the top and bottom burial depths of the lithologic segment.

6. The method for constructing a three-dimensional model of a lens according to claim 1, characterized in that: In step S3, a parabolic function is generated using the top elevation point and the bottom elevation point on each cross section on the horizontal maximum boundary, and an initial three-dimensional model of the lens is generated along the cross section.

7. The method for constructing a three-dimensional model of a lens according to claim 1, characterized in that: In step S4, the process of adjusting the size of the lens body is: providing an interface or script to adjust the horizontal boundary of the lens body, and recalculating and generating an adjusted initial three-dimensional model.

8. The method for constructing a three-dimensional model of a lens according to claim 1, characterized in that: In step S6, the process of screening and merging lithology segments is as follows: based on an interactive interface, using visualization tools to display lists, sections or three-dimensional views, and merging lithology segments by manually screening and selecting them.

9. The method for constructing a three-dimensional model of a lens according to claim 1, characterized in that: In step S7, the process of determining the interpolation function and function parameters is as follows: selecting the interpolation function according to the geological characteristics and lens morphology, and optimizing the parameters.

10. The method for constructing a three-dimensional model of a lens according to claim 1, characterized in that: In step S8, a number of splitting lines are extracted along the top and bottom surface boundaries as sweep paths; an interpolation function is calculated along each sweep path to obtain a side surface on the path; and all side surfaces are connected in order to form a complete side surface of the lens body.

Images