Geological profile map value evaluation method and system

Abstract

The invention discloses a geological profile map value evaluation method and system, and the method specifically comprises the steps: S1, obtaining a to-be-evaluated geological profile map, and carrying out the preprocessing of the geological profile map; s2, performing grid division on the geological profile map obtained by preprocessing; s3, constructing a stratum matrix based on position information and stratum information of all grids in the divided geological profile map, and coding the stratum information to obtain a stratum coding matrix; and S4, calculating the two-dimensional entropy of the stratum coding matrix, and substituting the calculated two-dimensional entropy into the constructed pricing function to obtain a pricing result of the geological profile map. According to the invention, a clear geological profile map pricing process and formula are given, and the method is used for pricing and trading the geological profile map, so that the conversion from qualitative pricing to quantitative pricing of the geological profile map is realized; the pricing strategy introduces the concept of information entropy, and the complexity of the stratum is quantitatively considered in the pricing process.

Description

Technical Field

[0001] The present invention relates to the technical field of engineering geology, and particularly relates to a method and system for evaluating the value of a geological section diagram. Background Art

[0002] The geological information of a certain section line is stored in the geological section diagram, which is the basic data for designers in the engineering field to carry out foundation and subgrade design. At present, with the development of the digital economy, the secondary trading of the achievement data of geological section diagrams has been gradually promoted. The production cost of traditional geological data can be directly obtained by calculating its engineering quantity and unit price, while the trading of the achievement of geological section diagrams focuses on the secondary utilization of existing geological section diagram achievements. Due to the replicability of electronic data and the non-destructiveness in the trading process, the production cost of existing geological section diagrams is almost zero. Therefore, how to reasonably determine the price of geological section diagrams is an important research direction. At present, in the process of trading geological achievement data, the price is mainly determined through negotiation. The negotiation-based pricing is highly subjective, which is not conducive to the development of geological data trading and is difficult to reflect the true value of geological section diagrams. Summary of the Invention

[0003] This application provides a method and system for evaluating the value of a geological section diagram to solve the problem that the existing negotiation-based pricing for the trading of geological section diagram achievements is highly subjective and difficult to reflect the true value of geological section diagrams.

[0004] According to a first aspect, in one embodiment, a method for evaluating the value of a geological section diagram is provided. The method includes:

[0005] Step S1, obtaining the geological section diagram to be evaluated and preprocessing the geological section diagram;

[0006] Step S2, dividing the preprocessed geological section diagram into grids;

[0007] Step S3, based on the position information and formation information of all grids in the divided geological section diagram, constructing a formation matrix and obtaining a formation coding matrix by encoding the formation information;

[0008] Step S4, calculating the two-dimensional entropy of the formation coding matrix, substituting the calculated two-dimensional entropy into the constructed pricing function, and obtaining the pricing result of the geological section diagram.

[0009] Further, the step S1 specifically includes:

[0010] Performing cleaning processing on the geological section diagram, only retaining the necessary formation stratification information on the longitudinal section, and deleting the unnecessary annotation information on the geological section diagram.

[0011] Further, the step S2 specifically includes:

[0012] Obtain the minimum bounding rectangle of the preprocessed geological profile, and divide the minimum bounding rectangle into small grids of equal spacing size.

[0013] Further, step S3 specifically includes:

[0014] Step S31, construct a formation matrix according to the divided grids. The initial formation matrix is:

[0015]

[0016] There are m rows and n columns in the divided grids, where diceng represents the formation name in each grid.

[0017] Further, step S3 specifically includes:

[0018] Step S32, encode all formations on the profile according to their categories. After encoding, the formation matrix is transformed into:

[0019]

[0020] There are m rows and n columns in total, where bianma represents the encoding of the formation category in each grid.

[0021] Further, step S32 specifically includes:

[0022] Use positive integers increasing from 0 to encode all formation categories on the geological profile. The null value within the grid range is encoded as 0, and the remaining formations are encoded sequentially starting from 1 in the order of the formation age from new to old, with each increment being 1.

[0023] Further, step S4 specifically includes:

[0024] Step S41, define i as any element in the formation encoding matrix B, J as a sub - matrix of B that contains element i, the number of rows and columns of the sub - matrix J is less than that of matrix B, j as the sum of all elements in the sub - matrix J, traverse the entire matrix B to obtain all (i, j) pairs, and calculate the two - dimensional entropy of the formation encoding matrix B:

[0025]

[0026] In the formula, bianma max is the maximum encoding value after formation encoding, and P(i, j) is the probability of the occurrence of (i, j).

