A method, apparatus, computer device and medium for generating a formation boundary surface

Abstract

Embodiments of the present application provide a method and device for generating a stratigraphic interface, computer equipment and a medium. The method interprets geological numbers, and partitions known points obtained by interpretation and to-be-interpolated points required for generating a stratigraphic interface, so that the known point blocks and the to-be-interpolated point blocks correspond to each other, converts global search in the interpolation process into local search, and then performs spatial interpolation on a block-by-block basis, thereby improving program running efficiency.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of strata, and particularly relates to a method and device for generating a strata interface, a computer device and a medium. BACKGROUND

[0002] In the field of seismic exploration, seismic interpretation data is an important basis for understanding the structure and distribution of underground strata. These data are usually distributed along the survey line and have the characteristics of large data volume and anisotropic data distribution. Each survey line carries rich geological information, which is crucial for subsequent geological exploration and development.

[0003] In actual engineering applications, in order to more intuitively understand the distribution of underground strata, it is often necessary to generate a strata interface based on seismic interpretation data. This process usually relies on spatial interpolation technology, which can predict the strata attribute of unknown points according to known seismic interpretation data points through a certain algorithm. However, directly importing the original seismic interpretation data into the interpolation algorithm for calculation often causes problems such as low program running efficiency and running out of memory.

[0004] For example, the seismic interpretation data is imported into the program, a number of interpolation points in the engineering site are selected, and distance inverse interpolation is used. In the interpolation process, it is often necessary to set an interpolation influence radius, that is, before predicting each interpolation point, a part of the known points within the interpolation influence radius needs to be searched from the entire engineering site interpretation. However, the data volume of the seismic interpretation data is large, and in this case, the program runs very time-consuming. SUMMARY

[0005] The embodiments of the present application provide a method and device for generating a strata interface, a computer device and a medium, aiming to solve the problem of low program running efficiency caused by directly importing the original seismic interpretation data into the interpolation algorithm for calculation in the existing interpolation technology.

[0006] In a first aspect, the embodiments of the present application provide a method for generating a strata interface using seismic interpretation data, comprising:

[0007] Step 1: obtaining geological data obtained by geological exploration of an engineering area through a survey line; wherein the geological exploration includes seismic exploration;

[0008] Step 2: comprehensively interpreting the geological data to obtain a series of data points of the strata interface distributed along the survey line;

[0009] Step 3: delineating the engineering area to obtain a delineated range of the engineering area, and taking a minimum rectangle containing the delineated range along the boundary of the delineated range to obtain a rectangular range;

[0010] Step 4: dividing the rectangular range to obtain a plurality of grids, taking each node of the grid as an interpolation point of the stratigraphic interface;

[0011] Step 5: selecting any data point of the stratigraphic interface as a known point, and dividing all the interpolation points and the known points by using a spatial search strategy to obtain a plurality of corresponding interpolation point blocks and known point blocks;

[0012] Step 6: interpolating the corresponding interpolation point block by using each known point block to obtain the current stratigraphic interface;

[0013] Step 7: repeating steps 5-6 to obtain all the stratigraphic interfaces in the engineering area.

[0014] In a second aspect, an embodiment of the present application provides a device for generating a stratigraphic interface by using seismic interpretation data, comprising:

[0015] a data acquisition unit configured to acquire geological data obtained by geological exploration of an engineering area by a survey line; wherein the geological exploration comprises seismic exploration;

[0016] an interpretation unit configured to comprehensively interpret the geological data to obtain a series of data points of stratigraphic interfaces distributed along the survey line;

[0017] a delineation unit configured to delineate the engineering area to obtain a delineated range of the engineering area, and take a minimum rectangle containing the delineated range along the boundary of the delineated range to obtain a rectangular range;

[0018] a grid division unit configured to divide the rectangular range to obtain a plurality of grids, and take each node of the grid as an interpolation point of the stratigraphic interface;

[0019] a block division unit configured to select any data point of the stratigraphic interface as a known point, and divide all the interpolation points and the known points by using a spatial search strategy to obtain a plurality of corresponding interpolation point blocks and known point blocks;

[0020] an interpolation unit configured to interpolate each interpolation point block and known point block to obtain the current stratigraphic interface;

[0021] a repeating unit configured to repeat the block division unit and the interpolation unit to obtain all the stratigraphic interfaces in the engineering area.

[0022] In a third aspect, an embodiment of the present application provides a computer device, comprising a memory, a processor, and a computer program stored in the memory and capable of running on the processor, and the processor implements the method for generating a stratigraphic interface by using seismic interpretation data as described above when executing the computer program.

[0023] In a fourth aspect, an embodiment of the present application provides a computer readable storage medium, wherein the computer readable storage medium stores a computer program, and the computer program, when executed by a processor, implements the method for generating a stratigraphic interface by using seismic interpretation data.

[0024] The method for generating a stratigraphic interface, the device, the computer device and the medium provided by the embodiment of the present application interpret geological numbers, and divide known points obtained by interpretation and to-be-interpolated points required for generating a stratigraphic interface into partitions, so that the known point blocks and the to-be-interpolated point blocks correspond to each other, convert global search in the interpolation process into local search, and then perform spatial interpolation on the partitions one by one, thereby improving program running efficiency. BRIEF DESCRIPTION OF DRAWINGS

