Method, device, computer equipment and storage medium for constructing geological layer map

Abstract

The present application discloses a method, device, computer equipment and storage medium for constructing a geological layer map, which belongs to the field of geophysical exploration technology. The method obtains multiple profile data of the target stratum, and determines multiple parallel profiles based on the multiple profile data; based on the predetermined parallel profile parameters, multiple interpolation positions are determined, and multiple first layer segment lines of the target stratum are determined; based on the three-dimensional coordinates of the first control point corresponding to each first layer segment line, the coordinates of multiple interpolation points of each interpolation position are determined; based on the coordinates of the multiple interpolation points of each interpolation position, multiple second layer segment lines are inserted, and a geological layer map is constructed based on the multiple first layer segment lines and multiple second layer segment lines of the target stratum. This method can construct a smooth and natural geological layer map when there is less existing profile data, thereby improving the effect of the constructed geological layer map.

Description

Technical Field

[0001] The present application relates to the field of geophysical exploration technology, and more particularly to a method, apparatus, computer equipment, and storage medium for constructing a geological layer map. Background Art

[0002] The seismic data acquisition phase primarily involves collecting profile data, which can be used to construct three-dimensional geological layer maps. Smooth geological layer maps are a prerequisite for forward modeling and illumination of three-dimensional geological models. Therefore, the quality of geological layer maps has a significant impact on the effectiveness of these two aspects.

[0003] In the related technology, a continuous layer triangulation network is first constructed through the collected profile data, and then the vertex coordinates of the layer triangulation network are determined using a surface fitting algorithm. The vertex coordinates of the layer triangulation network are then connected to obtain a geological layer map.

[0004] However, during the seismic data acquisition stage, there are usually only a few profile data. Due to the small amount of profile data, the layer triangulation constructed in the relevant technology is not smooth and natural enough, resulting in the constructed geological layer map not being smooth and natural enough, resulting in poor effect of the constructed geological layer map. Summary of the Invention

[0005] The present invention provides a method, apparatus, computer device, and storage medium for constructing a geological layer map, which can improve the effect of the constructed geological layer map. The specific technical solution is as follows:

[0006] In one aspect, an embodiment of the present application provides a method for constructing a geological layer map, the method comprising:

[0007] Acquire multiple cross-sectional data of a target stratum, and determine multiple parallel cross-sectional views based on the multiple cross-sectional data, wherein the cross-sectional data includes three-dimensional coordinates of a first control point, where the first control point is any position point in the target stratum;

[0008] Based on predetermined parallel section parameters, a plurality of interpolation positions are determined, wherein the interpolation positions are used to insert layer segment lines parallel to the plurality of parallel sections between the plurality of parallel sections;

[0009] Determine a plurality of first layer segment lines of the target stratum, wherein the first layer segment lines are formed by connecting a plurality of first control points;

[0010] determining coordinates of a plurality of interpolation points at each interpolation position based on the three-dimensional coordinates of the first control point corresponding to each first layer segment line;

[0011] inserting a plurality of second-layer segment lines based on the coordinates of the plurality of interpolation points at each interpolation position, wherein the second-layer segment lines are formed by connecting the plurality of interpolation points;

[0012] A geological layer map is constructed based on the multiple first layer segment lines and the multiple second layer segment lines of the target stratum.

[0013] In a possible implementation, determining the coordinates of multiple interpolation points at each interpolation position based on the three-dimensional coordinates of the first control point corresponding to each first layer segment line includes:

[0014] Determining a target layer segment line having the largest number of first control points based on the number of first control points corresponding to each first layer segment line;

[0015] Determining a plurality of position ratios based on the target layer segment line, wherein the position ratio is a ratio of a distance between a first first control point and each of the other first control points in the target layer segment line to a segment length of the target layer segment line, wherein an order of the first control points in the target layer segment line is determined according to a coordinate size of each first control point in a direction of a first coordinate axis, wherein the first coordinate axis is parallel to the first layer segment line;

[0016] Determining a plurality of second control points based on the plurality of position ratios and the segment lengths of other layer segment lines in the plurality of first layer segment lines, where the second control points are control points on the other layer segment lines corresponding to the plurality of position ratios;

[0017] determining the plurality of second control points as the plurality of first control points on the other layer segment lines;

[0018] The coordinates of the multiple interpolation points of each interpolation position are determined based on the re-determined three-dimensional coordinates of the first control point on the other layer segment line and the three-dimensional coordinates of the multiple first control points on the target layer segment line.

[0019] In another possible implementation, determining the coordinates of multiple interpolation points at each interpolation position based on the re-determined three-dimensional coordinates of the first control point on the other layer segment line and the three-dimensional coordinates of multiple first control points on the target layer segment line includes:

[0020] For each first control point of each position ratio, determine a relationship curve consisting of each adjacent preset number of first control points of the position ratio on the other layer segment lines and the target layer segment line, wherein the relationship curve is a curve with a first distance ratio as an independent variable, the first distance ratio is a ratio of a first distance to a second distance, the first distance is a distance between a third control point and a first first control point among a preset number of first control points, the third control point is a control point at any position on the relationship curve, and the second distance is a sum of distances between every two adjacent first control points among the preset number of first control points;

[0021] Determining a plurality of curve coefficients of the relationship curve based on the three-dimensional coordinates of each first control point constituting the relationship curve, wherein the curve coefficients include a constant term coefficient, a linear term coefficient, and a quadratic term coefficient, and each curve coefficient is a three-dimensional coordinate point, the constant term coefficient is a coefficient without an independent variable, the linear term coefficient is a coefficient with an independent variable exponent of one, and the quadratic term coefficient is a coefficient with an independent variable exponent of two;

[0022] determining a three-dimensional expression of the relationship curve based on the plurality of curve coefficients;

[0023] determining a sub-coefficient of each curve coefficient in the direction of a second coordinate axis, where the second coordinate axis is perpendicular to the first coordinate axis and is in the same horizontal plane;

[0024] For each interpolation position, based on the sub-coefficients of each curve coefficient, the three-dimensional expression of the relationship curve and the coordinates of the interpolation position on the second coordinate axis, the coordinates of the interpolation point corresponding to the interpolation position at the position scale are determined.

[0025] In another possible implementation, if the relationship curve where the interpolation position is located is a first relationship curve without overlap, determining the coordinates of the interpolation point corresponding to the interpolation position at the position scale based on the sub-coefficients of each curve coefficient, the three-dimensional expression of the relationship curve, and the coordinates of the interpolation position on the second coordinate axis includes:

[0026] determining a one-dimensional expression of the first relationship curve based on the sub-coefficients of each curve coefficient;

[0027] Correcting the constant term sub-coefficient in the one-dimensional expression based on the constant term sub-coefficient in the one-dimensional expression and the coordinate of the interpolation position in the direction of the second coordinate axis to obtain a corrected constant term sub-coefficient;

[0028] Determine a one-dimensional expression composed of the modified constant term sub-coefficient, the linear term sub-coefficient, and the quadratic term sub-coefficient in the one-dimensional expression;

[0029] determining a second distance ratio corresponding to when the value of the reconstructed one-dimensional expression is zero;

[0030] When determining that the independent variable of the first relationship curve is the second distance ratio, the value of the three-dimensional expression of the first relationship curve is used to obtain the coordinates of the interpolation point corresponding to the interpolation position at the position ratio.

[0031] In another possible implementation, if the relationship curve where the interpolation position is located is a first relationship curve obtained by overlapping a second relationship curve and a third relationship curve, determining the coordinates of the interpolation point corresponding to the interpolation position at the position scale based on the sub-coefficients of each curve coefficient, the three-dimensional expression of the relationship curve, and the coordinates of the interpolation position on the second coordinate axis includes:

[0032] determining a plurality of second curve coefficients of the second relationship curve and a plurality of third curve coefficients of the third relationship curve;

[0033] Determine a third distance ratio and a fourth distance ratio, wherein the third distance ratio is a ratio of a third distance to a fourth distance, wherein the third distance is the distance between a first first control point and a second first control point in the second relationship curve, and the fourth distance is the sum of the distances between every two adjacent first control points in the second relationship curve. The fourth distance ratio is a ratio of a fifth distance to a sixth distance, wherein the fifth distance is the distance between a first first control point and a second first control point in the third relationship curve, and the sixth distance is the sum of the distances between every two adjacent first control points in the third relationship curve.

[0034] determining, based on the plurality of second curve coefficients, the plurality of third curve coefficients, the third distance ratio, and the fourth distance ratio, a three-dimensional expression of a first relationship curve with a fifth distance ratio as an independent variable, the fifth distance ratio being a ratio of a seventh distance to an eighth distance, the seventh distance being a distance between an interpolation point at the interpolation position and a first control point at a beginning of the first relationship curve, and the eighth distance being a distance between the first control point at the beginning and the first control point at an end of the first relationship curve;

[0035] Determine a one-dimensional expression of the three-dimensional expression of the first relationship curve in the direction of the second coordinate axis;

[0036] Correcting the constant term sub-coefficient in the one-dimensional expression based on the constant term sub-coefficient in the one-dimensional expression, the coordinate of the interpolation position in the direction of the second coordinate axis, and the third distance ratio to obtain a corrected constant term sub-coefficient;

