Method, system and equipment for identifying top and bottom of sedimentary rock sequence based on gravity center analysis

Abstract

The invention discloses a sedimentary rock sequence top and bottom identification method, system and equipment based on gravity center analysis, and relates to the technical field of geology, photogrammetry and remote sensing, and the method comprises the steps: obtaining a target image, determining the vector boundary of each gravel in each sampling region in the target image, taking any gravel as the current gravel, and determining the current gravel; determining the section gravity center of the current gravel based on the vector boundary of the current gravel; determining a plurality of angle bisectors of the current gravel by taking the gravity center of the section of the current gravel as a circle center and a first preset angle as a step length, and determining intersection points of the angle bisectors and the vector boundary of the current gravel; determining the gravity direction of the current gravel based on all angle bisectors and all intersection points of the current gravel; and based on the gravity directions of all gravels in the to-be-identified area, determining the sequence top-bottom relationship of the sedimentary rock sample in the to-be-identified area. According to the method, the accuracy and applicability of top-bottom relationship identification of the sedimentary rock sequence are improved.

Description

Technical Field

[0001] The present application relates to the fields of geology, photogrammetry and remote sensing technology, and in particular to a method, system and equipment for identifying the top and bottom of a sedimentary rock sequence based on centroid analysis. Background Art

[0002] The investigation of sedimentary rock sequence is a common basic problem in geological work. The predecessors usually judged the top and bottom of sedimentary rocks through a small number of top-indicating structures, such as cross-bedding, wave-shaped ripples, mud cracks, etc. However, in geological practice, it is very difficult to find the above-mentioned specific top-indicating structures in the field, which is difficult to meet the needs of a wide range of sedimentary rock sequence top and bottom identification work and has poor applicability. In addition, affected by geological structural transformation, the occurrence and stratigraphic position of some sedimentary rocks have changed significantly. For example, the occurrence of sedimentary rocks in the field appears to be inverted, upright, inclined, or has obvious folds. Neither field investigations nor indoor studies can accurately identify sedimentary rock sequences and their top and bottom relationships.

[0003] At present, there is no reliable method for identifying the top and bottom relationship of sedimentary rock sequences. Therefore, there is an urgent need for a method for identifying the top and bottom of sedimentary rock sequences to solve the above problems. Summary of the invention

[0004] The purpose of this application is to provide a method, system and equipment for identifying the top and bottom of a sedimentary rock sequence based on centroid analysis, thereby improving the accuracy and applicability of the identification of the top and bottom relationship of a sedimentary rock sequence.

[0005] To achieve the above objectives, this application provides the following solutions:

[0006] In a first aspect, the present application provides a method for identifying the top and bottom of a sedimentary rock sequence based on centroid analysis, the method for identifying the top and bottom of a sedimentary rock sequence based on centroid analysis comprising:

[0007] Acquire a target image; the target image is an image of a sedimentary rock sample in the area to be identified; the area to be identified includes a plurality of sampling areas; each sampling area includes a plurality of gravels;

[0008] Determine the vector boundary of each gravel in each sampling area in the target image;

[0009] Set any gravel as the current gravel;

[0010] Based on the vector boundary of the current gravel, determine the center of gravity of the current gravel section;

[0011] Taking the center of gravity of the cross section of the current gravel as the center of the circle and the first preset angle as the step length, determining multiple angle bisectors of the current gravel, and determining the intersection of each angle bisector and the vector boundary of the current gravel;

[0012] Based on all angle bisectors and all intersections of the current gravel, determine the gravity direction of the current gravel; wherein the gravity direction is the direction pointing to the bottom of the rock layer of the sedimentary rock sample;

[0013] Based on the gravity directions of all gravels in the area to be identified, the top and bottom relationships of the sequences of sedimentary rock samples in the area to be identified are determined.

[0014] Optionally, obtaining a target image specifically includes:

[0015] Carry out cross-sectional photography of sedimentary rock samples in the area to be identified to obtain an initial sedimentary rock sample image;

[0016] The initial sedimentary rock sample image is preprocessed to obtain the target image.

[0017] Optionally, the preprocessing method includes image binarization, sharpening and contrast adjustment.

[0018] Optionally, based on all angle bisectors and all intersection points of the current gravel, the gravity direction of the current gravel is determined, specifically including:

[0019] Based on all angle bisectors and all intersections of the current gravel, draw a gravity direction analysis and judgment diagram of the current gravel;

[0020] Based on the current gravel gravity direction analysis and judgment diagram, the current gravel gravity direction is determined.

[0021] Optionally, based on all angle bisectors and all intersections of the current gravel, a gravity direction analysis and judgment diagram of the current gravel is drawn, specifically including:

[0022] With the angle as the horizontal axis, the first preset angle as the step length, and the distance between the center of the circle along the angle bisector of the current gravel and the corresponding intersection point as the vertical axis, a gravity direction analysis and judgment diagram of the current gravel is drawn.

