CN-122021115-A - Automatic generation method and device for multi-stratum continuous model for numerical simulation

CN122021115ACN 122021115 ACN122021115 ACN 122021115ACN-122021115-A

Abstract

The application provides an automatic generation method and device of a multi-stratum continuous model for numerical simulation, wherein the method comprises the steps of obtaining geological identification information of an area to be molded, extracting initial demarcation point sets of stratums corresponding to lithology in the area to be molded according to the geological identification information, processing the demarcation point sets of stratums corresponding to lithology by adopting a spatial interpolation method to obtain a prediction demarcation point set of stratums corresponding to lithology in the area to be molded, constructing a single stratum model of stratums corresponding to lithology by adopting a single stratum body topology processing method and the prediction demarcation point set of stratums corresponding to lithology in the area to be molded, and processing the single stratum model of stratums by adopting a multi-stratum body topology processing method to obtain the multi-stratum continuous model of the area to be molded. The technical scheme provided by the application realizes automatic, efficient and high-precision generation from geological data to a numerical simulation model.

Inventors

LIU YIQIN
LI HAITAO
DU WEISHENG
CAI SHAOYANG

Assignees

煤炭科学研究总院有限公司

Dates

Publication Date: 20260512
Application Date: 20251224

Claims (10)

1. An automatic generation method of a multi-stratum continuous model for numerical simulation, which is characterized by comprising the following steps: Obtaining geological identification information of an area to be modeled, and extracting an initial demarcation point set of each stratum corresponding to each lithology in the area to be modeled according to the geological identification information, wherein the geological identification information comprises geological age, coal seam or rock stratum information; processing the demarcation point set of each stratum corresponding to each lithology by adopting a spatial interpolation method to obtain a prediction demarcation point set of each stratum corresponding to each lithology in the region to be modeled; adopting a single stratum body topology processing method and a prediction demarcation point set of each stratum corresponding to each lithology in the region to be modeled to construct a single stratum model of each stratum corresponding to each lithology; and processing the single stratum model of each stratum corresponding to each lithology by adopting a multi-stratum topological processing method to obtain the multi-stratum continuous model of the region to be modeled.
2. The method of claim 1, wherein the extracting the initial set of demarcation points for each formation corresponding to each lithology within the region to be modeled based on the geological identification information comprises: Extracting plane coordinate system coordinate values and elevation values of each top surface boundary point of each stratum corresponding to each lithology in the region to be molded and plane coordinate system coordinate values and elevation values of each bottom surface boundary point based on drilling holes and the geological identification information, and constructing initial boundary point sets of each stratum corresponding to each lithology in the region to be molded based on plane coordinate system coordinate values and elevation values of each top surface boundary point of each stratum corresponding to each lithology in the region to be molded and plane coordinate system coordinate values and elevation values of each bottom surface boundary point, wherein the initial boundary point sets comprise a top surface initial boundary point set and a bottom surface initial boundary point set; Judging whether stratum tip vanishing points exist in each stratum corresponding to each lithology in the area to be molded according to drilling holes, and generating plane coordinate system coordinate values and elevation values of a pair of three-dimensional coordinate overlapped mark points at the position of the tip vanishing points when the stratum tip vanishing points exist; and respectively adding the plane coordinate system coordinate values and the elevation values of the pair of marking points to the top surface initial demarcation point set and the bottom surface initial demarcation point set of the stratum, and marking the marking points as point vanishing points.
3. The method of claim 2, wherein the processing the set of demarcation points for each formation corresponding to each lithology using the spatial interpolation method to obtain the set of predicted demarcation points for each formation corresponding to each lithology in the region to be modeled comprises: Performing interpolation calculation on the demarcation point set of each stratum corresponding to each lithology by adopting an anisotropic K neighbor sample point spatial interpolation method to obtain a prediction demarcation point set of each stratum corresponding to each lithology in the region to be modeled; The anisotropic K neighbor sample point spatial interpolation method comprises the following steps: dividing the space into a plurality of quadrants by taking a point to be interpolated as a center; Sorting from small to large Euclidean distances between each sample point in each quadrant and the point to be interpolated, and selecting sample points with K bits before sorting to form a sample point subset for interpolation calculation of the point to be interpolated; and carrying out interpolation calculation based on the sample point subset to be subjected to interpolation calculation.
4. A method as claimed in claim 3, wherein the method further comprises: Judging whether all the nearest neighboring 3 drilling stratum demarcation points at the periphery of the sample point to be interpolated are pinch-out points, and if so, marking the lattice point obtained by interpolation calculation of the sample point subset as the pinch-out lattice point, wherein the lattice point where the sample point is positioned in the triangle formed by the 3 drilling stratum demarcation points or on a boundary line.
5. The method of claim 4, wherein constructing a single formation model for each formation corresponding to each lithology using the single formation topology processing method and the set of predicted demarcation points for each formation corresponding to each lithology in the region to be modeled comprises: Triangulating a top surface prediction demarcation point set and a bottom surface prediction demarcation point set in the prediction demarcation point set of each stratum respectively to generate a top surface triangulation and a bottom surface triangulation of each stratum corresponding to each lithology in the region to be modeled; Acquiring a top triangle network boundary and a bottom triangle network boundary of each stratum; Constructing surrounding triangular nets of each stratum, which are connected with the top triangular net and the bottom triangular net, according to the top triangular net boundary and the bottom triangular net boundary of each stratum; and stitching the top triangular mesh, the bottom triangular mesh and the four surrounding triangular meshes of each stratum respectively to form a single stratum model of each stratum which is defined by a vertex set and a triangular patch set and is geometrically closed.
6. The method of claim 5, wherein the method further comprises: Traversing each triangle in the top triangle net and the bottom triangle net of each stratum, respectively judging whether each vertex of the triangle is marked as a pinch-out point or a pinch-out grid point, and if so, deleting the triangle from the triangle net corresponding to the triangle.
7. The method of claim 6, wherein the processing the single stratum model of each stratum corresponding to each lithology by using a multi-stratum topology processing method to obtain the multi-stratum continuous model of the region to be modeled comprises: f1, determining a top-down merging sequence according to a vertical superposition relation of single stratum models of all strata in space; And F2, sequentially carrying out one-to-one topological combination on the adjacent upper single stratum model and the lower single stratum model from top to bottom according to the combination sequence, wherein the one-to-one topological combination on the adjacent upper single stratum model and the lower single stratum model comprises the following steps: Identifying a vertex set used for representing a top interface of the lower single stratum model, establishing a one-to-one index mapping relation between each vertex in the vertex set and a first corresponding vertex with the same coordinates in the fused vertex set; Traversing all triangle patches of the lower single stratum model, for each triangle patch, inquiring the index mapping relation or a new storage position of each vertex in the fused vertex set according to each vertex index, determining a corrected index of each vertex in the fused vertex set, generating new triangle patch data based on the corrected index, and adding the new triangle patch data to the fused triangle set; Based on the fused vertex set and the fused triangle set, generating an intermediate continuous model after the current upper single stratum model and the lower single stratum model are combined, and taking the intermediate continuous model as an upper stratum model for the next combination; And F3, obtaining the multi-stratum continuous model of the region to be modeled after finishing iterative combination of all stratum.
8. The method of claim 7, wherein the method further comprises: calculating the area of each triangle in the multi-stratum continuous model, and marking the triangle with the area of zero or smaller than a first threshold value as a degraded triangle; geometrically merging vertices of the degenerate triangle into one vertex and updating a triangle index referencing the vertices of the degenerate triangle in the multi-stratigraphic continuous model.
9. An automatic generation device of a multi-stratum continuous model for numerical simulation, characterized in that the device comprises: The extraction module is used for acquiring geological identification information of the region to be molded and extracting an initial demarcation point set of each stratum corresponding to each lithology in the region to be molded according to the geological identification information, wherein the geological identification information comprises geological age, coal seam or rock stratum information; The interpolation module is used for processing the demarcation point set of each stratum corresponding to each lithology by adopting a spatial interpolation method to obtain the prediction demarcation point set of each stratum corresponding to each lithology in the region to be modeled; the first construction module is used for constructing a single stratum model of each stratum corresponding to each lithology by adopting a single stratum body topology processing method and a prediction demarcation point set of each stratum corresponding to each lithology in the region to be modeled; and the second construction module is used for processing the single stratum model of each stratum corresponding to each lithology by adopting a multi-stratum body topology processing method to obtain a multi-stratum continuous model of the region to be modeled.
10. A computer readable storage medium, on which a computer program is stored, characterized in that the program, when being executed by a processor, implements the method according to any of claims 1-8.