[0037] The coordinates of the interpolation point of the interpolation position at the position scale are determined based on the one-dimensional expression recomposed from the modified constant term sub-coefficient, the linear term sub-coefficient and the quadratic term sub-coefficient in the one-dimensional expression.

[0038] In another possible implementation, determining a plurality of interpolation positions based on predetermined parallel profile parameters includes:

[0039] Between two adjacent parallel sections, an interpolation position is determined for each distance corresponding to the parallel section parameter, thereby obtaining a plurality of interpolation positions.

[0040] In another possible implementation, inserting multiple second-layer segment lines based on the coordinates of the multiple interpolation points at each interpolation position includes:

[0041] For each interpolation position, multiple interpolation points corresponding to the interpolation position are sequentially connected to obtain a second-layer segment line corresponding to the interpolation position.

[0042] On the other hand, an embodiment of the present application provides a device for constructing a geological layer map, the device comprising:

[0043] an acquisition module, configured to acquire a plurality of cross-sectional data of a target stratum, and determine a plurality of parallel cross-sectional views based on the plurality of cross-sectional data, wherein the cross-sectional data includes a three-dimensional coordinate of a first control point, and the first control point is any position point in the target stratum;

[0044] A first determining module is configured to determine a plurality of interpolation positions based on predetermined parallel section parameters, wherein the interpolation positions are used to insert layer segment lines parallel to the plurality of parallel sections between the plurality of parallel sections;

[0045] A second determining module is configured to determine a plurality of first layer segment lines of the target stratum, wherein the first layer segment lines are formed by connecting a plurality of first control points;

[0046] a third determining module, configured to determine coordinates of a plurality of interpolation points at each interpolation position based on the three-dimensional coordinates of the first control point corresponding to each first layer segment line;

[0047] an interpolation module, configured to interpolate a plurality of second-layer segment lines based on the coordinates of the plurality of interpolation points at each interpolation position, wherein the second-layer segment lines are formed by connecting the plurality of interpolation points;

[0048] A construction module is used to construct a geological layer map based on multiple first layer segment lines and multiple second layer segment lines of the target stratum.

[0049] In one possible implementation, the third determination module is used to determine a target layer segment line with the largest number of first control points based on the number of first control points corresponding to each first layer segment line; determine multiple position ratios based on the target layer segment line, where the position ratio is the ratio of the distance between the first first control point and each other first control point in the target layer segment line to the line segment length of the target layer segment line, and the order of the first control points in the target layer segment line is determined according to the coordinate size of each first control point in the direction of the first coordinate axis, and the first coordinate axis is parallel to the first layer segment line; determine multiple second control points based on the multiple position ratios and the line segment lengths of other layer segment lines in the multiple first layer segment lines, where the second control points are control points on the other layer segment lines corresponding to the multiple position ratios; determine the multiple second control points as multiple first control points on the other layer segment lines; and determine the coordinates of multiple interpolation points at each interpolation position based on the re-determined three-dimensional coordinates of the first control points on the other layer segment lines and the three-dimensional coordinates of the multiple first control points on the target layer segment line.

[0050] In another possible implementation, the third determination module is used to determine, for each first control point of the position ratio, a relationship curve consisting of each adjacent preset number of first control points of the position ratio on the other layer segment lines and the target layer segment line, the relationship curve is a curve with a first distance ratio as an independent variable, the first distance ratio is the ratio of the first distance to the second distance, the first distance is the distance between the third control point and the first first control point among the preset number of first control points, the third control point is a control point at any position on the relationship curve, and the second distance is the sum of the distances between each adjacent two first control points among the preset number of first control points; based on the three-dimensional coordinates of each first control point constituting the relationship curve, the relationship curve is determined. multiple curve coefficients, the curve coefficients include a constant term coefficient, a linear term coefficient and a quadratic term coefficient, and each curve coefficient is a three-dimensional coordinate point, the constant term coefficient is a coefficient without an independent variable, the linear term coefficient is a coefficient with an independent variable exponent of one, and the quadratic term coefficient is a coefficient with an independent variable exponent of two; based on the multiple curve coefficients, determine the three-dimensional expression of the relationship curve; determine the sub-coefficient of each curve coefficient in the direction of the second coordinate axis, the second coordinate axis is perpendicular to the first coordinate axis and is in the same horizontal plane; for each interpolation position, based on the sub-coefficient of each curve coefficient, the three-dimensional expression of the relationship curve and the coordinates of the interpolation position on the second coordinate axis, determine the coordinates of the interpolation point corresponding to the interpolation position at the position scale.

[0051] In another possible implementation, if the relationship curve where the interpolation position is located is a first relationship curve without overlap, the third determination module is used to determine a one-dimensional expression of the first relationship curve based on the sub-coefficients of each curve coefficient; based on the constant term sub-coefficient in the one-dimensional expression and the coordinates of the interpolation position in the direction of the second coordinate axis, correct the constant term sub-coefficient in the one-dimensional expression to obtain a corrected constant term sub-coefficient; determine a one-dimensional expression recomposed of the corrected constant term sub-coefficient, the linear term sub-coefficient, and the quadratic term sub-coefficient in the one-dimensional expression; determine the second distance ratio corresponding to when the value of the recomposed one-dimensional expression is zero; determine the value of the three-dimensional expression of the first relationship curve when the independent variable of the first relationship curve is the second distance ratio, and obtain the coordinates of the interpolation point corresponding to the interpolation position at the position ratio.

[0052] In another possible implementation, if the relationship curve where the interpolation position is located is a first relationship curve obtained by overlapping the second relationship curve and the third relationship curve, the third determination module is used to determine multiple second curve coefficients of the second relationship curve and multiple third curve coefficients of the third relationship curve; determine a third distance ratio and a fourth distance ratio, the third distance ratio is the ratio of the third distance to the fourth distance, the third distance is the distance between the first first control point and the second first control point in the second relationship curve, the fourth distance is the sum of the distances between each two adjacent first control points in the second relationship curve, the fourth distance ratio is the ratio of the fifth distance to the sixth distance, the fifth distance is the distance between the first first control point and the second first control point in the third relationship curve, and the sixth distance is the sum of the distances between each two adjacent first control points in the third relationship curve; based on the multiple second curve coefficients, the multiple third curve coefficients, the The third distance ratio and the fourth distance ratio are used to determine a three-dimensional expression of the first relationship curve with the fifth distance ratio as the independent variable, where the fifth distance ratio is the ratio of the seventh distance to the eighth distance, the seventh distance being the distance between the interpolation point at the interpolation position and the first control point at the beginning of the first relationship curve, and the eighth distance being the distance between the first control point at the beginning and the first control point at the end of the first relationship curve; a one-dimensional expression of the three-dimensional expression of the first relationship curve in the direction of the second coordinate axis is determined; based on the constant term sub-coefficient in the one-dimensional expression, the coordinates of the interpolation position in the direction of the second coordinate axis, and the third distance ratio, the constant term sub-coefficient in the one-dimensional expression is corrected to obtain a corrected constant term sub-coefficient; and based on the one-dimensional expression recomposed from the corrected constant term sub-coefficient, the linear term sub-coefficient, and the quadratic term sub-coefficient in the one-dimensional expression, the coordinates of the interpolation point at the interpolation position at the position ratio are determined.

[0053] In another possible implementation, the first determining module is configured to determine an interpolation position between two adjacent parallel sections at each interval corresponding to the parallel section parameter, to obtain a plurality of interpolation positions.

[0054] In another possible implementation, the interpolation module is configured to sequentially connect, for each interpolation position, multiple interpolation points corresponding to the interpolation position to obtain a second-layer segment line corresponding to the interpolation position.

[0055] On the other hand, an embodiment of the present application provides a computer device, which includes a processor and a memory, wherein the memory stores at least one program code, and the at least one program code is loaded and executed by the processor to implement the method for constructing a geological layer map described in the embodiment of the present application.

[0056] On the other hand, an embodiment of the present application provides a computer-readable storage medium, in which at least one program code is stored. The at least one program code is loaded and executed by a processor to implement the method for constructing a geological layer map described in the embodiment of the present application.

[0057] On the other hand, an embodiment of the present application provides a computer program product, which stores at least one program code, and the at least one program code is loaded and executed by a processor to implement the method for constructing a geological layer map described in the embodiment of the present application.

[0058] The beneficial effects of the technical solution provided by the embodiments of the present application are:

[0059] An embodiment of the present application provides a method for constructing a geological layer map. The method interpolates multiple layer lines through existing profile data to obtain interpolated profile data, and then constructs a geological layer map based on the existing profile data and the interpolated profile data. Since the geological layer map is constructed with a large amount of profile data, the method can construct a smooth and natural geological layer map when there is less existing profile data, thereby improving the effect of the constructed geological layer map. BRIEF DESCRIPTION OF THE DRAWINGS

[0060] Figure 1 This is a flow chart of a method for constructing a geological layer map provided in an embodiment of the present application;

[0061] Figure 2 This is a schematic diagram of a plurality of parallel cross sections provided in an embodiment of the present application, all of which are parallel to the Y axis;

[0062] Figure 3 The embodiment of the present application provides 7 relationship curves composed of 3 adjacent first control points at the same position ratio;

