Stratum mechanics parameter finite element modeling method and terminal equipment based on drilling data

Abstract

The application discloses a stratum mechanics parameter finite element modeling method based on drilling data and a terminal device, and comprises the following steps: 1) establishing a finite element model without stratum parameters; 2) collecting drilling data in an engineering site; 3) discretizing the drilling data and giving corresponding parameters to the finite element model; 4) solving the parameters of units in non-drilling positions according to the interaction relationship between units; and 5) updating and classifying the stratum parameters in non-drilling positions and then giving the parameters to the finite element model. The application avoids problems such as grid distortion caused by directly dividing stratum geometry in the finite element model, can generate different stratum parameters in the grid-uniform finite element model according to the existing drilling data, reflects different stratum fluctuation change rules, and thus realizes finite element numerical simulation calculation of strata with inclination.

Description

Technical Field

[0001] This invention relates to the field of finite element modeling, and in particular to a finite element modeling method and terminal equipment for formation mechanical parameters based on borehole data. Background Technology

[0002] With the advancement and development of computer technology, the finite element method (FEM) is playing an increasingly important role in civil engineering design and calculation. Finite element calculation models for underground structures are mainly divided into two categories: structural mechanics models based on relaxation load theory and rock mechanics models based on rock bearing theory. Rock mechanics models, in particular, require the underground structure to be established together with the surrounding strata during finite element modeling. Accurately reflecting the actual characteristics of the strata during strata modeling in rock mechanics models is of great significance.

[0003] Existing modeling methods often select one or several unfavorable cross-sections from geological maps obtained from borehole data for modeling and calculation. Specifically, during modeling and calculation, the undulation characteristics of the strata are ignored, and the strata are constructed as horizontal based on the average thickness before the underground structure is designed. This design method is acceptable for strata with minimal undulation, but it has two drawbacks when applied to strata with significant undulation: 1) The bias effect caused by strata undulation cannot be considered in cross-sectional structural design, leading to inaccurate calculation results; 2) The longitudinal uneven load caused by strata undulation cannot be reflected in longitudinal section design calculations. For example, when calculating the thrust of a shield tunnel, if the strata are considered horizontal, the resulting thrust is a constant value, which is clearly inconsistent with reality. Therefore, in cases of significant strata undulation, finite element modeling calculations must consider the longitudinal uneven distribution of the strata.

[0004] There are two main methods for existing finite element modeling that consider the undulation characteristics of strata: one is to construct the shape of the undulations during geometric modeling, and the other is a method based on random field theory (e.g., CN113158315A) that considers the spatial variability of strata mechanical parameters. The first method often encounters situations where certain strata shapes have sharp points, which can easily lead to mesh distortion and other problems during mesh generation, resulting in significant errors in the calculation results. While the second method can characterize the heterogeneity of strata distribution, its implementation is more complex, requiring a large amount of computation to obtain stable statistical results. Summary of the Invention

[0005] The technical problem this invention aims to solve is to provide a finite element method and terminal device for modeling formation mechanical parameters based on borehole data, addressing the shortcomings of existing technologies and avoiding the massive computational burden of random field methods. It also avoids the mesh distortion problem caused by directly dividing the formation in the finite element model from a geometric modeling perspective.

[0006] To solve the above-mentioned technical problems, the technical solution adopted by the present invention is: a finite element modeling method for formation mechanical parameters based on borehole data, comprising the following steps:

[0007] S1. Establish a finite element model without stratum parameters; collect borehole data from the engineering site; the borehole data from the engineering site includes the location and number of each borehole, the burial depth of each soil layer in each borehole, and the mechanical parameters.

[0008] S2. Discretize the field drilling data and assign corresponding parameters to the finite element model;

[0009] S3. Calculate the parameters of the non-drilling location elements based on the interaction relationships between the elements in the finite element model.

[0010] S4. Update and classify the formation parameters at non-drilling locations, and assign the updated formation parameters at non-drilling locations to the finite element model.

[0011] This invention enables finite element formation modeling using only borehole data, eliminating the need for geological cross-sections and avoiding the massive computational burden of random field methods. It is computationally simple and easy to use. Furthermore, this invention achieves multi-stratum, multi-parameter modeling in finite element calculations, allowing for calculations of formation undulations and more accurately reflecting the impact of formation inhomogeneity on the results.

[0012] In step S1, the specific implementation process of establishing a finite element model without formation parameters includes:

[0013] Set the geometric dimensions of the finite element model, complete the geometric modeling of the overall dimensions of the finite element model, and obtain the overall geometric model of the finite element model;

[0014] Mesh the overall geometric model of the finite element method and number the resulting elements and nodes;

[0015] Derive the correspondence between the units and nodes, as well as the coordinate values ​​of each node.

[0016] This invention addresses the problem of mesh distortion caused by dividing the strata into distinct layers during geometric modeling, which is prone to mesh distortion. Instead, it performs geometric modeling on the overall model size, without considering the strata undulation characteristics during the geometric modeling stage. Instead, it divides the model into layers on the overall model and assigns different strata parameters to uniform meshes. This approach not only reflects the strata undulation characteristics but also solves the mesh distortion problem caused by geometric strata division, greatly reducing the error in the calculation results and improving the calculation accuracy.

[0017] Preferably, for ease of subsequent calculation, the correspondence between the units and nodes is stored in matrix A, where the first column of matrix A is the unit number and the remaining columns are the node numbers; the first number in each row of matrix A is the unit number, and the remaining numbers are the node numbers corresponding to that unit.

[0018] Preferably, for ease of subsequent calculation, the coordinate values ​​of each node are stored in matrix B, where the first column of matrix B is the node number and the remaining columns are the coordinate values; wherein the first number in each row of matrix B is the node number, and the remaining numbers are the coordinate values ​​corresponding to that node number.

[0019] In step S2, the specific implementation process of discretizing the borehole data and assigning corresponding parameters to the finite element model includes:

[0020] Calculate the coordinates of the center points of each element in the finite element model;

[0021] Based on the coordinates of the center points of each unit, first find the unit whose plane coordinates (other than the burial depth coordinates) are equal to the plane coordinates of each borehole. If there is no unit with the same plane coordinates as a certain borehole, then find the unit that is closest to the plane coordinates of that borehole and form a matrix E with all the units found.

[0022] By comparing the soil depth within the borehole with the unit's depth coordinates, the parameters corresponding to the borehole depth are assigned to the matrix E.

[0023] After step S2, the finite element model has been transformed from a parameterless state to a finite element model with drilling parameters, that is, the elements located at the drilling position in the model have been assigned parameters.

[0024] The specific implementation process of step S3 includes:

[0025] Store the coordinates of the center point of each unit in matrix D;

[0026] The difference between matrix D and matrix E is used to obtain the set of center coordinates F of the non-drilling position units.

[0027] In this invention, the specific implementation process of step S4 includes:

[0028] The i-th element parameter f at the non-drilling location ij Subtract the parameter of each element at the drilling location from the parameter of the element at the drilling location, select the element with the smallest difference, and use the element with the smallest difference as the i-th element parameter at the non-drilling location, replacing f. ij ;

[0029] After traversing all non-drilling location units, the parameters of all non-drilling location units are updated. The updated parameters of each non-drilling location unit are one of the parameters of the drilling location units, thus completing the unit parameter classification, i.e., realizing the formation classification.

[0030] The updated element parameters are assigned to the finite element model according to the element number to obtain a finite element model with full formation parameters.

[0031] in, e mj p is the j-th mechanical parameter of the m-th element in the borehole location element matrix E; n is the total number of the borehole location element matrix E; im To calculate the weight of the m-th element parameter in the borehole location element matrix E when the i-th element parameter is not in the borehole location, dx im dy im Let dx represent the horizontal and vertical distances between the m-th cell at the borehole location and the i-th cell at the non-drilling location, respectively. il dy il , respectively, represent the horizontal and vertical distances between the l-th unit at the borehole location and the i-th unit at the non-drilling location, k is the scaling factor, and n is the number of units in the borehole location unit matrix E.

[0032] Step S4 improves the inherent inverse distance weighted average method by introducing a proportionality coefficient k when calculating weights. The improved method can fully consider the correlation of formation mechanical parameters in different directions. The inverse distance weighted average method is a geological concept, belonging to interpolation methods, which calculates the elevation of an unknown point using the elevations of several known points. Because elevation information has no obvious directional correlation, the inverse distance weighted average method is reasonably effective. However, formation mechanical parameter information differs from elevation information; most strata are distributed horizontally and have obvious directional correlation. In this invention, by changing the value of k, the correlation of strata in different directions (i.e., whether the strata are mainly distributed horizontally or vertically) can be reflected. When k equals 1, it indicates that the influence of distance in the x and y directions is equal. A k value greater than 1 indicates that the horizontal distance between elements has a stronger effect on the parameters than the vertical distance; a k value less than 1 indicates that the vertical distance between elements has a greater effect on the parameters than the horizontal distance. Since the horizontal correlation of strata is generally greater than the vertical correlation (i.e., strata have strong horizontal continuity and are mostly distributed along the horizontal direction), the value of k is generally taken as a number greater than 1.

