A three-dimensional geological modeling method, device and electronic equipment

Abstract

The present disclosure relates to the field of computer technology, and provides a three-dimensional geological modeling method, device and electronic device. In order to solve the technical problem that the fault characteristics lead to inaccurate results of three-dimensional geological modeling, the method includes: obtaining first sampling point data and fault line data representing a preset area; dividing the preset area into at least two geological continuous blocks according to the fault line data; determining a first three-dimensional geological model of the preset area according to the first sampling point data of the first sampling point in each geological continuous block; the three-dimensional geological model represents the geological attribute values ​​of each grid point in the preset area; obtaining a second three-dimensional geological model of the preset area according to the first three-dimensional geological model corresponding to each geological continuous block; and obtaining a three-dimensional geological modeling result of the preset area according to the second three-dimensional geological model. In this way, when the preset area includes fault terrain, the geological attributes in the preset area can be accurately modeled in three dimensions.

Description

Technical Field

[0001] The present disclosure relates to the field of computer technology, and more specifically, to a three-dimensional geological modeling method, device, and electronic equipment. Background Art

[0002] Three-dimensional geological modeling is a technology used to study the field of earth science by using computer technology to combine tools such as spatial information management, geological interpretation, spatial analysis and prediction, earth science statistics, entity content analysis and graphic visualization.

[0003] At present, the existing geological 3D modeling technology commonly used in the industry is: geological 3D modeling based on a small number of drill holes.

[0004] Faults are weak points formed in natural mountain masses by tectonic movement. Intersections between fault planes and engineering excavation surfaces or other fault planes can negatively impact slope stability and are crucial for engineering survey and design. However, existing techniques often inaccurate 3D geological modeling results due to fault characteristics within the modeling area. Summary of the Invention

[0005] An object of the present disclosure is to provide a three-dimensional geological modeling method, apparatus, and electronic device.

[0006] According to a first aspect of the present disclosure, a three-dimensional geological modeling method is provided, comprising:

[0007] Acquire first sampling point data representing a geological attribute value of a first sampling point in a preset area and fault line data for describing at least one fault line in the preset area;

[0008] Dividing the preset area into at least two geologically continuous blocks according to the fault line data;

[0009] Determining a first three-dimensional geological model of the preset area based on first sampling point data of the first sampling point in each geological continuous block; the three-dimensional geological model represents geological attribute values ​​of each grid point in the preset area;

[0010] Obtaining a second three-dimensional geological model of the preset area according to the first three-dimensional geological model corresponding to each geological continuous block;

[0011] A three-dimensional geological modeling result of the preset area is obtained according to the second three-dimensional geological model.

[0012] Optionally, obtaining the second three-dimensional geological model of the preset area according to the first three-dimensional geological model corresponding to each geological continuous block includes:

[0013] Determining, based on the fault line data, a weight of each of the grid points in the preset area corresponding to each geological continuous block;

[0014] A second three-dimensional geological model of the preset area is obtained according to the weight of each grid point corresponding to each geological continuous block and the first three-dimensional geological model corresponding to each geological continuous block.

[0015] Optionally, determining the weight of each grid point corresponding to each geological continuous block in the preset area according to the fault line data includes:

[0016] Determining the minimum distance between a grid point in each geological continuous block and the fault line based on the fault line data;

[0017] Obtaining a weight of each grid point corresponding to each geological continuous block according to a minimum distance between the grid point in each geological continuous block and the fault line;

[0018] The weight of each grid point corresponding to each geological continuous block is normalized.

[0019] Optionally, obtaining the weight of each grid point corresponding to each geological continuous block according to the minimum distance from the grid point in each geological continuous block to the fault line includes:

[0020] When any grid point is within the first geological continuous block and not within the second geological continuous block, determining a weight of the any grid point corresponding to the second geological continuous block to be a set value;

[0021] Acquire first mapping data representing a mapping relationship between distance and weight;

[0022] A weight of the any grid point corresponding to the first geological continuous block is obtained according to the minimum distance from the any grid point to the fault line and the first mapping data.

[0023] Optionally, obtaining the second three-dimensional geological model of the preset area according to the weight of each grid point corresponding to each geological continuous block and the first three-dimensional geological model corresponding to each geological continuous block includes:

[0024] For each of the grid points, a second geological attribute value of the corresponding grid point is obtained according to the first geological attribute value of the corresponding grid point in each first three-dimensional geological model and the weight of each geological continuous block corresponding to the corresponding grid point;

[0025] The second three-dimensional geological model is obtained according to the second geological attribute value of each grid point in the preset area.

[0026] Optionally, dividing the preset area into at least two geologically continuous blocks according to the fault line data includes:

[0027] generating a fault line extension line from an endpoint of the fault line to a boundary of the preset area;

[0028] The preset area is divided into at least two geologically continuous blocks according to the fault line and the extended line of the fault line.

[0029] Optionally, obtaining a three-dimensional geological modeling result of the preset area according to the second three-dimensional geological model includes:

[0030] acquiring second sampling point data representing a geological attribute value of a second sampling point according to the second three-dimensional geological model, wherein the second sampling point is located on an extension line of the fault line;

[0031] Determining a third three-dimensional geological model of the preset area based on the first sampling point data and the second sampling point data of the first sampling point in each geological continuous block;

[0032] Obtaining a fourth three-dimensional geological model of the preset area according to the third three-dimensional geological model corresponding to each geological continuous block;

[0033] The fourth three-dimensional geological model is used as the three-dimensional geological modeling result of the preset area.

[0034] Optionally, determining the first three-dimensional geological model of the preset area according to the first sampling point data of the first sampling point in each geological continuous block includes:

[0035] Based on a difference algorithm corresponding to the attributes of the geological attribute values, a first three-dimensional geological model of the preset area is determined according to first sampling point data of the first sampling point in each geological continuous block.

[0036] According to a second aspect of the present disclosure, a three-dimensional geological modeling device is provided, comprising:

[0037] a data acquisition module, configured to acquire first sampling point data representing a geological attribute value of a first sampling point within a preset area and fault line data describing at least one fault line within the preset area;

[0038] A geological division module, configured to divide the preset area into at least two geologically continuous blocks according to the fault line data;