[0063] Figure 4 This embodiment of the present application provides a Figure 2 A schematic diagram of inserting 14 second-layer segment lines based on the 4 first-layer segment lines shown;

[0064] Figure 5 This embodiment of the present application provides a Figure 4 A schematic diagram of a geological layer map constructed based on multiple first layer segment lines and multiple second layer segment lines shown;

[0065] Figure 6 This is a schematic structural diagram of a device for constructing a geological layer map provided in an embodiment of the present application;

[0066] Figure 7 This is a structural block diagram of a computer device provided in an embodiment of the present application. DETAILED DESCRIPTION

[0067] In order to make the technical solutions and advantages of the present application clearer, the implementation methods of the present application are described in further detail below.

[0068] The present application embodiment provides a method for constructing a geological layer map, which is executed by a computer device, see Figure 1 , the method comprising:

[0069] Step 101: A computer device obtains a plurality of cross-sectional data of a target formation, and determines a plurality of parallel cross-sectional data based on the plurality of cross-sectional data.

[0070] The target stratum may include one stratum or multiple strata. If the target stratum includes multiple strata, the computer device acquires multiple cross-sectional data for the multiple strata. Based on the multiple cross-sectional data for the multiple strata, multiple parallel cross-sectional planes are determined. The cross-sectional data includes the three-dimensional coordinates of a first control point, which is any location point in the target stratum.

[0071] The computer device determines the position of each first control point in the three-dimensional coordinate system based on the three-dimensional coordinates of each first control point, and determines multiple parallel cross-sections based on the position of each first control point in the three-dimensional coordinate system. The multiple parallel cross-sections are parallel to a first coordinate axis in the three-dimensional coordinate system, where the first coordinate axis is either the X-axis or the Y-axis. In the embodiments of the present application, only the Y-axis is used as an example for description.

[0072] After the computer device determines multiple parallel sections, it can sort the parallel sections from small to large according to the coordinates of each parallel section in the X-axis direction to obtain sorted parallel sections. Figure 2 , Figure 2 There are 4 parallel sections in it, and these 4 parallel sections are arranged from small to large according to their coordinates in the X-axis direction.

[0073] It should be noted that before step 101, the computer device can pre-construct a three-dimensional geological model, which is composed of a three-dimensional coordinate system. For example, the scope of the three-dimensional geological model is: 8000m×6000m×5000m. The computer device can also obtain user-defined parallel section parameters, which are parameters on the X-axis or Y-axis. If the parallel section is parallel to the X-axis, then the parallel section parameter is a parameter on the Y-axis, for example, the parallel section parameter Dy is 400m. If the parallel section is parallel to the Y-axis, then the parallel section parameter is a parameter on the X-axis, for example, the parallel section parameter Dx is 400m. The computer device can also determine multiple parallel section positions parallel to the X-axis or Y-axis within the scope of the three-dimensional geological model, and subsequently determine the parallel section position corresponding to the section data after obtaining the section data.

[0074] Step 102: The computer device determines a plurality of interpolation positions based on predetermined parallel profile parameters.

[0075] The computer device can determine an interpolation position between two adjacent parallel sections at each distance corresponding to the parallel section parameters, thereby obtaining multiple interpolation positions.

[0076] For example, if the parallel section parameter Dx is 400, an interpolation position is determined every 400 m on the X-axis. The distance between two parallel sections in the X-axis direction is 1200 m. Two interpolation positions can be determined, and then two layer segment lines can be inserted subsequently.

[0077] Step 103: The computer device determines a plurality of first layer segment lines of the target formation.

[0078] The first layer segment line is composed of a plurality of first control points. The first layer segment line can be continuous or discontinuous, and there is no specific limitation on this. The computer device can identify the first layer segment line by the formation name and segment number.

[0079] If the target stratum includes a single stratum, the computer device determines multiple first-layer segment lines belonging to the stratum. If the target stratum includes multiple strata, in this step, the computer device determines multiple first-layer segment lines belonging to the same stratum. For each stratum, the computer device constructs a geological layer map based on the multiple first-layer segment lines of the stratum, performing the same steps. Therefore, in this embodiment of the application, only the example of a target stratum including a single stratum is used for description.

[0080] There can be one or more first layer segment lines in each parallel section. If there are multiple first layer segment lines, these multiple first layer segment lines are parallel to the first coordinate axis. Figure 2 , Figure 2There is a first layer segment line in each parallel section, and each first layer segment line is parallel to the Y axis.

[0081] Step 104 : The computer device determines a target layer segment line having the largest number of first control points based on the number of first control points corresponding to each first layer segment line.

[0082] Different first-layer segment lines may be composed of different numbers of first control points. In this step, the computer device first determines the number of first control points corresponding to each first-layer segment line, and then determines the first-layer segment line with the largest number of first control points as the target layer segment line.

[0083] For example, the target stratum corresponds to 4 layer segment lines, the first layer segment line is composed of 3 first control points, the second layer segment line is composed of 4 first control points, the third layer segment line is composed of 5 first control points, and the fourth layer segment line is composed of 7 first control points. Then the layer segment line composed of 7 first control points is determined as the target layer segment line.

[0084] Step 105: The computer device determines multiple position scales based on the target layer segment line.

[0085] The position ratio is the ratio of the distance between the first first control point in the target layer segment line and each other first control point to the line segment length of the target layer segment line. The order of the first control points in the target layer segment line is determined according to the coordinate size of each first control point in the direction of the first coordinate axis. The first coordinate axis is parallel to the first layer segment line.

[0086] For example, the target layer segment line is composed of 7 first control points, and these 7 first control points are sorted from small to large according to the coordinates in the X-axis direction. The computer device determines the ratio of the distance between the first first control point and the second first control point to the line segment length of the target layer segment line to obtain the first position ratio. The computer device determines the ratio of the distance between the first first control point and the third first control point to the line segment length of the target layer segment line to obtain the second position ratio. The computer device sequentially determines the ratio of the distance between the first first control point and the fourth first control point to the line segment length of the target layer segment line, the ratio of the distance between the first first control point and the fifth first control point to the line segment length of the target layer segment line, the ratio of the distance between the first first control point and the sixth first control point to the line segment length of the target layer segment line, and the ratio of the distance between the first first control point and the seventh first control point to the line segment length of the target layer segment line to obtain the third position ratio, the fourth position ratio, the fifth position ratio and the sixth position ratio respectively.

[0087] Step 106: The computer device determines a plurality of second control points based on the plurality of position ratios and the segment lengths of other layer segment lines in the plurality of first layer segment lines.

[0088] For each of the plurality of first layer segment lines other than the target layer segment line, the computer device determines the product of the segment length of the other layer segment line and each position ratio to obtain a plurality of first lengths. For each first length, the computer device determines the coordinates corresponding to the first control point of the other layer segment line in the direction of the first coordinate axis after the first length has been separated, and determines the point corresponding to the coordinates as a second control point.

[0089] Step 107: The computer device determines the plurality of second control points as the plurality of first control points on other layer segment lines.

[0090] For each other layer segment line, except for the first first control point whose position remains unchanged, the position of each of the remaining first control points is adjusted accordingly according to the multiple position ratios determined in step 105. After the adjustment, each first layer segment line has the same number of first control points with the same position ratio.

[0091] Step 108: The computer device determines the coordinates of multiple interpolation points at each interpolation position based on the three-dimensional coordinates of the first control points on the re-determined other layer segment lines and the three-dimensional coordinates of multiple first control points on the target layer segment line.

[0092] Step 108 can be implemented by following the steps (1) to (5), including:

[0093] (1) For each first control point of a position ratio, the computer device determines a relationship curve consisting of a preset number of adjacent first control points of the position ratio on other layer segment lines and the target layer segment line.

[0094] The relationship curve is a curve with a first distance ratio as an independent variable, the first distance ratio is the ratio of the first distance to the second distance, the first distance is the distance between the third control point and the first first control point among a preset number of first control points, the third control point is a control point at any position on the relationship curve, and the second distance is the sum of the distances between each two adjacent first control points among the preset number of first control points.

[0095] The preset number can be set and changed as needed, for example, the preset number is 3 or 4. In the embodiment of the present application, only the preset number of 3 is used for explanation. That is, when the preset number is 3, step (1) is for the computer device to determine the relationship curve composed of every three adjacent first control points of the position ratio on other layer segment lines and the target layer segment line.

[0096] Continue with Figure 2For example, the four first-layer segment lines in the figure are used to illustrate this. After position scaling, each first-layer segment line includes seven control points. For the first control point of each first-layer segment line, the computer device sequentially determines the relationship curves formed by each of the three adjacent first control points according to the coordinates of each first-layer segment line on the X-axis, obtaining two relationship curves, which are then overlapped to form one relationship curve. For the second control point of each first-layer segment line, the computer device sequentially determines the relationship curves formed by each of the three adjacent second control points in the X-axis direction, obtaining two relationship curves, which are then overlapped to form one relationship curve. For the third control point of each first-layer segment line, the computer device sequentially determines the relationship curves formed by each of the three adjacent third control points in the X-axis direction, obtaining two relationship curves, which are then overlapped to form one relationship curve. By analogy, the relationship curve composed of the last three adjacent control points in the X-axis direction is determined, and two relationship curves are obtained respectively. These two relationship curves are overlapped to form one relationship curve, and finally 7 relationship curves are obtained. These 7 relationship curves are all relationship curves in the X-axis direction. Figure 3 .