[0033] When the finite element model is a three-dimensional model, p im Calculated using the following formula:

[0034]

[0035] The above formula allows for the establishment of a three-dimensional finite element mechanical parameter model based on borehole data, thus expanding the applicability of the method of this invention.

[0036] As an inventive concept, the present invention also provides a terminal device, including a memory, a processor, and a computer program stored in the memory; the processor executes the computer program to implement the steps of the method described above.

[0037] As an inventive concept, the present invention also provides a computer-readable storage medium having a computer program / instructions stored thereon; when the computer program / instructions are executed by a processor, they implement the steps of the method described above.

[0038] Compared with the prior art, the beneficial effects of the present invention are as follows:

[0039] 1) This invention realizes multi-stratum and multi-parameter modeling in finite element calculation, which can calculate the working conditions of stratum undulation and more accurately reflect the influence of stratum heterogeneity on the calculation results.

[0040] 2) This invention can perform finite element formation modeling using only borehole data, without the need for geological cross-section maps, making it more convenient to use;

[0041] 3) The method of the present invention reflects the undulation characteristics of the formation by assigning different formation parameters to a uniform grid, which solves the problem of grid distortion caused by directly geometrically dividing the formation in the finite element model.

[0042] 4) The formation parameters obtained by the method of this invention are deterministic values, and the parameters are unique for a single unit. Therefore, a stable result can be obtained by calculating only once. Compared with the random field method, which requires a large number of calculations to obtain statistically stable values, the computational cost is low.

[0043] 5) This invention avoids problems such as mesh distortion caused by directly dividing the formation geometry in the finite element model. It can generate different formation parameters in the finite element model with uniform mesh based on existing borehole data, thereby reflecting the variation law of different formation undulations and realizing finite element numerical simulation calculation of inclined formations. Attached Figure Description

[0044] Figure 1 This is a flowchart of the method steps in Embodiment 1 of the present invention;

[0045] Figure 2 This is a schematic diagram of the finite element model meshing and element node numbering in Embodiment 1 of the present invention.

[0046] Figures 3(a) and 3(b) show the finite element model with drilling parameters in Embodiment 1 of the present invention. Figure 3(a) shows the distribution of elastic modulus parameters of the finite element model, and Figure 3(b) shows the distribution of Poisson's ratio parameters of the finite element model.

[0047] Figures 4(a) and 4(b) show the finite element model of the initial soil parameters in Embodiment 1 of the present invention. Figure 4(a) shows the distribution of the elastic modulus parameters of the finite element model, and Figure 4(b) shows the distribution of the Poisson's ratio parameters of the finite element model.

[0048] Figures 5(a) and 5(b) show the finite element model of the final soil parameters in Embodiment 1 of the present invention. Figure 5(a) shows the distribution of the elastic modulus parameters of the finite element model, and Figure 5(b) shows the distribution of the Poisson's ratio parameters of the finite element model. Detailed Implementation

[0049] The technical solution of the present invention will be described in detail below with reference to the accompanying drawings and examples. It should be noted that, for the sake of simplicity and convenience in explaining the present invention, a two-dimensional engineering example has been selected for illustration, which does not mean that the present invention can only be used for two-dimensional finite element modeling. The same method can be used in three-dimensional finite element modeling.

[0050] Example 1

[0051] For a specific project, a two-dimensional finite element model needs to be established with a length of 50m and a burial depth of 50m. Within this range, there are two borehole data points. Borehole 1 information is as follows: the borehole is located 9m long, and within a burial depth of 50m, the soil and rock mass is divided into three layers with depths from top to bottom of 20m, 10m, and 20m, respectively. The elastic moduli are 100MPa, 20MPa, and 100MPa, and the Poisson's ratios are 0.3, 0.2, and 0.3, respectively. Borehole 2 information is as follows: the borehole is located 43m long, and within a burial depth of 50m, the soil and rock mass is divided into three layers with depths from top to bottom of 10m, 30m, and 20m, respectively. The elastic moduli are 100MPa, 20MPa, and 100MPa, and the Poisson's ratios are 0.3, 0.2, and 0.3, respectively.