[0039] a first modeling module, configured to determine a first three-dimensional geological model of the preset area based on first sampling point data of a first sampling point in each geological continuous block; the three-dimensional geological model represents geological attribute values ​​of each grid point in the preset area;

[0040] A second modeling module is used to obtain a second three-dimensional geological model of the preset area according to the first three-dimensional geological model corresponding to each geological continuous block;

[0041] The result determination module is used to obtain the three-dimensional geological modeling result of the preset area according to the second three-dimensional geological model.

[0042] According to a third aspect of the present disclosure, there is provided an electronic device, including:

[0043] The device according to the second aspect of the present disclosure; or

[0044] A processor and a memory, wherein the memory is used to store a computer program, and the computer program is used to control the processor to execute the method according to the first aspect of the present disclosure.

[0045] According to a fourth aspect of the present disclosure, a computer-readable storage medium is provided, on which a computer program is stored. When the computer program is executed by a processor, the method according to the first aspect of the present disclosure is implemented.

[0046] Through this embodiment, the preset area is divided into at least two geologically continuous blocks based on fault line data used to describe at least one fault line in the preset area. A first three-dimensional geological model of the preset area is determined based on the first sampling point data of the first sampling point in each geologically continuous block. A second three-dimensional geological model of the preset area is obtained based on the first three-dimensional geological model corresponding to each geologically continuous block, and a three-dimensional geological modeling result of the preset area is obtained based on the second three-dimensional geological model. In this way, even if the preset area includes fault terrain, accurate three-dimensional modeling of the geological properties in the preset area can be performed.

[0047] Other features and advantages of the present disclosure will become apparent from the following detailed description of exemplary embodiments of the present disclosure with reference to the accompanying drawings. BRIEF DESCRIPTION OF THE DRAWINGS

[0048] The accompanying drawings, which are incorporated in and constitute a part of this specification, illustrate embodiments of the present disclosure and, together with the description, serve to explain the principles of the present disclosure.

[0049] Figure 1 A block diagram showing a hardware configuration of an electronic device that can implement an embodiment of the present disclosure;

[0050] Figure 2 A flow chart of a three-dimensional geological modeling method provided by one embodiment of the present disclosure is shown;

[0051] Figure 3 A schematic diagram showing a drawing provided by one embodiment of the present disclosure is shown;

[0052] Figure 4 A flow chart of a three-dimensional geological modeling method provided by another embodiment of the present disclosure is shown;

[0053] Figure 5 A block diagram of a three-dimensional geological modeling device provided by an embodiment of the present disclosure is shown;

[0054] Figure 6 A block diagram of an electronic device provided by an embodiment of the present disclosure is shown. DETAILED DESCRIPTION

[0055] Various exemplary embodiments of the present disclosure will now be described in detail with reference to the accompanying drawings. It should be noted that unless otherwise specifically stated, the relative arrangement of components and steps, numerical expressions and numerical values ​​set forth in these embodiments do not limit the scope of the present disclosure.

[0056] The following description of at least one exemplary embodiment is merely illustrative in nature and is in no way intended to limit the present disclosure, its application, or uses.

[0057] Technologies, methods, and equipment known to ordinary technicians in the relevant art may not be discussed in detail, but where appropriate, the technologies, methods, and equipment should be considered part of the specification.

[0058] In all examples shown and discussed herein, any specific values ​​should be interpreted as merely exemplary and not limiting. Therefore, other examples of the exemplary embodiments may have different values.

[0059] It should be noted that like reference numerals and letters refer to like items in the following figures, and therefore, once an item is defined in one figure, it need not be further discussed in subsequent figures.

[0060] <Hardware Configuration>

[0061] Figure 1 is a block diagram illustrating a hardware configuration of an electronic device 1000 that can implement an embodiment of the present disclosure.

[0062] The electronic device 1000 may be a portable computer, a desktop computer, a mobile phone, a tablet computer, etc. Figure 1As shown, electronic device 1000 may include a processor 1100, a memory 1200, an interface device 1300, a communication device 1400, a display device 1500, an input device 1600, a speaker 1700, a microphone 1800, and the like. Processor 1100 may be a CPU, a microprocessor MCU, or the like. Memory 1200 may include, for example, ROM (read-only memory), RAM (random access memory), or a non-volatile memory such as a hard disk. Interface device 1300 may include, for example, a USB port or a headphone jack. Communication device 1400 may be capable of wired or wireless communication, specifically, Wi-Fi, Bluetooth, or 2G / 3G / 4G / 5G communication. Display device 1500 may be, for example, an LCD display or a touchscreen display. Input device 1600 may include, for example, a touchscreen, a keyboard, or a motion sensor input device. Users may input and output voice information through speaker 1700 and microphone 1800.

[0063] Figure 1 The electronic device shown is merely illustrative and does not in any way limit the present disclosure, its application or use. In the embodiments of the present disclosure, the memory 1200 of the electronic device 1000 is used to store instructions, which are used to control the processor 1100 to operate to perform any of the three-dimensional geological modeling methods provided in the embodiments of the present disclosure. It should be understood by those skilled in the art that although Figure 1 While multiple devices are shown for electronic device 1000, the present disclosure may only relate to some of these devices. For example, electronic device 1000 may only relate to processor 1100 and storage device 1200. A skilled person can design instructions based on the solutions disclosed in this disclosure. How instructions control processor operations is well known in the art and will not be described in detail here.

[0064] <Method Example>

[0065] In this embodiment, a three-dimensional geological modeling method is provided. The three-dimensional geological modeling method can be implemented by an electronic device. The electronic device can be Figure 1 The electronic device 1000 is shown.

[0066] according to Figure 2 As shown, the 3D geological modeling method of this embodiment may include the following steps S2100 to S2500:

[0067] Step S2100 : obtaining first sampling point data representing a geological attribute value of a first sampling point in a preset area and fault line data for describing at least one fault line in the preset area.

[0068] The preset area of ​​this embodiment may be set in advance according to application scenarios or specific requirements.

[0069] In this embodiment, the first sampling point data of each first sampling point may include the geographic coordinates of the first sampling point and the collected geographic attribute values, wherein the geographic attribute values ​​may be any one or more of elevation values, soil biochemical properties, and mineral content.