[0097] (2) The computer device determines a plurality of curve coefficients of the relationship curve based on the three-dimensional coordinates of each first control point constituting the relationship curve.

[0098] The curve coefficients include the constant term coefficient, the linear term coefficient, and the quadratic term coefficient, and each coefficient is a three-dimensional coordinate point. Among them, the constant term coefficient is the coefficient without an independent variable, the linear term coefficient is the coefficient with the independent variable exponent equal to one, and the quadratic term coefficient is the coefficient with the independent variable exponent equal to two.

[0099] In step (1), it has been explained that the relationship curve is a quadratic curve. Then the quadratic curve can be expressed as: S(t) = A1 + A2t + A3t 2 , where t represents the first distance ratio and t∈[0,1], A1, A2 and A3 all represent curve coefficients, where A1 represents the constant term coefficient, A2 represents the linear term coefficient, A3 represents the quadratic term coefficient, and S(t) represents the coordinates of the interpolation point corresponding to the interpolation position when the first distance ratio is t.

[0100] For example, the first control points constituting the relationship curve are P1, P2, and P3, respectively. P1 is the first control point at the beginning of the relationship curve, P3 is the first control point at the end of the relationship curve, and P2 is the first control point between the beginning and the end. Based on the coordinates of the three first control points, it can be obtained that:

[0101]

[0102] in, d1=|P1-P2|, d2=|P2-P3|, d1 represents the first distance between points P1 and P2, d2 represents the seventh distance between points P2 and P3, t0 represents the ratio of the first distance to the second distance, and the second distance is the sum of the first distance and the seventh distance.

[0103] Since the coordinates of P1, P2 and P3 are known, the first distance and the seventh distance can be determined, thereby determining t0 at point P2, and then determining A1, A2 and A3 based on the coordinates of P1, P2 and P3 and t0.

[0104] (3) The computer device determines a three-dimensional expression of the relationship curve based on the multiple curve coefficients.

[0105] The computer device determines the functional relationship corresponding to the multiple curve coefficients to obtain a three-dimensional expression of the relationship curve.

[0106] (4) The computer device determines the sub-coefficient of each curve coefficient in the direction of the second coordinate axis.

[0107] Each curve coefficient is a three-dimensional coordinate point, and the computer device determines a sub-coefficient of each curve coefficient in the direction of a second coordinate axis, where the second coordinate axis is perpendicular to the first coordinate axis and is in the same horizontal plane.

[0108] For example, the first coordinate axis is the Y axis, and the second coordinate axis is the X axis. The computer device determines the sub-coefficient of each curve coefficient in the X axis direction, that is, determines the X coordinate of each curve coefficient as the sub-coefficient of the curve coefficient in the X axis direction.

[0109] (5) For each interpolation position, the computer device determines the coordinates of the interpolation point corresponding to the interpolation position at the position scale based on the sub-coefficients of each curve coefficient, the three-dimensional expression of the relationship curve, and the coordinates of the interpolation position on the second coordinate axis.

[0110] From step (1), it can be seen that for the first control point of the same position scale, the final relationship curve is formed by overlapping two relationship curves. Therefore, the relationship curves of some lengths overlap, while the relationship curves of other lengths do not overlap. If the relationship curve where the interpolation position is located is the first relationship curve without overlap, the coordinates of the interpolation point corresponding to the interpolation position at the position scale can be directly determined based on the sub-coefficients of each curve coefficient, the three-dimensional expression of the relationship curve, and the coordinates of the interpolation position on the second coordinate axis. This process can be seen in the following steps (B-1) to (B-5):

[0111] (B-1) The computer device determines a one-dimensional expression of the first relationship curve based on the sub-coefficients of each curve coefficient.

[0112] The computer device determines an expression composed of multiple sub-coefficients of the curve coefficients based on the sub-coefficients of each curve coefficient, and obtains a one-dimensional expression of the first relationship curve, which is an expression of the first relationship curve in the direction of the second coordinate axis.

[0113] For example, the second coordinate axis is the X axis, with A 1x represents the constant term coefficient, A 2x represents the coefficient of the first-order term, A 3x represents the coefficient of the quadratic term, then the one-dimensional expression of the first relationship curve in the X-axis direction is S(t)=A 1x +A 2x t+A 3x t 2 .

[0114] (B-2) The computer device corrects the constant term sub-coefficient in the one-dimensional expression based on the constant term sub-coefficient in the one-dimensional expression and the coordinates of the interpolation position in the direction of the second coordinate axis to obtain a corrected constant term sub-coefficient.

[0115] The computer device determines the value of the constant term sub-coefficient in the one-dimensional expression minus the coordinate of the interpolation position in the direction of the second coordinate axis to obtain a corrected constant term sub-coefficient.

[0116] For example, if the coordinate of the interpolation position in the X-axis direction is X1, then the corrected constant term coefficient is A 1x -X1.

[0117] (B-3) The computer device determines the one-dimensional expression composed of the corrected constant term sub-coefficient, the linear term sub-coefficient and the quadratic term sub-coefficient in the one-dimensional expression.

[0118] The computer device keeps the linear term sub-coefficient and the quadratic term sub-coefficient in the one-dimensional expression of the first relationship curve unchanged, and recomposes the corrected constant term sub-coefficient, linear term sub-coefficient and quadratic term sub-coefficient into a one-dimensional expression.

[0119] For example, the reconstructed one-dimensional expression is S(t)=(A 1x -X1)+A 2x t+A 3x t 2 .

[0120] (B-4) The computer device determines the second distance ratio corresponding to when the value of the reorganized one-dimensional expression is zero.

[0121] For example, computer equipment determines (A 1x -X1)+A 2x t+A 3x t 2 =0, and the second distance ratio is obtained.

[0122] (B-5) When the computer device determines that the independent variable of the first relationship curve is the second distance ratio, the value of the three-dimensional expression of the first relationship curve is obtained, and the coordinates of the interpolation point corresponding to the interpolation position at the position ratio are obtained.

[0123] The computer device substitutes the second distance ratio into the three-dimensional expression of the first relationship curve, determines the value corresponding to the three-dimensional expression, and obtains the interpolation point coordinates corresponding to the interpolation position at the position ratio.

[0124] For example, if the second distance ratio is t′, then S(t′)=A1+A2t′+A3t′ 2 Since A1, A2, A3 and t′ are all known, the value of S(t′) can be determined. The value of S(t′) is the coordinate of the interpolation point at the interpolation position at the position scale.

[0125] It should be noted that the coordinates of the interpolation point determined along the second coordinate axis are identical to the coordinates of the interpolation position along the second coordinate axis. The process of determining the coordinates of the interpolation point corresponding to an interpolation position at a certain position scale is essentially the process of determining the coordinates of the interpolation position along a third coordinate axis at that position scale. The third coordinate axis is perpendicular to the horizontal plane formed by the first and second coordinate axes. For example, if the first coordinate axis is the Y axis and the second coordinate axis is the X axis, then the third coordinate axis is the Z axis.

[0126] If the relationship curve at the interpolation position is the first relationship curve obtained by overlapping the second and third relationship curves, it is necessary to redefine the three-dimensional expression of the overlapping relationship curve based on the second and third relationship curves, and then determine the coordinates of the interpolation point corresponding to the interpolation position. This process can be seen in the following steps (C-1) to (C-6):

[0127] (C-1) The computer device determines a plurality of second curve coefficients of the second relationship curve and a plurality of third curve coefficients of the third relationship curve.

[0128] The computer device determines multiple second curve coefficients of the second relationship curve based on the three-dimensional coordinates of each first control point constituting the second relationship curve; and determines multiple third curve coefficients of the third relationship curve based on the three-dimensional coordinates of each first control point constituting the third relationship curve.

[0129] The process of determining the plurality of second curve coefficients and the plurality of third curve coefficients by the computer device is similar to the process of determining the plurality of curve coefficients in step (2), and will not be repeated here.

[0130] (C-2) The computer device determines a third distance ratio and a fourth distance ratio.

[0131] The third distance ratio is the ratio of the third distance to the fourth distance, the third distance is the distance between the first first control point and the second first control point in the second relationship curve, and the fourth distance is the sum of the distances between each two adjacent first control points in the second relationship curve.

[0132] The fourth distance ratio is the ratio of the fifth distance to the sixth distance, the fifth distance is the distance between the first first control point and the second first control point in the third relationship curve, and the sixth distance is the sum of the distances between each two adjacent first control points in the third relationship curve.

[0133] Because the second relationship curve and the third relationship curve overlap, the first and second first control points in the third relationship curve are also two of the first control points that make up the second relationship curve. For the second relationship curve, the range of values of its independent variable in the overlapping portion is greater than the third distance ratio and less than 1. For the third relationship curve, the range of values of its independent variable in the overlapping portion is greater than 0 and less than the fourth distance ratio.