[0027] Further, step S4 specifically includes:

[0028] Step S42, construct a pricing function, and the pricing function needs to meet the following conditions:

[0029] For any geological section diagram p1, its pricing is always greater than 0, which is expressed as:

[0030] For any two geological section diagrams p1 and p2, if the same-sized grid division is adopted and the two-dimensional entropy H(p2) of p2 is greater than or equal to the two-dimensional entropy H(p1) of p1, then the pricing of p2 is greater than or equal to the pricing of p1, which is expressed as:

[0031]

[0032] For any two geological section diagrams p1 and p2, for the new geological section diagram p1 + p2 formed by combining p1 and p2, the pricing of the new geological section diagram p1 + p2 should be less than or equal to the sum of the pricing of p1 and the pricing of p2, which is expressed as: price(p1 + p2) ≤ price(p1) + price(p2).

[0033] Further, step S4 specifically includes:

[0034] In step S43, the constructed pricing function has the form: price(p) = k * H(p), where k is a set coefficient, k > 0, and H(p) is the two-dimensional entropy of the geological section diagram p. The pricing result is obtained by substituting the calculated two-dimensional entropy value of the geological section diagram into the pricing function.

[0035] According to the second aspect, an embodiment provides a geological section diagram value evaluation system, which includes:

[0036] A preprocessing module, which is used to obtain the geological section diagram to be evaluated and preprocess the geological section diagram;

[0037] A grid division module, which is used to perform grid division on the geological section diagram obtained by preprocessing;

[0038] A matrix construction module, which is used to construct a stratum matrix based on the position information and stratum information of all grids in the divided geological section diagram, and obtain a stratum coding matrix by encoding the stratum information;

[0039] A pricing module, which is used to calculate the two-dimensional entropy of the stratum coding matrix, substitute the calculated two-dimensional entropy into the constructed pricing function, and obtain the pricing result of the geological section diagram.

[0040] This application provides a geological section diagram value evaluation method and system, which gives a clear geological section diagram pricing process and formula, and realizes the transformation from qualitative pricing to quantitative pricing of the geological section diagram; this pricing strategy introduces the concept of information entropy, and quantitatively considers the complexity of the stratum during the pricing process. It has the following beneficial effects:

[0041] (1) Improve the tradability of data: By valuing geological cross-sections, they can be transformed into quantifiable data assets and traded in the data market. This helps realize the value of geological data and promotes the circulation and sharing of data resources.

[0042] (2) Enhance the value of data assets: As part of data assets, the valuation of geological cross-sections can increase the value on the balance sheet and bring potential economic benefits to enterprises. This not only helps the performance of enterprises in the capital market but also serves as an important basis for corporate mergers and acquisitions and investment decisions.

[0043] (3) Promote the development of the digital economy: The valuation and digital management of geological cross-sections are important parts of the digital economy. Through digital means, the utilization efficiency of geological data can be improved, promoting the digital transformation of related industries and providing new impetus for economic growth. BRIEF DESCRIPTION OF THE DRAWINGS

[0044] Figure 1 It is a flowchart of a method for evaluating the value of a geological cross-section provided by an embodiment of the present invention;

[0045] Figure 2 It is the geological cross-section before processing in a method for evaluating the value of a geological cross-section provided by an embodiment of the present invention;

[0046] Figure 3 It is the geological cross-section after processing in a method for evaluating the value of a geological cross-section provided by an embodiment of the present invention;

[0047] Figure 4 It is the maximum and minimum value nodes of the horizontal and vertical coordinates of the geological cross-section in a method for evaluating the value of a geological cross-section provided by an embodiment of the present invention;

[0048] Figure 5 It is the minimum circumscribed rectangle area of the geological cross-section in a method for evaluating the value of a geological cross-section provided by an embodiment of the present invention;

[0049] Figure 6 It is a schematic diagram of the logical structure of a system for evaluating the value of a geological cross-section provided by an embodiment of the present invention. DETAILED DESCRIPTION OF THE INVENTION

[0050] The present invention will be further described in detail below in conjunction with the specific embodiments and the accompanying drawings. Similar elements in different embodiments are denoted by related similar element numbers. In the following embodiments, many detailed descriptions are provided to enable a better understanding of the present application. However, those skilled in the art can easily recognize that some of the features can be omitted in different situations, or can be replaced by other elements, materials, or methods. In some cases, some operations related to the present application are not shown or described in the specification in order to avoid the core part of the present application being overwhelmed by excessive descriptions. For those skilled in the art, it is not necessary to describe these related operations in detail, and they can fully understand the related operations based on the description in the specification and the general technical knowledge in the art.