[0025] In order to more clearly illustrate the technical solutions of the embodiments of the present application, the drawings required in the embodiment description will be briefly introduced. Obviously, the drawings in the following description are some embodiments of the present application, and other drawings can be obtained by those skilled in the art without creative labor.

[0026] Figure 1 A flowchart of a method for generating a stratigraphic interface by using seismic interpretation data provided by the embodiment of the present application;

[0027] Figure 2 A sub-flowchart of a method for generating a stratigraphic interface by using seismic interpretation data provided by the embodiment of the present application;

[0028] Figure 3 Another sub-flowchart of a method for generating a stratigraphic interface by using seismic interpretation data provided by the embodiment of the present application;

[0029] Figure 4 Another sub-flowchart of a method for generating a stratigraphic interface by using seismic interpretation data provided by the embodiment of the present application;

[0030] Figure 5 Another sub-flowchart of a method for generating a stratigraphic interface by using seismic interpretation data provided by the embodiment of the present application;

[0031] Figure 6 Another sub-flowchart of a method for generating a stratigraphic interface by using seismic interpretation data provided by the embodiment of the present application;

[0032] Figure 7 Another sub-flowchart of a method for generating a stratigraphic interface by using seismic interpretation data provided by the embodiment of the present application;

[0033] Figure 8 A schematic block diagram of a device for generating a stratigraphic interface using seismic interpretation data is provided for an embodiment of the present application;

[0034] Figure 9 A schematic diagram of a seismic interpretation data survey line and a delineated range of an engineering area;

[0035] Figure 10 A schematic diagram of a delineated range forming a rectangular range;

[0036] Figure 11 A schematic diagram of generating points to be interpolated according to the rectangular range;

[0037] Figure 12 A schematic diagram of dividing blocks of known points according to the rectangular range;

[0038] Figure 13 A schematic diagram of a one-to-one correspondence between the block division of known points and the block division of points to be interpolated;

[0039] Figure 14 A schematic diagram of establishing an overlap block between two blocks of known points and an overlap block of points to be interpolated;

[0040] Figure 15 A complete schematic diagram of a one-to-one correspondence between the overlap block of two blocks of known points and the overlap block of points to be interpolated;

[0041] Figure 16 A schematic diagram of establishing an overlap block between four blocks of known points and an overlap block of points to be interpolated;

[0042] Figure 17 A schematic diagram of a region not covered by the overlap block or the block of points to be interpolated;

[0043] Figure 18 A complete schematic diagram of all block divisions and the establishment of overlap blocks after the establishment of overlap blocks;

[0044] Figure 19 A schematic diagram of an actual engineering survey line, a wind field range, and a rectangular range;

[0045] Figure 20 A schematic diagram of the division of blocks of known points in an actual engineering;

[0046] Figure 21 A complete schematic diagram of all block divisions and the establishment of overlap blocks after the establishment of overlap blocks in an actual engineering;

[0047] Figure 22 A display effect diagram of a generated stratigraphic interface. DETAILED DESCRIPTION

[0048] With reference to the drawings of the embodiments of the present application, the technical solutions in the embodiments of the present application will be clearly and completely described, obviously, the described embodiments are a part of the embodiments of the present application, rather than all the embodiments. Based on the embodiments of the present application, all other embodiments obtained by a person of ordinary skill in the art without creative work belong to the protection scope of the present application.

[0049] It should be understood that the terms "comprising" and "including" as used in the specification and the appended claims indicate the presence of the described features, integers, steps, operations, elements, and / or components, but do not preclude the presence or addition of one or more other features, integers, steps, operations, elements, components, and / or groups thereof.

[0050] It should also be understood that the terms used in the present application specification are only for the purpose of describing specific embodiments and are not intended to limit the present application. As used in the present application specification and the appended claims, the singular forms "a", "an" and "the" are intended to include the plural forms, unless the context clearly indicates otherwise.

[0051] It should be further understood that the term "and / or" as used in the present application specification and the appended claims means any combination of one or more of the associated listed items and all possible combinations, and includes these combinations.

[0052] Please refer to Figure 1 The embodiments of the present application provide a method for generating a stratigraphic boundary surface by using seismic interpretation data, comprising S10-S70:

[0053] S10, obtaining geological data obtained by geological exploration of an engineering area through a survey line; wherein the geological exploration comprises seismic exploration;

[0054] In this step, a series of survey lines need to be arranged in the engineering area, the survey lines can be straight lines, broken lines or curves, and the specific shape and arrangement mode depend on the geological conditions, exploration target and exploration accuracy requirement, and then geological exploration is carried out along the survey lines to obtain geological data. The geological exploration includes seismic exploration, which is a kind of geophysical exploration method, which uses the characteristics of artificial excited seismic wave propagation in underground medium to infer the properties and morphology of underground rock stratum. In the process of seismic exploration, artificial seismic sources such as explosive explosion and controllable seismic source vibration are used to excite seismic wave, and then precise instruments are used to record the propagation of seismic wave on the ground and underground. Through the seismic exploration instrument, the data of the propagation of seismic wave in the underground medium are collected, including the amplitude, frequency, phase and other parameters of the seismic wave. These parameters are used as seismic exploration data for subsequent data processing and analysis.