[0023] Optionally, based on the gravity directions of all gravels in the area to be identified, determining the top-bottom relationship of the sequence of sedimentary rock samples in the area to be identified specifically includes:

[0024] Based on the gravity direction of all gravels in each sampling area, a rosette diagram of gravity advantage direction of the corresponding sampling area is drawn;

[0025] Based on the gravity advantage pointing rose diagram of each sampling area, the gravity direction of the corresponding sampling area is determined;

[0026] Based on the gravity direction of all sampling areas, the top and bottom relationships of the sedimentary rock samples in the area to be identified are determined.

[0027] Optionally, based on the gravity direction of all gravels in each sampling area, a rose diagram of gravity dominance pointing to the corresponding sampling area is drawn, specifically including:

[0028] Count the frequencies of different gravity directions of gravels in each sampling area;

[0029] Based on the frequency of occurrence of different gravity directions of gravel in each sampling area, a rose diagram of the gravity advantage direction of the corresponding sampling area was drawn.

[0030] Optionally, based on the gravity directions of all sampling areas, determining the top and bottom relationships of the sequences of sedimentary rock samples in the area to be identified specifically includes:

[0031] Based on the gravity directions of all sampling areas, a gravity direction point map of the area to be identified is drawn; the gravity direction point map includes multiple gravity direction points; one gravity direction point corresponds to the gravity direction of one sampling area;

[0032] Taking any two gravity direction points as the point groups to be screened, multiple point groups to be screened are obtained;

[0033] Determine any point group to be screened as the current point group, and determine the number of gravity direction points whose gravity directions are between the gravity directions corresponding to two gravity direction points in the current point group as the current number;

[0034] Determine whether the current point group meets the preset conditions; the preset conditions are that the difference in gravity direction of the sampling area corresponding to the two gravity direction points in the current point group does not exceed 30 degrees, and the current number is greater than all non-current numbers; the non-current number is the number of gravity direction points whose gravity directions are between the gravity directions corresponding to the two gravity direction points in the non-current point group; the non-current point group is all the point groups to be screened except the current point group in all the point groups to be screened;

[0035] If so, then the gravity directions of the two gravity direction points in the current point group and the gravity direction points whose gravity directions are between the gravity directions corresponding to the two gravity direction points in the current point group are averaged to obtain the top and bottom relationship of the sequence of the sedimentary rock samples in the area to be identified.

[0036] In a second aspect, the present application provides a sedimentary rock sequence top and bottom identification system based on centroid analysis, the sedimentary rock sequence top and bottom identification system based on centroid analysis is used to implement the sedimentary rock sequence top and bottom identification method based on centroid analysis, and the sedimentary rock sequence top and bottom identification system based on centroid analysis includes:

[0037] The target image acquisition unit is used to acquire a target image; the target image is an image of a sedimentary rock sample in the area to be identified; the area to be identified includes a plurality of sampling areas; each sampling area includes a plurality of gravels;

[0038] A vector boundary determination unit, used for determining the vector boundary of each gravel in each sampling area in the target image;

[0039] A current gravel determination unit, used for taking any gravel as the current gravel;

[0040] A section center determination unit, used for determining the section center of the current gravel based on the vector boundary of the current gravel;

[0041] An angle bisector and intersection determination unit, used to determine multiple angle bisectors of the current gravel with the center of gravity of the section of the current gravel as the center of the circle and the first preset angle as the step length, and determine the intersection of each angle bisector and the vector boundary of the current gravel;

[0042] The unit for determining the gravity direction of the gravel is used to determine the gravity direction of the current gravel based on all angle bisectors and all intersections of the current gravel; wherein the gravity direction is the direction pointing to the bottom of the rock layer of the sedimentary rock sample;

[0043] The stratigraphic top-bottom relationship determination unit is used to determine the stratigraphic top-bottom relationship of the sedimentary rock samples in the area to be identified based on the gravity directions of all gravels in the area to be identified.

[0044] In a third aspect, the present application provides a computer device comprising: a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the computer program to implement any of the above-mentioned methods for identifying the top and bottom of a sedimentary rock sequence based on centroid analysis.

[0045] According to the specific embodiments provided in this application, this application has the following technical effects:

[0046] The present application discloses a method, system and equipment for identifying the top and bottom of a sedimentary rock sequence based on centroid analysis. By analyzing the stable natural characteristics of gravel during the sedimentation process, calculating the centroid of the gravel section, determining the gravity advantage direction of the gravel, and thus determining the gravity direction of the gravel, the top and bottom relationship of the sequence of the sedimentary rock sample in the area to be identified can be finally identified. The top and bottom relationship of the sequence of the sedimentary rock sample can be accurately determined, and a series of difficult problems encountered in the top and bottom identification of the traditional sedimentary rock sequence can be effectively solved, especially in the case of complex geological structures and rare or unclear top-indicating structures, which provides a new solution for the top and bottom identification of sedimentary rock sequences. At the same time, this method is applicable to a wider range of geological environments and sedimentary rock types, and enhances the applicability and flexibility of the top and bottom identification of sedimentary sequences. BRIEF DESCRIPTION OF THE DRAWINGS

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

