Method and device for determining three-dimensional outcrop model

Abstract

The present application provides a method and device for determining a three-dimensional outcrop model, the method comprising: obtaining a cross-sectional image and an irregular triangulated network model corresponding to a target carbonate outcrop, and binarizing the cross-sectional image to obtain a binary grid image; based on the height and width of the cross-sectional image, obtaining the two-dimensional texture coordinates corresponding to each of the fracture vertices of the fracture object in the binary grid image; applying the irregular triangulated network model and the two-dimensional texture coordinates to obtain the three-dimensional global coordinates of the intersection between the irregular triangulated network model and the fracture object, and the three-dimensional global coordinates of each of the fracture vertices; determining the three-dimensional outcrop model corresponding to the target carbonate outcrop according to the three-dimensional global coordinates of the intersection, the fracture vertices and the irregular triangulated network model. The present application can improve the reliability of the three-dimensional outcrop model, and thus can improve the reliability of using the three-dimensional outcrop model to determine underground mineral reservoirs and reflect geological movements.

Description

Technical Field

[0001] The present application relates to the technical field of three-dimensional mapping, and in particular to a method and device for determining a three-dimensional outcrop model. Background Art

[0002] Outcrops are geological bodies exposed on the surface, which can be used as approximate analogs of underground reservoirs. As favorable evidence for exploring underground mineral reservoirs and reflecting geological movements, outcrops play a key role in reservoir modeling. At present, two-dimensional outcrop modeling technology has many disadvantages in terms of modeling efficiency, accuracy and visualization. Therefore, three-dimensional digital outcrop representation has received more and more attention.

[0003] In the existing technology, a three-dimensional outcrop model can be constructed by applying technologies such as drone oblique photography, ground-based lidar and computer virtualization to three-dimensional digital outcrop representation. However, the current three-dimensional outcrop model is mainly based on the display of outcrop profiles, lacking vivid interactive query functions and in-depth analysis methods, resulting in low reliability of the three-dimensional outcrop model. Summary of the invention

[0004] In response to the problems in the prior art, the present application proposes a three-dimensional outcrop model determination method and device, which can improve the reliability of the three-dimensional outcrop model, and further improve the reliability of using the three-dimensional outcrop model to determine underground mineral reservoirs and reflect geological movements.

[0005] In order to solve the above technical problems, this application provides the following technical solutions:

[0006] In a first aspect, the present application provides a method for determining a three-dimensional outcrop model, comprising:

[0007] Obtain the cross-sectional image and irregular triangulated network model corresponding to the target carbonate rock outcrop, and obtain a binary raster image after binarization processing of the cross-sectional image;

[0008] Based on the height and width of the cross-section image, obtain the two-dimensional texture coordinates corresponding to each of the vertices of the hole object in the binary raster image;

[0009] Applying the irregular triangulated network model and the two-dimensional texture coordinates, obtaining the three-dimensional global coordinates of the intersection between the irregular triangulated network model and the seam hole object, and the three-dimensional global coordinates of each of the seam hole vertices;

[0010] According to the three-dimensional global coordinates of the intersection points and the fracture-cavity vertices and the irregular triangulated network model, a three-dimensional outcrop model corresponding to the target carbonate outcrop is determined.

[0011] Furthermore, the obtaining of the two-dimensional texture coordinates corresponding to each of the seam hole vertices of the seam hole object in the binary raster image based on the height and width of the cross-sectional image includes:

[0012] Obtaining raster data of the seam hole object in the binary raster image;

[0013] Using the grid data, the two-dimensional image coordinates of each of the vertices of the crack hole are obtained;

[0014] Based on the two-dimensional image coordinates, the height and width of the cross-sectional image, the two-dimensional texture coordinates corresponding to each of the vertices of the seam hole are obtained.

[0015] Furthermore, the applying of the irregular triangulated network model and the two-dimensional texture coordinates to obtain the three-dimensional global coordinates of the intersection between the irregular triangulated network model and the seam hole object, and the three-dimensional global coordinates of each of the seam hole vertices, includes:

[0016] Acquire triangles corresponding to the vertices of the seam holes from the irregular triangulated network model;

[0017] Applying the triangles and two-dimensional texture coordinates corresponding to the vertices of each of the seam holes to obtain the three-dimensional global coordinates of the vertices of each of the seam holes;

[0018] The three-dimensional global coordinates of the intersection points between the irregular triangulated network model and each of the hole objects are obtained by using the two-dimensional texture coordinates corresponding to each of the hole vertices.

[0019] Furthermore, the obtaining of triangles corresponding to the vertices of the seam holes from the irregular triangulated network model includes:

[0020] Determine a triangle set corresponding to the seam hole object according to the minimum circumscribed rectangle of the seam hole object and an irregular triangulated network model, wherein the irregular triangulated network model includes a plurality of triangles;

[0021] From the triangle set, triangles corresponding to the respective seam hole vertices of the corresponding seam hole object are obtained.

[0022] Furthermore, the applying of the triangles and two-dimensional texture coordinates corresponding to the respective vertices of the seam holes to obtain the three-dimensional global coordinates of the vertices of the seam holes includes:

[0023] Obtaining three-dimensional global coordinates and two-dimensional texture coordinates corresponding to the triangle vertices of the triangle;

[0024] According to the three-dimensional global coordinates corresponding to the triangle vertices of the triangle, a rotation matrix from the three-dimensional global coordinate system corresponding to the irregular triangulated network model to the three-dimensional local coordinate system corresponding to the triangle is obtained;

[0025] Based on the three-dimensional global coordinates corresponding to the triangle vertices and the rotation matrix, obtain the three-dimensional local coordinates corresponding to the triangle vertices;

[0026] Based on the three-dimensional local coordinates and the two-dimensional texture coordinates corresponding to the triangle vertices, obtaining a target affine transformation matrix corresponding to the triangle;

[0027] Applying the target affine transformation matrix corresponding to the triangle and the two-dimensional texture coordinates of the seam hole vertex corresponding to the triangle to obtain the three-dimensional local coordinates of the seam hole vertex;

[0028] Based on the three-dimensional local coordinates of the vertex of the crack hole and the rotation matrix, the three-dimensional global coordinates of the vertex of the crack hole are obtained.

[0029] Furthermore, the applying of the grid data to obtain the two-dimensional image coordinates of each of the vertices of the seam hole includes:

[0030] According to the preset unmarked point determination order, the target unmarked point in the binary raster image is determined. If the crack hole object corresponding to the target unmarked point is a crack linear object, the traversal step is performed:

[0031] Determine whether there is an unmarked point adjacent to the target unmarked point in the binary raster image, and if so, record the direction from the target unmarked point as the starting point to the adjacent unmarked point, and determine whether the direction is consistent with the direction recorded last time;

[0032] If the direction is inconsistent with the direction recorded last time, the grid data corresponding to the target unmarked point is converted into two-dimensional image coordinates and saved, and the target unmarked point is modified into a marked point;

[0033] The traversal step is performed again using the adjacent unmarked point as the target unmarked point until there is no adjacent unmarked point to the target unmarked point;

[0034] According to the preset unmarked point determination order, the next target unmarked point in the binary raster image is determined, and if the crack hole object corresponding to the next target unmarked point is a crack linear object, the traversal step is performed again until all the unmarked points in the binary raster image are traversed;

[0035] The unmarked points are pixels corresponding to the seam hole objects in the binary raster image and the pixels are not marked.

[0036] Furthermore, the three-dimensional outcrop model determination method further includes:

[0037] If the direction is consistent with the direction recorded last time, the traversal step is performed again with the adjacent unmarked point as the target unmarked point until there is no adjacent marked point to the target unmarked point.

[0038] Furthermore, the determining whether there is an unmarked point adjacent to the target unmarked point in the binary raster image includes:

[0039] It is determined whether there is an unmarked point adjacent to the target unmarked point in the direction recorded last time; if not, it is determined whether there is an unmarked point adjacent to the target unmarked point in the direction rotated by a preset angle.

[0040] Furthermore, the three-dimensional outcrop model determination method further includes:

[0041] If the hole object corresponding to the unmarked point of the target is a hole surface object, the contour detection step is performed:

[0042] The cvFindContours function is used to perform contour detection on the hole surface object, and the two-dimensional image coordinates corresponding to the grid data of the crack hole vertices of the hole surface object are obtained;

[0043] According to the preset unmarked point determination order, the next target unmarked point in the binary raster image is determined, and if the crack hole object corresponding to the next target unmarked point is a crack linear object, the traversal step is performed again until all the unmarked points in the binary raster image are traversed;

[0044] The unmarked points are pixels corresponding to the seam hole objects in the binary raster image and the pixels are not marked.

[0045] Furthermore, the three-dimensional outcrop model determination method further includes:

[0046] If the seam hole object corresponding to the unmarked point of the next target is a hole surface object, the contour detection step is performed again.

[0047] In a second aspect, the present application provides a three-dimensional outcrop model determination device, comprising:

[0048] An acquisition module is used to acquire a cross-sectional image and an irregular triangulated network model corresponding to a target carbonate rock outcrop, and obtain a binary raster image after performing a binarization process on the cross-sectional image;

[0049] A texture coordinate acquisition module, used for obtaining two-dimensional texture coordinates corresponding to each seam hole vertex of the seam hole object in the binary raster image based on the height and width of the cross-section image;

[0050] A global coordinate acquisition module, used to apply the irregular triangulated network model and the two-dimensional texture coordinates to obtain the three-dimensional global coordinates of the intersection between the irregular triangulated network model and the seam hole object, and the three-dimensional global coordinates of each of the seam hole vertices;

[0051] A determination module is used to determine the three-dimensional outcrop model corresponding to the target carbonate rock outcrop according to the three-dimensional global coordinates of the intersection points and the fracture-hole vertices and the irregular triangulated network model.

[0052] Furthermore, the texture coordinate acquisition module includes:

[0053] A raster data acquisition unit, used for acquiring raster data of seam hole objects in the binary raster image;

[0054] An image coordinate acquisition unit, used to obtain the two-dimensional image coordinates of each of the vertices of the seam hole by applying the grid data;

[0055] The texture coordinate acquisition unit is used to obtain the two-dimensional texture coordinates corresponding to each of the vertices of the seam hole based on the two-dimensional image coordinates and the height and width of the cross-sectional image.

[0056] Furthermore, the global coordinate acquisition module includes:

[0057] A triangle acquisition unit, used for acquiring triangles corresponding to the vertices of the seam holes from the irregular triangulated network model;

[0058] A first global coordinate acquisition unit, used for obtaining the three-dimensional global coordinates of each of the vertices of the seam hole by applying the triangles and two-dimensional texture coordinates corresponding to each of the vertices of the seam hole;

[0059] The second global coordinate acquisition unit is used to apply the two-dimensional texture coordinates corresponding to each of the seam hole vertices to obtain the three-dimensional global coordinates of the intersection between the irregular triangulated network model and each of the seam hole objects.

[0060] Furthermore, the triangle acquisition unit is used to perform the following:

[0061] Determine a triangle set corresponding to the seam hole object according to the minimum circumscribed rectangle of the seam hole object and an irregular triangulated network model, wherein the irregular triangulated network model includes a plurality of triangles;

[0062] From the triangle set, triangles corresponding to the respective seam hole vertices of the corresponding seam hole object are obtained.

[0063] Furthermore, the first global coordinate acquisition unit is used to perform the following:

[0064] Obtaining three-dimensional global coordinates and two-dimensional texture coordinates corresponding to the triangle vertices of the triangle;