[0051] In addition, the features, operations, or characteristics described in the specification can be combined in any appropriate manner to form various embodiments. At the same time, the steps or actions in the method description can also be reordered or adjusted in a manner obvious to those skilled in the art. Therefore, the various sequences in the specification and the drawings are only for clearly describing a certain embodiment and do not mean that they are the necessary sequences, unless it is stated that a certain sequence must be followed.

[0052] The present invention aims to provide a quantitative pricing strategy for geological cross-sections to promote the trading of the achievement data of geological cross-sections. A method for evaluating the value of a geological cross-section provided in the first embodiment of the present invention will be described in detail below in conjunction with Figure 1 this.

[0053] As Figure 1 shown, in step S1, a geological cross-section to be evaluated is obtained, and the geological cross-section is preprocessed.

[0054] The above steps specifically include: performing a cleaning process on the geological cross-section, only retaining the necessary formation stratification information on the longitudinal section, and deleting the unnecessary annotation information on the geological cross-section.

[0055] Specifically, for the existing geological cross-section, the annotations such as borehole text annotations, borehole lines, in-situ test data, sampling symbols, water level lines, weathering lines, and descriptions are deleted. Only the formation stratification information on the longitudinal section is retained, such as formation patterns and formation lines. For example Figure 2 is the geological cross-section before processing, and the elements on the drawing include formation name standards, borehole number annotations, stratification information annotations, sampling symbols, water level lines, formation lines, lithology pattern fills, etc. After the cleaning process, Figure 3 only includes formation lines and formation patterns.

[0056] As Figure 1 shown, in step S2, the preprocessed geological cross-section is divided into grids.

[0057] The above steps specifically include: obtaining the minimum bounding rectangle of the preprocessed geological profile, and dividing the minimum bounding rectangle into small grids of equal spacing size.

[0058] First, obtain the minimum bounding rectangle of the effective data of the geological profile. The effective data refers to the stratigraphic data of the geological profile. The maximum and minimum values of the abscissa and the maximum and minimum values of the ordinate on the geological profile form the four corner point coordinates of the minimum bounding rectangle. The four corner point coordinates of the minimum bounding rectangle are (maximum abscissa, minimum ordinate), (maximum abscissa, maximum ordinate), (minimum abscissa, minimum ordinate), and (minimum abscissa, maximum ordinate). As Figure 4 shown, the maximum and minimum values of the abscissa and ordinate of the geological profile can be obtained from nodes 1, 2, and 3 in the schematic diagram. The minimum bounding rectangle of the geological profile formed by the maximum and minimum values of the abscissa and ordinate is as Figure 5 shown.

[0059] Then, divide the minimum bounding rectangle of the geological profile into small grids of equal spacing size. The actual length represented by the side length of the grid is not greater than 1 cm, and the length and width of the grid may not be equal. Figure 5 In, the length of the minimum bounding rectangle is 78 m, the height is 35 m, and the side lengths in the length and height directions of the divided small grids are both 1 cm. Then the total number of divided grids is 780 * 350 = 273000.

[0060] As Figure 1 shown, in step S3, based on the position information and stratigraphic information of all grids in the divided geological profile, construct a stratigraphic matrix, and obtain a stratigraphic coding matrix by encoding the stratigraphic information.

[0061] The above steps specifically include:

[0062] Step S31, construct a stratigraphic matrix according to the divided grids. The initial stratigraphic matrix is:

[0063]

[0064] The divided grids have m rows and n columns, where diceng represents the stratigraphic name in each grid.

[0065] Step S32, encode all the strata on the profile. After encoding, the stratigraphic matrix is transformed into:

[0066]

[0067] It has m rows and n columns, where bianma represents the code of the strata in each grid.

[0068] In this embodiment, positive integers increasing from 0 are used to encode all formation categories on the profile section. The null value within the rectangle is 0, and for the remaining formations, they are encoded incrementally starting from 1 in the order of the formation time of the formations from the youngest to the oldest, with an increment value of 1 each time.

[0069] For example, for a certain geological profile section, after rasterization, the formation matrix:

[0070]

[0071] Among them, the null value is numbered with 0. The formation time of the miscellaneous fill is later than that of the loess. The miscellaneous fill is numbered 1, and the loess is numbered 2. After encoding, the formation matrix is:

[0072]

[0073] As Figure 1 shown, in step S4, calculate the two-dimensional entropy of the formation encoding matrix, substitute the calculated two-dimensional entropy into the constructed pricing function, and obtain the pricing result of the geological profile section.