[0052] The formation mechanical parameters finite element modeling method based on borehole data according to Embodiment 1 of the present invention is processed according to the following steps:

[0053] Step 1) Establish a finite element model without formation parameters:

[0054] A geometric model with a length of 50m (x-direction) and a burial depth of 50m (y-direction) was established based on the research problem. The geometric model was then meshed, and the elements (mesh elements) and nodes were numbered. To simply present the method of Embodiment 1 of this invention, a mesh size of 2.5m was chosen. The meshing and element node numbering are as follows... Figure 2 As shown, a total of 400 grids were divided, resulting in a finite element model without stratum parameters.

[0055] Export the element and node information, specifically including the correspondence between elements and nodes and the coordinate values ​​of each node;

[0056] To facilitate subsequent data processing, the correspondence between cells and nodes is stored in matrix A, where the first column of matrix A is the cell number and the remaining columns are the node numbers. In each row of matrix A, the first number is the cell number, and the remaining numbers are the node numbers corresponding to that cell. In this embodiment, matrix A is as follows:

[0057]

[0058] Preferably, to facilitate subsequent data processing, the coordinate values ​​of each node are stored in matrix B, where the first column of matrix B is the node number and the remaining columns are the coordinate values; the first number in each row of matrix B is the node number, and the remaining numbers are the coordinate values ​​corresponding to that node number.

[0059]

[0060] The main contents of the on-site data collection in step S2 include:

[0061] Location and number of each borehole, depth of each soil layer and mechanical parameters within each borehole;

[0062] Preferably, to facilitate subsequent batch operations, the burial depth and mechanical parameter information of each soil layer in each borehole can be represented by matrix C, where the first column of matrix C is the burial depth and thickness of each soil layer, and the remaining columns of matrix C are the mechanical parameters of the corresponding soil layers, such as unit weight, elastic modulus, Poisson's ratio, cohesion, internal friction angle, etc.

[0063] This embodiment has two boreholes, so the information is stored using two matrices C1 and C2 respectively, as follows:

[0064]

[0065]

[0066] The specific implementation process of step S3 includes:

[0067] The coordinate values ​​of the center points of each element in the finite element model in step S1 are calculated using matrices A and B as described in step S1. Specifically, the node number corresponding to each element is extracted from matrix A, and then the coordinate value of each node is found in matrix B based on the node number. The average of the x, y, and z coordinate values ​​of all nodes corresponding to each element is calculated (the average of the x and y coordinate values ​​in the two-dimensional model) to obtain the coordinate values ​​of the center points of each element. The element number and coordinate value information are stored in matrix D.

[0068] Based on the center point coordinate information (matrix D) of each element obtained above, firstly, find elements whose plane coordinates (excluding burial depth coordinates) are equal to the plane coordinates of each borehole. If no element has the same plane coordinates as a certain borehole, then find the element closest to that borehole plane coordinate, and form matrix E with the found elements. Subsequently, by comparing the soil depth within the borehole with the burial depth coordinates of the elements, assign the parameters of the corresponding borehole depth to the above-mentioned borehole location element matrix E, thereby completing the assignment of formation mechanical parameters to the borehole location elements in the meshed finite element model described in step S1.

[0069] After step S3, the finite element model has been transformed from the original parameterless state into a finite element model with drilling parameters. That is, the elements located at the drilling position in the model have been assigned parameters, as shown in Figures 3(a) and 3(b).

[0070] The specific implementation process of step S4 includes:

[0071] By taking the difference between the center coordinates of all cells (matrix D) in step S3 and the center coordinates of the borehole location cells (matrix E), we obtain the set of center coordinates F of the non-drillhole location cells.

[0072] The j-th mechanical parameter f of the i-th element in the set of non-drilling location center coordinates F ij Obtained by the following formula:

[0073]

[0074] In the above formula, e mj This refers to the j-th mechanical parameter of the m-th element in the borehole location element matrix E;

[0075] n is the total number of element matrix E.

[0076] p im The weight of the m-th element parameter in the borehole position element matrix E for the i-th element parameter at the non-drilling position is calculated using the following formula:

[0077]

[0078] In the above formula:

[0079] The subscript i represents the i-th element in the set of non-drilled location elements;

[0080] The subscripts m and l represent the m-th and l-th elements of the borehole location element matrix E, respectively;

[0081] dx and dy represent the horizontal and vertical distances between the two units, respectively.