[0065] According to the three-dimensional global coordinates corresponding to the triangle vertices of the triangle, a rotation matrix from the three-dimensional global coordinate system corresponding to the irregular triangulated network model to the three-dimensional local coordinate system corresponding to the triangle is obtained;

[0066] Based on the three-dimensional global coordinates corresponding to the triangle vertices and the rotation matrix, obtain the three-dimensional local coordinates corresponding to the triangle vertices;

[0067] Based on the three-dimensional local coordinates and the two-dimensional texture coordinates corresponding to the triangle vertices, obtaining a target affine transformation matrix corresponding to the triangle;

[0068] Applying the target affine transformation matrix corresponding to the triangle and the two-dimensional texture coordinates of the seam hole vertex corresponding to the triangle to obtain the three-dimensional local coordinates of the seam hole vertex;

[0069] Based on the three-dimensional local coordinates of the vertex of the crack hole and the rotation matrix, the three-dimensional global coordinates of the vertex of the crack hole are obtained.

[0070] Furthermore, the image coordinate acquisition unit is used to perform the following:

[0071] According to the preset unmarked point determination order, the target unmarked point in the binary raster image is determined. If the crack hole object corresponding to the target unmarked point is a crack linear object, the traversal step is performed:

[0072] Determine whether there is an unmarked point adjacent to the target unmarked point in the binary raster image, and if so, record the direction from the target unmarked point as the starting point to the adjacent unmarked point, and determine whether the direction is consistent with the direction recorded last time;

[0073] If the direction is inconsistent with the direction recorded last time, the grid data corresponding to the target unmarked point is converted into two-dimensional image coordinates and saved, and the target unmarked point is modified into a marked point;

[0074] The traversal step is performed again using the adjacent unmarked point as the target unmarked point until there is no adjacent unmarked point to the target unmarked point;

[0075] According to the preset unmarked point determination order, the next target unmarked point in the binary raster image is determined, and if the crack hole object corresponding to the next target unmarked point is a crack linear object, the traversal step is performed again until all the unmarked points in the binary raster image are traversed; wherein the unmarked point is a pixel point corresponding to the crack hole object in the binary raster image and the pixel point is not marked.

[0076] Furthermore, the image coordinate acquisition unit is further used to execute the following:

[0077] If the direction is consistent with the direction recorded last time, the traversal step is performed again with the adjacent unmarked point as the target unmarked point until there is no adjacent marked point to the target unmarked point.

[0078] Furthermore, the image coordinate acquisition unit is further used to execute the following:

[0079] If the hole object corresponding to the unmarked point of the target is a hole surface object, the contour detection step is performed:

[0080] The cvFindContours function is used to perform contour detection on the hole surface object, and the two-dimensional image coordinates corresponding to the grid data of the crack hole vertices of the hole surface object are obtained;

[0081] According to the preset unmarked point determination order, the next target unmarked point in the binary raster image is determined, and if the crack hole object corresponding to the next target unmarked point is a crack linear object, the traversal step is performed again until all the unmarked points in the binary raster image are traversed; wherein the unmarked point is a pixel point corresponding to the crack hole object in the binary raster image and the pixel point is not marked.

[0082] Furthermore, the image coordinate acquisition unit is further configured to execute the contour detection step again if the seam hole object corresponding to the unmarked point of the next target is a hole surface object.

[0083] In a third aspect, the present application provides an electronic device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor implements the three-dimensional outcrop model determination method when executing the program.

[0084] In a fourth aspect, the present application provides a computer-readable storage medium having computer instructions stored thereon, which implement the three-dimensional outcrop model determination method when the instructions are executed.

[0085] As can be seen from the above technical solution, the present application provides a method and device for determining a three-dimensional outcrop model. The method includes: obtaining a cross-sectional image and an irregular triangulated network model corresponding to a target carbonate outcrop, and binarizing the cross-sectional image to obtain a binary grid image; based on the height and width of the cross-sectional image, obtaining the two-dimensional texture coordinates corresponding to each of the fracture-hole vertices of the fracture-hole object in the binary grid image; applying the irregular triangulated network model and the two-dimensional texture coordinates to obtain the three-dimensional global coordinates of the intersection between the irregular triangulated network model and the fracture-hole object, and the three-dimensional global coordinates of each of the fracture-hole vertices; determining the three-dimensional outcrop model corresponding to the target carbonate outcrop according to the three-dimensional global coordinates of the intersection, the fracture-hole vertices and the irregular triangulated network model, which can improve the reliability of the three-dimensional outcrop model, and further improve the reliability of using the three-dimensional outcrop model to determine underground mineral reservoirs and reflect geological movements; specifically, it can also improve the accuracy and efficiency of determining the three-dimensional outcrop model, save memory resources, and improve the accuracy of simulating, depicting and reproducing the structural characteristics of field outcrops. BRIEF DESCRIPTION OF THE DRAWINGS

[0086] In order to more clearly illustrate the embodiments of the present application or the technical solutions in the prior art, the drawings required for use in the embodiments or the description of the prior art will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present application. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative work.

[0087] Figure 1 This is a first flow chart of a method for determining a three-dimensional outcrop model in an embodiment of the present application;

[0088] Figure 2 2 is a schematic diagram of a second process of the method for determining a three-dimensional outcrop model in an embodiment of the present application;

[0089] Figure 3 is a third flow chart of the method for determining a three-dimensional outcrop model in an embodiment of the present application;

[0090] Figure 4 is a fourth flow chart of the method for determining a three-dimensional outcrop model in an embodiment of the present application;

[0091] Figure 5 It is a flowchart of steps 501 to 506 in the method for determining a three-dimensional outcrop model in an embodiment of the present application;

[0092] Figure 6 It is a flowchart of steps 601 to 603 in the method for determining a three-dimensional outcrop model in an embodiment of the present application;

[0093] Figure 7is a flow chart of step 202 of the method for determining a three-dimensional outcrop model in the application example of the present application;

[0094] Figure 8 It is a flow chart of steps 801 to 803 in the method for determining a three-dimensional outcrop model in an embodiment of the present application;

[0095] Fig. 9 It is a logical diagram of the customized eight directions in the application example of this application;

[0096] Fig.10 It is a flowchart of the process of traversing a raster image using the improved eight-direction depth traversal algorithm in the application example of this application;

[0097] Fig.11 It is a comparative diagram of the image coordinate system and the texture coordinate system;

[0098] Fig.12 It is a schematic diagram of the rotation of the three-dimensional coordinate system in a specific application example of the present application;

[0099] Fig.13 It is a flowchart of the process of converting image coordinates and global coordinates of a three-dimensional model in a method for determining a three-dimensional outcrop model in a specific application example of the present application;

[0100] Fig.14 This is a cross-section photo of a carbonate outcrop as an example;

[0101] Fig.15 This is an example of a carbonate rock outcrop profile image with an outcrop hole polygonal layer loaded;

[0102] Fig.16 This is an example of a carbonate rock outcrop profile image with an outcrop fracture linear polygon layer loaded;

[0103] Fig.17 It is a first structural schematic diagram of a three-dimensional outcrop model determination device in an embodiment of the present application;

[0104] Fig.18 is a second structural schematic diagram of the three-dimensional outcrop model determination device in an embodiment of the present application;

[0105] Fig.19 3. It is a third structural schematic diagram of the device for determining a three-dimensional outcrop model in an embodiment of the present application;

[0106] Fig. 20 This is a schematic block diagram of the system structure of the electronic device 9600 according to an embodiment of the present application. DETAILED DESCRIPTION

[0107] In order to enable those skilled in the art to better understand the technical solutions in this specification, the technical solutions in the embodiments of this application will be clearly and completely described below in conjunction with the drawings in the embodiments of this application. Obviously, the described embodiments are only part of the embodiments of this application, not all of them. Based on the embodiments in this application, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of this application.

[0108] In the technical fields of three-dimensional geographic information system, three-dimensional mapping and virtual reality, it is mainly to build a virtual geographic environment based on the comprehensive DEM digital elevation model data and vector basic geographic information data, especially to use vector data to qualitatively and quantitatively describe the characteristic elements such as points, lines and surfaces on the three-dimensional terrain surface, and further perform measurement and spatial analysis on the characteristic elements within the specified area, so as to achieve the effect of efficient discovery and transmission of various thematic map information. However, the above methods are always based on the acquired ground control point coordinates and elevation information to carry out three-dimensional terrain modeling and the acquisition and expression of thematic information of points, lines and surfaces. Based on this, in order to solve the problem of seamless superposition of thematic elements such as fractures and holes in the field of three-dimensional digital outcrop representation and models, this application considers establishing an association between the two-dimensional image coordinate system and the global coordinate system of the three-dimensional model based on the recognition of the two-dimensional image coordinates of thematic elements such as fractures and holes, so as to realize the three-dimensional digital display of the outcrop model and integrate the geological thematic information of the outcrop surface into the outcrop model intelligently, in real time and seamlessly. Among them, the integration of geological thematic information and models is inseparable from the extraction of thematic information. Regarding the extraction of characteristic elements such as cracks and holes on the outcrop surface, the two-dimensional image target automatic recognition technology based on deep learning has been quite mature. Through training, it is entirely possible to quickly identify and extract cracks and holes on the outcrop surface on two-dimensional images.

[0109] Based on this, in order to improve the reliability of the three-dimensional outcrop model, and further improve the accuracy of using the three-dimensional outcrop model to determine the underground mineral reservoir and reflect the geological movement, the embodiment of the present application provides a three-dimensional outcrop model determination device, which can be a server or client device, and the client device can include a smart phone, a tablet electronic device, a network set-top box, a portable computer, a desktop computer, a personal digital assistant (PDA), a vehicle-mounted device, and a smart wearable device, etc. Among them, the smart wearable device can include smart glasses, smart watches, and smart bracelets, etc.

[0110] In practical applications, the part of determining the three-dimensional outcrop model can be performed on the server side as described above, or all operations can be completed in the client device. The selection can be made based on the processing capability of the client device and the limitations of the user's usage scenario. This application does not limit this. If all operations are completed in the client device, the client device may also include a processor.

[0111] The client device may have a communication module (i.e., a communication unit) that can communicate with a remote server to achieve data transmission with the server. The server may include a server on the task scheduling center side, and other implementation scenarios may also include a server on an intermediate platform, such as a server on a third-party server platform that has a communication link with the task scheduling center server. The server may include a single computer device, or a server cluster consisting of multiple servers, or a server structure of a distributed device.

[0112] The server and the client device may communicate using any suitable network protocol, including network protocols that have not yet been developed on the date of filing this application. The network protocol may include, for example, TCP / IP protocol, UDP / IP protocol, HTTP protocol, HTTPS protocol, etc. Of course, the network protocol may also include, for example, RPC protocol (Remote Procedure Call Protocol) and REST protocol (Representational State Transfer) used on top of the above protocols.

[0113] The details are described through the following embodiments.

[0114] In order to improve the reliability of the three-dimensional outcrop model, and further improve the accuracy of using the three-dimensional outcrop model to determine underground mineral reservoirs and reflect geological movements, this embodiment provides a three-dimensional outcrop model determination method in which the execution subject is a three-dimensional outcrop model determination device, such as Figure 1 As shown, the method specifically includes the following contents:

[0115] Step 101: Obtain a cross-sectional image and an irregular triangulated network model corresponding to a target carbonate rock outcrop, and perform binarization processing on the cross-sectional image to obtain a binary raster image.