[0074] The above steps specifically include:

[0075] Step S41, define i as any element in the formation encoding matrix B, J contains the element i, J is a sub-matrix of B, the number of rows and columns of the sub-matrix J is less than the number of rows and columns of the matrix B, j is the sum of all elements in the sub-matrix J, traverse the entire matrix B to obtain all (i, j) pairs, and calculate the two-dimensional entropy of the formation encoding matrix B:

[0076]

[0077] In the formula, bianma max is the maximum encoding value after formation encoding, and P(i, j) is the probability of the occurrence of (i, j). Using the above formula to characterize the complexity of the profile section not only includes the information of the formation area ratio but also reflects the relationship between adjacent formations.

[0078] For example, for the matrix If the size of the sub-matrix J is 2*2, then after traversing the matrix B, the distributions of J, j, and (i, j) are shown in the following table:

[0079]

[0080]

[0081] Statistical probability of the occurrence of each (i, j) pair, the occurrence probability of (0, 3) is 1 / 9, the occurrence probability of (1, 4) is 2 / 9, the occurrence probability of (1, 6) is 3 / 9 = 1 / 3, and the occurrence probability of (2, 8) is 3 / 9 = 1 / 3.

[0082] Calculate the corrected two-dimensional entropy. As can be seen from the two-dimensional entropy formula, the calculation range is non-zero formation codes. Then

[0083]

[0084] Step S42: Construct a pricing function. The pricing function needs to meet the following conditions:

[0085] (1) For any geological profile p1, its pricing is always greater than 0, which is expressed as:

[0086]

[0087] (2) For any two geological profiles p1 and p2, with the same grid division size, if the two-dimensional entropy H(p2) of p2 is greater than or equal to the two-dimensional entropy H(p1) of p1, then the pricing of p2 is greater than or equal to the pricing of p1, which is expressed as:

[0088] (3) For any two geological profiles p1 and p2, the new geological profile p1 + p2 formed by combining p1 and p2, the pricing of the new geological profile p1 + p2 should be less than or equal to the sum of the pricing of p1 and the pricing of p2, which is expressed as: price(p1 + p2) ≤ price(p1) + price(p2).

[0089] Step S43: The form of the pricing function constructed in this embodiment is: price(p) = k * H(p), where k is a set coefficient, k > 0, and H(p) is the two-dimensional entropy of the geological profile p. The pricing result is obtained by substituting the calculated two-dimensional entropy value of the geological profile into the pricing function. For example, the two-dimensional entropy H of the profile calculated in Step S41 is 1.54. Assuming k is taken as 100, then the pricing of this profile is price = 100 * 1.54 = 154.

[0090] Corresponding to the above-disclosed method for evaluating the value of a geological profile, the embodiment of the present invention also discloses a system for evaluating the value of a geological profile, as Figure 6 shown, which specifically includes:

[0091] A preprocessing module for obtaining the geological profile to be evaluated and preprocessing the geological profile;

[0092] A grid division module for dividing the geological profile obtained by preprocessing into grids;

[0093] A matrix construction module for constructing a formation matrix based on the position information and formation information of all grids in the divided geological profile, and obtaining a formation coding matrix by encoding the formation information;

[0094] A pricing module is configured to calculate the two-dimensional entropy of a formation coding matrix, substitute the calculated two-dimensional entropy into a constructed pricing function, and obtain a pricing result of a geological profile.

[0095] It should be noted that for a detailed description of a geological profile value evaluation system provided in an embodiment of the present invention, reference can be made to the relevant description of a geological profile value evaluation method provided in an embodiment of the present application, which will not be elaborated here.

[0096] Those skilled in the art can understand that all or part of the functions of the above methods can be implemented in a hardware manner or in a computer program manner. When all or part of the functions in the above embodiments are implemented in a computer program manner, the program can be stored in a computer-readable storage medium. The storage medium may include: read-only memory, random access memory, magnetic disk, optical disk, hard disk, etc. The above functions can be realized by a computer executing the program. For example, the program is stored in the memory of a device, and when the processor executes the program in the memory, the above all or part of the functions can be realized. In addition, when all or part of the functions in the above embodiments are implemented in a computer program manner, the program can also be stored in a storage medium such as a server, another computer, magnetic disk, optical disk, flash drive or mobile hard disk, downloaded or copied and saved to the memory of a local device, or the system of the local device is updated. When the processor executes the program in the memory, the above all or part of the functions in the embodiments can be realized.