[0048] Figure 1 This is an application environment diagram of a sedimentary rock sequence top and bottom identification method based on centroid analysis in one embodiment of the present application;

[0049] Figure 2 A schematic flow chart of a method for identifying the top and bottom of a sedimentary rock sequence based on centroid analysis provided in one embodiment of the present application;

[0050] Figure 3 A schematic diagram of a target image provided by an embodiment of the present application;

[0051] Figure 4 A schematic diagram of the vector boundary and the center of gravity of the section of gravel provided in an embodiment of the present application;

[0052] Figure 5 A schematic diagram of multiple angle bisectors of a current gravel and intersections of each angle bisector with a vector boundary of the current gravel provided in an embodiment of the present application;

[0053] Figure 6 A diagram for analyzing and determining the gravity direction of gravel 1 provided in one embodiment of the present application;

[0054] Figure 7 A diagram for analyzing and determining the gravity direction of gravel 2 provided in one embodiment of the present application;

[0055] Figure 8 A diagram for analyzing and determining the gravity direction of gravel 3 provided in one embodiment of the present application;

[0056] Fig. 9 A rosette diagram is provided for the gravity advantage provided in one embodiment of the application;

[0057] Fig.10 A schematic diagram of functional modules of a sedimentary rock sequence top and bottom identification system based on centroid analysis provided in one embodiment of the present application;

[0058] Fig.11 A schematic diagram of the structure of a computer device provided in one embodiment of the present application.

[0059] Reference numerals:

[0060] Terminal 102, server 104, target image acquisition unit 1, vector boundary determination unit 2, current gravel determination unit 3, section center of gravity determination unit 4, angle bisector and intersection determination unit 5, gravel gravity direction determination unit 6, stratigraphic top and bottom relationship determination unit 7. DETAILED DESCRIPTION

[0061] The following will be combined with the drawings in the embodiments of the present application to clearly and completely describe the technical solutions in the embodiments of the present application. Obviously, the described embodiments are only part of the embodiments of the present application, not all of the embodiments. Based on the embodiments in the present application, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of this application.

[0062] In order to make the above-mentioned objects, features and advantages of the present application more obvious and easy to understand, the present application is further described in detail below in conjunction with the accompanying drawings and specific implementation methods.

[0063] The method for identifying the top and bottom of a sedimentary rock sequence based on centroid analysis provided in the embodiments of the present application can be applied to Figure 1 In the application environment shown. Among them, the terminal 102 communicates with the server 104 through the network. The data storage system can store the data that the server 104 needs to process. The data storage system can be set up separately, or it can be integrated on the server 104, or it can be placed on the cloud or other servers. The terminal 102 can send the target image to the server 104. After the server 104 receives the target image, the server 104 performs the top and bottom identification of the sedimentary rock sequence based on the center of gravity analysis for the target image. The server 104 can feed back the top and bottom identification of the sequence of the sedimentary rock samples in the area to be identified to the terminal 102. In addition, in some embodiments, the top and bottom identification method of the sedimentary rock sequence based on the center of gravity analysis can also be implemented separately by the server 104 or the terminal 102, such as the terminal 102 can directly perform the top and bottom identification of the sedimentary rock sequence based on the center of gravity analysis for the target image, or the server 104 can obtain the target image from the data storage system, and perform the top and bottom identification of the sedimentary rock sequence based on the center of gravity analysis for the target image.

[0064] The terminal 102 may be, but is not limited to, various desktop computers, laptop computers, smart phones, tablet computers, IoT devices, and portable wearable devices. The IoT devices may be smart speakers, smart TVs, smart air conditioners, smart vehicle-mounted devices, etc. The portable wearable devices may be smart watches, smart bracelets, head-mounted devices, etc. The server 104 may be implemented as an independent server or a server cluster consisting of multiple servers, or may be a cloud server.

[0065] In an exemplary embodiment, Figure 2As shown, a method for identifying the top and bottom of a sedimentary rock sequence based on centroid analysis is provided. The method is executed by a computer device, and can be executed by a computer device such as a terminal or a server alone, or by a terminal and a server together. In the embodiment of the present application, the method is applied to Figure 1 The server 104 in the example is used for explanation, and the steps include the following steps S1 to S7.

[0066] in:

[0067] Step S1, obtaining a target image; the target image is an image of a sedimentary rock sample in the area to be identified; the target image schematic diagram is as follows Figure 3 As shown in the figure, the area to be identified includes multiple sampling areas; each sampling area includes multiple gravels. Among them, at least three sampling ranges are randomly determined in the area to be identified, marked as Q(1), Q(2), Q(3)...Q(n), as sampling areas. The schematic diagram of the sampling area is shown in the figure Figure 3 The area circled by Q(1), Q(2), and Q(3).