Description

Automatic generation method and device for multi-stratum continuous model for numerical simulation Technical Field The application relates to the technical fields of computer aided design, mathematical geology and numerical simulation, in particular to an automatic generation method and device of a multi-stratum continuous model for numerical simulation. Background In the field of underground engineering such as coal exploitation, petroleum and natural gas exploration and development, carbon dioxide geological storage and the like, numerical simulation is an important means for predicting geomechanical behaviors. The convergence of the numerical simulation calculation process and the accuracy of the results depend on the grid quality of the geologic model. In order to accurately characterize the heterogeneity of a subsurface space, computing a mesh model requires strict adherence to geometrical fluctuations and contact relationships of sedimentary formations, constructing a multi-formation continuous model that can reflect the true formation structure and geologic structure. However, existing geologic modeling techniques still have limitations in processing complex formations and geologic structures. When the stratum pinch-out is processed, a human-computer interaction is needed to encircle the pinch-out range, and automatic processing is difficult to realize. The grid data of the stratum model is generated by focusing on visual presentation, the quality requirement of the calculated grid model data of numerical simulation is not met, and a large amount of man-machine interaction adjustment is needed for the model. Disclosure of Invention The application provides an automatic generation method and device for a multi-stratum continuous model for numerical simulation, which at least solve the problem that the existing geological modeling technology generates errors or has low efficiency when processing stratum pinch-out, thin stratum and model grid topology, and realize efficient generation of a three-dimensional geological model meeting the technical requirements of numerical simulation. An embodiment of a first aspect of the present application provides a method for automatically generating a multi-stratum continuous model for numerical simulation, the method comprising: Obtaining geological identification information of an area to be modeled, and extracting an initial demarcation point set of each stratum corresponding to each lithology in the area to be modeled according to the geological identification information, wherein the geological identification information comprises geological age, coal seam or rock stratum information; processing the demarcation point set of each stratum corresponding to each lithology by adopting a spatial interpolation method to obtain a prediction demarcation point set of each stratum corresponding to each lithology in the region to be modeled; adopting a single stratum body topology processing method and a prediction demarcation point set of each stratum corresponding to each lithology in the region to be modeled to construct a single stratum model of each stratum corresponding to each lithology; and processing the single stratum model of each stratum corresponding to each lithology by adopting a multi-stratum topological processing method to obtain the multi-stratum continuous model of the region to be modeled. Preferably, the extracting, according to the geological identifier information, the initial demarcation point set of each stratum corresponding to each lithology in the area to be modeled includes: Extracting plane coordinate system coordinate values and elevation values of each top surface boundary point of each stratum corresponding to each lithology in the region to be molded and plane coordinate system coordinate values and elevation values of each bottom surface boundary point based on drilling holes and the geological identification information, and constructing initial boundary point sets of each stratum corresponding to each lithology in the region to be molded based on plane coordinate system coordinate values and elevation values of each top surface boundary point of each stratum corresponding to each lithology in the region to be molded and plane coordinate system coordinate values and elevation values of each bottom surface boundary point, wherein the initial boundary point sets comprise a top surface initial boundary point set and a bottom surface initial boundary point set; Judging whether stratum tip vanishing points exist in each stratum corresponding to each lithology in the area to be molded according to drilling holes, and generating plane coordinate system coordinate values and elevation values of a pair of three-dimensional coordinate overlapped mark points at the position of the tip vanishing points when the stratum tip vanishing points exist; and respectively adding the plane coordinate system coordinate values and the elevation values of the p