[0070] In this embodiment, the fault line can be a polyline. The fault line data for a fault line can include multiple geographic coordinates describing the fault line. The multiple geographic coordinates describing the fault line can include the geographic coordinates of the fault line's endpoints and inflection points. Based on this, the fault line can be a polyline formed by sequentially connecting multiple geographic coordinates.

[0071] In this embodiment, the first sampling point data and the fault line data of the first sampling point may be obtained through measurement.

[0072] In one example, the first sampling point can be Figure 3 As shown in the points in part (a), the fault line can be as follows Figure 3 As shown by the broken line in part (a).

[0073] Step S2200: Divide the preset area into at least two geological continuous blocks according to the fault line data.

[0074] In one embodiment of the present disclosure, the endpoints of the fault line may be located on the boundary of a preset area or outside the preset area, and the preset area may be directly divided into at least two geologically continuous blocks based on the fault line data.

[0075] In another embodiment of the present disclosure, the plurality of geographic coordinates used to describe the fault line are all located within a preset area, and the preset area cannot be directly divided into at least two geologically continuous blocks. Then, dividing the preset area into at least two geologically continuous blocks based on the fault line data may include steps S2210 to S2220 as follows:

[0076] Step S2210 : generating a fault line extension line from the endpoint of the fault line to the boundary of a preset area.

[0077] In this embodiment, starting from the endpoint of the fault line and following the direction from the inflection point adjacent to the endpoint of the fault line to the endpoint, a line segment from the endpoint to the boundary of the preset area is generated, which is the fault line extension line.

[0078] For a fault line, if one endpoint of the fault line is located within a preset area, a fault line extension line is generated with the endpoint as the starting point; if both endpoints of the fault line are located within the preset area, two fault line extension lines are generated with the two endpoints as the starting points respectively.

[0079] In one example, the fault line extension line can be as follows Figure 3 Indicated by the dotted line in part (b).

[0080] Step S2220: Divide the preset area into at least two geological continuous blocks according to the fault line and the fault line extension line.

[0081] In this embodiment, the boundary of any geological continuous block may be composed of at least one of a fault line, an extension of the fault line, and a boundary of a preset area. In other words, any geological continuous block does not have a fault line or an extension of the fault line.

[0082] In this embodiment, the geological volume is homogeneous within each geological contiguous block.

[0083] In such Figure 3 In the example shown in (b), the preset area can be divided into two geological continuous blocks based on a fault line and the corresponding fault line extension line.

[0084] Step S2300 : determining a first three-dimensional geological model of a preset area according to first sampling point data of a first sampling point in each geological continuous block.

[0085] The three-dimensional geological model represents the geological attribute values ​​of each grid point in the preset area. The first three-dimensional geological model represents the first geological attribute values ​​of each grid point in the preset area.

[0086] In this embodiment, evenly distributed grid points may be pre-set in a preset area, and the geographical location of each grid point may be obtained.

[0087] In one embodiment of the present disclosure, for each geological continuous block, a first three-dimensional geological model corresponding to the geological continuous block may be determined based on the geographical location of the grid points within the geological continuous block and the first sampling point data of the first sampling point within the geological continuous block.