[0068] As an optional implementation, step S1 specifically includes:

[0069] Step S11, taking cross-sectional photographs of the sedimentary rock samples in the area to be identified to obtain an initial sedimentary rock sample image.

[0070] Specifically, sedimentary rock samples include typical field sections, geological drill cores, and typical sedimentary rock samples; typical field sections include geological sections and soil sections; sedimentary rock samples should avoid the influence of strong deformation or metamorphism of tectonic belts, and gravel or sediments should not be significantly affected by tectonic transformation; sedimentary rock samples should meet the conditions of natural or artificial exposure in the field and meet certain scale or quantity conditions; sedimentary rock samples include field sampling and indoor identification samples, and this type of operation method should record and be able to restore the field occurrence of sedimentary rocks. Sedimentary rock samples cover traditional sedimentary rocks and undiagenerated sediments in the field of earth sciences.

[0071] The initial sedimentary rock sample image obtained by profile photography is a digital image or image; profile photography includes field conglomerate profile photography, conglomerate cutting profile photography, borehole core cutting profile photography, and rock section photography under a microscope; among them, field conglomerate profile photography refers to the section to be identified in field geological work, which is large in volume and cannot be moved. The section is cleaned by natural or artificial methods to make it clear that the boundary relationship between gravel and filling materials can be clearly seen, and the profile image of the section is obtained by real-life three-dimensional mode and close-up photography; conglomerate cutting profile photography refers to the field collection of typical conglomerate samples to be identified, marking the spatial position and occurrence of the conglomerate during field sampling, cutting the conglomerate indoors by a cutting machine to form a cutting profile, and obtaining the conglomerate cutting profile image of the conglomerate to be identified by photography; borehole core cutting profile photography refers to obtaining gravel cores during drilling work, and obtaining borehole core cutting profile photography by cutting the cores. Microscopic rock section photography refers to the process of identifying sedimentary rocks such as limestone, sandstone, and gravel limestone, in which sections are obtained by cutting rocks, cutting cores, etc. Since sandstone or gravel has a small diameter, its microscopic rock section photography is obtained by microscopic means.

[0072] Profile photography develops different digital imaging or image acquisition methods for different application environments, including: 1) For larger-scale geological profiles, the initial sedimentary rock sample image is obtained by close-range photography using a digital camera; 2) For geological drilling cores, a cutter is used to cut the rock to be inspected, and the initial sedimentary rock sample image is obtained through a scanner or digital camera; 3) For typical sedimentary rock samples, such as psammitic limestone and conglomerate limestone, a cutter is used to cut the limestone to obtain a typical cross-section surface, and the initial sedimentary rock sample image is obtained through a digital camera or a digital microscope.

[0073] Step S12, preprocessing the initial sedimentary rock sample image to obtain a target image.

[0074] Specifically, the target image is a typical digital photographic image; a typical digital photographic image refers to an image with a resolution greater than 300 DPI, and can clearly distinguish the boundaries between gravel, sand, gravel and other sedimentary rock samples and interstitial objects. The typical digital photographic image can be obtained by processing the digital enhancement method of the boundary between gravel and interstitial objects. The typical digital photographic image is rotated, mirrored, etc. to make the target image correspond to the upper, lower, left and right relationship of the sedimentary rock sample section.

[0075] Preprocessing methods include image binarization, sharpening and contrast adjustment, which are used to improve the boundary between gravel and matrix in sedimentary rocks, facilitating subsequent identification.

[0076] Step S2, determining the vector boundary of each gravel in each sampling area in the target image.

[0077] Specifically, typical gravels L{1, 2, 3, n|n>40} are randomly selected in the random range Q(1), Q(2), Q(3)...Q(n), and the number of typical gravels selected is not less than 40. The vector boundaries of the gravels are obtained manually or automatically using GIS software, or the vector boundaries of the gravels are obtained by visual method. The schematic diagram of the vector boundaries of the gravels is shown in Figure 4 shown.

[0078] Step S3: taking any gravel as the current gravel.

[0079] Step S4, based on the vector boundary of the current gravel, determine the center of gravity of the current gravel section. Assuming that the density of the gravel is approximately uniform, the center of gravity of the section Z{...} is the center of gravity of the boundary range of the two-dimensional gravel image vector. The center of gravity of the gravel section is as follows: Figure 4 shown.

[0080] Step S5, taking the center of gravity of the cross section of the current gravel as the center of the circle and the first preset angle as the step length, determining multiple angle bisectors of the current gravel, and determining the intersection of each angle bisector and the vector boundary of the current gravel.

[0081] Specifically, the angle bisector is based on the coordinate system of the target image, with the 0° direction as the starting point, the center of gravity of the section as the center of the circle, and the rays are used to divide the full circle where the gravel is located into equal angles. The equal angle step size is set to 5°, which has a better effect; the rays intersect with the vector boundary of the current gravel to form an intersection point set P{1, 2, 3...n}, and the schematic diagram of multiple angle bisectors of the current gravel and the intersection points of each angle bisector and the vector boundary of the current gravel (i.e., the equal angle radius change diagram) is as follows: Figure 5 shown.