[0134] For example, there are four first control points, namely P i 、P i+1 、P i+2 and P i+3 , every three adjacent first control points form a relationship curve, then P i 、P i+1 and P i+2 Composing the second relationship curve, P i+1 、P i+2 and P i+3 The third relationship curve is formed, and the first control point common to the second and third relationship curves is P i+1 and P i+2 The first relationship curve obtained by overlapping the second relationship curve and the third relationship curve is P i+1 and P i+2 Correspondingly, with r1 representing the third distance ratio and r2 representing the fourth distance ratio, the third distance ratio and the fourth distance ratio can be expressed as follows:

[0135] and

[0136] Among them, |P i -P i+1 | indicates P i and P i+1 The distance between |P i+1 -P i+2 | indicates P i+1 and P i+2 The distance between |P i+2 -P i+3| indicates P i+2 and P i+3 The distance between them.

[0137] For the second curve S1(t1)=A4+A5t1+A6t1 2 For example, the value range of t1 in the overlapping part satisfies t1∈[r1,1]; for the third curve S2(t2)=A7+A8t2+A9t2 2 For example, the value range of t2 in the overlapping part satisfies t2∈[0,r2].

[0138] (C-3) The computer device determines a three-dimensional expression of the first relationship curve with the fifth distance ratio as an independent variable based on the plurality of second curve coefficients, the plurality of third curve coefficients, the third distance ratio, and the fourth distance ratio.

[0139] The fifth distance ratio is the ratio of the seventh distance to the eighth distance, the seventh distance is the distance between the interpolation point at the interpolation position and the first control point at the beginning of the first relationship curve, and the eighth distance is the distance between the first control point at the beginning and the first control point at the end of the first curve.

[0140] In this step, the computer device determines the first weight of the second relationship curve and the second weight of the third relationship curve based on the position of the interpolation position in the first relationship curve, and then based on the first weight and the second weight, performs weighted summation on the three-dimensional expression of the second relationship curve and the three-dimensional expression of the third relationship curve, respectively, to obtain the three-dimensional expression of the first relationship curve.

[0141] The closer the interpolation position is to the first control point at the beginning of the first relationship curve, the greater the first weight of the second relationship curve. The closer the interpolation position is to the first control point at the end of the first relationship curve, the greater the second weight of the third relationship curve. Based on this, if t1 is used to represent the first weight k and the second weight m, then If t2 is used to represent the first weight k and the second weight m, then but According to k or m, it can be determined

[0142] Based on the third distance ratio and the fourth distance ratio, the computer device represents the independent variable of either the second or third relationship curve using the independent variable of the other relationship curve, and then represents the three-dimensional expression of the other relationship curve using the independent variable of the other relationship curve. For example, the computer device represents the independent variable of the third relationship curve using the independent variable of the second relationship curve, and then represents the three-dimensional expression of the third relationship curve using the independent variable of the second relationship curve; or, based on the third distance ratio and the fourth distance ratio, represents the independent variable of the second relationship curve using the independent variable of the third relationship curve, and then represents the three-dimensional expression of the second relationship curve using the independent variable of the third relationship curve.

[0143] The computer device performs a weighted summation of the three-dimensional expression of the second relationship curve and the three-dimensional expression of the third relationship curve based on the first weight and the second weight to obtain the three-dimensional expression of the first relationship curve, wherein the three-dimensional expression of the second relationship curve is expressed using the independent variable of the third relationship curve, or the three-dimensional expression of the third relationship curve is expressed using the independent variable of the second relationship curve.

[0144] For example, Here we only take t1 to represent t2 as an example. like Then t2 = a(t1-r1). Expressing S2(t2) by t1, then

[0145] Let t1 represent the first weight k and the second weight m, To simplify the description, if Then k = b(1-t1), m = b(t1-r1), the computer device performs weighted summation of S1(t1) and S2(t2) based on the first weight and the second weight, then

[0146] According to the above content, S2(t2)=A7+A8a(t1-r1)+A9a 2 (t1-r1) 2 ,but By splitting and merging, the three-dimensional expression of the first relationship curve can be obtained as follows:

[0147] t1=t, then t is replaced by t, and the three-dimensional expression of the first relationship curve is:

[0148]

[0149] (C-4) The computer device determines a one-dimensional expression of the three-dimensional expression of the first relationship curve in the direction of the second coordinate axis.

[0150] The computer device determines a sub-coefficient of each curve coefficient in the direction of the second coordinate axis in the three-dimensional expression of the first relationship curve, and composes a one-dimensional expression of the first relationship curve based on the sub-coefficient of each curve coefficient.

[0151] For example, if the second coordinate axis is the X axis, the computer device determines the X coordinates of A4, A5, A6, A7, A8, and A9, and then composes a one-dimensional expression of the first relationship curve.

[0152] (C-5) The computer device corrects the constant term sub-coefficient in the one-dimensional expression based on the constant term sub-coefficient in the one-dimensional expression, the coordinates of the interpolation position in the direction of the second coordinate axis, and the third distance ratio to obtain a corrected constant term sub-coefficient.

[0153] The computer device first determines 1 minus the value of the third distance ratio to obtain a first difference value, then determines the product of the coordinate of the interpolation position in the direction of the second coordinate axis and the first difference value to obtain a first product value, and finally determines the value of the constant term sub-coefficient minus the first product value to obtain a corrected constant term sub-coefficient.

[0154] For example, the constant term coefficient is (A4+A8ar1 2 -A9a 2 r1 3 -A7r1), the coordinate of the interpolation position in the direction of the second coordinate axis is X2, then the computer device determines (A4+A8ar1 2 -A9a 2 r1 3 -A7r1)-(X2×(1-r1)) to obtain the corrected constant term coefficient.

[0155] (C-6) The computer device determines the coordinates of the interpolation point at the interpolation position at the position scale based on the one-dimensional expression recomposed from the corrected constant term sub-coefficient, the linear term sub-coefficient and the quadratic term sub-coefficient in the one-dimensional expression.

[0156] The computer device maintains the linear sub-coefficient and the quadratic sub-coefficient in the one-dimensional expression of the first curve unchanged, and reconstructs the modified constant sub-coefficient, the linear sub-coefficient, and the quadratic sub-coefficient into a one-dimensional expression. The computer device determines the sixth distance ratio corresponding to a value of zero in the one-dimensional expression, determines the value of the three-dimensional expression of the first relationship curve when the independent variable of the first curve is the sixth distance ratio, and obtains the coordinates of the interpolation point corresponding to the interpolation position at the position ratio.

[0157] This process is similar to steps (B-3) to (B-5) and will not be described again here.

[0158] Step 109: The computer device interpolates a plurality of second-layer segment lines based on the coordinates of the plurality of interpolation points at each interpolation position.

[0159] The second-layer segment lines are composed of multiple interpolation points. For each interpolation position, the computer device sequentially connects the multiple interpolation points corresponding to the interpolation position to obtain the second-layer segment lines corresponding to the interpolation position, and ultimately obtains multiple second-layer segment lines.

[0160] See also Figure 4 , Figure 4 For Figure 2 FIG. 1 is a schematic diagram of 14 second-layer segment lines inserted by the above method based on the 4 first-layer segment lines shown, and the 14 second-layer segment lines are all parallel to the Y axis.

[0161] Step 110: The computer device constructs a geological layer map based on the multiple first layer segment lines and the multiple second layer segment lines of the target stratum.

[0162] Based on the multiple first-layer segment lines and multiple second-layer segment lines of the target stratum, the computer constructs a continuous triangulated mesh. Using a surface fitting algorithm, the vertex coordinates of the triangulated mesh are determined and connected to create a three-dimensional geological layer map. This process is similar to methods used in related technologies and will not be further described here.

[0163] See also Figure 5 , Figure 5 Based on Figure 4 The geological layer map constructed by the multiple first layer segment lines and the multiple second layer segment lines shown in FIG. Figure 5 It can be seen from the figure that the geological layer map is smooth and flat, so the effect of the constructed geological layer map is good.

[0164] An embodiment of the present application provides a method for constructing a geological layer map. The method interpolates multiple layer lines through existing profile data to obtain interpolated profile data, and then constructs a geological layer map based on the existing profile data and the interpolated profile data. Since the geological layer map is constructed with a large amount of profile data, the method can construct a smooth and natural geological layer map when there is less existing profile data, thereby improving the effect of the constructed geological layer map.

[0165] The present application embodiment provides a device for constructing a geological layer map, see Figure 6 , the device comprises:

[0166] An acquisition module 601 is configured to acquire a plurality of cross-sectional data of a target stratum and determine a plurality of parallel cross-sectional views based on the plurality of cross-sectional data, wherein the cross-sectional data includes a three-dimensional coordinate of a first control point, which is any position point in the target stratum.

[0167] A first determining module 602 is configured to determine a plurality of interpolation positions based on predetermined parallel section parameters, wherein the interpolation positions are used to insert layer segment lines parallel to the plurality of parallel sections between the plurality of parallel sections;

[0168] A second determining module 603 is used to determine a plurality of first layer segment lines of the target stratum, where the first layer segment lines are formed by connecting a plurality of first control points;

[0169] a third determining module 604 for determining coordinates of a plurality of interpolation points at each interpolation position based on the three-dimensional coordinates of the first control point corresponding to each first layer segment line;

[0170] an interpolation module 605 for interpolating a plurality of second-layer segment lines based on the coordinates of a plurality of interpolation points at each interpolation position, wherein the second-layer segment lines are formed by connecting the plurality of interpolation points;