[0097] The above uses specific examples to illustrate the present invention, which is only for helping to understand the present invention and is not intended to limit the present invention. For those skilled in the art of the present invention, according to the idea of the present invention, several simple deductions, deformations or substitutions can also be made.

Claims

1. A method for evaluating the value of a geological profile, characterized in that: The method comprises: Step S1, obtaining a geological profile to be evaluated, and preprocessing the geological profile; Step S2, gridding the geological profile obtained by preprocessing; Step S3, constructing a stratigraphic matrix based on the position information and stratigraphic information of all grids in the divided geological profile, and obtaining a stratigraphic coding matrix by encoding the stratigraphic information; Step S4, calculate the two-dimensional entropy of the stratigraphic coding matrix, substitute the calculated two-dimensional entropy into the constructed pricing function, and obtain the pricing result of the geological profile.

2. A method for evaluating the value of a geological profile according to claim 1, characterized in that: The step S1 specifically includes: The geological profile is cleaned to retain only the necessary stratigraphic information on the longitudinal section and delete unnecessary annotation information on the geological profile.

3. A method for evaluating the value of a geological profile according to claim 1, characterized in that: The step S2 specifically includes: The minimum bounding rectangle of the preprocessed geological profile is obtained, and the minimum bounding rectangle is divided into fine grids of equal spacing.

4. A method for evaluating the value of a geological profile according to claim 1, characterized in that: The step S3 specifically includes: Step S31, constructing a stratigraphic matrix according to the divided grids, the initial stratigraphic matrix is: The divided grid has m rows and n columns, where diceng represents the name of the stratum in each grid.

5. A method for evaluating the value of a geological profile as claimed in claim 4, characterized in that: The step S3 specifically includes: Step S32, encode all the strata on the cross-section diagram, and the stratum matrix after encoding is converted into: There are m rows and n columns in total, where bianma represents the code of the stratum in each grid.

6. A method for evaluating the value of a geological profile according to claim 5, characterized in that: The step S32 specifically includes: All stratigraphic categories on the geological profile are coded using positive integers increasing from 0. The null value within the grid range is coded as 0, and the remaining strata are coded in ascending order from the newest to the oldest, with each increment being 1.

7. A method for evaluating the value of a geological profile according to claim 5, characterized in that: The step S4 specifically includes: Step S41, define i as any element in the stratigraphic coding matrix B, J as a submatrix of B, J contains element i, the number of rows and columns of submatrix J is less than the number of rows and columns of matrix B, j is the sum of all elements in submatrix J, traverse the entire matrix B, obtain all (i, j) pairs, and calculate the two-dimensional entropy of stratigraphic coding matrix B: In the formula, bianma max is the maximum coding value after the stratigraphic category is coded, and P(i,j) is the probability of (i,j) occurring.

8. A method for evaluating the value of a geological profile according to claim 7, characterized in that: The step S4 specifically includes: Step S42: construct a pricing function, which must meet the following conditions: For any geological profile p1, its price is always greater than 0, expressed as: For any geological profiles p1 and p2, using the same size of grid division, if the two-dimensional entropy H(p2) of p2 is greater than or equal to the two-dimensional entropy H(p1) of p1, then the pricing of p2 is greater than or equal to the pricing of p1, expressed as: price(p2)≥price(p1); For any geological profiles p1 and p2, a new geological profile p1+p2 is formed by merging p1 and p2. The pricing of the new geological profile p1+p2 should be less than or equal to the sum of the pricing of p1 and the pricing of p2, expressed as: price(p1+p2)≤price(p1)+price(p2).

9. A method for evaluating the value of a geological profile as claimed in claim 8, characterized in that: The step S4 specifically includes: Step S43, the pricing function constructed is in the form of: price(p)=k*H(p), k is the set coefficient, k>0, H(p) is the two-dimensional entropy of the geological profile p, and the pricing result is obtained by substituting the calculated two-dimensional entropy value of the geological profile into the pricing function.

10. A geological profile value assessment system, characterized in that: The system comprises: A preprocessing module, used for obtaining a geological profile to be evaluated and preprocessing the geological profile; A grid division module is used to perform grid division on the geological profile obtained through preprocessing; A matrix construction module is used to construct a stratigraphic matrix based on the position information and stratigraphic information of all grids in the divided geological profile, and obtain a stratigraphic coding matrix by encoding the stratigraphic information; The pricing module is used to calculate the two-dimensional entropy of the stratigraphic coding matrix, substitute the calculated two-dimensional entropy into the constructed pricing function, and obtain the pricing result of the geological profile.

Images