[0082] Step S6, determining the gravity direction of the current gravel based on all angle bisectors and all intersections of the current gravel; wherein the gravity direction is the direction pointing to the bottom of the rock layer of the sedimentary rock sample.

[0083] As an optional implementation, step S6 specifically includes:

[0084] Step S61, based on all angle bisectors and all intersection points of the current gravel, a gravity direction analysis and judgment diagram of the current gravel is drawn.

[0085] The gravity direction analysis and judgment diagram of the current gravel is drawn with the angle as the horizontal axis, the first preset angle as the step length, and the distance between the center of the circle and the angle bisector of the current gravel to the corresponding intersection point as the vertical axis. The gravity direction analysis and judgment diagram is drawn one by one in the sampling areas Q(1), Q(2), Q(3)...Q(n); the schematic diagram of the gravity direction analysis and judgment diagram of the drawn gravel is as follows: Figure 6-Figure 8 shown.

[0086] Step S62, based on the current gravity direction analysis and judgment diagram of the gravel, determine the current gravity direction of the gravel.

[0087] Specifically, the gravity advantage direction is the natural tendency of the center of gravity area to be stable under the influence of flowing water and gravity during the gravel deposition and diagenesis process under natural sedimentation conditions. The gravity direction and horizontal direction of the gravel are determined according to the gravity direction analysis and judgment diagram. Figure 6-Figure 8 As shown in the figure, the gravity direction analysis and judgment diagram is generally in a typical "W" shape. In general, the lower the center of gravity is, the more stable it is. The two points A and B with longer radii represent the direction of water flow, and the directions C and D with shorter radii represent the direction parallel to gravity. The gravity advantage direction is selected from the shorter radii C and D, and the shorter radius represents the gravity direction ZL{...}. Therefore, based on the current gravel gravity direction analysis and judgment diagram, the gravity advantage direction of the gravel can be judged, the gravity advantage direction of each gravel can be judged and recorded, and included in ZL{}.

[0088] like Figure 6-Figure 8 As shown in the gravity direction analysis and judgment diagram, there are two high-value directions, such as Figure 6-Figure 8 Among points A and B, points A and B are far away from the center of gravity, and the AB line represents the horizontal direction; points C and D are close to the center of gravity, and the CD line represents parallel to the direction of gravity; under natural conditions, the lower the center of gravity of gravel, the more stable the deposition, so continue to judge the size of C and D, and identify the value of points C and D that is obviously smaller as the direction of gravity.

[0089] Specifically, Figure 6 In gravel 1, A and B values ​​are the largest and approximately symmetrical, representing the horizontal direction; C and D are short and approximately symmetrical, representing the parallel gravity direction, and D is obviously smaller than C, so D represents the gravity direction of gravel 1.

[0090] In gravel 2, the maximum values ​​of A and B represent the horizontal direction, the short positions of C and D represent the parallel gravity direction, and D is obviously smaller than C, so D represents the gravity direction of gravel 2.

[0091] In gravel 3, the maximum A and B values ​​represent the horizontal direction, the short positions of C and D represent the parallel gravity direction, and C and D are basically equal. Therefore, the gravity advantage of gravel 3 is not obvious, and its gravity direction cannot be determined, so it is eliminated.

[0092] Step S7, determining the top and bottom relationship of the sequence of the sedimentary rock samples in the area to be identified based on the gravity directions of all gravels in the area to be identified.

[0093] As an optional implementation, step S7 specifically includes:

[0094] Step S71: based on the gravity directions of all gravels in each sampling area, a rose diagram of gravity dominance pointing to the corresponding sampling area is drawn. The rose diagram of gravity dominance pointing to the corresponding sampling area is shown in FIG. Fig. 9 As shown. The gravel gravity directions are obtained within the range of Q(1), Q(2), Q(3)...Q(n), and a gravity advantage direction rose diagram is drawn. The gravity advantage direction rose diagram is drawn using the gravel gravity direction recorded by ZL{}. The gravity advantage direction is the direction in which the gravity direction is particularly obvious in the gravity advantage direction rose diagram; the gravity advantage direction is based on a comprehensive analysis of the sedimentary environment, excluding some rapid sedimentation environments such as glaciers and debris flows where the gravity effect is not obvious. The gravity advantage direction rose diagram is drawn with equal angle spacing-frequency as the horizontal axis and vertical axis.

[0095] Step S72, based on the gravity advantage pointing rose diagram of each sampling area, determine the gravity direction of the corresponding sampling area.

[0096] Step S73, based on the gravity directions of all sampling areas, determine the top and bottom relationships of the sedimentary rock samples in the area to be identified.