[0082] k is a proportionality coefficient, reflecting the ratio of the influence intensity of the horizontal distance between elements to the vertical distance on soil parameters. When k equals 1, it indicates that the influence intensity of the distance in the x and y directions is the same. Since the correlation of the strata in the area where this example is located is much greater in the horizontal direction than in the vertical direction, k = 10 is taken.

[0083] Preferably, when the finite element model is a three-dimensional model, p im Calculated using the following formula:

[0084]

[0085] In the above formula, dx and dy represent the horizontal distance between two elements, and dz represents the vertical distance between two elements.

[0086] After step S4, the elements in the non-drilled locations of the finite element model have preliminary soil parameters. These preliminary parameters are generally different from the original parameters of the borehole. If the initial parameters are directly assigned to the finite element model, the results shown in Figure 4(a) and Figure 4(b) will be generated.

[0087] The specific implementation process of step S5 includes the following steps:

[0088] The i-th element parameter f at the non-drilling location ij Subtract the parameter from the parameter of each element at the borehole location, and select the element with the smallest absolute difference. Use the parameter of that borehole location as the parameter of the i-th element at the non-drilling location to replace f. ij This method iterates through all non-drilled location units, thus modifying the parameters of all non-drilled location units. The modified parameters of each non-drilled location unit are then one type of those for drilled location units, thereby completing unit parameter classification and achieving formation classification.

[0089] The updated element parameters are assigned to the finite element model according to the element number to realize the finite element model of all formation parameters. The finite element model of all formation parameters is shown in Figure 5(a) and Figure 5(b).

[0090] Example 2

[0091] Embodiment 3 of the present invention provides a terminal device corresponding to Embodiment 1 above. The terminal device can be a processing device for a client, such as a mobile phone, a laptop, a tablet computer, a desktop computer, etc., to execute the method of the above embodiments.

[0092] The terminal device in this embodiment includes a memory, a processor, and a computer program stored in the memory; the processor executes the computer program in the memory to implement the steps of the method in Embodiment 1 described above.

[0093] In some implementations, the memory may be high-speed random access memory (RAM), and may also include non-volatile memory, such as at least one disk storage device.

[0094] In other implementations, the processor can be any type of general-purpose processor, such as a central processing unit (CPU) or a digital signal processor (DSP), and there is no limitation here.

[0095] Example 3

[0096] Embodiment 3 of the present invention provides a computer-readable storage medium corresponding to the above embodiments, wherein a computer program / instructions are stored thereon. When the computer program / instructions are executed by a processor, they implement the steps of the method of Embodiment 1 described above.

[0097] A computer-readable storage medium can be a tangible device that holds and stores instructions for use by an instruction execution device. A computer-readable storage medium can be, for example, but not limited to, an electrical storage device, a magnetic storage device, an optical storage device, an electromagnetic storage device, a semiconductor storage device, or any combination thereof.

[0098] Those skilled in the art will understand that embodiments of this application can be provided as methods, systems, or computer program products. Therefore, this application can take the form of a completely hardware embodiment, a completely software embodiment, or an embodiment combining software and hardware aspects. Furthermore, this application can take the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) containing computer-usable program code. The solutions in the embodiments of this application can be implemented in various computer languages, such as the object-oriented programming language Java and the interpreted scripting language JavaScript.

[0099] This application is described with reference to flowchart illustrations and / or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of this application. It will be understood that each block of the flowchart illustrations and / or block diagrams, and combinations of blocks in the flowchart illustrations and / or block diagrams, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, special-purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, generate instructions for implementing the flowchart... Figure 1 One or more processes and / or boxes Figure 1 A device that provides the functions specified in one or more boxes.

[0100] These computer program instructions may also be loaded onto a computer or other programmable data processing equipment to cause a series of operational steps to be performed on the computer or other programmable equipment to produce a computer-implemented process, thereby providing instructions that execute on the computer or other programmable equipment for implementing the process. Figure 1 One or more processes and / or boxes Figure 1 The steps of the function specified in one or more boxes.

[0101] Although preferred embodiments of this application have been described, those skilled in the art, upon learning the basic inventive concept, can make other changes and modifications to these embodiments. Therefore, the appended claims are intended to be interpreted as including the preferred embodiments as well as all changes and modifications falling within the scope of this application.

[0102] Obviously, those skilled in the art can make various modifications and variations to this application without departing from the spirit and scope of this application. Therefore, if such modifications and variations fall within the scope of the claims of this application and their equivalents, this application also intends to include such modifications and variations.