[0055] In a specific implementation, the geological exploration further includes borehole exploration or static force exploration. The borehole exploration is to obtain a direct sample of an underground rock layer through drilling. The static force exploration (CPT) is an in-situ testing technology. A probe is pressed into the soil at a certain rate by using a pressure device. The soil layer properties are evaluated by measuring various data of the probe in the penetration process. In the implementation process, the seismic exploration data and the borehole data (or CPT data) are required as the geological data.

[0056] S20, comprehensively interpreting the geological data to obtain a series of data points of a stratigraphic interface distributed along the survey line;

[0057] As can be seen from the above steps, the geological data includes the seismic exploration data and the borehole data (or CPT data). For these data, interpretation is required. For the interpretation of the seismic exploration data, different strata are identified by observing the reflection waveform, amplitude and continuity and other characteristics in the seismic profile. The depth and thickness of each stratum are calculated according to the trend and velocity information of the reflection wave. The physical properties of the stratum are analyzed through the propagation velocity of the seismic wave. For the interpretation of the borehole data (or CPT data), the type, composition and structure of the stratum are analyzed through the core description and physical and mechanical property test data of the borehole sample (or the penetration resistance and side wall friction resistance of the stratum are analyzed by using the CPT test data, and the strength and deformation characteristics of the stratum are further understood). According to the analysis results of the borehole data / CPT data, the position of the stratigraphic interface is determined. The interpretation results of the seismic exploration data are compared and verified. Specifically, the seismic exploration data and the borehole data (or CPT data) are matched and aligned in space and time to ensure the consistency of the two. The specific position of the stratigraphic interface is determined by combining the interpretation results of the seismic exploration data and the borehole data / CPT data, so as to obtain a series of data points of the stratigraphic interface distributed along the survey line.

[0058] S30, delineating the engineering area to obtain a delineated range of the engineering area, and taking a minimum rectangle containing the delineated range along the boundary of the delineated range to obtain a rectangular range;

[0059] In this step, the engineering area is first delineated so as to have a clear delineated range. In a specific implementation, the delineated range can be irregular or regular. Therefore, the delineated range can be simplified as a minimum circumscribed rectangle of the delineated range, so as to obtain the rectangular range. The rectangular range can improve the processing efficiency in the subsequent process.

[0060] In an embodiment, as shown in FIG. 3, S30 includes: Figure 2

[0061] S31, inputting a plurality of vertex coordinates and connecting the vertex coordinates to obtain the delineated range of the engineering area; ​

[0062] S32, X values and Y values of all vertex coordinates in the range are taken out to form X set and Y set;

[0063] S33, maximum value and minimum value in the X set and the Y set are calculated respectively, and the maximum value of the X set is combined with the maximum value of the Y set, the maximum value of the X set is combined with the minimum value of the Y set, the minimum value of the X set is combined with the maximum value of the Y set, and the minimum value of the X set is combined with the minimum value of the Y set to obtain four coordinate points;

[0064] S34, the four coordinate points are connected to obtain a rectangular range.