[0097] The gravity direction must be the same or similar in the sampling area of ​​Q(1), Q(2), Q(3)...Q(n), and the difference between the gravity directions should not exceed 30°. The average value of the gravity direction angles is taken as the final gravity direction. The gravity dominant direction is the gravity direction in the rose diagram where the gravity is more concentrated or has a significant advantage. If the gravity dominant direction is not obvious, re-slice and re-identify.

[0098] As an optional implementation, step S71 specifically includes:

[0099] Step S711, counting the frequency of occurrence of different gravity directions of gravels in each sampling area.

[0100] Step S712: based on the frequency of occurrence of different gravity directions of gravels in each sampling area, a rose diagram of gravity dominance direction of the corresponding sampling area is drawn.

[0101] As an optional implementation, step S73 specifically includes:

[0102] Step S731, based on the gravity directions of all sampling areas, draw a gravity direction point map of the area to be identified; the gravity direction point map includes multiple gravity direction points; one gravity direction point corresponds to the gravity direction of one sampling area.

[0103] Step S732: taking any two gravity direction points as the point groups to be screened, and obtaining a plurality of point groups to be screened.

[0104] Step S733: determine any point group to be screened as the current point group, and determine the number of gravity direction points whose gravity directions are between the gravity directions corresponding to two gravity direction points in the current point group as the current number.

[0105] Step S734, determine whether the current point group meets the preset conditions; the preset conditions are that the difference in gravity direction of the sampling area corresponding to two gravity direction points in the current point group does not exceed 30 degrees, and the current number is greater than all non-current numbers; the non-current number is the number of gravity direction points whose gravity directions are between the gravity directions corresponding to two gravity direction points in the non-current point group; the non-current point group is all the point groups to be screened except the current point group.

[0106] Step S735, if yes, then average the gravity directions of the two gravity direction points in the current point group and the gravity direction points whose gravity directions are between the gravity directions corresponding to the two gravity direction points in the current point group to obtain the top and bottom relationship of the stratigraphic sequence of the sedimentary rock samples in the area to be identified.

[0107] Step S736: if not, reselect the sedimentary rock sample slice and return to step S1; or draw a conclusion: the top-bottom relationship is unclear.

[0108] like Fig. 9 As shown, point H and point F are located in the symmetric direction, and point E and point G are located in the symmetric direction. The values ​​corresponding to point H and point F are similar, the values ​​corresponding to point E and point G are quite different, and the direction of point G has obvious advantages. The direction of point G is identified as the direction of gravity advantage, and the angle corresponding to G is recorded; when the values ​​of point H and point F are not similar and the difference between the numbers of point E and point G is not large, the gravity advantage is not obvious, and re-sectioning is required for identification; when the difference between the numbers of point E and point G is small, the gravity advantage is not obvious.

[0109] Based on the same inventive concept, the embodiment of the present application also provides a sedimentary rock sequence top and bottom identification system based on gravity analysis for implementing the sedimentary rock sequence top and bottom identification method based on gravity analysis involved above. The implementation scheme for solving the problem provided by the system is similar to the implementation scheme recorded in the above method, so the specific limitations in one or more embodiments of the sedimentary rock sequence top and bottom identification system based on gravity analysis provided below can refer to the limitations of the sedimentary rock sequence top and bottom identification method based on gravity analysis above, and will not be repeated here.

[0110] In an exemplary embodiment, Fig.10 As shown, a system for identifying the top and bottom of a sedimentary rock sequence based on centroid analysis is provided, including:

[0111] The target image acquisition unit 1 is used to acquire a target image; the target image is an image of a sedimentary rock sample in an area to be identified; the area to be identified includes a plurality of sampling areas; each sampling area includes a plurality of gravels.

[0112] The vector boundary determination unit 2 is used to determine the vector boundary of each gravel in each sampling area in the target image.

[0113] The current gravel determining unit 3 is used to take any gravel as the current gravel.

[0114] The section center of gravity determining unit 4 is used to determine the section center of gravity of the current gravel based on the vector boundary of the current gravel.

[0115] The angle bisector and intersection determination unit 5 is used to determine multiple angle bisectors of the current gravel with the center of gravity of the cross section of the current gravel as the center of the circle and the first preset angle as the step length, and determine the intersection of each angle bisector and the vector boundary of the current gravel.

[0116] The gravel gravity direction determination unit 6 is used to determine the gravity direction of the current gravel based on all angle bisectors and all intersections of the current gravel; wherein the gravity direction is the direction pointing to the bottom of the rock layer of the sedimentary rock sample.

[0117] The stratigraphic top-bottom relationship determining unit 7 is used to determine the stratigraphic top-bottom relationship of the sedimentary rock samples in the area to be identified based on the gravity directions of all gravels in the area to be identified.

[0118] In an exemplary embodiment, a computer device is provided, including: a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the computer program to implement a method for identifying the top and bottom of a sedimentary rock sequence based on centroid analysis.