[0116] The triangulated irregular network (TIN) model can use a series of connected triangles to fit the ground surface or other irregular surfaces to construct a digital ground model, especially a digital elevation model. The Delaunay subdivision method can be used to obtain the triangulated irregular network model corresponding to the target carbonate rock outcrop.

[0117] Furthermore, in order to improve the efficiency of determining the three-dimensional outcrop model, the binary raster image can be divided into n×m sub-binary raster images of the same specifications, and each sub-binary raster image can be used as a binary raster image to execute the three-dimensional outcrop model determination process provided in this embodiment.

[0118] Step 102: Based on the height and width of the cross-sectional image, obtain the two-dimensional texture coordinates corresponding to each of the seam hole vertices of the seam hole object in the binary raster image.

[0119] Specifically, the crack object may be a crack linear object or a hole surface object, the crack linear object may be a crack linear polygon image in a binary raster image, and the hole surface object may be a hole surface polygon image in a binary raster image; the binary raster image may contain multiple crack objects, and the multiple crack objects may include crack linear objects and / or hole surface objects. The height unit and width unit of the cross-sectional image are pixels. The crack vertex is the vertex of the crack object.

[0120] Step 103: Apply the irregular triangulated network model and the two-dimensional texture coordinates to obtain the three-dimensional global coordinates of the intersection between the irregular triangulated network model and the seam hole object, and the three-dimensional global coordinates of each of the seam hole vertices.

[0121] Specifically, the irregular triangulated network model and the two-dimensional texture coordinates corresponding to each of the vertices of the seam hole are applied to obtain the three-dimensional global coordinates of each of the vertices of the seam hole, wherein each of the vertices of the seam hole corresponds to a two-dimensional texture coordinate and a three-dimensional global coordinate. The irregular triangulated network model and the two-dimensional texture coordinates corresponding to each of the vertices of the seam hole are applied to obtain the three-dimensional global coordinates of the intersection between the irregular triangulated network model and each of the seam hole objects.

[0122] Step 104: Determine a three-dimensional outcrop model corresponding to the target carbonate outcrop according to the three-dimensional global coordinates of the intersection points and the fracture-cavity vertices and the irregular triangulated network model.

[0123] Specifically, the three-dimensional global coordinates of the fracture vertices and the three-dimensional global coordinates of the intersection points can be used as feature points, and the cross-sectional image can be mapped onto the irregular triangulated network model to generate a three-dimensional outcrop model corresponding to the target carbonate rock outcrop containing fracture graphic features.

[0124] As can be seen from the above description, the 3D outcrop model determination method provided in this embodiment can determine the 3D global coordinates of each vertex of the fracture and cave object in the target carbonate outcrop, as well as the 3D global coordinates of the intersection with the irregular triangulated network model, by acquiring the cross-sectional image and irregular triangulated network model of the target carbonate outcrop. It can achieve the maximum simulation, characterization and reproduction of the structural characteristics of the field outcrop, improve the reliability and efficiency of determining the 3D outcrop model, and further improve the reliability of using the 3D outcrop model to determine underground mineral reservoirs and reflect geological movements.

[0125] In order to improve the accuracy of obtaining the two-dimensional texture coordinates of the fracture hole vertices, and thus improve the accuracy of subsequently determining the three-dimensional outcrop model using reliable two-dimensional texture coordinates, in one embodiment of the present application, see Figure 2 , step 102 comprises:

[0126] Step 201: Obtain raster data of the seam hole object in the binary raster image.

[0127] Step 202: Apply the grid data to obtain the two-dimensional image coordinates of each of the vertices of the seam hole.

[0128] Specifically, grid data of the seam hole object is obtained from the binary grid image, and the grid data of the seam hole object is applied to obtain the two-dimensional image coordinates corresponding to the seam hole vertices of the seam hole object.

[0129] Step 203: Based on the two-dimensional image coordinates, the height and width of the cross-sectional image, the two-dimensional texture coordinates corresponding to each of the vertices of the seam hole are obtained.

[0130] Specifically, the two-dimensional texture coordinates corresponding to the apex of the crack hole can be obtained based on the height and width of the cross-sectional image and the two-dimensional image coordinates of the apex of the crack hole. Substitute the two-dimensional image coordinates (i, j) of the apex of the crack hole P into the conversion formula:

[0131]

[0132] Among them, H and W represent the height and width of the cross-section image, respectively. In the vertex pixel coordinate system, the origin is the upper left corner of the image, i and j represent the number of rows and columns, respectively. The number of columns increases horizontally to the right, and the number of rows increases vertically downward, and the two-dimensional texture coordinates (u, v) of the seam vertex P are obtained.

[0133] In order to improve the accuracy of obtaining the three-dimensional global coordinates of the fracture vertices and the three-dimensional global coordinates of the intersections between the irregular triangulated network model and each of the fracture objects, and thus improve the accuracy of the subsequent application of reliable three-dimensional global coordinates to determine the three-dimensional outcrop model, in one embodiment of the present application, see Figure 3 , step 103 comprises:

[0134] Step 301: Obtain triangles corresponding to the vertices of the seam holes from the irregular triangulated network model.

[0135] Step 302: Apply the triangles and two-dimensional texture coordinates corresponding to each of the seam hole vertices to obtain the three-dimensional global coordinates of each of the seam hole vertices.

[0136] Specifically, the triangle and the two-dimensional texture coordinates corresponding to any of the vertices of the seam hole may be used to obtain the three-dimensional global coordinates of the vertex of the seam hole.

[0137] Step 303: applying the two-dimensional texture coordinates corresponding to each of the seam hole vertices to obtain the three-dimensional global coordinates of the intersection points between the irregular triangulated network model and each of the seam hole objects.

[0138] Specifically, the 2D texture coordinates of any of the vertices of the seam hole can be used to obtain the 3D global coordinates of the intersection between the irregular triangulated network model and the seam hole object. The texture coordinates of the intersection between the seam hole object and the irregular triangulated network model are calculated based on the 2D texture coordinates corresponding to each seam hole vertex of the seam hole object, and the 3D global coordinates of the intersection are calculated by linear interpolation.

[0139] In order to improve the efficiency of determining the triangles corresponding to the fracture-hole vertices, and thus improve the efficiency of determining the three-dimensional outcrop model, in one embodiment of the present application, see Figure 4 , step 301 comprises:

[0140] Step 401: Determine a triangle set corresponding to the seam hole object according to a minimum circumscribed rectangle of the seam hole object and an irregular triangulated network model, wherein the irregular triangulated network model includes a plurality of triangles.

[0141] Specifically, the minimum bounding rectangle of the seam hole object is determined in the binary raster image, each of the seam hole objects has a corresponding minimum bounding rectangle, the irregular triangulated network model contains multiple triangles, and the minimum bounding rectangle is applied to the irregular triangulated network model to filter triangles to obtain a triangle set S. The screening rules are as follows: based on the two-dimensional texture coordinates of each vertex of the triangle and the minimum bounding rectangle, ① if the minimum value of the horizontal coordinates of the three vertices of the triangle is greater than the maximum value of the horizontal coordinates of each point on the minimum bounding rectangle; ② the maximum value of the horizontal coordinates of the three vertices of the triangle is less than the minimum value of the horizontal coordinates of each point on the minimum bounding rectangle; ③ the minimum value of the vertical coordinates of the three vertices of the triangle is greater than the maximum value of the vertical coordinates of each point on the minimum bounding rectangle; ④ the maximum value of the vertical coordinates of the three vertices of the triangle is less than the minimum value of the vertical coordinates of each point on the minimum bounding rectangle, triangles that do not meet the above four conditions will be retained to form a triangle set S.

[0142] Step 402: From the triangle set, obtain the triangles corresponding to the respective seam hole vertices of the corresponding seam hole object.

[0143] Specifically, from the triangle set, triangles corresponding to each seam vertex of the seam object corresponding to the triangle set are obtained; taking each seam vertex of the seam object P as a unit, in the triangle set S, the triangle T to which a vertex a' belongs is judged, and the judging steps are as follows: the triangles containing the seam vertex are screened out in the triangle set S; and the triangle T containing the vertex a' is determined among the triangles, wherein the inclusion of the vertex a' in the triangle T indicates that the vertex a' is inside or on the triangle T.

[0144] In order to obtain reliable 3D local coordinates and 3D global coordinates of the fracture-cavity vertices, and then use the reliable 3D local coordinates and 3D global coordinates to improve the reliability of the subsequent determination of the 3D outcrop model, in one embodiment of the present application, see Figure 5 , step 302 comprises:

[0145] Step 501: Obtain the three-dimensional global coordinates and two-dimensional texture coordinates corresponding to the triangle vertices of the triangle.

[0146] Step 502: According to the three-dimensional global coordinates corresponding to the triangle vertices of the triangle, a rotation matrix from the three-dimensional global coordinate system corresponding to the irregular triangulated network model to the three-dimensional local coordinate system corresponding to the triangle is obtained.

[0147] Specifically, the y-axis of the three-dimensional global coordinate system may be rotated to a direction consistent with the normal vector of the triangle to obtain the rotation matrix.

[0148] The three-dimensional global coordinate system corresponding to the irregular triangulated network model may be a world coordinate system.

[0149] Step 503: Based on the three-dimensional global coordinates corresponding to the triangle vertices and the rotation matrix, obtain the three-dimensional local coordinates corresponding to the triangle vertices.

[0150] Specifically, the rotation matrix can be obtained according to the three-dimensional global coordinates of each triangle vertex of the triangle; the three-dimensional local coordinates of each triangle vertex of the triangle can be obtained according to the rotation matrix and the three-dimensional global coordinates corresponding to the triangle vertex; each triangle vertex has corresponding three-dimensional global coordinates, three-dimensional local coordinates and two-dimensional texture coordinates, and the three-dimensional global coordinates and two-dimensional texture coordinates corresponding to the triangle vertex can be determined according to the irregular triangle model, and the triangle vertex is the vertex of the triangle. Assuming that the three triangle vertices of the triangle are P1, P2 and P3 respectively, the three unit vectors corresponding to the triangle are obtained according to the three-dimensional global coordinates of P1, P2 and P3 and the following formulas

[0151]

[0152] Assumptions according to and The rotation matrix R can be determined:

[0153]

[0154] The three-dimensional local coordinates of the triangle vertices can be obtained according to the following formula:

[0155]