[0065] In this embodiment, by inputting multiple vertex coordinates, the input can be clockwise or counterclockwise, and the vertex coordinates are sequentially connected to form a polygon, which is the range of the engineering area (see the yellow frame in Figure 9 The X values and Y values of the vertex coordinates are taken out to form {X} set: {x1, x2, …, x n} and {Y} set: {y1, y2, …, y n}, the minimum value x min and the maximum value x max are calculated from the {X} set, y min and y max are calculated from the {Y} set, and they are combined to obtain coordinate points Point1: (x min, y min ), Point 2: ( xmin ,y max ), Point 3: (x max ,y max ), and Point4: (x max, y min ), and the four coordinate points are connected to form a rectangular range (see the red frame in Figure 10 ).

[0066] S40, the rectangular range is divided to obtain multiple grids, and each grid node is taken as a to-be-interpolated point of the stratigraphic interface;

[0067] In this step, the rectangular range is divided into multiple grids, and after division, each grid node (i.e., grid point) is taken as a to-be-interpolated point of the stratigraphic interface, thereby improving the accuracy of subsequent interpolation.

[0068] In an embodiment, S40 includes:

[0069] The rectangular range is divided along the x direction according to a first predetermined length, and the rectangular range is divided along the y direction according to a second predetermined length to obtain multiple grids, and the nodes of all the grids are used as the points to be interpolated of the stratum interface.

[0070] In this embodiment, the first predetermined length and the second predetermined length can be the same length or different lengths. If they are the same length, the grid is a square grid, and the distance between each node and the adjacent node is the same. If they are different lengths, the grid is a rectangular grid, and the distance between each node and the adjacent node is different. In a specific implementation, in order to improve the accuracy of subsequent interpolation, it is preferred to set the first predetermined length (g x ) and a second predetermined length (g y ) are set to the same length, the last grid length along the x direction is g xf =(x max -x min )%g x , the length of the last grid segment along the y direction is g yf =(y max -y min )%g y , so that the distance between each interpolation point is the same and the distribution is more uniform. Among them, % means remainder. The range of the divided rectangle can be referred to Figure 11 .

[0071] S50, selecting any data point of the stratum interface as a known point, and dividing all the to-be-interpolated points and known points by using a spatial search strategy to obtain a plurality of corresponding to-be-interpolated point blocks and known point blocks;

[0072] In this step, data points distributed along the survey line on one of the stratigraphic interfaces are selected, and these data points are used as spatial interpolation to generate known points of this stratigraphic interface. These known points are used for subsequent interpolation of the points to be interpolated. The interpolation principle is: use the coordinates of the point to be interpolated as the center of the circle and r as the radius to divide a range to be interpolated, and use the known points within the range to interpolate the points to be interpolated. Before interpolation, it is necessary to use a spatial search strategy to partition the interpolation points and known points respectively, and obtain multiple known point blocks and to-be-interpolated point blocks, and realize a one-to-one correspondence between the two, that is, one known point block corresponds to one to-be-interpolated point block, so that the corresponding to-be-interpolated point block is interpolated through the known point block.

[0073] In one embodiment, if Figure 3 As shown, S50 includes:

[0074] S51, dividing the rectangular range along the x-direction by a third predetermined length and along the y-direction by a fourth predetermined length from the lower left corner thereof to obtain a plurality of known point blocks;

[0075] S52, each known point block is reduced by a first predetermined distance towards the center point corresponding to the known point block, to obtain a corresponding interpolation point block.

[0076] In this embodiment, the known point blocks are first divided, that is, by defining the length of the known point block in the X direction (third predetermined length l x ) and the length in the Y direction (fourth predetermined length l y ), starting from the lower left corner of the rectangular range, until the rectangular range is fully covered, thereby obtaining a plurality of known point blocks, and the length of the last block in the x direction is still l x , part of which exceeds the rectangular range, and the length of the last block in the y direction is still l y , part of which exceeds the rectangular range, and the known point block division schematic can be referred to Figure 12 . Then each divided known point block is reduced by a first predetermined distance towards the center point of the respective known point block, thereby ensuring that each known point block obtains a corresponding interpolation point block, that is, the number of known point blocks and interpolation point blocks is the same, and the interpolation point block division schematic can be referred to Figure 13 . The third predetermined length and the fourth predetermined length can be the same or different, and can be set according to the specific application scenario, and the first predetermined distance is less than the third predetermined length and the fourth predetermined length, which can be set according to the specific scene. The corresponding value.

[0077] In an embodiment, as shown in Figure 4 , S50 further comprises:

[0078] S53, respectively, all interpolation point blocks and known point blocks are established to establish overlapping blocks;

[0079] S54, obtaining all known point blocks not covered by the overlapping blocks or interpolation point blocks;

[0080] S55, the interpolation points are divided into the nearest overlapping blocks or interpolation point blocks.

[0081] In this embodiment, when interpolating the to-be-interpolated points by using the known points, the known points within a certain influence radius r (i.e., the to-be-interpolated range with the to-be-interpolated point as the center and the radius r) around the to-be-interpolated points are selected, and the to-be-interpolated points are interpolated by using the known points in the to-be-interpolated range. Since the to-be-interpolated point blocks are smaller than the known point blocks by a third predetermined distance in the up, down, left and right four directions when being divided, the plurality of to-be-interpolated point blocks that are divided cannot cover all the to-be-interpolated points. At this time, the overlapping regions between the known point blocks and the overlapping regions between the to-be-interpolated point blocks need to be established, and then the to-be-interpolated points that are not covered by the overlapping blocks or the to-be-interpolated point blocks in all the known point blocks are obtained. These to-be-interpolated points are located in the regions that are not covered by the overlapping blocks or the to-be-interpolated point blocks (see Figure 17 ), and then it is judged which overlapping block or to-be-interpolated point block the to-be-interpolated points are close to, and the to-be-interpolated points are divided into the overlapping block or to-be-interpolated point block closest to them (see Figure 18 ). At this time, all the to-be-interpolated points are preliminarily processed.

[0082] In an embodiment, as shown in Figure 5 , S53 comprises:

[0083] S531, the number of blocks of the known point blocks is selected;

[0084] S532, if the selected number of blocks is a first value, the first known point overlapping blocks are established according to a first establishment mode;

[0085] S533, if the selected number of blocks is a second value, the second known point overlapping blocks are established according to a second establishment mode;

[0086] S534, the first known point overlapping blocks are reduced by a second predetermined distance towards the center points corresponding to the first known point overlapping blocks, to obtain the first to-be-interpolated point overlapping blocks;

[0087] S535, the second known point overlapping blocks are reduced by a third predetermined distance towards the center points corresponding to the second known point overlapping blocks, to obtain the second to-be-interpolated point overlapping blocks.

[0088] In this embodiment, the overlapping blocks between the known point blocks are established first. According to the different numbers of blocks of the known point blocks, there are two establishment modes, so that two kinds of known point overlapping blocks are obtained, and then two kinds of to-be-interpolated point overlapping blocks are established according to the two kinds of known point overlapping blocks. The establishment mode is to reduce the corresponding predetermined distances r, i.e., the second predetermined distance and the third predetermined distance, towards the center points corresponding to the known point overlapping blocks, respectively. The two predetermined distances can be set to be the same or different.

[0089] In specific implementation, the first value is 2, and the first establishment mode is that

[0090] As shown in Figure 14 and Figure 15 , a rectangular overlap block is established according to the overlap block length and the overlap block height, and serves as a first known point overlap block, with the overlap block length being set along a direction perpendicular to a symmetry axis of the intersection line of two adjacent known point blocks, and the overlap block height being set along a direction parallel to the symmetry axis. The overlap block length can be set as 4r (r can be self-determined), and the overlap block height can be set as the length of the intersection line.

[0091] In a specific implementation, the second number is 4, and the second establishment manner is as follows:

[0092] As shown in Figure 16 , a square overlap block is established according to the overlap block length and the overlap block height, and serves as a second known point overlap block, with the overlap block length and the overlap block height being set along directions perpendicular to a symmetry axis of the intersection point of four adjacent known point blocks, with the intersection point as a center point. The overlap block length and the overlap block height can both be set as 4r (r can be self-determined).

[0093] S60, interpolating each known point block to a corresponding to-be-interpolated point block to obtain a current stratigraphic boundary surface;

[0094] In this step, each known point block is used to interpolate a corresponding to-be-interpolated point block one by one, so that all to-be-interpolated point blocks are interpolated, thereby obtaining the entire stratigraphic boundary surface (i.e., the current stratigraphic boundary surface). This interpolation manner can improve the program running efficiency.

[0095] In an embodiment, as shown in Figure 6 , S60 includes:

[0096] S61, establishing a corresponding relationship for all blocks in the following manner: each known point block corresponds to a to-be-interpolated point block, and each known point overlap block corresponds to a to-be-interpolated point overlap block.

[0097] S62, predicting all to-be-interpolated points of each to-be-interpolated point block one by one, and if the prediction is successful, performing spatial interpolation on all to-be-interpolated points.

[0098] In this embodiment, in order to better interpolate all the to-be-interpolated blocks, all the blocks are preferably associated, which can be established by numbering, such as that the known point block 1 corresponds to the to-be-interpolated point block 1, the known point overlap block 1 corresponds to the to-be-interpolated point overlap block 1, and according to the numbering order, the interpolation points of all the to-be-interpolated point blocks are predicted one by one to determine whether the interpolation point can be spatially interpolated, if the prediction is successful, the interpolation point is spatially interpolated, otherwise, it is reserved first and then processed subsequently.

[0099] In one embodiment, as shown in FIG. 6, Figure 7 S62 includes:

[0100] S621, any of the to-be-interpolated point blocks is selected as a current to-be-interpolated point block, and the first to-be-interpolated point in the current to-be-interpolated point block is selected as a current to-be-interpolated point, and it is determined whether the current to-be-interpolated point is within the circumscribed range, if not, the current to-be-interpolated point is skipped, if yes, it is determined that the prediction is successful, and the current to-be-interpolated point is spatially interpolated by the known points within the predetermined range, until all the to-be-interpolated points in the current block are predicted and the to-be-interpolated points that are successfully predicted are spatially interpolated;

[0101] S622, the memory of the current to-be-interpolated point block is released, and then all the to-be-interpolated points in the remaining to-be-interpolated point blocks are predicted one by one and the to-be-interpolated points that are successfully predicted are spatially interpolated.

[0102] In this embodiment, the known point block and the to-be-interpolated point block numbered 1 can be selected, and the first to-be-interpolated point in the to-be-interpolated point block is selected, and it is determined whether the to-be-interpolated point is within the circumscribed range of the engineering area, which can be determined by the ray method, if not, the to-be-interpolated point is skipped, if yes, the known points within the known point block 1 are searched by the interpolation influence radius r to obtain a to-be-interpolated range, and the to-be-interpolated point is spatially interpolated by the known points within the to-be-interpolated range. The above process is repeated until all the to-be-interpolated points in the to-be-interpolated point block 1 are predicted, and then the memory of the known points and the to-be-interpolated points in the block is released, and then the spatial interpolation of the blocks 2, 3, …, the overlap block 1, the overlap block 2, … is performed, until all the blocks are predicted and interpolated.

[0103] S70, S50-S60 are repeated to obtain all the stratigraphic boundaries in the engineering area.

[0104] The stratigraphic boundaries are generated in the manner of the foregoing steps S50-S60 until all the stratigraphic boundaries in the engineering area are generated.

[0105] The embodiment of the present application converts global search in the interpolation process into local search by interpreting geological data, partitioning known points and interpolation points needed for generating stratigraphic boundaries, and making known point blocks and interpolation point blocks one-to-one correspondence, and then performing spatial interpolation on each partition, and releasing the memory of the data points stored in the partition after the interpolation of each partition is completed, thereby improving the program running efficiency and avoiding the situation of exceeding the memory.

[0106] In actual application, as shown in Figures 19-21 , there will also be a wind field range in the engineering area, which will affect the distribution of the survey lines in the engineering area, so that the delineation range of the engineering area is an irregular range, and then the aforementioned steps are used to simplify the delineation range into a minimum circumscribed rectangle, so as to obtain a rectangular range, and then the rectangular range is divided to obtain a plurality of known point blocks, and then interpolation point division and various overlapping block establishment are performed, and the division method and the establishment method can refer to the aforementioned steps, which will not be repeated here, and then the stratigraphic boundary shown in Figure 22 is generated, wherein, Figure 22 Z represents the depth value of the interpolation point.

[0107] The embodiment of the present application also provides a device for generating a stratigraphic boundary by using seismic interpretation data, which is used to execute any one of the aforementioned embodiments of the method for generating a stratigraphic boundary by using seismic interpretation data. Specifically, please refer to Figure 8 , Figure 8 is a schematic block diagram of a device for generating a stratigraphic boundary by using seismic interpretation data provided by the embodiment of the present application. The device 800 for generating a stratigraphic boundary by using seismic interpretation data comprises:

[0108] A data acquisition unit 810 is configured to acquire geological data obtained by geological exploration of an engineering area by a survey line, wherein the geological exploration comprises seismic exploration;

[0109] An interpretation unit 820 is configured to comprehensively interpret the geological data to obtain a series of data points of stratigraphic boundaries distributed along the survey line;

[0110] A delineation unit 830 is configured to delineate the engineering area to obtain a delineation range of the engineering area, and obtain a rectangular range by taking a minimum rectangle containing the delineation range along the boundary of the delineation range;

[0111] A grid division unit 840 is configured to divide the rectangular range to obtain a plurality of grids, and take each node of the grid as an interpolation point of the stratigraphic boundary;

[0112] The block division unit 850 is configured to select a data point of any of the formation interfaces as a known point, divide all the to-be-interpolated points and the known points respectively by using a spatial search strategy, and obtain a plurality of corresponding to-be-interpolated point blocks and known point blocks.

[0113] The interpolation unit 860 is configured to perform interpolation processing on each of the to-be-interpolated point blocks and the known point blocks, and obtain a current formation interface.

[0114] The repeating unit 870 is configured to repeat the block division unit and the interpolation unit to obtain all the formation interfaces in the engineering area.

[0115] In an embodiment, the delineation unit 830 includes:

[0116] The input unit is configured to input a plurality of vertex coordinates, and connect the vertex coordinates to obtain a delineated range of the engineering area.

[0117] The value taking unit is configured to take X values and Y values of all the vertex coordinates in the delineated range to form an X set and a Y set.

[0118] The combination unit is configured to calculate maximum values and minimum values in the X set and the Y set respectively, and combine the maximum value of the X set with the maximum value of the Y set, the maximum value of the X set with the minimum value of the Y set, the minimum value of the X set with the maximum value of the Y set, and the minimum value of the X set with the minimum value of the Y set to obtain four coordinate points.

[0119] The connection unit is configured to connect the four coordinate points to obtain a rectangular range.

[0120] In an embodiment, the grid division unit 840 includes:

[0121] The first division unit is configured to divide the rectangular range by a first predetermined length in an x direction, and divide the rectangular range by a second predetermined length in a y direction to obtain a plurality of grids, and take nodes of all the grids as to-be-interpolated points of the formation interfaces.

[0122] In an embodiment, the block division unit 850 includes:

[0123] The second division unit is configured to divide the rectangular range by a third predetermined length in the x direction from a lower left corner of the rectangular range, and divide the rectangular range by a fourth predetermined length in the y direction to obtain a plurality of known point blocks.

[0124] The first reduction unit is configured to reduce each known point block by a first predetermined distance towards a center point corresponding to the known point block to obtain a corresponding to-be-interpolated point block.

[0125] In an embodiment, the block division unit 850 further comprises:

[0126] an overlapping block establishing unit, configured to establish overlapping blocks for all the to-be-interpolated point blocks and the known point blocks respectively;

[0127] a to-be-interpolated point obtaining unit, configured to obtain to-be-interpolated points in all the known point blocks which are not covered by the overlapping blocks or the to-be-interpolated point blocks;

[0128] a third division unit, configured to divide the to-be-interpolated points to the nearest overlapping block or to-be-interpolated point block.

[0129] In an embodiment, the overlapping block establishing unit comprises:

[0130] a number selecting unit, configured to select a number of blocks of the known point blocks;

[0131] a first establishing unit, configured to, if the selected number of blocks is a first value, establish first known point overlapping blocks for all the known point blocks according to a first establishing manner;

[0132] a second establishing unit, configured to, if the selected number of blocks is a second value, establish second known point overlapping blocks for all the known point blocks according to a second establishing manner;

[0133] a second reducing unit, configured to reduce the first known point overlapping blocks by a second predetermined distance towards a center point of the corresponding first known point overlapping block, to obtain first to-be-interpolated point overlapping blocks;

[0134] a third reducing unit, configured to reduce the second known point overlapping blocks by a third predetermined distance towards a center point of the corresponding second known point overlapping block, to obtain second to-be-interpolated point overlapping blocks.

[0135] In an embodiment, the first establishing unit comprises:

[0136] the first value is 2, the first establishing unit sets an overlapping block length along a direction perpendicular to a symmetry axis of a boundary line between two adjacent known point blocks, and sets an overlapping block height along a direction parallel to the symmetry axis, and establishes a rectangular overlapping block according to the overlapping block length and the overlapping block height, and takes the rectangular overlapping block as the first known point overlapping block.

[0137] In an embodiment, the second establishing unit comprises:

[0138] The second number is 4, and a square overlap block is established according to an overlap block length and an overlap block height, the center point of which is the intersection point of four adjacent known point blocks, the overlap block length is set along the direction perpendicular to the symmetry axis, and the overlap block height is set along the direction parallel to the symmetry axis.

[0139] In an embodiment, the interpolation unit 860 comprises:

[0140] A third establishing unit is configured to establish a corresponding relationship for all blocks in the following manner: each known point block corresponds to a to-be-interpolated point block, and each known point overlap block corresponds to a to-be-interpolated point overlap block.

[0141] A predicting unit is configured to predict all to-be-interpolated points of each to-be-interpolated point block one by one, and if the prediction is successful, perform spatial interpolation on all to-be-interpolated points.

[0142] In an embodiment, the predicting unit comprises:

[0143] An interpolation subunit is configured to select any to-be-interpolated point block as a current to-be-interpolated point block, select a first to-be-interpolated point in the current to-be-interpolated point block as a current to-be-interpolated point, and determine whether the current to-be-interpolated point is within the circumscribed range, if not, skip the current to-be-interpolated point, if yes, determine that the prediction is successful, and perform spatial interpolation on the current to-be-interpolated point by using known points within a predetermined range, until all to-be-interpolated points in the current block are predicted and spatial interpolation is performed on the to-be-interpolated points whose prediction is successful.

[0144] A releasing unit is configured to release the memory of the current to-be-interpolated point block, and then predict all to-be-interpolated points of the remaining to-be-interpolated point blocks one by one and perform spatial interpolation on the to-be-interpolated points whose prediction is successful.

[0145] The embodiment of the present application provides a computer device, which comprises a memory, a processor, and a computer program stored in the memory and executable on the processor, and the processor implements the method for generating a stratigraphic interface by using seismic interpretation data when executing the computer program.

[0146] The embodiment of the present application provides a computer readable storage medium, which stores a computer program, and the computer program is executable on a processor to implement the method for generating a stratigraphic interface by using seismic interpretation data.

[0147] The above merely illustrates the specific embodiments of the present application, but the protection scope of the present application is not limited thereto, and any skilled person in the art can easily think of various equivalent modifications or replacements within the technical range disclosed by the present application, and these modifications or replacements shall be covered within the protection scope of the present application. Therefore, the protection scope of the present application shall be subject to the protection scope of the claims.

Claims

1. A method for generating stratigraphic interfaces using seismic interpretation data, characterized in that: include: Step 1: Acquire geological data after geological exploration of the project area through survey lines; wherein the geological exploration includes seismic exploration; Step 2: Comprehensively interpreting the geological data to obtain a series of data points of stratigraphic interfaces distributed along the survey line; Step 3: Delineate the project area to obtain a delineated range of the project area, and take a minimum rectangle containing the delineated range along the boundary of the delineated range to obtain a rectangular range; Step 4: Divide the rectangular range into multiple grids, and use the nodes of each grid as the interpolation points of the stratum interface; Step 5: Select any data point of the stratigraphic interface as a known point, and use a spatial search strategy to divide all the points to be interpolated and the known points to obtain a plurality of corresponding blocks of points to be interpolated and blocks of known points; Step 6: Use each known point block to interpolate the corresponding point block to obtain the current stratigraphic interface; Step 7: Repeat steps 5-6 to obtain all stratum interfaces in the project area; Step 5 includes: dividing the rectangular range from the lower left corner along the x-direction by a third predetermined length and dividing it along the y-direction by a fourth predetermined length to obtain a plurality of known point blocks; reducing each known point block by a first predetermined distance toward the center point of the corresponding known point block to obtain a corresponding to-be-interpolated point block; Establish overlapping blocks for all to-be-interpolated point blocks and known point blocks respectively; obtain to-be-interpolated points in all known point blocks that are not covered by the overlapping blocks or to-be-interpolated point blocks; and divide the to-be-interpolated points into the overlapping blocks or to-be-interpolated point blocks with the closest distance.

2. The method for generating stratum interfaces using seismic interpretation data according to claim 1, characterized in that: The step 3 includes: Input multiple vertex coordinates and connect them to get the scope of the project area; Taking out the X and Y values ​​of all vertex coordinates within the circled range to form an X set and a Y set; Calculate the maximum and minimum values ​​in the X set and the Y set respectively, and combine the maximum value of the X set with the maximum value of the Y set, the maximum value of the X set with the minimum value of the Y set, the minimum value of the X set with the maximum value of the Y set, and the minimum value of the X set with the minimum value of the Y set to obtain four coordinate points; Connect the four coordinate points to get a rectangular range.

3. The method for generating stratum interfaces using seismic interpretation data according to claim 1, wherein: The step 4 comprises: The rectangular range is divided along the x direction according to a first predetermined length, and the rectangular range is divided along the y direction according to a second predetermined length to obtain multiple grids, and the nodes of all the grids are used as the points to be interpolated of the stratum interface.

4. The method for generating stratum interfaces using seismic interpretation data according to claim 1, wherein: The step of establishing overlapping blocks for all the blocks of points to be interpolated and the blocks of known points respectively includes: Selecting the number of blocks of the known point block; If the selected number of blocks is a first value, a first known point overlapping block is established for all known point blocks according to a first establishment method; If the selected number of blocks is a second value, then a second known point overlapping block is established for all known point blocks according to a second establishment method; Shrinking the first known point overlapping block by a second predetermined distance toward the center point of the corresponding first known point overlapping block to obtain a first to-be-interpolated point overlapping block; The second known point overlapping block is reduced by a third predetermined distance toward the center point of the corresponding second known point overlapping block to obtain a second to-be-interpolated point overlapping block.

5. The method for generating stratum interfaces using seismic interpretation data according to claim 4, characterized in that: The first value is 2, and the first establishment method is: Taking the boundary line of two adjacent known point blocks as the axis of symmetry, setting the overlapping block length in a direction perpendicular to the axis of symmetry, and setting the overlapping block height in a direction parallel to the axis of symmetry, establishing a rectangular overlapping block according to the overlapping block length and the overlapping block height, and serving as the first known point overlapping block.

6. The method for generating stratum interfaces using seismic interpretation data according to claim 5, characterized in that: The second value is 4, and the second establishment method is: Taking the intersection point of four adjacent known point blocks as the center point, setting the overlapping block length in a direction perpendicular to the axis of symmetry, and setting the overlapping block height in a direction parallel to the axis of symmetry, establishing a square overlapping block according to the overlapping block length and the overlapping block height, and serving as the second known point overlapping block.

7. The method for generating stratum interfaces using seismic interpretation data according to claim 6, characterized in that: The step 6 comprises: Establish a correspondence between all blocks as follows: each known point block corresponds to a block of points to be interpolated, and each overlapping block of known points corresponds to an overlapping block of points to be interpolated; All the points to be interpolated in each of the blocks of points to be interpolated are predicted one by one. If the prediction is successful, spatial interpolation is performed on all the points to be interpolated.

8. The method for generating stratum interfaces using seismic interpretation data according to claim 7, characterized in that: The step of predicting all the points to be interpolated in each of the blocks of points to be interpolated one by one, and if the prediction is successful, performing spatial interpolation on all the points to be interpolated, includes: Select any of the blocks of points to be interpolated as the current block of points to be interpolated, select the first point to be interpolated from the current block of points to be interpolated as the current point to be interpolated, and determine whether the current point to be interpolated is within the circled range. If not, skip the current point to be interpolated. If it is, determine that the prediction is successful, and perform spatial interpolation on the current point to be interpolated using known points within the predetermined range until all points to be interpolated in the current block are predicted and spatial interpolation is performed on the successfully predicted points to be interpolated. Release the memory of the current block of points to be interpolated, then predict all the points to be interpolated in the remaining blocks of points to be interpolated one by one and perform spatial interpolation on the successfully predicted points to be interpolated.

9. A device for generating stratum interfaces using seismic interpretation data, used to implement the method for generating stratum interfaces using seismic interpretation data according to any one of claims 1 to 8, characterized in that: include: A data acquisition unit, configured to acquire geological data obtained through geological exploration of the project area through survey lines; wherein the geological exploration includes seismic exploration; An interpretation unit, configured to perform comprehensive interpretation on the geological data to obtain a series of data points of stratigraphic interfaces distributed along the survey line; a delineating unit, configured to delineate the engineering area to obtain a delineated range of the engineering area, and to obtain a rectangular range by taking a minimum rectangle containing the delineated range along a boundary of the delineated range; A grid division unit, configured to divide the rectangular range into a plurality of grids, and use the nodes of each grid as the points to be interpolated on the stratum interface; A block division unit is used to select any data point of the stratum interface as a known point, and divide all the points to be interpolated and the known points respectively by using a spatial search strategy to obtain a plurality of corresponding blocks of points to be interpolated and blocks of known points; An interpolation unit is used to interpolate the corresponding block of points to be interpolated using each known point block to obtain the current stratigraphic interface; Repeating unit, used for repeating block division unit and interpolation unit to obtain all stratigraphic interfaces in the project area; a block division unit, specifically configured to divide the rectangular range from the lower left corner along the x-direction by a third predetermined length and along the y-direction by a fourth predetermined length to obtain a plurality of known point blocks; and reduce each known point block by a first predetermined distance toward the center point of the corresponding known point block to obtain a corresponding block of points to be interpolated; Establish overlapping blocks for all to-be-interpolated point blocks and known point blocks respectively; obtain to-be-interpolated points in all known point blocks that are not covered by the overlapping blocks or to-be-interpolated point blocks; and divide the to-be-interpolated points into the overlapping blocks or to-be-interpolated point blocks with the closest distance.

10. A computer device, characterized in that: The method comprises a memory, a processor and a computer program stored in the memory and executable on the processor, wherein when the processor executes the computer program, the method for generating stratigraphic interfaces using seismic interpretation data as described in any one of claims 1 to 8 is implemented.

11. A computer-readable storage medium, characterized in that The computer-readable storage medium stores a computer program, which, when executed by a processor, implements the method for generating stratigraphic interfaces using seismic interpretation data according to any one of claims 1 to 8.

Images