[0088] Specifically, in the division In the case of geological continuous blocks, it can be for the kth (k is less than or equal to (a positive integer) geological continuous blocks, and according to the geographical locations of the grid points in the kth geological continuous block and the first sampling point data of the first sampling point in the kth geological continuous block, obtain the first geological attribute value of each grid point in the preset area corresponding to the kth geological continuous block.

[0089] In one embodiment of the present disclosure, determining a first three-dimensional geological model of a preset area based on first sampling point data of a first sampling point in each geological continuous block may include: determining the first three-dimensional geological model of the preset area based on first sampling point data of a first sampling point in each geological continuous block based on a difference algorithm corresponding to an attribute of a geological attribute value.

[0090] The difference algorithm may include any one or more of the least squares distance weighted interpolation method, the linear interpolation method, and the Kriging difference method.

[0091] When the geological attribute value is an elevation value, the corresponding difference algorithm can be the Kriging difference method.

[0092] Step S2400 : obtaining a second three-dimensional geological model of a preset area according to the first three-dimensional geological model corresponding to each geological continuous block.

[0093] In one embodiment of the present disclosure, obtaining a second three-dimensional geological model of a preset area based on the first three-dimensional geological model corresponding to each geological continuous block may include: determining an average value of the first three-dimensional geological model corresponding to each geological continuous block to obtain the second three-dimensional geological model.

[0094] Specifically, for each grid point in the preset area, the average value of the first geological attribute value corresponding to each geological continuous block of the grid point can be determined as the second geological attribute value of the grid point; and the second three-dimensional geological model can be obtained based on the second attribute values ​​of all grid points in the preset area.

[0095] In another embodiment of the present disclosure, obtaining a second three-dimensional geological model of a preset area according to the first three-dimensional geological model corresponding to each geological continuous block includes the following steps S2410 to S2420:

[0096] Step S2410: Determine the weight of each grid point corresponding to each geological continuous block in a preset area based on the fault line data.

[0097] In one embodiment of the present disclosure, determining the weight of each grid point corresponding to each geological continuous block in a preset area based on fault line data may include the following steps S2411 to S2412:

[0098] Step S2411: Determine the minimum distance from a grid point to a fault line in each geological continuous block based on the fault line data.

[0099] In this embodiment, when the preset area includes at least two fault lines, the minimum distance from each grid point to each fault line may be determined for each grid point, and the smallest distance is used as the minimum distance from the grid point to the fault line.

[0100] Step S2412: Obtain the weight of each grid point corresponding to each geological continuous block based on the minimum distance from the grid point in each geological continuous block to the fault line.

[0101] In one embodiment of the present disclosure, obtaining the weight of each grid point corresponding to each geological continuous block based on the minimum distance from the grid point to the fault line in each geological continuous block may include: when any grid point is within the first geological continuous block and not within the second geological continuous block, determining the weight of any grid point corresponding to the second geological continuous block to be a set value; obtaining first mapping data representing a mapping relationship between distance and weight; and obtaining the weight of any grid point corresponding to the first geological continuous block based on the minimum distance from any grid point to the fault line and the first mapping data.

[0102] In this embodiment, for each geological continuous block, the weight of each grid point corresponding to the geological continuous block may be determined.

[0103] Specifically, in the division In the case of a geological continuous block, where the preset area includes m grid points, it can be for the kth (k is less than or equal to For each of the i-th (i is a positive integer less than or equal to m) geologically continuous blocks, when the i-th grid point is in the k-th geologically continuous block, the weight of the i-th grid point corresponding to the k-th geologically continuous block is determined based on the minimum distance from the i-th grid point to the fault line and the first mapping data. The weight of the i-th grid point corresponding to geologically continuous blocks other than the k-th geologically continuous block is a set value. The set value can be pre-set based on the application scenario or specific requirements. For example, the set value can be 1 or 0.9.

[0104] The first mapping data may be a first mapping function, or a first comparison table, etc., which is not limited here.

[0105] For the first mapping function, the dependent variable of the first mapping function is the weight, and the independent variable is the distance. In this way, by substituting the minimum distance from the i-th grid point to the fault line into the first mapping function, the weight corresponding to the minimum distance from the i-th grid point to the fault line can be obtained, which is used as the weight of the i-th grid point corresponding to the k-th geological continuous block.

[0106] For the first lookup table, the weight corresponding to the minimum distance from the i-th grid point to the fault line can be found in the first lookup table and used as the weight for the i-th grid point corresponding to the k-th geological contiguous block. If the minimum distance from the i-th grid point to the fault line cannot be directly found in the first lookup table, two values ​​adjacent to the minimum distance from the i-th grid point to the fault line can be found. Based on these two values ​​and the weights corresponding to these two values, the weight corresponding to the minimum distance from the i-th grid point to the fault line can be obtained by interpolation and used as the weight for the i-th grid point corresponding to the k-th geological contiguous block.

[0107] In one embodiment, Figure 3 (c) in the figure can represent the grid points in the preset area. Figure 3 (b) Weight distribution diagram of a geological continuous block, Figure 3 (d) in the figure can represent the grid points in the preset area. Figure 3 (b) Weight distribution diagram for another geological contiguous block, where white represents a weight of 1.0 and black represents a weight of 0.0. It can be seen that the weight values ​​jump on both sides of the fault line and the extended fault line, while the weight transitions smoothly in areas without faults.

[0108] In one embodiment of the present disclosure, based on the fault line data, the weight of each grid point corresponding to each geological continuous block in the preset area is determined. On the basis of the aforementioned steps S2411 to S2412, step S2413 as shown below can also be included to normalize the weight of each grid point corresponding to each geological continuous block.

[0109] Step S2420 : obtaining a second three-dimensional geological model of the preset area according to the weight of each grid point corresponding to each geological continuous block and the first three-dimensional geological model corresponding to each geological continuous block.

[0110] In one embodiment of the present disclosure, obtaining a second three-dimensional geological model of a preset area according to the weight of each grid point corresponding to each geological continuous block and the first three-dimensional geological model corresponding to each geological continuous block may include the following steps S2421 to S2422:

[0111] Step S2421 : for each grid point, obtain the second geological attribute value of the corresponding grid point according to the first geological attribute value of the corresponding grid point in each first three-dimensional geological model and the weight of each geological continuous block corresponding to the corresponding grid point.

[0112] In an embodiment in which the weight of each grid point corresponding to each geological continuous block is normalized, for each grid point, a second geological attribute value of the corresponding grid point is obtained based on the first geological attribute value of the corresponding grid point in each first three-dimensional geological model and the weight of the corresponding grid point corresponding to each geological continuous block. This can be achieved by performing a weighted sum of the first geological attribute values ​​of each grid point corresponding to each geological continuous block for each grid point in a preset area to serve as the second geological attribute value of the grid point; and obtaining a second three-dimensional geological model based on the second attribute values ​​of all grid points in the preset area.

[0113] In an embodiment where the weight of each grid point corresponding to each geological continuous block is not normalized, for each grid point, the second geological attribute value of the corresponding grid point is obtained based on the first geological attribute value of the corresponding grid point in each first three-dimensional geological model and the weight of the corresponding grid point corresponding to each geological continuous block. This can be done by determining, for each grid point in a preset area, a weighted average of the first geological attribute values ​​of the grid point corresponding to each geological continuous block as the second geological attribute value of the grid point; and obtaining a second three-dimensional geological model based on the second attribute values ​​of all grid points in the preset area.

[0114] Step S2422: Obtain a second three-dimensional geological model based on the second geological attribute value of each grid point in the preset area.

[0115] In one embodiment, Figure 3 (e) in FIG. 5 may represent a second three-dimensional geological model. Figure 3 The color or grayscale of (e) may represent the second geological attribute value corresponding to the geographic location.

[0116] Step S2500: Obtain a three-dimensional geological modeling result of a preset area based on the second three-dimensional geological model.

[0117] Through this embodiment, the preset area is divided into at least two geologically continuous blocks based on fault line data used to describe at least one fault line in the preset area. A first three-dimensional geological model of the preset area is determined based on the first sampling point data of the first sampling point in each geologically continuous block. A second three-dimensional geological model of the preset area is obtained based on the first three-dimensional geological model corresponding to each geologically continuous block, and a three-dimensional geological modeling result of the preset area is obtained based on the second three-dimensional geological model. In this way, even if the preset area includes fault terrain, accurate three-dimensional modeling of the geological properties in the preset area can be performed.

[0118] In one embodiment of the present disclosure, the second three-dimensional geological model may be used as the three-dimensional geological modeling result of the preset area.

[0119] In another embodiment of the present disclosure, obtaining a 3D geological modeling result of a preset area according to the second 3D geological model may include steps S2510 to S2540 as follows:

[0120] Step S2510: acquiring second sampling point data representing a geological attribute value of a second sampling point according to the second three-dimensional geological model, wherein the second sampling point is located on an extension line of the fault line.

[0121] In this embodiment, the second geological attribute value of the second sampling point in the second three-dimensional geological model may be used as the second sampling point data of the second sampling point.

[0122] In one embodiment, Figure 3 The point in part (f) may represent the second sampling point.

[0123] In one embodiment, the number of second sampling points on each fault line extension line may be the same, for example, may all be the first number.

[0124] In another embodiment, sampling may be performed at set intervals to obtain at least one second sampling point on each extended line of the fault line.

[0125] Step S2520 : determining a third three-dimensional geological model of the preset area according to the first sampling point data and the second sampling point data of the first sampling point in each geological continuous block.

[0126] In this embodiment, the third three-dimensional geological model represents the third geological attribute value of each grid point in the preset area.

[0127] In one embodiment of the present disclosure, for each geological continuous block, a third three-dimensional geological model corresponding to the geological continuous block may be determined based on the geographical location of the grid points within the geological continuous block, the first sampling point data of the first sampling point within the geological continuous block, and the second sampling point data of all the second sampling points within the geological continuous block.

[0128] Specifically, in the division In the case of geological continuous blocks, it can be for the kth (k is less than or equal to (a positive integer) geological continuous blocks, and obtaining, based on the geographical locations of the grid points in the kth geological continuous block, the first sampling point data of the first sampling point in the kth geological continuous block, and the second sampling point data of all the second sampling points, a third geological attribute value corresponding to the kth geological continuous block for each grid point in the preset area is obtained.

[0129] In this embodiment, the difference algorithm used to determine the third three-dimensional geological model may be the same as the difference algorithm used to determine the first three-dimensional geological model.

[0130] Step S2530: Obtain a fourth three-dimensional geological model of the preset area based on the third three-dimensional geological model corresponding to each geological continuous block.

[0131] In one embodiment of the present disclosure, obtaining a fourth three-dimensional geological model of a preset area based on the third three-dimensional geological model corresponding to each geological continuous block may include: determining an average value of the third three-dimensional geological model corresponding to each geological continuous block to obtain the fourth three-dimensional geological model.

[0132] Specifically, for each grid point in the preset area, the average value of the third geological attribute value corresponding to each geological continuous block of the grid point can be determined as the fourth geological attribute value of the grid point; and the fourth three-dimensional geological model can be obtained based on the fourth attribute values ​​of all grid points in the preset area.

[0133] In another embodiment of the present disclosure, obtaining a fourth three-dimensional geological model of a preset area according to the third three-dimensional geological model corresponding to each geological continuous block includes the following steps S2531 to S2532:

[0134] Step S2531: Determine the weight of each grid point corresponding to each geological continuous block in a preset area based on the fault line data.

[0135] In one embodiment of the present disclosure, determining the weight of each grid point corresponding to each geological continuous block in a preset area based on fault line data may include the following steps S25311 to S25312:

[0136] Step S25311, based on the fault line data, determine the minimum distance from the grid point to the fault line in each geological continuous block.

[0137] In this embodiment, when the preset area includes at least two fault lines, the minimum distance from each grid point to each fault line may be determined for each grid point, and the smallest distance is used as the minimum distance from the grid point to the fault line.

[0138] Step S25312: Obtain the weight of each grid point corresponding to each geological continuous block based on the minimum distance from the grid point in each geological continuous block to the fault line.

[0139] In one embodiment of the present disclosure, obtaining the weight of each grid point corresponding to each geological continuous block based on the minimum distance from the grid point to the fault line in each geological continuous block may include: when any grid point is within the first geological continuous block and not within the second geological continuous block, determining the weight of any grid point corresponding to the second geological continuous block to be a set value; obtaining first mapping data representing a mapping relationship between distance and weight; and obtaining the weight of any grid point corresponding to the first geological continuous block based on the minimum distance from any grid point to the fault line and the first mapping data.

[0140] In this embodiment, for each geological continuous block, the weight of each grid point corresponding to the geological continuous block may be determined.

[0141] Specifically, in the division In the case of a geological continuous block, where the preset area includes m grid points, it can be for the kth (k is less than or equal to For each of the i-th (i is a positive integer less than or equal to m) geologically continuous blocks, when the i-th grid point is in the k-th geologically continuous block, the weight of the i-th grid point corresponding to the k-th geologically continuous block is determined based on the minimum distance from the i-th grid point to the fault line and the first mapping data. The weight of the i-th grid point corresponding to geologically continuous blocks other than the k-th geologically continuous block is a set value. The set value can be pre-set based on the application scenario or specific requirements. For example, the set value can be 1 or 0.9.

[0142] Step S2532 : obtaining a fourth three-dimensional geological model of the preset area according to the weight of each grid point corresponding to each geological continuous block and the third three-dimensional geological model corresponding to each geological continuous block.

[0143] In one embodiment of the present disclosure, obtaining a fourth three-dimensional geological model of a preset area based on the weight of each grid point corresponding to each geological continuous block and the third three-dimensional geological model corresponding to each geological continuous block may include the following steps S25321-S25322:

[0144] Step S25321: For each grid point, the fourth geological attribute value of the corresponding grid point is obtained according to the third geological attribute value of the corresponding grid point in each third three-dimensional geological model and the weight of each geological continuous block corresponding to the corresponding grid point.

[0145] In an embodiment in which the weight of each grid point corresponding to each geological continuous block is normalized, for each grid point, a fourth geological attribute value of the corresponding grid point is obtained based on the third geological attribute value of the corresponding grid point in each third three-dimensional geological model and the weight of the corresponding grid point corresponding to each geological continuous block. This can be achieved by performing a weighted summation of the third geological attribute values ​​of each grid point corresponding to each geological continuous block for each grid point in a preset area to serve as the fourth geological attribute value of the grid point; and obtaining a fourth three-dimensional geological model based on the second attribute values ​​of all grid points in the preset area.

[0146] In an embodiment where the weight of each grid point corresponding to each geological continuous block is not normalized, for each grid point, the fourth geological attribute value of the corresponding grid point is obtained based on the third geological attribute value of the corresponding grid point in each third three-dimensional geological model and the weight of the corresponding grid point corresponding to each geological continuous block. The fourth geological attribute value of the grid point can be obtained by determining the weighted average of the third geological attribute values ​​of the grid point corresponding to each geological continuous block for each grid point in the preset area; and the fourth three-dimensional geological model is obtained based on the second attribute values ​​of all grid points in the preset area.

[0147] Step S25322: Obtain a fourth three-dimensional geological model based on the fourth geological attribute value of each grid point in the preset area.

[0148] In one embodiment, Figure 3 (g) in the figure may represent a fourth three-dimensional geological model. Figure 3 The color or grayscale of (g) in FIG. 5 may represent the fourth geological attribute value of the corresponding geographic location.

[0149] Step S2540: Using the fourth three-dimensional geological model as the three-dimensional modeling result of the preset area.

[0150] Through this embodiment, the accuracy of the three-dimensional modeling result of the preset area can be improved.

[0151] Figure 4 A flowchart of a three-dimensional geological modeling method provided by another embodiment of the present disclosure is shown.

[0152] like Figure 4 As shown, the method may include the following steps:

[0153] Step S4001: Acquire first sampling point data representing a geological attribute value of a first sampling point in a preset area and fault line data for describing at least one fault line in the preset area.

[0154] Step S4002: Generate a fault line extension line from the endpoint of the fault line to the boundary of a preset area.

[0155] Step S4003: Divide the preset area into at least two geological continuous blocks according to the fault line and the fault line extension line.

[0156] Step S4004: determining a first three-dimensional geological model of a preset area according to the first sampling point data of the first sampling point in each geological continuous block.

[0157] Step S4005: Determine the weight of each grid point corresponding to each geological continuous block in the preset area based on the fault line data.

[0158] Step S4006 : obtaining a second three-dimensional geological model of the preset area according to the weight of each grid point corresponding to each geological continuous block and the first three-dimensional geological model corresponding to each geological continuous block.

[0159] Step S4007: acquiring second sampling point data representing a geological attribute value of a second sampling point according to the second three-dimensional geological model, wherein the second sampling point is located on an extension line of the fault line.

[0160] Step S4008 : determining a third three-dimensional geological model of the preset area according to the first sampling point data and the second sampling point data of the first sampling point in each geological continuous block.

[0161] Step S4009, obtaining a fourth three-dimensional geological model of the preset area according to the weight of each grid point corresponding to each geological continuous block and the third three-dimensional geological model corresponding to each geological continuous block;

[0162] Step S4010: Using the fourth three-dimensional geological model as the three-dimensional modeling result of the preset area.

[0163] <Device Example>

[0164] The present disclosure provides a three-dimensional geological modeling device 5000, such as Figure 5 As shown, the three-dimensional geological modeling device 5000 includes a data acquisition module 5100 , a geological division module 5200 , a first modeling module 5300 , a second modeling module 5400 and a result determination module 5500 .

[0165] The data acquisition module 5100 is used to acquire first sampling point data representing a geological attribute value of a first sampling point in a preset area and fault line data for describing at least one fault line in the preset area.

[0166] The geological division module 5200 is used to divide the preset area into at least two geological continuous blocks according to the fault line data.

[0167] The first modeling module 5300 is used to determine a first three-dimensional geological model of the preset area based on the first sampling point data of the first sampling point in each geological continuous block; the three-dimensional geological model represents the geological attribute values ​​of each grid point in the preset area.

[0168] The second modeling module 5400 is used to obtain a second three-dimensional geological model of the preset area according to the first three-dimensional geological model corresponding to each geological continuous block.

[0169] The result determination module 5500 is used to obtain the three-dimensional geological modeling result of the preset area according to the second three-dimensional geological model.

[0170] In one embodiment of the present disclosure, the second modeling module 5400 is used to:

[0171] Determining, based on the fault line data, a weight of each of the grid points in the preset area corresponding to each geological continuous block;

[0172] A second three-dimensional geological model of the preset area is obtained according to the weight of each grid point corresponding to each geological continuous block and the first three-dimensional geological model corresponding to each geological continuous block.

[0173] In one embodiment of the present disclosure, determining the weight of each grid point corresponding to each geological continuous block in the preset area based on the fault line data includes:

[0174] Determining the minimum distance between a grid point in each geological continuous block and the fault line based on the fault line data;

[0175] Obtaining a weight of each grid point corresponding to each geological continuous block according to a minimum distance between the grid point in each geological continuous block and the fault line;

[0176] The weight of each grid point corresponding to each geological continuous block is normalized.

[0177] In one embodiment of the present disclosure, obtaining the weight of each grid point corresponding to each geological continuous block according to the minimum distance between the grid point in each geological continuous block and the fault line includes:

[0178] When any grid point is within the first geological continuous block and not within the second geological continuous block, determining a weight of the any grid point corresponding to the second geological continuous block to be a set value;

[0179] Acquire first mapping data representing a mapping relationship between distance and weight;

[0180] A weight of the any grid point corresponding to the first geological continuous block is obtained according to the minimum distance from the any grid point to the fault line and the first mapping data.

[0181] In one embodiment of the present disclosure, obtaining the second three-dimensional geological model of the preset area according to the weight of each grid point corresponding to each geological continuous block and the first three-dimensional geological model corresponding to each geological continuous block includes:

[0182] For each of the grid points, a second geological attribute value of the corresponding grid point is obtained according to the first geological attribute value of the corresponding grid point in each first three-dimensional geological model and the weight of each geological continuous block corresponding to the corresponding grid point;

[0183] The second three-dimensional geological model is obtained according to the second geological attribute value of each grid point in the preset area.

[0184] In one embodiment of the present disclosure, the geological partitioning module 5200 is used to:

[0185] generating a fault line extension line from an endpoint of the fault line to a boundary of the preset area;

[0186] The preset area is divided into at least two geologically continuous blocks according to the fault line and the extended line of the fault line.

[0187] In one embodiment of the present disclosure, the result determination module 5500 is used to:

[0188] acquiring second sampling point data representing a geological attribute value of a second sampling point according to the second three-dimensional geological model, wherein the second sampling point is located on an extension line of the fault line;

[0189] Determining a third three-dimensional geological model of the preset area based on the first sampling point data and the second sampling point data of the first sampling point in each geological continuous block;

[0190] Obtaining a fourth three-dimensional geological model of the preset area according to the third three-dimensional geological model corresponding to each geological continuous block;

[0191] The fourth three-dimensional geological model is used as the three-dimensional geological modeling result of the preset area.

[0192] In one embodiment of the present disclosure, the first modeling module 5300 is used to determine the first three-dimensional geological model of the preset area based on the first sampling point data of the first sampling point in each geological continuous block based on a difference algorithm corresponding to the attributes of the geological attribute value.

[0193] <Electronic Equipment Example>

[0194] In this embodiment, an electronic device 6000 is also provided. The electronic device 6000 may be Figure 1 The electronic device 1000 is shown.

[0195] In one aspect, the electronic device 6000 may include the aforementioned 3D geological modeling device 6000 , for implementing the 3D geological modeling method of any embodiment of the present disclosure.

[0196] On the other hand, Figure 6 As shown, the electronic device 6000 may further include a processor 6100 and a memory 6200, wherein the memory 6200 is used to store executable instructions; the processor 6100 is used to control the operation of the electronic device 6000 under the control of the instructions to execute the three-dimensional geological modeling method according to any embodiment of the present disclosure.

[0197] In this embodiment, the electronic device 6000 may be a mobile phone, tablet computer, PDA, desktop computer, laptop computer, etc. For example, the electronic device 6000 may be an electronic product with a 3D geological modeling function.

[0198] <Computer-readable storage medium embodiment>

[0199] In this embodiment, a computer-readable storage medium is further provided, on which a computer program is stored. When the computer program is executed by a processor, the three-dimensional geological modeling method of any embodiment of the present disclosure is implemented.

[0200] The present disclosure may be a system, method and / or computer program product. The computer program product may include a computer-readable storage medium carrying computer-readable program instructions for causing a processor to implement various aspects of the present disclosure.

[0201] A computer-readable storage medium can be a tangible device that can hold and store instructions for use by an instruction execution device. A computer-readable storage medium can be, for example, but not limited to, an electrical storage device, a magnetic storage device, an optical storage device, an electromagnetic storage device, a semiconductor storage device, or any suitable combination thereof. More specific examples (a non-exhaustive list) of computer-readable storage media include: a portable computer disk, a hard disk, random access memory (RAM), read-only memory (ROM), erasable programmable read-only memory (EPROM or flash memory), static random access memory (SRAM), a portable compact disc read-only memory (CD-ROM), a digital versatile disk (DVD), a memory stick, a floppy disk, a mechanical encoding device, such as a punch card or raised-in-groove structure on which instructions are stored, and any suitable combination thereof. As used herein, a computer-readable storage medium is not to be construed as a transient signal per se, such as a radio wave or other freely propagating electromagnetic wave, an electromagnetic wave propagating through a waveguide or other transmission medium (e.g., a light pulse through a fiber optic cable), or an electrical signal transmitted through an electrical wire.

[0202] The computer-readable program instructions described herein can be downloaded from a computer-readable storage medium to each computing / processing device, or downloaded to an external computer or external storage device via a network, such as the Internet, a local area network, a wide area network, and / or a wireless network. The network can include copper transmission cables, fiber optic transmission, wireless transmission, routers, firewalls, switches, gateway computers, and / or edge servers. The network adapter card or network interface in each computing / processing device receives the computer-readable program instructions from the network and forwards the computer-readable program instructions to be stored in the computer-readable storage medium in each computing / processing device.

[0203] The computer program instructions for performing the operations of the present disclosure may be assembly instructions, instruction set architecture (ISA) instructions, machine instructions, machine-dependent instructions, microcode, firmware instructions, state setting data, or source code or object code written in any combination of one or more programming languages, including object-oriented programming languages ​​such as Smalltalk, C++, and conventional procedural programming languages ​​such as "C" or similar programming languages. The computer-readable program instructions may be executed entirely on the user's computer, partially on the user's computer, as a stand-alone software package, partially on the user's computer and partially on a remote computer, or entirely on a remote computer or server. In the case of a remote computer, the remote computer may be connected to the user's computer via any type of network, including a local area network (LAN) or a wide area network (WAN), or may be connected to an external computer (e.g., via the Internet using an Internet service provider). In some embodiments, the state information of the computer-readable program instructions is used to personalize an electronic circuit, such as a programmable logic circuit, a field programmable gate array (FPGA), or a programmable logic array (PLA), so that the electronic circuit can execute the computer-readable program instructions, thereby implementing various aspects of the present disclosure.

[0204] Various aspects of the present disclosure are described herein with reference to flowcharts and / or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the present disclosure. It should be understood that each block of the flowcharts and / or block diagrams, and combinations of blocks in the flowcharts and / or block diagrams, can be implemented by computer-readable program instructions.

[0205] These computer-readable program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, or other programmable data processing device, thereby producing a machine, so that when these instructions are executed by the processor of the computer or other programmable data processing device, a device is generated that implements the functions / actions specified in one or more blocks in the flowchart and / or block diagram. These computer-readable program instructions can also be stored in a computer-readable storage medium, where these instructions cause the computer, programmable data processing device, and / or other device to operate in a specific manner. Thus, the computer-readable medium storing the instructions comprises an article of manufacture that includes instructions for implementing various aspects of the functions / actions specified in one or more blocks in the flowchart and / or block diagram.

[0206] Computer-readable program instructions may also be loaded onto a computer, other programmable data processing apparatus, or other device so that a series of operational steps are performed on the computer, other programmable data processing apparatus, or other device to produce a computer-implemented process, thereby causing the instructions executed on the computer, other programmable data processing apparatus, or other device to implement the functions / actions specified in one or more blocks in the flowchart and / or block diagram.

[0207] The flowcharts and block diagrams in the accompanying drawings show the possible implementation architecture, functions and operations of the systems, methods and computer program products according to multiple embodiments of the present disclosure. In this regard, each box in the flowchart or block diagram can represent a module, program segment or part of an instruction, and the module, program segment or part of the instruction contains one or more executable instructions for implementing the specified logical function. In some alternative implementations, the functions marked in the box can also occur in an order different from that marked in the accompanying drawings. For example, two consecutive boxes can actually be executed substantially in parallel, and they can sometimes be executed in the opposite order, depending on the functions involved. It should also be noted that each box in the block diagram and / or flowchart, and the combination of boxes in the block diagram and / or flowchart, can be implemented by a dedicated hardware-based system that performs the specified function or action, or can be implemented by a combination of dedicated hardware and computer instructions. It is well known to those skilled in the art that implementation by hardware, implementation by software, and implementation by a combination of software and hardware are all equivalent.

[0208] The embodiments of the present disclosure have been described above. The above description is exemplary, not exhaustive, and is not limited to the disclosed embodiments. Many modifications and variations will be apparent to those skilled in the art without departing from the scope and spirit of the described embodiments. The terms used herein are selected to best explain the principles of the embodiments, their practical applications, or technical improvements in the marketplace, or to enable other persons skilled in the art to understand the embodiments disclosed herein. The scope of the present disclosure is defined by the appended claims.

Claims

1. A three-dimensional geological modeling method, characterized in that: include: Acquire first sampling point data representing a geological attribute value of a first sampling point in a preset area and fault line data for describing at least one fault line in the preset area; Dividing the preset area into at least two geologically continuous blocks according to the fault line data; Determining a first three-dimensional geological model of the preset area according to first sampling point data of the first sampling point in each geological continuous block; the three-dimensional geological model represents the geological attribute values ​​of each grid point in the preset area; Obtaining a second three-dimensional geological model of the preset area according to the first three-dimensional geological model corresponding to each geological continuous block; Obtaining a three-dimensional geological modeling result of the preset area according to the second three-dimensional geological model; The step of obtaining the second three-dimensional geological model of the preset area according to the first three-dimensional geological model corresponding to each geological continuous block includes: Determining, according to the fault line data, the weight of each of the grid points in the preset area corresponding to each geological continuous block; Obtaining a second three-dimensional geological model of the preset area according to the weight of each of the grid points corresponding to each of the geological continuous blocks and the first three-dimensional geological model corresponding to each of the geological continuous blocks; The second three-dimensional geological model of the preset area is obtained according to the weight of each of the grid points corresponding to each geological continuous block and the first three-dimensional geological model corresponding to each geological continuous block, including: For each of the grid points, a second geological attribute value of the corresponding grid point is obtained according to the first geological attribute value of the corresponding grid point in each first three-dimensional geological model and the weight of each geological continuous block corresponding to the corresponding grid point; The second three-dimensional geological model is obtained according to the second geological attribute value of each grid point in the preset area.

2. The method according to claim 1, characterized in that Determining the weight of each of the grid points in the preset area corresponding to each geological continuous block according to the fault line data includes: Determining the minimum distance from a grid point in each geological continuous block to the fault line according to the fault line data; According to the minimum distance between the grid points in each geological continuous block and the fault line, the weight of each grid point corresponding to each geological continuous block is obtained; The weight of each of the grid points corresponding to each geological continuous block is normalized.

3. The method according to claim 1, characterized in that The step of obtaining the weight of each grid point corresponding to each geological continuous block according to the minimum distance from the grid point in each geological continuous block to the fault line includes: When any grid point is within the first geological continuous block and not within the second geological continuous block, determining the weight of the any grid point corresponding to the second geological continuous block to be a set value; Acquire first mapping data representing a mapping relationship between distance and weight; The weight of any grid point corresponding to the first geological continuous block is obtained according to the minimum distance from any grid point to the fault line and the first mapping data.

4. The method according to claim 1, characterized in that: The step of dividing the preset area into at least two geologically continuous blocks according to the fault line data comprises: generating a fault line extension line from an endpoint of the fault line to a boundary of the preset area; According to the fault line and the extended line of the fault line, the preset area is divided into at least two geological continuous blocks.

5. The method according to claim 4, characterized in that The step of obtaining the three-dimensional geological modeling result of the preset area according to the second three-dimensional geological model includes: According to the second three-dimensional geological model, obtaining second sampling point data representing a geological attribute value of a second sampling point, wherein the second sampling point is located on an extension line of the fault line; Determine a third three-dimensional geological model of the preset area according to the first sampling point data and the second sampling point data of the first sampling point in each geological continuous block respectively; Obtaining a fourth three-dimensional geological model of the preset area according to the third three-dimensional geological model corresponding to each geological continuous block; The fourth three-dimensional geological model is used as the three-dimensional geological modeling result of the preset area.

6. The method according to claim 1, characterized in that The step of determining the first three-dimensional geological model of the preset area according to the first sampling point data of the first sampling point in each geological continuous block comprises: Based on a difference algorithm corresponding to the attributes of the geological attribute values, a first three-dimensional geological model of the preset area is determined according to first sampling point data of the first sampling point in each geological continuous block.

7. A three-dimensional geological modeling device, characterized in that: include: A data acquisition module, used to acquire first sampling point data representing a geological attribute value of a first sampling point in a preset area and fault line data for describing at least one fault line in the preset area; A geological division module, used for dividing the preset area into at least two geological continuous blocks according to the fault line data; A first modeling module is used to determine a first three-dimensional geological model of the preset area according to first sampling point data of a first sampling point in each geological continuous block; the three-dimensional geological model represents the geological attribute values ​​of each grid point in the preset area; A second modeling module, used for obtaining a second three-dimensional geological model of the preset area according to the first three-dimensional geological model corresponding to each geological continuous block; A result determination module, used for obtaining a three-dimensional geological modeling result of the preset area according to the second three-dimensional geological model; The step of obtaining the second three-dimensional geological model of the preset area according to the first three-dimensional geological model corresponding to each geological continuous block includes: Determining, according to the fault line data, the weight of each of the grid points in the preset area corresponding to each geological continuous block; Obtaining a second three-dimensional geological model of the preset area according to the weight of each of the grid points corresponding to each of the geological continuous blocks and the first three-dimensional geological model corresponding to each of the geological continuous blocks; The second three-dimensional geological model of the preset area is obtained according to the weight of each of the grid points corresponding to each geological continuous block and the first three-dimensional geological model corresponding to each geological continuous block, including: For each of the grid points, a second geological attribute value of the corresponding grid point is obtained according to the first geological attribute value of the corresponding grid point in each first three-dimensional geological model and the weight of each geological continuous block corresponding to the corresponding grid point; The second three-dimensional geological model is obtained according to the second geological attribute value of each grid point in the preset area.

8. An electronic device, characterized in that: include: The device according to claim 7; or, A processor and a memory, wherein the memory is used to store a computer program, wherein the computer program is used to control the processor to execute the method according to any one of claims 1 to 6.

Images