[0156] Where (x', y', z') are the three-dimensional local coordinates of the triangle, and y' is a constant on the entire triangle surface in the three-dimensional local coordinate system.

[0157] Step 504: Based on the three-dimensional local coordinates and two-dimensional texture coordinates corresponding to the triangle vertices, obtain the target affine transformation matrix corresponding to the triangle.

[0158] Specifically, assuming the target affine transformation matrix Substituting the coordinate groups corresponding to the three triangle vertices of the triangle into the following formulas respectively, the target affine transformation matrix M corresponding to the triangle can be obtained. Each coordinate group includes: three-dimensional local coordinates and two-dimensional texture coordinates corresponding to a unique triangle vertex.

[0159]

[0160] Among them, (u', v') is the two-dimensional texture coordinate corresponding to the triangle vertex, (x', y', z') is the three-dimensional local coordinate of the triangle vertex, and y' on the entire triangle surface is a constant in the three-dimensional local coordinate system, so we only need to consider the conversion relationship between (u', v') and (y', z').

[0161] Step 505: Apply the target affine transformation matrix corresponding to the triangle and the two-dimensional texture coordinates of the seam hole vertex corresponding to the triangle to obtain the three-dimensional local coordinates of the seam hole vertex.

[0162] Specifically, the three-dimensional local coordinates of the vertices of the crack holes can be obtained according to the following formula:

[0163]

[0164] Among them, M is the target affine transformation matrix, (u, v) is the two-dimensional texture coordinate of the hole vertex, (x1, y1, z1) is the three-dimensional local coordinate of the hole vertex, and since y' on the entire triangle surface in the local coordinate system is a constant, y1=y' is directly obtained.

[0165] Step 506: Based on the three-dimensional local coordinates of the vertex of the crack hole and the rotation matrix, obtain the three-dimensional global coordinates of the vertex of the crack hole.

[0166] Specifically, the inverse matrix of the rotation matrix can be obtained first, and the three-dimensional global coordinates of the vertex of the crack hole can be obtained according to the three-dimensional local coordinates of the vertex of the crack hole and the inverse matrix of the rotation matrix; the three-dimensional global coordinates of the vertex of the crack hole can be obtained according to the following formula:

[0167]

[0168] Among them, R -1 is the inverse matrix of the rotation matrix, (x1, y1, z1) is the three-dimensional local coordinate of the vertex of the crack hole, and (x, y, z) is the three-dimensional global coordinate of the vertex of the crack hole.

[0169] In order to improve the reliability and efficiency of obtaining the two-dimensional image coordinates corresponding to each of the vertices of the seam hole and save memory usage, in one embodiment of the present application, see Figure 6 , step 202 comprises:

[0170] Step 601: according to the preset unmarked point determination order, determine the target unmarked point from multiple unmarked points in the binary raster image. If the crack object corresponding to the target unmarked point is a crack linear object, perform a traversal step: determine whether there is an unmarked point adjacent to the target unmarked point in the binary raster image. If so, record the direction from the target unmarked point as the starting point to the adjacent unmarked point, and determine whether the direction is consistent with the direction recorded last time; if the direction is inconsistent with the direction recorded last time, convert the raster data corresponding to the target unmarked point into two-dimensional image coordinates and save them, and modify the target unmarked point to a marked point.

[0171] Specifically, the order of determining the unmarked points can be set according to actual needs. For example, each pixel point in the binary raster image can be searched in the order from left to right and then from top to bottom. If an unmarked point appears for the first time, it is determined as a target unmarked point; the unmarked point is a pixel point corresponding to a seam hole object in the binary raster image and the pixel point is not marked, and each unmarked point corresponds to a pixel point. The direction can be one of the upper left, lower left, left, upper right, lower right, right, upper and lower. The direction recorded last time can be the direction recorded when the traversal step was performed last time and the corresponding seam hole object is the same as the seam hole object corresponding to the traversal step performed this time; it can be understood that if there is no direction recorded last time, that is, the current target unmarked point is the starting point of the corresponding seam hole object, the raster data corresponding to the target unmarked point is converted into a two-dimensional image coordinate and saved, and the adjacent unmarked point is used as the target unmarked point to perform the traversal step again.

[0172] Among them, the unmarked points adjacent to the target unmarked point refer to the unmarked points adjacent to the target unmarked point except the previous target unmarked point corresponding to the target unmarked point. After the raster data corresponding to the target unmarked point is converted into two-dimensional image coordinates, the target unmarked point can be modified to a marked point, which can avoid repeated traversal to improve efficiency.

[0173] In order to further improve the efficiency of obtaining the two-dimensional image coordinates corresponding to each of the vertices of the seam hole, the determination of whether there are unmarked points adjacent to the target unmarked point in the binary raster image may include: determining whether there are unmarked points adjacent to the target unmarked point in the direction of the previous record, and if not, determining whether there are unmarked points adjacent to the target unmarked point in the direction after rotating the preset angle. The preset angle can be set according to actual needs, and this application does not impose any restrictions on this. Each time an adjacent point is searched, the direction will be recorded. When searching for adjacent points next time, the direction of the next point will be searched first. If it is not found, it will be searched clockwise in sequence, which can greatly reduce the number of searches and improve efficiency.

[0174] Step 602: Execute the traversal step again with the adjacent unmarked point as the target unmarked point until there is no adjacent unmarked point to the target unmarked point.

[0175] It is understandable that if the target unmarked point has no adjacent unmarked point, it means that the two-dimensional image coordinates of each vertex of the current seam hole object have been determined, and the operation of the current seam hole object can be ended and the operation of the next seam hole object can be performed.

[0176] Step 603: Determine the order according to the preset unmarked points, determine the next target unmarked point in the binary raster image, if the crack object corresponding to the next target unmarked point is a crack linear object, execute the traversal step again until all the unmarked points in the binary raster image are traversed, wherein the unmarked point is a pixel point corresponding to the crack object in the binary raster image and the pixel point is not marked.

[0177] In order to further reduce the occupation of memory resources, if the direction is consistent with the direction recorded last time, the adjacent unmarked point is used as the target unmarked point to perform the traversal step again until there is no adjacent marked point to the target unmarked point.

[0178] From the above description, it can be seen that the three-dimensional outcrop model determination method provided in this embodiment automatically simplifies line elements during the traversal process, retaining only the first and last points of the same line segment, retaining the characteristic points (inflection points) of the curve to the greatest extent, and achieving the effect of compressing the fracture vector layer data.

[0179] To further illustrate this solution, in an application example of this application, see Figure 7 , step 202 comprises:

[0180] Step 701: Determine the target unmarked points in the binary raster image according to a preset unmarked point determination order.

[0181] Step 702: If the crack hole object corresponding to the target unmarked point is a crack linear object, a traversal step, namely step 703, is executed.

[0182] Step 7021: If the hole object corresponding to the target unmarked point is a hole surface object, perform a contour detection step: use the cvFindContours function to perform contour detection on the hole surface object to obtain the two-dimensional image coordinates corresponding to the grid data of the hole vertex of the hole surface object. After step 7021 is completed, execute step 708.

[0183] Step 703: Determine whether there is an unmarked point adjacent to the target unmarked point in the binary raster image. If yes, execute step 704; otherwise, execute step 708.

[0184] Step 704: Record the direction from the target unmarked point as the starting point to the adjacent unmarked point.

[0185] Step 705: Determine whether the direction is consistent with the direction recorded last time and whether the target unmarked point is not the starting point of the seam hole object. If not, execute step 706; if yes, directly execute step 706.

[0186] Specifically, the starting point of the seam hole object may be a target unmarked point when the traversal step corresponding to the seam hole object is executed for the first time.

[0187] Step 706: convert the grid data corresponding to the target unmarked point into two-dimensional image coordinates and save them, and modify the target unmarked point to a marked point, and modify the target unmarked point to a marked point.

[0188] Step 707: Take the adjacent unmarked point as the target unmarked point and execute step 703 again.

[0189] Step 708: Determine the order according to the preset unmarked points and judge whether there is a next target unmarked point in the binary raster image; if so, execute step 709; if not, end the current operation.

[0190] Step 709: Determine whether the crack hole object to which the next target unmarked point belongs is a crack linear object. If yes, take the next target unmarked point as the target unmarked point and execute step 703 again; if no, take the next target unmarked point as the target unmarked point and execute step 7021 again.

[0191] In order to further improve the reliability and efficiency of obtaining the two-dimensional image coordinates corresponding to each of the vertices of the seam hole, in one embodiment of the present application, see Figure 8 , the three-dimensional outcrop model determination method further includes:

[0192] Step 801: If the hole object corresponding to the target unmarked point is a hole surface object, execute the contour detection step: apply the cvFindContours function to perform contour detection on the hole surface object to obtain the two-dimensional image coordinates corresponding to the grid data of the hole vertices of the hole surface object.

[0193] The cvFindContours function is provided by the OPENCV image processing library.

[0194] Step 802: Determine the order according to the preset unmarked points, determine the next target unmarked point in the binary raster image, and if the crack object corresponding to the next target unmarked point is a crack linear object, execute the traversal step again until all the unmarked points in the binary raster image are traversed; wherein the unmarked point is a pixel point corresponding to the crack object in the binary raster image and the pixel point is not marked.

[0195] Further, see Figure 8 After step 802, the method may further include step 803: if the seam hole object corresponding to the unmarked point of the next target is a hole surface object, the contour detection step is performed again.

[0196] In order to further illustrate the scheme, the present application also provides an application example of a three-dimensional outcrop model determination method, which specifically includes the following contents:

[0197] S100: Input binary image raster data of the extraction result of two-dimensional lines and surface objects of crack-hole features on carbonate rock outcrop profile photos, and cut the crack-hole line and surface raster images into n×m raster images of the same specification by customizing the number of rows (n) and columns (m). In an example, the carbonate rock outcrop profile photos can be as follows: Fig.14 shown.

[0198] S200: using the improved eight-direction depth traversal algorithm, for the n×m grid images obtained in step S100, converting the grid data of the crack linear object into vector data, and saving the two-dimensional image coordinates of each vertex of the crack linear object in a .txt file in units of cracks. S200 specifically includes:

[0199] S201: Find the first black point in the raster image.

[0200] S202: Start from this point and extend the search curve in eight directions. Each time a point is found, it is set as point A, and the search direction is recorded and saved.

[0201] S203: Determine whether this direction is consistent with the direction saved last time. If it is consistent, and point A is not the starting point of the crack linear object, then point A is not saved in the deque. <point>Double queue, otherwise save the 2D image coordinates (i, j) of point A into deque <point>Double queue, then whiten to prevent duplicate searches.

[0202] S204: When a point without adjacent points is found, the point without adjacent points is saved as a vertex in the deque. <point>In the double queue, it indicates that the search is completed and the crack line object is saved.

[0203] S205: Search for another crack linear object until all traversals are completed.

[0204] Specifically, for the conversion of outcrop surface fracture raster data into vector data, the improved eight-direction depth traversal algorithm uses deque <point>The double queue is used to describe the curve, and the input processing image is required to be a thinned binary image, where the curve width is 1 pixel and the color is black, such as Fig. 9 As shown, the custom eight directions can be direction 1 to direction 8, and the user can adjust according to the actual situation. The improved eight-direction depth traversal algorithm is specifically described as follows:

[0205] 1. Find the first black point in the raster image; 2. Start from this point and extend the search curve in eight directions. Every time a point is found, record the search direction and save it; 3. Determine whether this direction is consistent with the direction saved last time. If it is consistent, and this point is not the starting point of the line object, do not save this point into the deque <point>Double queue, otherwise save the 2D image coordinates of this point into deque <point>Double queue, and then turn white (to prevent repeated search); 4. Search for a point without adjacent points, indicating that the search is completed, and save this crack linear object; 5. Search for another crack linear object until all traversals are completed.

[0206] See also Fig.10 , the process of traversing the raster image using the eight-direction depth traversal algorithm specifically includes:

[0207] S1: Initialize vector <deque <cv::point>>lines: Initializes a vector (array) for storing all line (crack) objects. This vector stores all line objects; S2: Initializes deque <cv::point>line: initialize a bidirectional queue for storing a line (crack), which stores each vertex on the line; set the search direction d1 = 1; S3: determine whether the first black point P1 is found, if so, execute step S4, otherwise end the current operation; S4: add P1 to line; S5: search in eight directions starting from the direction d1 of point P1; S6: determine whether the black point p is found, if so, execute step S7, otherwise execute step S9; S7: record direction d, d = {1, 2, 3, 4, 5, 6, 7, 8}; set P1 = p; turn point p white; S8: determine d1 == d, if so, execute step S5, otherwise, execute step S4; S9: add line to lines.

[0208] The improved eight-direction deep traversal algorithm has the following two outstanding advantages: First, each time an adjacent point is searched, the direction will be recorded. When searching for adjacent points next time, the direction of the next point will be searched first. If it is not found, it will be searched clockwise in sequence. This greatly reduces the number of searches and improves efficiency. Second, the algorithm automatically simplifies line features during the traversal process, retaining only the first and last points of the same line segment, retaining the characteristic points (inflection points) of the curve to the greatest extent, and achieving the effect of compressing the crack vector layer data.

[0209] S300: Using the function cvFindContours provided by the OPENCV image processing library, perform contour detection on the hole surface object for the n×m grid images obtained in step S100, obtain the two-dimensional image coordinates (i, j) of the vertices of the hole surface object, and save them in a .txt file in units of holes.

[0210] S400: Read the .txt file, taking the crack hole polygon as the unit. The crack hole polygon can be the above-mentioned crack linear object or hole surface object. The crack linear object is the crack linear polygon, and the hole surface object is the hole surface polygon. Fig.15 This is an example of a carbonate rock outcrop profile image with an outcrop hole polygon layer loaded. a1 to a3 are all hole polygons. Fig.16 This is an example of a carbonate rock outcrop profile image with an outcrop fracture linear polygon layer loaded, where b1 to b3 are all fracture linear polygons. For a certain fracture-hole polygon P, the vertex two-dimensional image coordinates (i, j) of the fracture-hole polygon P obtained in step S200 or S300 are substituted into the conversion formula: Where H and W represent the height and width of the image, respectively, and the two-dimensional texture coordinates (u, v) are obtained.

[0211] Specifically, in this application example, the coordinates of the vertices in each coordinate system are described uniformly: the point coordinates in the image coordinate system are (i, j); the two-dimensional texture coordinates of the seam hole vertices in the texture coordinate system are (u, v), and the two-dimensional texture coordinates of the triangle vertices are (u', v'); the three-dimensional local coordinates of the seam hole vertices in the three-dimensional local coordinate system are (x1, y1, z1), and the three-dimensional local coordinates of the triangle vertices are (x', y', z'); the three-dimensional global coordinates of the seam hole vertices in the three-dimensional global coordinate system are (x, y, z), and the three-dimensional global coordinates of the triangle vertices are (x2, y2, z2).

[0212] In the process of converting the image coordinates (i, j) of the vertex of the seam hole vector polygon to the texture coordinates (u, v), see Fig.11 In the image coordinate system, i.e., the vertex pixel coordinate system, the origin is the upper left corner of the image, i and j represent the row and column numbers respectively, the number of columns increases horizontally to the right, and the number of rows increases vertically downwards, while in the texture coordinate system, the origin of the coordinate system is the lower left corner of the image, u and v represent the horizontal and vertical coordinates respectively, u increases horizontally to the right, and v increases vertically upwards, and the value range is (0,0) to (1,1). Therefore, the conversion formula is: Where H and W represent the height and width of the image respectively.

[0213] S500: Find the enclosing rectangle of the fracture-hole polygon P, and use the enclosing rectangle to screen triangles on the carbonate rock outcrop TIN, that is, on the above irregular triangulated network model, to obtain a triangle set S.

[0214] S600: Taking each vertex of the seam hole polygon P as a unit, determine the triangle T to which a vertex a' belongs in the triangle set S, and calculate its rotation matrix R and its inverse matrix R based on the three-dimensional global coordinates (x2, y2, z2) of the vertex of triangle T. -1 , use the rotation matrix R to rotate the three vertices of triangle T to obtain the three-dimensional local coordinates of the vertices (x', y', z'), where y' is the same, record and save, and combine with the vertex texture coordinates (u', v') of triangle T, and bring it into the formula: Obtain the affine transformation matrix M.

[0215] Specifically, the process of establishing an association (affine transformation) between the texture coordinate system and the three-dimensional local coordinate system includes: through the above conversion, the three-dimensional local coordinates and texture coordinates of the three vertices of the known outcrop model TIN triangle surface are known. This application example assumes that the triangle surface in the model and the triangle surface in the corresponding texture space are in an affine transformation relationship. Therefore, the affine transformation relationship can be established based on the three-dimensional local coordinates (x', y', z') of the triangle vertex and the texture coordinates (u', v'), and the affine transformation matrix M is solved. Affine transformation matrix: In order to find the affine transformation matrix M, a matrix operation formula can be established: At this time, after the above 3D global coordinates are transformed into 3D local coordinates, the y' of each point on the same triangle is the same, which is the constant yc, so it does not need to be considered. Next, the 3D local coordinates and texture coordinates of the three vertices of the triangle are substituted into the formula, and the parameters a, b, c, d, m, and n are solved to obtain the affine transformation matrix M.

[0216] S700: Using the affine transformation matrix M obtained in step S600, the two-dimensional texture coordinates (u, v) of the vertex a' of the seam hole polygon P are substituted into the formula: Perform affine transformation to obtain its three-dimensional local coordinates (x1, y1, z1); then combine the inverse matrix R of the rotation matrix R obtained in step S600 -1 , put it into the formula:

[0217] The three-dimensional global coordinates (x, y, z) of vertex a' are solved by rotation transformation and saved as a .txt file; return to step S600 and continue to traverse the next vertex of the seam hole polygon P until all vertices are traversed.

[0218] Specifically, in this application example, the process of establishing an association between the three-dimensional global coordinate system and the three-dimensional local coordinate system (rotation of the three-dimensional coordinate system) includes: the three-dimensional global coordinates (x2, y2, z2) and texture coordinates (u', v') of the three vertices of each triangle in the known outcrop model TIN, assuming that the triangular faces in the model and the triangular faces in the corresponding texture space are in an affine transformation relationship, so in order to keep the three-dimensional local coordinate system matching the texture coordinate system, it is first necessary to rotate the three-dimensional coordinate system and calculate the rotation matrix R.

[0219] This application example rotates the coordinates to the direction where the Y axis is consistent with the normal vector of the triangle surface, customizes the directions of the X axis and Z axis, and obtains the coordinate values ​​of the three vertices in the three-dimensional local coordinates according to the rotation transformation. Fig.12 As shown, the UZ direction pointed by the side of triangle P1P3 is selected as the positive direction of the local coordinate system Z axis, and the UY direction pointed by the normal vector N of the triangle plane is selected as the positive direction of the local coordinate system Y axis. Therefore, the UX direction perpendicular to the plane formed by the Z axis and the Y axis is selected as the positive direction of the local coordinate system X axis.

[0220] Then the rotation matrix R can be directly expressed as: See also Fig.12 , given the three-dimensional global coordinates of the vertices P1, P2, and P3 of the triangle, find the three unit vectors of the new coordinate axis XYZ in the original coordinate system but It is expressed as:

[0221]

[0222] At this time, according to the matrix operation formula: A coordinate transformation relationship is established between the three-dimensional global coordinate system and the three-dimensional local coordinate system. At this time, the triangular face in the local coordinate system is parallel to the plane formed by the coordinate axes X and Z, so the y' value of each point on the triangular face is the same and is a constant.

[0223] S800: For the fracture polygon P, according to the two-dimensional texture coordinates (u, v) of the vertices of the fracture polygon P calculated in step S400, the texture coordinates (u1, v1) of the intersection of the polygon P and the carbonate outcrop TIN are calculated, and the three-dimensional global coordinates (x3, y3, z3) of the intersection are calculated by linear interpolation of the texture coordinates on the adjacent triangulated mesh vertices, and saved in a .txt file.

[0224] Specifically, this application example establishes associations between the point image coordinate system, texture coordinate system, 3D local coordinate system and 3D global coordinate system on the outcrop model TIN. For the point that needs to be mapped, when the image coordinates of the point are input, its texture coordinates can be easily calculated. Then, by determining the outcrop model TIN triangle to which the point belongs, the rotation matrix R and the affine transformation matrix M are calculated according to the 3D global and texture coordinates of the triangle vertices. Finally, the 3D global coordinates of the point can be obtained by simply performing the inverse transformation of the affine transformation and the rotation transformation on the texture coordinates of the point.

[0225] S900: Return to step S400, traverse the next seam hole polygon, and continue to the next step until all seam hole polygons are traversed.

[0226] S110: Read the .txt file and generate a .shp file according to the three-dimensional global coordinates of the vertices and intersections of the seam hole polygons.

[0227] Specifically, the three-dimensional outcrop model corresponding to the carbonate outcrop can be determined according to the .shp file.

[0228] From the above description, it can be seen that the three-dimensional outcrop model determination method provided by this application example, based on the cutting of two-dimensional fracture hole line and surface grid images, uses the improved eight-direction depth traversal algorithm and cvFindContours function to convert the fracture hole polygon grid data into vector data, which can greatly improve the efficiency of polygon grid data feature point retrieval and conversion; uses affine transformation and three-dimensional coordinate system rotation methods to establish a mapping relationship between the two-dimensional image and the three-dimensional model, converts the image coordinates of the two-dimensional fracture hole vector graphic feature points to texture coordinates, and then to three-dimensional local coordinates, and finally successfully converts to three-dimensional global coordinates and generates a fracture hole three-dimensional .shp layer file, so as to achieve seamless superposition of vector fracture hole features on the three-dimensional digital outcrop; in order to achieve accurate description of fracture hole vector objects, not only the vertices of the fracture hole polygon are used as feature points, but also the three-dimensional model global coordinates of the intersection of the fracture hole line, surface polygon and carbonate rock outcrop TIN are calculated by linear interpolation method and mapped to the carbonate rock outcrop model surface as feature points, which can achieve the maximum simulation, characterization and reproduction of the field outcrop structural characteristics.

[0229] In order to further illustrate the scheme, the present application also provides a specific application example of a three-dimensional outcrop model determination method, including:

[0230] S10: First, the seam hole line and surface raster binary images are cut into n×m raster images of the same specifications by customizing the number of rows (n) and columns (m). Considering that the original seam hole line and surface raster images have many pixels, for example, the image size in this specific application example is 7264 pixels wide and 5440 pixels high, if the image elements are to be traversed and retrieved one by one, the efficiency of the algorithm is very low. Therefore, the original image is cut into 32 rows and 32 columns, that is, cut into n×m images of the same specifications of 227 pixels wide and 170 pixels high for storage.

[0231] S20: Convert the crack linear object raster data to vector data for the 32×32 binary images obtained in step S10, and obtain the two-dimensional image coordinates (i, j) of the vertices of the crack linear object. Use the improved eight-direction depth traversal method to traverse the crack linear object, and finally obtain the two-dimensional image coordinates of the vertices of the crack linear object, then splice the image coordinates, and store them in units of crack hole polygons and save them in a .txt file.

[0232] S30: Convert the raster data of the hole surface object to vector data of the 32×32 binary images obtained in step S10, and obtain the two-dimensional image coordinates (i, j) of the vertex of the hole surface object. Use the cvFindContours function provided by the OPENCV image processing library to detect the contour of the object, obtain the two-dimensional image coordinates of the vertex of the hole surface element, and save it in a .txt file.

[0233] S40: Read the .txt file to obtain the two-dimensional image coordinates of n hole polygons. Taking the hole polygon as a unit, for a certain polygon P, the two-dimensional image coordinates (i, j) of the vertex of the hole polygon obtained in steps S20 and S30 are substituted into the conversion formula: In this specific application example, the image width W=7264 pixels, the height H=5440 pixels, and the two-dimensional texture coordinates (u, v) of the polygon vertices are obtained by calculation.

[0234] S50: Find the enclosing rectangle of polygon P, and use the enclosing rectangle to filter triangles on the carbonate outcrop TIN to obtain a triangle set S. The screening rules are as follows: ① If the minimum value of the horizontal coordinate u' of the three vertices of the triangle is greater than the maximum value of the horizontal coordinates of each point on the enclosing rectangle; ② The maximum value of the horizontal coordinate u' of the three vertices of the triangle is less than the minimum value of the horizontal coordinates of each point on the enclosing rectangle; ③ The minimum value of the vertical coordinate v' of the three vertices of the triangle is greater than the maximum value of the vertical coordinates of each point on the enclosing rectangle; ④ The maximum value of the vertical coordinate v' of the three vertices of the triangle is less than the minimum value of the vertical coordinates of each point on the enclosing rectangle, the triangles that do not meet the above 4 conditions will be retained to form a set S.

[0235] S60: Taking each vertex of polygon P as a unit, determine in set S which triangle T a vertex a' belongs to. The determination steps are as follows: ① Filter out the triangles whose bounding rectangles contain the vertex in set S; ② Determine whether the vertex is inside or on the triangle. For triangle T, based on the three-dimensional global coordinates (x2, y2, z2) of the vertices of triangle T, calculate its rotation matrix R and the inverse matrix of R, and use the rotation matrix R to rotate and transform the three vertices of triangle T to obtain the three-dimensional local coordinates of the vertex (x', y', z'), where y' is the same, record and save, and combine with the texture coordinates (u', v') of the vertex of triangle T and substitute into the formula: Solve the affine transformation matrix M.

[0236] S70: Using the affine transformation matrix M solved in step S60, the texture coordinates (u, v) of the vertex a' of the polygon P are substituted into the formula: Perform affine transformation to obtain its three-dimensional local coordinates (x1, y1, z1); then combine the inverse matrix R of the rotation matrix R obtained in step S60 -1 , put it into the formula: The three-dimensional global coordinates (x, y, z) of vertex a' are solved by rotation transformation and saved as a .txt file; return to step S60 and continue to traverse the next vertex of polygon P until all vertices are traversed.

[0237] S80: For polygon P, according to the two-dimensional texture coordinates (u, v) of the vertices of polygon P calculated in step S40, calculate the texture coordinates (u1, v1) of the intersection of polygon P and carbonate outcrop TIN, and calculate the three-dimensional global coordinates (x3, y3, z3) of the intersection on the adjacent triangulated mesh vertices according to the texture coordinates, and save them to a .txt file.

[0238] S90: Return to step S40, traverse the next polygon, and continue to the next step until all polygons are traversed.

[0239] S91: Read the .txt file and generate a .shp file according to the three-dimensional global coordinates of the vertices and intersections of the seam hole polygons.

[0240] S92: Determine, according to the .shp file, a three-dimensional model corresponding to the graphic features of the fractures and caves in the carbonate outcrop.

[0241] From the above description, it can be seen that the three-dimensional outcrop model determination method provided in this application example can establish an association between the two-dimensional image coordinate system and the three-dimensional model global coordinate system, integrate various geological thematic information into the model, and realize the image and three-dimensional model coordinate mapping based on affine transformation. It can improve the reliability of the three-dimensional outcrop model, and then improve the accuracy of using the three-dimensional outcrop model to determine underground mineral reservoirs and reflect geological movements.

[0242] In order to further illustrate the present scheme, the present application also provides a specific application example of the process of converting image coordinates and global coordinates of a three-dimensional model in a three-dimensional outcrop model determination method, see Fig.13 , the specific application examples include:

[0243] S01: Read the hole polygon vertex .txt file to obtain n polygons.

[0244] S02: Convert the image coordinates of the mth polygon vertex into texture coordinates.

[0245] S03: Find the enclosing rectangle of the mth polygon, and use the enclosing rectangle to filter triangles on the TIN to obtain a triangle set S.

[0246] S04: Determine the triangle L to which the k-th vertex of the m-th polygon belongs: 1. Filter out the triangles whose outer rectangles contain the vertex in the set S; 2. Determine whether the vertex is inside or on the triangle.

[0247] S05: Calculate the rotation matrix R and affine transformation matrix M of triangle L.

[0248] S06: Convert polygon vertex texture coordinates to 3D global coordinates.

[0249] S07: Find the intersection point between the hole polygon and TIN, and perform linear interpolation to obtain the 3D global coordinates of the intersection point.

[0250] S08: The seam hole polygons are saved as .shp files.

[0251] From the above description, it can be seen that the image coordinate and three-dimensional model global coordinate conversion process provided by this specific application example, based on affine transformation and three-dimensional coordinate system rotation, establishes a two-dimensional and three-dimensional transformation relationship for the fracture vector polygon that needs to be mapped according to the two-dimensional image coordinates and three-dimensional model coordinates of the three vertices of the outcrop TIN triangle where the polygon vertex is located, and then maps it to the vertex in the polygon. The overall process involves: conversion of vertex image coordinates to texture coordinates; then conversion from texture coordinates to three-dimensional local coordinates; and finally conversion from three-dimensional local coordinates to three-dimensional model global coordinates.

[0252] From the software level, in order to improve the reliability of the three-dimensional outcrop model, and thus improve the accuracy of using the three-dimensional outcrop model to determine underground mineral reservoirs and reflect geological movements, the present application provides an embodiment of a three-dimensional outcrop model determination device for implementing all or part of the content of the three-dimensional outcrop model determination method, see Fig.17 The three-dimensional outcrop model determination device specifically includes the following contents:

[0253] The acquisition module 10 is used to acquire the cross-sectional image and irregular triangulated network model corresponding to the target carbonate rock outcrop, and obtain a binary raster image after performing binarization processing on the cross-sectional image.

[0254] The texture coordinate acquisition module 20 is used to obtain the two-dimensional texture coordinates corresponding to each of the seam hole vertices of the seam hole object in the binary raster image based on the height and width of the cross-section image.

[0255] The global coordinate acquisition module 30 is used to apply the irregular triangulated network model and the two-dimensional texture coordinates to obtain the three-dimensional global coordinates of the intersection between the irregular triangulated network model and the seam hole object, and the three-dimensional global coordinates of each of the seam hole vertices.

[0256] The determination module 40 is used to determine the three-dimensional outcrop model corresponding to the target carbonate rock outcrop according to the three-dimensional global coordinates of the intersection points and the fracture-cavity vertices and the irregular triangulated network model.

[0257] In one embodiment of the present application, see Fig.18 , the texture coordinate acquisition module 20 comprises:

[0258] The raster data acquisition unit 21 is used to acquire the raster data of the seam hole object in the binary raster image.

[0259] The image coordinate acquisition unit 22 is used to apply the grid data to obtain the two-dimensional image coordinates of each of the vertices of the seam hole.

[0260] The texture coordinate acquisition unit 23 is used to obtain the two-dimensional texture coordinates corresponding to each of the vertices of the seam hole based on the two-dimensional image coordinates and the height and width of the cross-sectional image.

[0261] In one embodiment of the present application, see Fig.19 , the global coordinate acquisition module 30 comprises:

[0262] A triangle acquisition unit 31 is used to acquire triangles corresponding to the vertices of the seam holes from the irregular triangulated network model;

[0263] The first global coordinate acquisition unit 32 is used to obtain the three-dimensional global coordinates of each of the vertices of the seam hole by applying the triangle and the two-dimensional texture coordinates corresponding to each of the vertices of the seam hole.

[0264] The second global coordinate acquisition unit 33 is used to obtain the three-dimensional global coordinates of the intersection points between the irregular triangulated network model and each of the hole objects by applying the two-dimensional texture coordinates corresponding to each of the hole vertices.

[0265] In one embodiment of the present application, the triangle acquisition unit is used to perform the following:

[0266] A triangle set corresponding to the seam hole object is determined according to a minimum circumscribed rectangle of the seam hole object and an irregular triangulated network model, wherein the irregular triangulated network model includes a plurality of triangles.

[0267] From the triangle set, triangles corresponding to the respective seam hole vertices of the corresponding seam hole object are obtained.

[0268] In one embodiment of the present application, the first global coordinate acquisition unit is used to perform the following:

[0269] Obtain three-dimensional global coordinates and two-dimensional texture coordinates corresponding to the triangle vertices of the triangle.

[0270] According to the three-dimensional global coordinates corresponding to the triangle vertices of the triangle, a rotation matrix from the three-dimensional global coordinate system corresponding to the irregular triangulated network model to the three-dimensional local coordinate system corresponding to the triangle is obtained.

[0271] Based on the three-dimensional global coordinates corresponding to the triangle vertices and the rotation matrix, the three-dimensional local coordinates corresponding to the triangle vertices are obtained.

[0272] Based on the three-dimensional local coordinates and the two-dimensional texture coordinates corresponding to the triangle vertices, a target affine transformation matrix corresponding to the triangle is obtained.

[0273] The target affine transformation matrix corresponding to the triangle and the two-dimensional texture coordinates of the seam hole vertex corresponding to the triangle are applied to obtain the three-dimensional local coordinates of the seam hole vertex.

[0274] Based on the three-dimensional local coordinates of the vertex of the crack hole and the rotation matrix, the three-dimensional global coordinates of the vertex of the crack hole are obtained.

[0275] In one embodiment of the present application, the image coordinate acquisition unit is used to perform the following:

[0276] According to the preset unmarked point determination order, the target unmarked point in the binary raster image is determined. If the crack object corresponding to the target unmarked point is a crack linear object, the traversal step is performed: it is determined whether there is an unmarked point adjacent to the target unmarked point in the binary raster image. If so, the direction from the target unmarked point as the starting point to the adjacent unmarked point is recorded, and it is determined whether the direction is consistent with the direction recorded last time; if the direction is inconsistent with the direction recorded last time, the raster data corresponding to the target unmarked point is converted into two-dimensional image coordinates and saved, and the target unmarked point is modified to a marked point.

[0277] The traversal step is performed again using the adjacent unmarked point as the target unmarked point until there is no adjacent unmarked point to the target unmarked point.

[0278] According to the preset unmarked point determination order, the next target unmarked point in the binary raster image is determined, and if the crack hole object corresponding to the next target unmarked point is a crack linear object, the traversal step is performed again until all the unmarked points in the binary raster image are traversed; wherein the unmarked point is a pixel point corresponding to the crack hole object in the binary raster image and the pixel point is not marked.

[0279] In one embodiment of the present application, the image coordinate acquisition unit is further configured to execute the following:

[0280] If the direction is consistent with the direction recorded last time, the traversal step is performed again with the adjacent unmarked point as the target unmarked point until there is no adjacent marked point to the target unmarked point.

[0281] In one embodiment of the present application, the image coordinate acquisition unit is further configured to execute the following:

[0282] If the hole object corresponding to the unmarked point of the target is a hole surface object, the contour detection step is performed:

[0283] The cvFindContours function is used to perform contour detection on the hole surface object to obtain the two-dimensional image coordinates corresponding to the grid data of the seam vertices of the hole surface object.

[0284] According to the preset unmarked point determination order, the next target unmarked point in the binary raster image is determined, and if the crack hole object corresponding to the next target unmarked point is a crack linear object, the traversal step is performed again until all the unmarked points in the binary raster image are traversed; wherein the unmarked point is a pixel point corresponding to the crack hole object in the binary raster image and the pixel point is not marked.

[0285] In one embodiment of the present application, the image coordinate acquisition unit is further used to execute the contour detection step again if the seam hole object corresponding to the unmarked point of the next target is a hole surface object.

[0286] The embodiment of the three-dimensional outcrop model determination device provided in this specification can be specifically used to execute the processing flow of the embodiment of the above-mentioned three-dimensional outcrop model determination method. Its functions are not repeated here, and reference can be made to the detailed description of the embodiment of the above-mentioned three-dimensional outcrop model determination method.

[0287] From the above description, it can be seen that the three-dimensional outcrop model determination method and device provided in the present application can improve the reliability of the three-dimensional outcrop model, and thus can improve the reliability of using the three-dimensional outcrop model to determine underground mineral reservoirs and reflect geological movements; specifically, it can also improve the accuracy and efficiency of determining the three-dimensional outcrop model, save memory resources, and can improve the accuracy of simulation, characterization and reproduction of field outcrop structural characteristics.

[0288] From the hardware level, in order to improve the efficiency and flexibility of determining the three-dimensional outcrop model, and thus improve the user experience, the present application provides an embodiment of an electronic device for implementing all or part of the content of the three-dimensional outcrop model determination method. The electronic device specifically includes the following content:

[0289] Processor, memory, communication interface and bus; wherein the processor, memory and communication interface communicate with each other through the bus; the communication interface is used to realize information transmission between the three-dimensional outcrop model determination device and related devices such as user terminals; the electronic device can be a desktop computer, a tablet computer and a mobile terminal, etc., but the present embodiment is not limited thereto. In the present embodiment, the electronic device can be implemented with reference to the embodiment for realizing the three-dimensional outcrop model determination method and the embodiment for realizing the three-dimensional outcrop model determination device, and the contents are incorporated herein, and the repeated parts are not repeated.

[0290] Fig. 20 FIG. 9 is a schematic block diagram of the system structure of the electronic device 9600 according to an embodiment of the present application. Fig. 20 As shown, the electronic device 9600 may include a central processor 9100 and a memory 9140; the memory 9140 is coupled to the central processor 9100. It is worth noting that Fig. 20 is exemplary; other types of structures may also be used to supplement or replace this structure to implement telecommunication functions or other functions.

[0291] In one or more embodiments of the present application, the 3D outcrop model determination function may be integrated into the central processor 9100. The central processor 9100 may be configured to perform the following control:

[0292] Step 101: Obtain a cross-sectional image and an irregular triangulated network model corresponding to a target carbonate rock outcrop, and perform binarization processing on the cross-sectional image to obtain a binary raster image.

[0293] Step 102: Based on the height and width of the cross-sectional image, obtain the two-dimensional texture coordinates corresponding to each of the seam hole vertices of the seam hole object in the binary raster image.

[0294] Step 103: Apply the irregular triangulated network model and the two-dimensional texture coordinates to obtain the three-dimensional global coordinates of the intersection between the irregular triangulated network model and the seam hole object, and the three-dimensional global coordinates of each of the seam hole vertices.

[0295] Step 104: Determine a three-dimensional outcrop model corresponding to the target carbonate outcrop according to the three-dimensional global coordinates of the intersection points and the fracture-cavity vertices and the irregular triangulated network model.

[0296] From the above description, it can be seen that the electronic device provided by the embodiment of the present application can improve the reliability of the three-dimensional outcrop model, and further improve the reliability of using the three-dimensional outcrop model to determine underground mineral reservoirs and reflect geological movements.

[0297] In another embodiment, the three-dimensional outcrop model determination device can be configured separately from the central processing unit 9100. For example, the three-dimensional outcrop model determination device can be configured as a chip connected to the central processing unit 9100, and the three-dimensional outcrop model determination function is realized through the control of the central processing unit.

[0298] like Fig. 20 As shown, the electronic device 9600 may also include: a communication module 9110, an input unit 9120, an audio processor 9130, a display 9160, and a power supply 9170. It is worth noting that the electronic device 9600 does not necessarily have to include Fig. 20 In addition, the electronic device 9600 may also include Fig. 20 For components not shown, reference may be made to the prior art.

[0299] like Fig. 20 As shown, the central processing unit 9100 is sometimes also referred to as a controller or an operation control, and may include a microprocessor or other processor device and / or logic device. The central processing unit 9100 receives input and controls the operation of various components of the electronic device 9600.

[0300] The memory 9140 may be, for example, one or more of a cache, a flash memory, a hard drive, a removable medium, a volatile memory, a non-volatile memory or other suitable devices. The above-mentioned information related to the failure may be stored, and a program for executing the relevant information may also be stored. The CPU 9100 may execute the program stored in the memory 9140 to implement information storage or processing, etc.

[0301] The input unit 9120 provides input to the central processing unit 9100. The input unit 9120 is, for example, a key or a touch input device. The power supply 9170 is used to provide power to the electronic device 9600. The display 9160 is used to display display objects such as images and texts. The display may be, for example, an LCD display, but is not limited thereto.

[0302] The memory 9140 may be a solid-state memory, such as a read-only memory (ROM), a random access memory (RAM), a SIM card, etc. It may also be a memory that saves information even when the power is off, can be selectively erased, and is provided with more data, examples of which are sometimes referred to as EPROMs, etc. The memory 9140 may also be some other type of device. The memory 9140 includes a buffer memory 9141 (sometimes referred to as a buffer). The memory 9140 may include an application / function storage unit 9142, which is used to store application programs and function programs or processes for executing the operation of the electronic device 9600 through the central processor 9100.

[0303] The memory 9140 may also include a data storage unit 9143 for storing data, such as contacts, digital data, pictures, sounds, and / or any other data used by the electronic device. The driver storage unit 9144 of the memory 9140 may include various drivers for communication functions of the electronic device and / or for executing other functions of the electronic device (such as messaging applications, address book applications, etc.).

[0304] The communication module 9110 is a transmitter / receiver 9110 that sends and receives signals via an antenna 9111. The communication module (transmitter / receiver) 9110 is coupled to the central processor 9100 to provide input signals and receive output signals, which may be the same as the case of a conventional mobile communication terminal.

[0305] Based on different communication technologies, multiple communication modules 9110 may be provided in the same electronic device, such as a cellular network module, a Bluetooth module and / or a wireless LAN module, etc. The communication module (transmitter / receiver) 9110 is also coupled to a speaker 9131 and a microphone 9132 via an audio processor 9130 to provide an audio output via the speaker 9131 and receive an audio input from the microphone 9132, thereby realizing a common telecommunication function. The audio processor 9130 may include any suitable buffer, decoder, amplifier, etc. In addition, the audio processor 9130 is also coupled to the central processor 9100, so that recording can be performed on the local machine through the microphone 9132, and the sound stored on the local machine can be played through the speaker 9131.

[0306] It can be seen from the above description that the electronic device provided by the embodiment of the present application can improve the reliability of the three-dimensional outcrop model, and further can improve the reliability of using the three-dimensional outcrop model to determine underground mineral reservoirs and reflect geological movements.

[0307] The embodiments of the present application also provide a computer-readable storage medium capable of implementing all the steps in the method for determining a three-dimensional outcrop model in the above embodiments. The computer-readable storage medium stores a computer program. When the computer program is executed by a processor, all the steps in the method for determining a three-dimensional outcrop model in the above embodiments are implemented. For example, when the processor executes the computer program, the following steps are implemented:

[0308] Step 101: Obtain a cross-sectional image and an irregular triangulated network model corresponding to a target carbonate rock outcrop, and perform binarization processing on the cross-sectional image to obtain a binary raster image.

[0309] Step 102: Based on the height and width of the cross-sectional image, obtain the two-dimensional texture coordinates corresponding to each of the seam hole vertices of the seam hole object in the binary raster image.

[0310] Step 103: Apply the irregular triangulated network model and the two-dimensional texture coordinates to obtain the three-dimensional global coordinates of the intersection between the irregular triangulated network model and the seam hole object, and the three-dimensional global coordinates of each of the seam hole vertices.

[0311] Step 104: Determine a three-dimensional outcrop model corresponding to the target carbonate outcrop according to the three-dimensional global coordinates of the intersection points and the fracture-cavity vertices and the irregular triangulated network model.

[0312] From the above description, it can be seen that the computer-readable storage medium provided in the embodiment of the present application can improve the reliability of the three-dimensional outcrop model, and further can improve the reliability of using the three-dimensional outcrop model to determine underground mineral reservoirs and reflect geological movements.

[0313] In this application, each embodiment of the above method is described in a progressive manner, and the same or similar parts between the embodiments can be referred to each other, and each embodiment focuses on the differences from other embodiments. For relevant parts, refer to the partial description of the method embodiment.

[0314] Those skilled in the art will appreciate that the embodiments of the present application may be provided as methods, systems, or computer program products. Therefore, the present application may adopt the form of a complete hardware embodiment, a complete software embodiment, or an embodiment in combination with software and hardware. Moreover, the present application may adopt the form of a computer program product implemented in one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) that contain computer-usable program code.

[0315] The present application is described with reference to the flowcharts and / or block diagrams of the methods, devices (systems), and computer program products according to the embodiments of the present application. It should be understood that each process and / or box in the flowchart and / or block diagram, as well as the combination of the processes and / or boxes in the flowchart and / or block diagram, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing device to generate a machine, so that the instructions executed by the processor of the computer or other programmable data processing device generate instructions for implementing the processes in the flowchart and / or block diagram. Figure 1 A process or multiple processes and / or boxes Figure 1 A device that provides the functions specified in a block or multiple blocks.

[0316] These computer program instructions may also be stored in a computer-readable memory capable of directing a computer or other programmable data processing device to operate in a specific manner, so that the instructions stored in the computer-readable memory produce an article of manufacture comprising an instruction device, which implements the process Figure 1 A process or multiple processes and / or boxes Figure 1 A function specified in one or more boxes.

[0317] These computer program instructions can also be loaded onto a computer or other programmable data processing device so that a series of operating steps are executed on the computer or other programmable device to produce a computer-implemented process, thereby providing instructions for implementing the process. Figure 1 A process or multiple processes and / or boxes Figure 1 The steps for the functions specified in one or more boxes.

[0318] Specific embodiments are used in this application to illustrate the principles and implementation methods of this application. The description of the above embodiments is only used to help understand the method and core idea of ​​this application. At the same time, for those skilled in the art, according to the idea of ​​this application, there will be changes in the specific implementation methods and application scope. In summary, the content of this specification should not be understood as a limitation on this application.< / cv::point> < / cv::point> < / point> < / point> < / point> < / point> < / point> < / point>

Claims

1. A method for determining a three-dimensional outcrop model, characterized in that: include: Obtain the cross-sectional image and irregular triangulated network model corresponding to the target carbonate rock outcrop, and obtain a binary raster image after binarization processing of the cross-sectional image; Based on the height and width of the cross-section image, obtain the two-dimensional texture coordinates corresponding to each of the vertices of the hole object in the binary raster image; Applying the irregular triangulated network model and the two-dimensional texture coordinates, obtaining the three-dimensional global coordinates of the intersection between the irregular triangulated network model and the seam hole object, and the three-dimensional global coordinates of each of the seam hole vertices; According to the three-dimensional global coordinates of the intersection points and the fracture-cavity vertices and the irregular triangulated network model, a three-dimensional outcrop model corresponding to the target carbonate outcrop is determined.

2. The three-dimensional outcrop model determination method according to claim 1, characterized in that: The step of obtaining the two-dimensional texture coordinates corresponding to each of the vertices of the hole object in the binary raster image based on the height and width of the cross-section image comprises: Obtaining raster data of the seam hole object in the binary raster image; Using the grid data, the two-dimensional image coordinates of each of the vertices of the crack hole are obtained; Based on the two-dimensional image coordinates, the height and width of the cross-sectional image, the two-dimensional texture coordinates corresponding to each of the vertices of the seam hole are obtained.

3. The three-dimensional outcrop model determination method according to claim 1, characterized in that: The applying of the irregular triangulated network model and the two-dimensional texture coordinates to obtain the three-dimensional global coordinates of the intersection between the irregular triangulated network model and the seam hole object, and the three-dimensional global coordinates of each of the seam hole vertices, includes: Acquire triangles corresponding to the vertices of the seam holes from the irregular triangulated network model; Applying the triangles and two-dimensional texture coordinates corresponding to the vertices of each of the seam holes to obtain the three-dimensional global coordinates of the vertices of each of the seam holes; The three-dimensional global coordinates of the intersection points between the irregular triangulated network model and each of the hole objects are obtained by using the two-dimensional texture coordinates corresponding to each of the hole vertices.

4. The three-dimensional outcrop model determination method according to claim 3, characterized in that: The step of obtaining triangles corresponding to the vertices of the seam holes from the irregular triangulated network model includes: Determine a triangle set corresponding to the seam hole object according to the minimum circumscribed rectangle of the seam hole object and an irregular triangulated network model, wherein the irregular triangulated network model includes a plurality of triangles; From the triangle set, triangles corresponding to the respective seam hole vertices of the corresponding seam hole object are obtained.

5. The method for determining a three-dimensional outcrop model according to claim 3, characterized in that: The applying the triangles and two-dimensional texture coordinates corresponding to each of the vertices of the seam hole to obtain the three-dimensional global coordinates of each of the vertices of the seam hole comprises: Obtaining three-dimensional global coordinates and two-dimensional texture coordinates corresponding to the triangle vertices of the triangle; According to the three-dimensional global coordinates corresponding to the triangle vertices of the triangle, a rotation matrix from the three-dimensional global coordinate system corresponding to the irregular triangulated network model to the three-dimensional local coordinate system corresponding to the triangle is obtained; Based on the three-dimensional global coordinates corresponding to the triangle vertices and the rotation matrix, obtain the three-dimensional local coordinates corresponding to the triangle vertices; Based on the three-dimensional local coordinates and the two-dimensional texture coordinates corresponding to the triangle vertices, obtaining a target affine transformation matrix corresponding to the triangle; Applying the target affine transformation matrix corresponding to the triangle and the two-dimensional texture coordinates of the seam hole vertex corresponding to the triangle to obtain the three-dimensional local coordinates of the seam hole vertex; Based on the three-dimensional local coordinates of the vertex of the crack hole and the rotation matrix, the three-dimensional global coordinates of the vertex of the crack hole are obtained.

6. The method for determining a three-dimensional outcrop model according to claim 2, characterized in that: The applying of the grid data to obtain the two-dimensional image coordinates of each of the vertices of the seam hole includes: According to the preset unmarked point determination order, a target unmarked point is determined from multiple unmarked points in the binary raster image. If the crack hole object corresponding to the target unmarked point is a crack linear object, the traversal step is performed: Determine whether there is an unmarked point adjacent to the target unmarked point in the binary raster image, and if so, record the direction from the target unmarked point as the starting point to the adjacent unmarked point, and determine whether the direction is consistent with the direction recorded last time; If the direction is inconsistent with the direction recorded last time, the grid data corresponding to the target unmarked point is converted into two-dimensional image coordinates and saved, and the target unmarked point is modified into a marked point; The traversal step is performed again using the adjacent unmarked point as the target unmarked point until there is no adjacent unmarked point to the target unmarked point; According to the preset unmarked point determination order, the next target unmarked point in the binary raster image is determined, and if the crack hole object corresponding to the next target unmarked point is a crack linear object, the traversal step is performed again until all the unmarked points in the binary raster image are traversed; The unmarked points are pixels corresponding to the seam hole objects in the binary raster image and the pixels are not marked.

7. The method for determining a three-dimensional outcrop model according to claim 6, characterized in that: Also includes: If the direction is consistent with the direction recorded last time, the traversal step is performed again with the adjacent unmarked point as the target unmarked point until there is no adjacent marked point to the target unmarked point.

8. The method for determining a three-dimensional outcrop model according to claim 6, characterized in that: The determining whether there is an unmarked point adjacent to the target unmarked point in the binary raster image includes: It is determined whether there is an unmarked point adjacent to the target unmarked point in the direction recorded last time; if not, it is determined whether there is an unmarked point adjacent to the target unmarked point in the direction rotated by a preset angle.

9. The method for determining a three-dimensional outcrop model according to claim 6, characterized in that: Also includes: If the hole object corresponding to the unmarked point of the target is a hole surface object, the contour detection step is performed: The cvFindContours function is used to perform contour detection on the hole surface object, and the two-dimensional image coordinates corresponding to the grid data of the crack hole vertices of the hole surface object are obtained; According to the preset unmarked point determination order, the next target unmarked point in the binary raster image is determined, and if the crack hole object corresponding to the next target unmarked point is a crack linear object, the traversal step is performed again until all the unmarked points in the binary raster image are traversed; The unmarked points are pixels corresponding to the seam hole objects in the binary raster image and the pixels are not marked.

10. The method for determining a three-dimensional outcrop model according to claim 9, characterized in that: Also includes: If the seam hole object corresponding to the unmarked point of the next target is a hole surface object, the contour detection step is performed again.

11. A three-dimensional outcrop model determination device, characterized in that: include: An acquisition module is used to acquire a cross-sectional image and an irregular triangulated network model corresponding to a target carbonate rock outcrop, and obtain a binary raster image after performing a binarization process on the cross-sectional image; A texture coordinate acquisition module, used for obtaining two-dimensional texture coordinates corresponding to each seam hole vertex of the seam hole object in the binary raster image based on the height and width of the cross-section image; A global coordinate acquisition module, used to apply the irregular triangulated network model and the two-dimensional texture coordinates to obtain the three-dimensional global coordinates of the intersection between the irregular triangulated network model and the seam hole object, and the three-dimensional global coordinates of each of the seam hole vertices; A determination module is used to determine the three-dimensional outcrop model corresponding to the target carbonate rock outcrop according to the three-dimensional global coordinates of the intersection points and the fracture-hole vertices and the irregular triangulated network model.

12. The three-dimensional outcrop model determination device according to claim 11, characterized in that: The texture coordinate acquisition module comprises: A raster data acquisition unit, used for acquiring raster data of seam hole objects in the binary raster image; An image coordinate acquisition unit, used to obtain the two-dimensional image coordinates of each of the vertices of the seam hole by applying the grid data; The texture coordinate acquisition unit is used to obtain the two-dimensional texture coordinates corresponding to each of the vertices of the seam hole based on the two-dimensional image coordinates and the height and width of the cross-sectional image.

13. The three-dimensional outcrop model determination device according to claim 11, characterized in that: The global coordinate acquisition module comprises: A triangle acquisition unit, used for acquiring triangles corresponding to the vertices of the seam holes from the irregular triangulated network model; A first global coordinate acquisition unit, used for obtaining the three-dimensional global coordinates of each of the vertices of the seam hole by applying the triangles and two-dimensional texture coordinates corresponding to each of the vertices of the seam hole; The second global coordinate acquisition unit is used to apply the two-dimensional texture coordinates corresponding to each of the seam hole vertices to obtain the three-dimensional global coordinates of the intersection between the irregular triangulated network model and each of the seam hole objects.

14. The three-dimensional outcrop model determination device according to claim 13, characterized in that: The triangle acquisition unit is used to perform the following: Determine a triangle set corresponding to the seam hole object according to the minimum circumscribed rectangle of the seam hole object and an irregular triangulated network model, wherein the irregular triangulated network model includes a plurality of triangles; From the triangle set, triangles corresponding to the respective seam hole vertices of the corresponding seam hole object are obtained.

15. The three-dimensional outcrop model determination device according to claim 13, characterized in that: The first global coordinate acquisition unit is used to execute the following: Obtaining three-dimensional global coordinates and two-dimensional texture coordinates corresponding to the triangle vertices of the triangle; According to the three-dimensional global coordinates corresponding to the triangle vertices of the triangle, a rotation matrix from the three-dimensional global coordinate system corresponding to the irregular triangulated network model to the three-dimensional local coordinate system corresponding to the triangle is obtained; Based on the three-dimensional global coordinates corresponding to the triangle vertices and the rotation matrix, obtain the three-dimensional local coordinates corresponding to the triangle vertices; Based on the three-dimensional local coordinates and the two-dimensional texture coordinates corresponding to the triangle vertices, obtaining a target affine transformation matrix corresponding to the triangle; Applying the target affine transformation matrix corresponding to the triangle and the two-dimensional texture coordinates of the seam hole vertex corresponding to the triangle to obtain the three-dimensional local coordinates of the seam hole vertex; Based on the three-dimensional local coordinates of the vertex of the crack hole and the rotation matrix, the three-dimensional global coordinates of the vertex of the crack hole are obtained.

16. The three-dimensional outcrop model determination device according to claim 12, characterized in that: The image coordinate acquisition unit is used to perform the following: According to the preset unmarked point determination order, the target unmarked point in the binary raster image is determined. If the crack hole object corresponding to the target unmarked point is a crack linear object, the traversal step is performed: Determine whether there is an unmarked point adjacent to the target unmarked point in the binary raster image, and if so, record the direction from the target unmarked point as the starting point to the adjacent unmarked point, and determine whether the direction is consistent with the direction recorded last time; If the direction is inconsistent with the direction recorded last time, the grid data corresponding to the target unmarked point is converted into two-dimensional image coordinates and saved, and the target unmarked point is modified into a marked point; The traversal step is performed again using the adjacent unmarked point as the target unmarked point until there is no adjacent unmarked point to the target unmarked point; According to the preset unmarked point determination order, the next target unmarked point in the binary raster image is determined, and if the crack hole object corresponding to the next target unmarked point is a crack linear object, the traversal step is performed again until all the unmarked points in the binary raster image are traversed; The unmarked points are pixels corresponding to the seam hole objects in the binary raster image and the pixels are not marked.

17. The three-dimensional outcrop model determination device according to claim 16, characterized in that: The image coordinate acquisition unit is further used to execute the following: If the direction is consistent with the direction recorded last time, the traversal step is performed again with the adjacent unmarked point as the target unmarked point until there is no adjacent marked point to the target unmarked point.

18. The three-dimensional outcrop model determination device according to claim 16, characterized in that: The image coordinate acquisition unit is further used to execute the following: If the hole object corresponding to the unmarked point of the target is a hole surface object, the contour detection step is performed: The cvFindContours function is used to perform contour detection on the hole surface object, and the two-dimensional image coordinates corresponding to the grid data of the crack hole vertices of the hole surface object are obtained; According to the preset unmarked point determination order, the next target unmarked point in the binary raster image is determined, and if the crack hole object corresponding to the next target unmarked point is a crack linear object, the traversal step is performed again until all the unmarked points in the binary raster image are traversed; The unmarked points are pixels corresponding to the seam hole objects in the binary raster image and the pixels are not marked.

19. The three-dimensional outcrop model determination device according to claim 18, characterized in that: The image coordinate acquisition unit is further used to execute the contour detection step again if the seam hole object corresponding to the unmarked point of the next target is a hole surface object.

20. An electronic device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, characterized in that: When the processor executes the program, the three-dimensional outcrop model determination method according to any one of claims 1 to 10 is implemented.

21. A computer-readable storage medium having computer instructions stored thereon, characterized in that: When the instructions are executed, the three-dimensional outcrop model determination method according to any one of claims 1 to 10 is implemented.

Images