[0119] In an exemplary embodiment, a computer device is provided. The computer device may be a server or a terminal. The internal structure diagram thereof may be as follows: Fig.11 As shown. The computer device includes a processor, a memory, an input / output interface (Input / Output, referred to as I / O) and a communication interface. The processor, the memory and the input / output interface are connected through a system bus, and the communication interface is connected to the system bus through the input / output interface. The processor of the computer device is used to provide computing and control capabilities. The memory of the computer device includes a non-volatile storage medium and an internal memory. The non-volatile storage medium stores an operating system, a computer program and a database. The internal memory provides an environment for the operation of the operating system and the computer program in the non-volatile storage medium. The input / output interface of the computer device is used to exchange information between the processor and an external device. The communication interface of the computer device is used to communicate with an external terminal through a network connection. When the computer program is executed by the processor, a method for identifying the top and bottom of a sedimentary rock sequence based on centroid analysis is implemented.

[0120] Those skilled in the art will understand that Fig.11The structure shown in the figure is only a block diagram of a part of the structure related to the solution of the present application, and does not constitute a limitation on the computer device to which the solution of the present application is applied. The specific computer device may include more or fewer components than those shown in the figure, or combine certain components, or have a different arrangement of components.

[0121] It should be noted that the user information (including but not limited to user device information, user personal information, etc.) and data (including but not limited to data used for analysis, stored data, displayed data, etc.) involved in this application are all information and data authorized by the user or fully authorized by all parties, and the collection, use and processing of relevant data must comply with relevant regulations.

[0122] Those skilled in the art can understand that all or part of the processes in the above-mentioned embodiment methods can be completed by instructing the relevant hardware through a computer program, and the computer program can be stored in a non-volatile computer-readable storage medium. When the computer program is executed, it can include the processes of the embodiments of the above-mentioned methods. Among them, any reference to the memory, database or other medium used in the embodiments provided in the present application can include at least one of non-volatile and volatile memory. Non-volatile memory can include read-only memory (ROM), magnetic tape, floppy disk, flash memory, optical memory, high-density embedded non-volatile memory, resistive random access memory (ReRAM), magnetic random access memory (MRAM), ferroelectric random access memory (FRAM), phase change memory (PCM), graphene memory, etc. Volatile memory can include random access memory (RAM) or external cache memory, etc. By way of illustration and not limitation, RAM may be in various forms, such as static random access memory (SRAM) or dynamic random access memory (DRAM).

[0123] The database involved in each embodiment provided in this application may include at least one of a relational database and a non-relational database. The non-relational database may include a distributed database based on blockchain, etc., but is not limited thereto. The processor involved in each embodiment provided in this application may be a general-purpose processor, a central processing unit, a graphics processor, a digital signal processor, a programmable logic device, a data processing logic device based on quantum computing, etc., but is not limited thereto.

[0124] The technical features of the above embodiments may be combined arbitrarily. To make the description concise, not all possible combinations of the technical features in the above embodiments are described. However, as long as there is no contradiction in the combination of these technical features, they should be considered to be within the scope of this specification.

[0125] This article uses specific examples to illustrate the principles and implementation methods of this application. The description of the above embodiments is only used to help understand the method and core ideas of this application. At the same time, for those skilled in the art, according to the ideas of this application, there will be changes in the specific implementation methods and application scope. In summary, the content of this specification should not be understood as limiting this application.

Claims

1. A method for identifying the top and bottom of a sedimentary rock sequence based on centroid analysis, characterized in that: The method for identifying the top and bottom of a sedimentary rock sequence based on centroid analysis includes: Acquire a target image; the target image is an image of a sedimentary rock sample in the area to be identified; the area to be identified includes a plurality of sampling areas; each sampling area includes a plurality of gravels; Determine the vector boundary of each gravel in each sampling area in the target image; Set any gravel as the current gravel; Based on the vector boundary of the current gravel, determine the center of gravity of the current gravel section; Taking the center of gravity of the cross section of the current gravel as the center of the circle and the first preset angle as the step length, determining multiple angle bisectors of the current gravel, and determining the intersection of each angle bisector and the vector boundary of the current gravel; Based on all angle bisectors and all intersections of the current gravel, determine the gravity direction of the current gravel; wherein the gravity direction is the direction pointing to the bottom of the rock layer of the sedimentary rock sample; Based on the gravity directions of all gravels in the area to be identified, the top and bottom relationships of the sequences of sedimentary rock samples in the area to be identified are determined.

2. The method for identifying the top and bottom of a sedimentary rock sequence based on centroid analysis according to claim 1, characterized in that: Get the target image, including: Carry out cross-sectional photography of sedimentary rock samples in the area to be identified to obtain an initial sedimentary rock sample image; The initial sedimentary rock sample image is preprocessed to obtain the target image.

3. The method for identifying the top and bottom of a sedimentary rock sequence based on centroid analysis according to claim 2, characterized in that: Preprocessing methods include image binarization, sharpening, and contrast adjustment.