[0171] The construction module 606 is configured to construct a geological layer map based on the plurality of first layer segment lines and the plurality of second layer segment lines of the target stratum.

[0172] In one possible implementation, the third determination module 604 is used to determine the target layer segment line with the largest number of first control points based on the number of first control points corresponding to each first layer segment line; determine multiple position ratios based on the target layer segment line, where the position ratio is the ratio of the distance between the first first control point and each other first control point in the target layer segment line to the line segment length of the target layer segment line, and the order of the first control points in the target layer segment line is determined according to the coordinate size of each first control point in the direction of the first coordinate axis, and the first coordinate axis is parallel to the first layer segment line; determine multiple second control points based on the multiple position ratios and the line segment lengths of other layer segment lines in the multiple first layer segment lines, where the second control points are control points on other layer segment lines corresponding to the multiple position ratios; determine the multiple second control points as multiple first control points on other layer segment lines; and determine the coordinates of multiple interpolation points at each interpolation position based on the three-dimensional coordinates of the re-determined first control points on other layer segment lines and the three-dimensional coordinates of the multiple first control points on the target layer segment line.

[0173] In another possible implementation, the third determination module 604 is used to determine, for each first control point of each position ratio, a relationship curve composed of each adjacent preset number of first control points of the position ratio on other layer segment lines and the target layer segment line, the relationship curve is a curve with the first distance ratio as the independent variable, the first distance ratio is the ratio of the first distance to the second distance, the first distance is the distance between the third control point and the first first control point among the preset number of first control points, the third control point is a control point at any position on the relationship curve, and the second distance is the sum of the distances between each adjacent two first control points among the preset number of first control points; based on the three-dimensional coordinates of each first control point constituting the relationship curve, determine the relationship curve. A plurality of curve coefficients of a relationship curve, the curve coefficients include a constant term coefficient, a linear term coefficient and a quadratic term coefficient, and each curve coefficient is a three-dimensional coordinate point, the constant term coefficient is a coefficient without an independent variable, the linear term coefficient is a coefficient with an independent variable exponent of one, and the quadratic term coefficient is a coefficient with an independent variable exponent of two; based on the plurality of curve coefficients, a three-dimensional expression of the relationship curve is determined; a sub-coefficient of each curve coefficient in the direction of a second coordinate axis is determined, where the second coordinate axis is perpendicular to the first coordinate axis and is in the same horizontal plane; for each interpolation position, the coordinates of the interpolation point corresponding to the interpolation position under the position scale are determined based on the sub-coefficients of each curve coefficient, the three-dimensional expression of the relationship curve and the coordinates of the interpolation position on the second coordinate axis.

[0174] In another possible implementation, if the relationship curve where the interpolation position is located is a first relationship curve without overlap, the third determination module 604 is used to determine a one-dimensional expression of the first relationship curve based on the sub-coefficients of each curve coefficient; based on the constant term sub-coefficients in the one-dimensional expression and the coordinates of the interpolation position in the direction of the second coordinate axis, correct the constant term sub-coefficients in the one-dimensional expression to obtain corrected constant term sub-coefficients; determine a one-dimensional expression recomposed of the corrected constant term sub-coefficients, the linear term sub-coefficients, and the quadratic term sub-coefficients in the one-dimensional expression; determine the second distance ratio corresponding to when the value of the recomposed one-dimensional expression is zero; determine the value of the three-dimensional expression of the first relationship curve when the independent variable of the first relationship curve is the second distance ratio, and obtain the coordinates of the interpolation point corresponding to the interpolation position under the position ratio.

[0175] In another possible implementation, if the relationship curve where the interpolation position is located is the first relationship curve obtained by overlapping the second relationship curve and the third relationship curve, the third determination module 604 is used to determine multiple second curve coefficients of the second relationship curve and multiple third curve coefficients of the third relationship curve; determine a third distance ratio and a fourth distance ratio, the third distance ratio is the ratio of the third distance to the fourth distance, the third distance is the distance between the first first control point and the second first control point in the second relationship curve, the fourth distance is the sum of the distances between every two adjacent first control points in the second relationship curve, the fourth distance ratio is the ratio of the fifth distance to the sixth distance, the fifth distance is the distance between the first first control point and the second first control point in the third relationship curve, and the sixth distance is the sum of the distances between every two adjacent first control points in the third relationship curve; based on the multiple second curve coefficients, the multiple third The curve coefficient, the third distance ratio and the fourth distance ratio are used to determine a three-dimensional expression of the first relationship curve with the fifth distance ratio as the independent variable, the fifth distance ratio is the ratio of the seventh distance to the eighth distance, the seventh distance is the distance between the interpolation point at the interpolation position and the first control point at the beginning of the first relationship curve, and the eighth distance is the distance between the first control point at the beginning and the first control point at the end of the first relationship curve; a one-dimensional expression of the three-dimensional expression of the first relationship curve in the direction of the second coordinate axis is determined; based on the constant term sub-coefficient in the one-dimensional expression, the coordinates of the interpolation position in the direction of the second coordinate axis and the third distance ratio, the constant term sub-coefficient in the one-dimensional expression is corrected to obtain the corrected constant term sub-coefficient; based on the one-dimensional expression recomposed from the corrected constant term sub-coefficient, the linear term sub-coefficient and the quadratic term sub-coefficient in the one-dimensional expression, the coordinates of the interpolation point at the interpolation position under the position ratio are determined.

[0176] In another possible implementation, the first determining module 602 is configured to determine an interpolation position between two adjacent parallel sections at every distance corresponding to the parallel section parameter, to obtain a plurality of interpolation positions.

[0177] In another possible implementation, the interpolation module 605 is configured to sequentially connect multiple interpolation points corresponding to each interpolation position to obtain a second-layer segment line corresponding to the interpolation position.

[0178] An embodiment of the present application provides a device for constructing a geological layer map. The device interpolates multiple layer lines through existing profile data to obtain interpolated profile data, and then constructs a geological layer map based on the existing profile data and the interpolated profile data. Since the geological layer map is constructed with a large amount of profile data, the device can construct a smooth and natural geological layer map with less existing profile data, thereby improving the effect of the constructed geological layer map.

[0179] It should be noted that the geological layer map construction device provided in the above embodiment only uses the division of the above-mentioned functional modules as an example to illustrate the construction of the geological layer map. In actual applications, the above-mentioned functions can be assigned to different functional modules as needed, that is, the internal structure of the device can be divided into different functional modules to complete all or part of the functions described above. In addition, the geological layer map construction device provided in the above embodiment and the geological layer map construction method embodiment are based on the same concept. The specific implementation process is detailed in the method embodiment and will not be repeated here.

[0180] Figure 7 The following is a block diagram of a computer device 700 according to an exemplary embodiment of the present application. The computer device 700 may be a portable mobile computer device, such as a smartphone, a tablet computer, an MP3 player (Moving Picture Experts Group Audio Layer III), an MP4 player (Moving Picture Experts Group Audio Layer IV), a laptop computer, or a desktop computer. The computer device 700 may also be referred to as a user device, a portable computer device, a laptop computer device, a desktop computer device, or other similar names.

[0181] Typically, the computer device 700 includes a processor 701 and a memory 702 .

[0182] The processor 701 may include one or more processing cores, such as a 4-core processor, an 8-core processor, etc. The processor 701 may be implemented in at least one hardware form of DSP (Digital Signal Processing), FPGA (Field-Programmable Gate Array), or PLA (Programmable Logic Array). The processor 701 may also include a main processor and a coprocessor. The main processor is a processor for processing data in the awake state, also known as a CPU (Central Processing Unit); the coprocessor is a low-power processor for processing data in the standby state. In some embodiments, the processor 701 may be integrated with a GPU (Graphics Processing Unit), which is responsible for rendering and drawing the content to be displayed on the display screen. In some embodiments, the processor 701 may also include an AI (Artificial Intelligence) processor, which is used to process computing operations related to machine learning.

[0183] Memory 702 may include one or more computer-readable storage media, which may be non-transitory. Memory 702 may also include high-speed random access memory and non-volatile memory, such as one or more magnetic disk storage devices or flash memory storage devices. In some embodiments, the non-transitory computer-readable storage medium in memory 702 is used to store at least one program code, which is executed by processor 701 to implement the method for constructing a geological layer map provided in the method embodiment of the present application.

[0184] In some embodiments, computer device 700 may optionally include a peripheral device interface 703 and at least one peripheral device. Processor 701, memory 702, and peripheral device interface 703 may be connected via a bus or signal lines. Each peripheral device may be connected to peripheral device interface 703 via a bus, signal lines, or circuit boards. Specifically, the peripheral device may include at least one of a radio frequency circuit 704, a display screen 705, a camera assembly 706, an audio circuit 707, a positioning assembly 708, and a power supply 709.

[0185] The peripheral device interface 703 can be used to connect at least one I / O (Input / Output)-related peripheral device to the processor 701 and the memory 702. In some embodiments, the processor 701, the memory 702, and the peripheral device interface 703 are integrated on the same chip or circuit board; in some other embodiments, any one or two of the processor 701, the memory 702, and the peripheral device interface 703 can be implemented on separate chips or circuit boards, which is not limited in this embodiment.