Claims

1. A finite element method for modeling formation mechanical parameters based on borehole data, characterized in that, Includes the following steps: S1. Establish a finite element model without stratum parameters; collect borehole data from the engineering site; the borehole data from the engineering site includes the location and number of each borehole, the burial depth of each soil layer in each borehole, and the mechanical parameters. S2. Discretize the field drilling data and assign corresponding parameters to the finite element model; S3. Calculate the parameters of the non-drilling location elements based on the interaction relationships between the elements in the finite element model. S4. Update and classify the formation parameters at non-drilling locations, and assign the updated formation parameters at non-drilling locations to the finite element model. The specific implementation process of step S4 includes: The i-th element parameter at the non-drilling location Subtract the parameter of each element at the drilling location from the parameter of the element at that location, select the element with the smallest difference, and use the element with the smallest difference at the drilling location as the parameter of the i-th element at the non-drilling location, replacing the parameter of that element. ; After traversing all non-drilling location units, the parameters of all non-drilling location units are updated. The updated parameters of each non-drilling location unit are one of the parameters of the drilling location units, thus completing the unit parameter classification, i.e., realizing the formation classification. The updated element parameters are assigned to the finite element model according to the element number to obtain a finite element model with full formation parameters. in, ; This refers to the j-th mechanical parameter of the m-th element in the borehole location element matrix E; This represents the total number of borehole location element matrix E; To calculate the weight of the m-th element parameter in the borehole location element matrix E when the i-th element parameter is not in the borehole location, dx im dy im Let dx represent the horizontal and vertical distances between the m-th cell at the borehole location and the i-th cell at the non-drilling location, respectively. il dy il These represent the horizontal and vertical distances between the l-th unit at the borehole location and the i-th unit outside the borehole location, respectively. is the scaling factor, and n is the number of cells in the borehole location cell matrix E.

2. The finite element modeling method for formation mechanical parameters based on borehole data according to claim 1, characterized in that, In step S1, the specific implementation process of establishing a finite element model without formation parameters includes: Set the geometric dimensions of the finite element model, complete the geometric modeling of the overall dimensions of the finite element model, and obtain the overall geometric model of the finite element model; Mesh the overall geometric model of the finite element method and number the resulting elements and nodes; Derive the correspondence between the units and nodes, as well as the coordinate values ​​of each node.

3. The finite element modeling method for formation mechanical parameters based on borehole data according to claim 2, characterized in that, The correspondence between the units and nodes is stored in matrix A. The first column of matrix A is the unit number, and the remaining columns are the node numbers. The first number in each row of matrix A is the unit number, and the remaining numbers are the node numbers corresponding to that unit.

4. The finite element modeling method for formation mechanical parameters based on borehole data according to claim 2, characterized in that, The coordinate values ​​of each node are stored in matrix B. The first column of matrix B is the node number, and the remaining columns are the coordinate values. In each row of matrix B, the first number is the node number, and the remaining numbers are the coordinate values ​​corresponding to that node number.

5. The finite element modeling method for formation mechanical parameters based on borehole data according to claim 1, characterized in that, In step S2, the specific implementation process of discretizing the borehole data and assigning corresponding parameters to the finite element model includes: Calculate the coordinates of the center points of each element in the finite element model; Based on the coordinates of the center points of each unit, find the unit whose planar coordinates are equal to the planar coordinates of each borehole. If there is no unit with the same planar coordinates as a certain borehole, find the unit that is closest to the planar coordinates of that borehole and form a matrix E with all the units found. By comparing the soil depth within the borehole with the unit's depth coordinates, the parameters corresponding to the borehole depth are assigned to the matrix E.

6. The finite element modeling method for formation mechanical parameters based on borehole data according to claim 5, characterized in that, The specific implementation process of step S3 includes: Store the coordinates of the center point of each unit in matrix D; The difference between matrix D and matrix E is used to obtain the set of center coordinates F of the non-drilling position units.

7. The finite element modeling method for formation mechanical parameters based on borehole data according to claim 1, characterized in that, When the finite element model is a three-dimensional model Calculated using the following formula: 。 8. A terminal device, comprising a memory, a processor, and a computer program stored in the memory; characterized in that, The processor executes the computer program to implement the steps of the method according to any one of claims 1 to 7.

9. A computer-readable storage medium having a computer program / instructions stored thereon; characterized in that, When the computer program / instructions are executed by the processor, they implement the steps of the method according to any one of claims 1 to 7.

Images