4. The method for identifying the top and bottom of a sedimentary rock sequence based on centroid analysis according to claim 1, characterized in that: Based on all angle bisectors and all intersections of the current gravel, determine the gravity direction of the current gravel, including: Based on all angle bisectors and all intersections of the current gravel, draw a gravity direction analysis and judgment diagram of the current gravel; Based on the current gravel gravity direction analysis and judgment diagram, the current gravel gravity direction is determined.

5. The method for identifying the top and bottom of a sedimentary rock sequence based on centroid analysis according to claim 4, characterized in that: Based on all the angle bisectors and all the intersections of the current gravel, draw a gravity direction analysis and judgment diagram of the current gravel, including: With the angle as the horizontal axis, the first preset angle as the step length, and the distance between the center of the circle along the angle bisector of the current gravel and the corresponding intersection point as the vertical axis, a gravity direction analysis and judgment diagram of the current gravel is drawn.

6. The method for identifying the top and bottom of a sedimentary rock sequence based on centroid analysis according to claim 1, characterized in that: Based on the gravity direction of all gravels in the area to be identified, determine the top and bottom relationship of the sequence of sedimentary rock samples in the area to be identified, including: Based on the gravity direction of all gravels in each sampling area, a rosette diagram of gravity advantage direction of the corresponding sampling area is drawn; Based on the gravity advantage pointing rose diagram of each sampling area, the gravity direction of the corresponding sampling area is determined; Based on the gravity direction of all sampling areas, the top and bottom relationships of the sedimentary rock samples in the area to be identified are determined.

7. The method for identifying the top and bottom of a sedimentary rock sequence based on centroid analysis according to claim 6, characterized in that: Based on the gravity direction of all gravels in each sampling area, a rose diagram of the gravity advantage direction of the corresponding sampling area is drawn, including: Count the frequencies of different gravity directions of gravels in each sampling area; Based on the frequency of occurrence of different gravity directions of gravel in each sampling area, a rose diagram of the gravity advantage direction of the corresponding sampling area was drawn.

8. The method for identifying the top and bottom of a sedimentary rock sequence based on centroid analysis according to claim 6, characterized in that: Based on the gravity direction of all sampling areas, determine the top and bottom relationship of the sequence of sedimentary rock samples in the area to be identified, including: Based on the gravity directions of all sampling areas, a gravity direction point map of the area to be identified is drawn; the gravity direction point map includes multiple gravity direction points; one gravity direction point corresponds to the gravity direction of one sampling area; Taking any two gravity direction points as the point groups to be screened, multiple point groups to be screened are obtained; Determine any point group to be screened as the current point group, and determine the number of gravity direction points whose gravity directions are between the gravity directions corresponding to two gravity direction points in the current point group as the current number; Determine whether the current point group meets the preset conditions; the preset conditions are that the difference in gravity direction of the sampling area corresponding to the two gravity direction points in the current point group does not exceed 30 degrees, and the current number is greater than all non-current numbers; the non-current number is the number of gravity direction points whose gravity directions are between the gravity directions corresponding to the two gravity direction points in the non-current point group; the non-current point group is all the point groups to be screened except the current point group in all the point groups to be screened; If so, then the gravity directions of the two gravity direction points in the current point group and the gravity direction points whose gravity directions are between the gravity directions corresponding to the two gravity direction points in the current point group are averaged to obtain the top and bottom relationship of the sequence of the sedimentary rock samples in the area to be identified.

9. A system for identifying the top and bottom of a sedimentary rock sequence based on centroid analysis, characterized in that: The top and bottom identification system of a sedimentary rock sequence based on centroid analysis is used to implement the top and bottom identification method of a sedimentary rock sequence based on centroid analysis according to any one of claims 1 to 8, and the top and bottom identification system of a sedimentary rock sequence based on centroid analysis comprises: The target image acquisition unit is used to acquire a target image; the target image is an image of a sedimentary rock sample in the area to be identified; the area to be identified includes a plurality of sampling areas; each sampling area includes a plurality of gravels; A vector boundary determination unit, used for determining the vector boundary of each gravel in each sampling area in the target image; A current gravel determination unit, used for taking any gravel as the current gravel; A section center determination unit, used for determining the section center of the current gravel based on the vector boundary of the current gravel; An angle bisector and intersection determination unit, used to determine multiple angle bisectors of the current gravel with the center of gravity of the section of the current gravel as the center of the circle and the first preset angle as the step length, and determine the intersection of each angle bisector and the vector boundary of the current gravel; The unit for determining the gravity direction of the gravel is used to determine the gravity direction of the current gravel based on all angle bisectors and all intersections of the current gravel; wherein the gravity direction is the direction pointing to the bottom of the rock layer of the sedimentary rock sample; The stratigraphic top-bottom relationship determination unit is used to determine the stratigraphic top-bottom relationship of the sedimentary rock samples in the area to be identified based on the gravity directions of all gravels in the area to be identified.

10. A computer device comprising: A memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the computer program to implement the method for identifying the top and bottom of a sedimentary rock sequence based on centroid analysis as described in any one of claims 1 to 8.

Images