[0186] The RF circuit 704 is used to receive and transmit RF (Radio Frequency) signals, also known as electromagnetic signals. The RF circuit 704 communicates with communication networks and other communication devices via electromagnetic signals. The RF circuit 704 converts electrical signals into electromagnetic signals for transmission, or converts received electromagnetic signals into electrical signals. Optionally, the RF circuit 704 includes an antenna system, an RF transceiver, one or more amplifiers, a tuner, an oscillator, a digital signal processor, a codec chipset, a user identity module card, and the like. The RF circuit 704 can communicate with other computer devices via at least one wireless communication protocol. Such wireless communication protocols include, but are not limited to, the World Wide Web, metropolitan area networks, intranets, various generations of mobile communication networks (2G, 3G, 4G, and 5G), wireless local area networks, and / or WiFi (Wireless Fidelity) networks. In some embodiments, the RF circuit 704 may also include circuits related to NFC (Near Field Communication), which is not limited in this application.

[0187] Display screen 705 is used to display a user interface (UI). This UI may include graphics, text, icons, videos, or any combination thereof. When display screen 705 is a touchscreen display, it is also capable of collecting touch signals on or above the surface of display screen 705. These touch signals can be input as control signals to processor 701 for processing. Display screen 705 can also be used to provide virtual buttons and / or a virtual keyboard, also known as soft buttons and / or a soft keyboard. In some embodiments, there can be one display screen 705, located on the front panel of computer device 700. In other embodiments, there can be at least two display screens 705, located on different surfaces of computer device 700 or in a foldable design. In other embodiments, display screen 705 can be a flexible display, located on a curved or foldable surface of computer device 700. Display screen 705 can also be configured as a non-rectangular, irregular shape, also known as a special-shaped screen. Display screen 705 can be made of materials such as LCD (Liquid Crystal Display) and OLED (Organic Light-Emitting Diode).

[0188] The camera assembly 706 is used to capture images or videos. Optionally, the camera assembly 706 includes a front camera and a rear camera. Typically, the front camera is set on the front panel of the computer device, and the rear camera is set on the back of the computer device. In some embodiments, there are at least two rear cameras, which are any one of a main camera, a depth of field camera, a wide-angle camera, and a telephoto camera, so as to realize the fusion of the main camera and the depth of field camera to realize the background blur function, the fusion of the main camera and the wide-angle camera to realize panoramic shooting and VR (Virtual Reality) shooting function or other fusion shooting functions. In some embodiments, the camera assembly 706 may also include a flash. The flash can be a single-color temperature flash or a dual-color temperature flash. A dual-color temperature flash refers to a combination of a warm light flash and a cold light flash, which can be used for light compensation at different color temperatures.

[0189] The audio circuit 707 may include a microphone and a speaker. The microphone is used to collect sound waves from the user and the environment, and convert the sound waves into electrical signals to be input into the processor 701 for processing, or input into the radio frequency circuit 704 to achieve voice communication. For the purpose of stereo sound collection or noise reduction, there can be multiple microphones, which are respectively arranged in different parts of the computer device 700. The microphone can also be an array microphone or an omnidirectional collection microphone. The speaker is used to convert the electrical signals from the processor 701 or the radio frequency circuit 704 into sound waves. The speaker can be a traditional thin film speaker or a piezoelectric ceramic speaker. When the speaker is a piezoelectric ceramic speaker, it can not only convert the electrical signals into sound waves audible to humans, but also convert the electrical signals into sound waves inaudible to humans for purposes such as ranging. In some embodiments, the audio circuit 707 may also include a headphone jack.

[0190] Positioning component 708 is used to locate the current geographic location of computer device 700 to implement navigation or LBS (Location Based Service). Positioning component 708 can be a positioning component based on the US GPS (Global Positioning System), China's Beidou system, or Russia's Galileo system.

[0191] Power supply 709 is used to power the various components of computer device 700. Power supply 709 can be AC power, DC power, a disposable battery, or a rechargeable battery. When power supply 709 includes a rechargeable battery, the rechargeable battery can be a wired rechargeable battery or a wireless rechargeable battery. A wired rechargeable battery is a battery that is charged via a wired line, while a wireless rechargeable battery is a battery that is charged via a wireless coil. The rechargeable battery can also be used to support fast charging technology.

[0192] In some embodiments, the computer device 700 further includes one or more sensors 710 , including but not limited to: an acceleration sensor 711 , a gyroscope sensor 712 , a pressure sensor 713 , a fingerprint sensor 714 , an optical sensor 715 , and a proximity sensor 716 .

[0193] The accelerometer 711 can detect the magnitude of acceleration along the three coordinate axes of the coordinate system established by the computer device 700. For example, the accelerometer 711 can be used to detect the components of gravity acceleration along the three coordinate axes. The processor 701 can control the display screen 705 to display the user interface in a landscape or portrait view based on the gravity acceleration signal collected by the accelerometer 711. The accelerometer 711 can also be used to collect game or user motion data.

[0194] The gyroscope sensor 712 can detect the orientation and rotation angle of the computer device 700. It can also work with the accelerometer 711 to collect 3D motions of the user on the computer device 700. Based on the data collected by the gyroscope sensor 712, the processor 701 can implement the following functions: motion sensing (e.g., changing the UI based on the user's tilt), image stabilization during shooting, game control, and inertial navigation.

[0195] The pressure sensor 713 can be installed on the side frame of the computer device 700 and / or below the display screen 705. When the pressure sensor 713 is installed on the side frame of the computer device 700, it can detect the user's grip signal of the computer device 700. The processor 701 can perform left and right hand recognition or shortcut operations based on the grip signal collected by the pressure sensor 713. When the pressure sensor 713 is installed below the display screen 705, the processor 701 controls the operational controls on the UI interface based on the user's pressure operation on the display screen 705. The operational controls include at least one of a button control, a scroll bar control, an icon control, and a menu control.

[0196] The fingerprint sensor 714 is used to collect the user's fingerprint. The processor 701 identifies the user's identity based on the fingerprint collected by the fingerprint sensor 714, or the fingerprint sensor 714 identifies the user's identity based on the collected fingerprint. When the user's identity is recognized as a trusted identity, the processor 701 authorizes the user to perform relevant sensitive operations, such as unlocking the screen, viewing encrypted information, downloading software, making payments, and changing settings. The fingerprint sensor 714 can be set on the front, back, or side of the computer device 700. When a physical button or manufacturer logo is set on the computer device 700, the fingerprint sensor 714 can be integrated with the physical button or manufacturer logo.

[0197] The optical sensor 715 is used to detect ambient light intensity. In one embodiment, the processor 701 can control the display brightness of the display screen 705 based on the ambient light intensity detected by the optical sensor 715. Specifically, when the ambient light intensity is high, the display brightness of the display screen 705 is increased; when the ambient light intensity is low, the display brightness of the display screen 705 is decreased. In another embodiment, the processor 701 can also dynamically adjust the shooting parameters of the camera assembly 706 based on the ambient light intensity detected by the optical sensor 715.

[0198] Proximity sensor 716, also known as a distance sensor, is typically located on the front panel of computer device 700. Proximity sensor 716 is used to detect the distance between the user and the front of computer device 700. In one embodiment, when proximity sensor 716 detects that the distance between the user and the front of computer device 700 is gradually decreasing, processor 701 controls display screen 705 to switch from the screen-on state to the screen-off state. When proximity sensor 716 detects that the distance between the user and the front of computer device 700 is gradually increasing, processor 701 controls display screen 705 to switch from the screen-off state to the screen-on state.

[0199] Those skilled in the art will understand that Figure 7 The structure shown in the figure does not constitute a limitation on the computer device 700, and the computer device 700 may include more or fewer components than shown in the figure, or combine some components, or adopt a different component arrangement.

[0200] An embodiment of the present application also provides a computer-readable storage medium, in which at least one program code is stored. The at least one program code is loaded and executed by a processor to implement the method for constructing a geological layer map in the embodiment of the present application.

[0201] An embodiment of the present application also provides a computer program product, which stores at least one program code, and the at least one program code is loaded and executed by a processor to implement the method for constructing a geological layer map in the embodiment of the present application.

[0202] In some embodiments, the computer program involved in the embodiments of the present application may be deployed and executed on a computer device, or on multiple computer devices located at one location, or on multiple computer devices distributed at multiple locations and interconnected through a communication network. Multiple computer devices distributed at multiple locations and interconnected through a communication network may constitute a blockchain system.

[0203] The above description is only for the purpose of facilitating those skilled in the art to understand the technical solution of this application and is not intended to limit this application. Any modifications, equivalent replacements, improvements, etc. made within the spirit and principles of this application shall be included in the scope of protection of this application.

Claims

1. A method for constructing a geological layer map, characterized in that: The method comprises: Acquire multiple cross-sectional data of a target stratum, and determine multiple parallel cross-sectional views based on the multiple cross-sectional data, wherein the cross-sectional data includes three-dimensional coordinates of a first control point, where the first control point is any position point in the target stratum; Based on predetermined parallel section parameters, a plurality of interpolation positions are determined, wherein the interpolation positions are used to insert layer segment lines parallel to the plurality of parallel sections between the plurality of parallel sections; Determine a plurality of first layer segment lines of the target stratum, wherein the first layer segment lines are formed by connecting a plurality of first control points; Determining a target layer segment line having the largest number of first control points based on the number of first control points corresponding to each first layer segment line; Determining a plurality of position ratios based on the target layer segment line, wherein the position ratio is a ratio of a distance between a first first control point and each of the other first control points in the target layer segment line to a segment length of the target layer segment line, wherein an order of the first control points in the target layer segment line is determined according to a coordinate size of each first control point in a direction of a first coordinate axis, wherein the first coordinate axis is parallel to the first layer segment line; Determining a plurality of second control points based on the plurality of position ratios and the segment lengths of other layer segment lines in the plurality of first layer segment lines, where the second control points are control points on the other layer segment lines corresponding to the plurality of position ratios; determining the plurality of second control points as the plurality of first control points on the other layer segment lines; determining coordinates of a plurality of interpolation points at each interpolation position based on the re-determined three-dimensional coordinates of the first control point on the other layer segment line and the three-dimensional coordinates of the plurality of first control points on the target layer segment line; inserting a plurality of second-layer segment lines based on the coordinates of the plurality of interpolation points at each interpolation position, wherein the second-layer segment lines are formed by connecting the plurality of interpolation points; A geological layer map is constructed based on the multiple first layer segment lines and the multiple second layer segment lines of the target stratum.

2. The method according to claim 1, characterized in that The step of determining the coordinates of multiple interpolation points at each interpolation position based on the re-determined three-dimensional coordinates of the first control point on the other layer segment line and the three-dimensional coordinates of multiple first control points on the target layer segment line includes: For each first control point of each position ratio, determine a relationship curve consisting of each adjacent preset number of first control points of the position ratio on the other layer segment lines and the target layer segment line, wherein the relationship curve is a curve with a first distance ratio as an independent variable, the first distance ratio is a ratio of a first distance to a second distance, the first distance is a distance between a third control point and a first first control point among a preset number of first control points, the third control point is a control point at any position on the relationship curve, and the second distance is a sum of distances between every two adjacent first control points among the preset number of first control points; Determining a plurality of curve coefficients of the relationship curve based on the three-dimensional coordinates of each first control point constituting the relationship curve, wherein the curve coefficients include a constant term coefficient, a linear term coefficient, and a quadratic term coefficient, and each curve coefficient is a three-dimensional coordinate point, the constant term coefficient is a coefficient without an independent variable, the linear term coefficient is a coefficient with an independent variable exponent of one, and the quadratic term coefficient is a coefficient with an independent variable exponent of two; determining a three-dimensional expression of the relationship curve based on the plurality of curve coefficients; determining a sub-coefficient of each curve coefficient in the direction of a second coordinate axis, where the second coordinate axis is perpendicular to the first coordinate axis and is in the same horizontal plane; For each interpolation position, based on the sub-coefficients of each curve coefficient, the three-dimensional expression of the relationship curve and the coordinates of the interpolation position on the second coordinate axis, the coordinates of the interpolation point corresponding to the interpolation position at the position scale are determined.

3. The method according to claim 2, characterized in that If the relationship curve where the interpolation position is located is a first relationship curve without overlap, determining the coordinates of the interpolation point corresponding to the interpolation position at the position scale based on the sub-coefficients of each curve coefficient, the three-dimensional expression of the relationship curve, and the coordinates of the interpolation position on the second coordinate axis includes: determining a one-dimensional expression of the first relationship curve based on the sub-coefficients of each curve coefficient; Correcting the constant term sub-coefficient in the one-dimensional expression based on the constant term sub-coefficient in the one-dimensional expression and the coordinate of the interpolation position in the direction of the second coordinate axis to obtain a corrected constant term sub-coefficient; Determine a one-dimensional expression composed of the modified constant term sub-coefficient, the linear term sub-coefficient, and the quadratic term sub-coefficient in the one-dimensional expression; determining a second distance ratio corresponding to when the value of the reconstructed one-dimensional expression is zero; When determining that the independent variable of the first relationship curve is the second distance ratio, the value of the three-dimensional expression of the first relationship curve is used to obtain the coordinates of the interpolation point corresponding to the interpolation position at the position ratio.

4. The method according to claim 2, characterized in that If the relationship curve where the interpolation position is located is a first relationship curve obtained by overlapping a second relationship curve and a third relationship curve, determining the coordinates of the interpolation point corresponding to the interpolation position at the position scale based on the sub-coefficients of each curve coefficient, the three-dimensional expression of the relationship curve, and the coordinates of the interpolation position on the second coordinate axis includes: determining a plurality of second curve coefficients of the second relationship curve and a plurality of third curve coefficients of the third relationship curve; Determine a third distance ratio and a fourth distance ratio, wherein the third distance ratio is a ratio of a third distance to a fourth distance, wherein the third distance is the distance between a first first control point and a second first control point in the second relationship curve, and the fourth distance is the sum of the distances between every two adjacent first control points in the second relationship curve. The fourth distance ratio is a ratio of a fifth distance to a sixth distance, wherein the fifth distance is the distance between a first first control point and a second first control point in the third relationship curve, and the sixth distance is the sum of the distances between every two adjacent first control points in the third relationship curve. determining, based on the plurality of second curve coefficients, the plurality of third curve coefficients, the third distance ratio, and the fourth distance ratio, a three-dimensional expression of a first relationship curve with a fifth distance ratio as an independent variable, the fifth distance ratio being a ratio of a seventh distance to an eighth distance, the seventh distance being a distance between an interpolation point at the interpolation position and a first control point at a beginning of the first relationship curve, and the eighth distance being a distance between the first control point at the beginning and the first control point at an end of the first relationship curve; Determine a one-dimensional expression of the three-dimensional expression of the first relationship curve in the direction of the second coordinate axis; Correcting the constant term sub-coefficient in the one-dimensional expression based on the constant term sub-coefficient in the one-dimensional expression, the coordinate of the interpolation position in the direction of the second coordinate axis, and the third distance ratio to obtain a corrected constant term sub-coefficient; The coordinates of the interpolation point of the interpolation position at the position scale are determined based on the one-dimensional expression recomposed from the modified constant term sub-coefficient, the linear term sub-coefficient and the quadratic term sub-coefficient in the one-dimensional expression.

5. The method according to claim 1, wherein The step of determining a plurality of interpolation positions based on predetermined parallel profile parameters includes: Between two adjacent parallel sections, an interpolation position is determined for each distance corresponding to the parallel section parameter, thereby obtaining a plurality of interpolation positions.

6. The method according to claim 1, characterized in that The inserting of a plurality of second-layer segment lines based on the coordinates of the plurality of interpolation points at each interpolation position comprises: For each interpolation position, multiple interpolation points corresponding to the interpolation position are sequentially connected to obtain a second-layer segment line corresponding to the interpolation position.

7. A device for constructing a geological layer map, characterized in that: The device comprises: an acquisition module, configured to acquire a plurality of cross-sectional data of a target stratum, and determine a plurality of parallel cross-sectional views based on the plurality of cross-sectional data, wherein the cross-sectional data includes a three-dimensional coordinate of a first control point, and the first control point is any position point in the target stratum; A first determining module is configured to determine a plurality of interpolation positions based on predetermined parallel section parameters, wherein the interpolation positions are used to insert layer segment lines parallel to the plurality of parallel sections between the plurality of parallel sections; A second determining module is configured to determine a plurality of first layer segment lines of the target stratum, wherein the first layer segment lines are formed by connecting a plurality of first control points; A third determination module is configured to determine a target layer segment line with the largest number of first control points based on the number of first control points corresponding to each first layer segment line; determine a plurality of position ratios based on the target layer segment line, wherein the position ratio is the ratio of the distance between the first first control point and each other first control point in the target layer segment line to the segment length of the target layer segment line, and the order of the first control points in the target layer segment line is determined according to the coordinate size of each first control point in the direction of a first coordinate axis, wherein the first coordinate axis is parallel to the first layer segment line; determine a plurality of second control points based on the plurality of position ratios and the segment lengths of other layer segment lines in the plurality of first layer segment lines, wherein the second control points are control points on the other layer segment lines corresponding to the plurality of position ratios; determine the plurality of second control points as a plurality of first control points on the other layer segment lines; determine the coordinates of a plurality of interpolation points at each interpolation position based on the three-dimensional coordinates of the re-determined first control points on the other layer segment lines and the three-dimensional coordinates of the plurality of first control points on the target layer segment line; an interpolation module, configured to interpolate a plurality of second-layer segment lines based on the coordinates of the plurality of interpolation points at each interpolation position, wherein the second-layer segment lines are formed by connecting the plurality of interpolation points; A construction module is used to construct a geological layer map based on multiple first layer segment lines and multiple second layer segment lines of the target stratum.

8. A computer device, characterized in that: The computer device includes a processor and a memory, wherein the memory stores at least one program code, and the at least one program code is loaded and executed by the processor to implement the method for constructing a geological layer map according to any one of claims 1 to 6.

9. A computer-readable storage medium, characterized in that The computer-readable storage medium stores at least one program code, and the at least one program code is loaded and executed by the processor to implement the method for constructing a geological layer map according to any one of claims 1 to 6.

Images