Method and device for generating stratigraphic interface, computer equipment and medium

Abstract

The embodiment of the invention provides a method and device for generating a stratum interface, computer equipment and a medium, and the method comprises the steps: explaining a geological number, carrying out the partitioning of a known point obtained through the explaining and a to-be-interpolated point needed by the generation of the stratum interface, enabling a known point block to be in one-to-one correspondence with a to-be-interpolated point block, and enabling the to-be-interpolated point block to be in one-to-one correspondence with the known point block; global search in the interpolation process is converted into local search, then spatial interpolation is carried out on partitions one by one, and the program running efficiency is improved.

Description

Technical Field

[0001] The present invention relates to the field of stratum technology, and in particular to a method, device, computer equipment and medium for generating stratum interfaces. Background Art

[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, with significant 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 stratigraphic interfaces based on seismic interpretation data. This process usually relies on spatial interpolation technology, which can predict the stratigraphic properties of unknown points through certain algorithms based on known seismic interpretation data points. However, directly importing the original seismic interpretation data into the interpolation algorithm for calculation often leads to problems such as low program running efficiency and running out of memory.

[0004] For example, when seismic interpretation data is imported into the program, several points to be interpolated in the engineering site are selected, and inverse distance interpolation is used, during the interpolation process, it is often necessary to set the interpolation influence radius. That is to say, before predicting each point to be interpolated, it is necessary to search for a part within the interpolation influence radius among the known points obtained from the interpretation of the entire engineering site. However, the amount of seismic interpretation data is large. In this case, the operation of the program is very time-consuming. Summary of the invention

[0005] The embodiments of the present invention provide a method, device, computer equipment and medium for generating stratigraphic interfaces, aiming to solve the problem that the existing interpolation technology directly imports the original seismic interpretation data into the interpolation algorithm for calculation, resulting in low program running efficiency.

[0006] In a first aspect, an embodiment of the present invention provides a method for generating a stratigraphic interface using seismic interpretation data, comprising:

[0007] Step 1: Acquire geological data after geological exploration of the project area through survey lines; wherein the geological exploration includes seismic exploration;

[0008] Step 2: Comprehensively interpret the geological data to obtain a series of data points of stratigraphic interfaces distributed along the survey line;

[0009] Step 3: delineate the engineering area to obtain a delineated range of the engineering area, and take a minimum rectangle including the delineated range along the boundary of the delineated range to obtain a rectangular range;

[0010] Step 4: Divide the rectangular range into multiple grids, and use the nodes of each grid as the interpolation points of the stratum interface;

[0011] 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;

[0012] Step 6: Use each known point block to interpolate the corresponding point block to be interpolated to obtain the current stratigraphic interface;

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

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

[0015] A data acquisition unit, used to acquire geological data after geological exploration of the engineering area through a survey line; wherein the geological exploration includes seismic exploration;

[0016] An interpretation unit, used for comprehensively interpreting the geological data to obtain a series of data points of stratigraphic interfaces distributed along the survey line;

[0017] A delineation unit is used 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 including the delineated range along a boundary of the delineated range;

[0018] A grid division unit, used for dividing the rectangular range into a plurality of grids, and using the nodes of each grid as the interpolation points of the stratum interface;

[0019] A block division unit is used to select any data point of the stratigraphic interface as a known point, and to 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;

[0020] An interpolation unit is used to perform interpolation processing on each block of points to be interpolated and each block of known points to obtain the current stratigraphic interface;

[0021] Repeating unit is used to repeat the block division unit and interpolation unit to obtain all the stratigraphic interfaces in the project area.

[0022] In a third aspect, an embodiment of the present invention provides a computer device, comprising 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 above is implemented.

[0023] In a fourth aspect, an embodiment of the present invention provides a computer-readable storage medium having a computer program stored thereon, and when the computer program is executed by a processor, the method for generating stratigraphic interfaces using seismic interpretation data as described above is implemented.

[0024] The embodiment of the present invention provides a method, device, computer equipment and medium for generating stratigraphic interfaces. The method interprets geological data and partitions the known points obtained by the interpretation and the interpolation points required for generating stratigraphic interfaces, so that the known point blocks and the interpolation point blocks correspond one to one, converts the global search in the interpolation process into a local search, and then performs spatial interpolation on each partition, thereby improving the program running efficiency. BRIEF DESCRIPTION OF THE DRAWINGS

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

[0026] Figure 1 A schematic flow chart of a method for generating a stratigraphic interface using seismic interpretation data provided by an embodiment of the present invention;

[0027] Figure 2 A schematic diagram of a sub-process of a method for generating a stratum interface using seismic interpretation data provided by an embodiment of the present invention;

[0028] Figure 3 A schematic diagram of another sub-process of a method for generating a stratum interface using seismic interpretation data provided by an embodiment of the present invention;

[0029] Figure 4 A schematic diagram of another sub-process of a method for generating a stratum interface using seismic interpretation data provided by an embodiment of the present invention;

[0030] Figure 5 A schematic diagram of another sub-process of a method for generating a stratum interface using seismic interpretation data provided by an embodiment of the present invention;

[0031] Figure 6 A schematic diagram of another sub-process of a method for generating a stratum interface using seismic interpretation data provided by an embodiment of the present invention;

[0032] Figure 7 A schematic diagram of another sub-process of a method for generating a stratum interface using seismic interpretation data provided by an embodiment of the present invention;

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

[0034] Fig. 9 A schematic diagram of the delineation of the seismic interpretation data survey line and the engineering area;

[0035] Fig.10 A schematic diagram of a rectangular range formed for encircling the range;

[0036] Fig.11 A schematic diagram for generating points to be interpolated based on a rectangular range;

[0037] Fig.12 A schematic diagram of dividing known points into blocks according to rectangular range;

[0038] Fig.13 A schematic diagram showing a one-to-one correspondence between the known point block division and the to-be-interpolated point block division;

[0039] Fig.14 A schematic diagram for establishing overlapping blocks between two known point blocks and overlapping blocks of points to be interpolated;

[0040] Fig.15 It is a complete schematic diagram of the one-to-one correspondence between the overlapping blocks of two known point blocks and the overlapping blocks of the points to be interpolated;

[0041] Fig.16 A schematic diagram for establishing overlapping blocks between four known point blocks and overlapping blocks of points to be interpolated;

[0042] Fig.17 A schematic diagram of an area not covered by overlapping blocks or blocks to be interpolated;

[0043] Fig.18 A complete schematic diagram after dividing all blocks and establishing overlapping blocks;

[0044] Fig.19 It is a schematic diagram of the actual engineering survey line, wind field range and rectangular range;

[0045] Fig. 20 This is a schematic diagram of the division of known point blocks in actual engineering;

[0046] Fig.21 A complete schematic diagram after dividing all blocks in the actual project and establishing overlapping blocks;

[0047] Fig. 22 This is the display effect diagram of the generated stratigraphic interface. DETAILED DESCRIPTION

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

[0049] It should be understood that when used in this specification and the appended claims, the terms "include" and "comprises" indicate the presence of described features, integers, steps, operations, elements and / or components, but do not exclude the presence or addition of one or more other features, integers, steps, operations, elements, components and / or combinations thereof.

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

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

[0052] See also Figure 1 The embodiment of the present invention provides a method for generating a stratigraphic interface using seismic interpretation data, including S10-S70:

[0053] S10, obtaining geological data after geological exploration of the project area through survey lines; wherein the geological exploration includes seismic exploration;

[0054] In this step, a series of survey lines need to be arranged in the project area. The survey lines can be straight lines, broken lines or curves. The specific shape and arrangement method depends on the geological conditions, exploration targets and exploration accuracy requirements. Then, geological exploration is carried out along the survey lines to obtain geological data. Among them, geological exploration includes seismic exploration. Seismic exploration is a geophysical exploration method that uses the characteristics of artificially excited seismic waves propagating in underground media to infer the properties and morphology of underground rock formations. In the process of seismic exploration, artificial seismic sources such as explosive explosions and controllable vibrators are used to excite seismic waves, and then precision instruments are used to record the propagation of seismic waves on the surface and underground. Through seismic exploration instruments, data on the propagation of seismic waves in underground media are collected, including parameters such as the amplitude, frequency, and phase of seismic waves. These parameters are used as seismic exploration data for subsequent data processing and analysis.

[0055] In specific implementation, geological exploration also includes borehole exploration or static exploration. Borehole exploration is to obtain direct samples of underground rock formations through drilling. Static exploration (CPT) is an in-situ testing technology that uses a pressure device to press the probe into the soil at a certain rate, and evaluates the properties of the soil layer by measuring various data during the penetration of the probe. In the implementation process, seismic exploration data and borehole data (or CPT data) are needed as geological data.

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

[0057] From the above steps, it can be seen that geological data includes seismic exploration data and drilling data (or CPT data). All of these data need to be interpreted. For the interpretation of seismic exploration data, different strata are identified by observing the reflection wave morphology, amplitude and continuity in the seismic profile. The depth and thickness of each stratum are calculated based on the trend and velocity information of the reflection wave, and the physical properties of the strata are analyzed through the propagation velocity of seismic waves. For the interpretation of drilling data (or CPT data), the type, composition and structure of the strata are analyzed through the core description of the drilling samples and the physical and mechanical property test data (or the penetration resistance and sidewall friction resistance of the strata are analyzed using the CPT test data to further understand the strength and deformation characteristics of the strata). According to the analysis results of the drilling data / CPT data, the position of the stratum interface is determined and compared and verified with the interpretation results of the seismic exploration data. Specifically, the seismic exploration data and the borehole data (or CPT data) are matched and aligned in space and time to ensure their consistency. The interpretation of the seismic exploration data and the results of the borehole data / CPT data interpretation are combined to determine the specific location of the stratigraphic interface, thereby obtaining a series of stratigraphic interface data points 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 including the delineated range along the boundary of the delineated range to obtain a rectangular range;

[0059] In this step, the project area must be delineated first so that the project area has a clear delineation range. In the specific implementation, the delineation range may be an irregular shape or a regular shape. Therefore, the delineation range can be simplified to the minimum circumscribed rectangle of the delineation range, thereby obtaining a rectangular range. The rectangular range can improve the subsequent processing efficiency.

[0060] In one embodiment, if Figure 2 As shown, S30 includes:

[0061] S31, input multiple vertex coordinates, and connect the vertex coordinates to obtain the delineated range of the project area;

[0062] S32, taking out the X values ​​and Y values ​​of all vertex coordinates within the circled range to form an X set and a Y set;

[0063] S33, respectively calculating the maximum value and the minimum value in the X set and the Y set, and combining 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;

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

[0065] In this embodiment, multiple vertex coordinates can be input clockwise or counterclockwise, and these vertex coordinates are connected in sequence to form a polygon, which is the demarcation range of the project area (see Fig. 9 The yellow box), then take out the X and Y values ​​of these vertex coordinates to form the {X} set: {x1,x2,…,x n} and {Y} set: {y1,y2,…,y n}, calculate the minimum value x from the set {X} min and the maximum value x max , y is calculated from the set {Y} min and max , and combine them to get the coordinate point Point1: (x min, y min ), Point 2: ( xmin ,y max ), Point 3: (x max ,y max ), Point4: (x max, y min ), connect these four coordinate points to form a rectangular range (refer to Fig.10 red frame).

[0066] S40, dividing the rectangular range into a plurality of grids, and using the nodes of each grid as the interpolation points of the stratum interface;

[0067] In this step, the rectangular range is divided into a plurality of grids. After the division, each grid node (ie, grid point) is used as a point to be interpolated at the stratum interface, thereby improving the accuracy of subsequent interpolation.

[0068] In one 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 a plurality of grids, and the nodes of all the grids are used as the interpolation points of the stratum interface.

[0070] In this embodiment, the first predetermined length and the second predetermined length may be the same length or different lengths. If they are the same lengths, 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 of 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 Fig.11 .

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

[0072] In this step, one of the data points distributed along the survey line of the stratigraphic interface is 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 in the range to interpolate the points to be interpolated. Before interpolation, it is necessary to use a spatial search strategy to partition the points to be interpolated and the 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 according to a third predetermined length and along the y direction according to a fourth predetermined length from the lower left corner thereof to obtain a plurality of known point blocks;

[0075] S52, shrinking each known point block by a first predetermined distance toward the center point of the corresponding known point block to obtain a corresponding point block to be interpolated.

[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 (the third predetermined length l x ) and the length in the Y direction (the fourth predetermined length l y ), starting from the lower left corner of the rectangular range, until the entire rectangular range is covered, thereby obtaining multiple known point blocks, and the length of the last block along the x direction is still l x , part of the length exceeds the rectangular range, and the length of the last block along the y direction is still l y , some of the lengths are beyond the rectangular range, the schematic diagram of the known point block division can be referred to Fig.12 Then, each of the divided known point blocks is reduced toward the center point of the respective known point blocks according to the first predetermined distance, so as to ensure that each known point block obtains a corresponding block of points to be interpolated, that is, the number of known point blocks and blocks of points to be interpolated is the same. The schematic diagram of the division of blocks of points to be interpolated can be referred to Fig.13 The third predetermined length and the fourth predetermined length may be set to be the same or different, and may be set according to specific application scenarios. The first predetermined distance is smaller than the third predetermined length and the fourth predetermined length, and corresponding values ​​may be set according to specific scenarios.

[0077] In one embodiment, if Figure 4 As shown, S50 also includes:

[0078] S53, respectively establishing overlapping blocks for all the to-be-interpolated point blocks and the known point blocks;

[0079] S54, obtaining the to-be-interpolated points in all known point blocks that are not covered by the overlapping blocks or the to-be-interpolated point blocks;

[0080] S55: Divide the to-be-interpolated points into the overlapping blocks or to-be-interpolated point blocks with the closest distance.

[0081] In this embodiment, when using known points to interpolate the points to be interpolated, it is necessary to select known points within a certain influence radius r (i.e., a range to be interpolated with a radius of r and the point to be interpolated as the center of the circle) around the points to be interpolated, and use the known points within the range to be interpolated to interpolate the points to be interpolated. Since the blocks of points to be interpolated are smaller than the blocks of known points by a third predetermined distance in the four directions of up, down, left, and right when divided, the divided multiple blocks of points to be interpolated cannot cover all the points to be interpolated. At this time, it is necessary to establish overlapping areas between the blocks of known points and overlapping areas between the blocks of points to be interpolated, and then obtain the points to be interpolated that are not covered by the overlapping blocks or the blocks of points to be interpolated in all the known point blocks. These points to be interpolated are located in the areas not covered by the overlapping blocks or the blocks to be interpolated (refer to Fig.17 ), then determine which overlapping block or interpolation point block these interpolation points are close to, and then divide the interpolation points into the overlapping block or interpolation point block closest to it (refer to Fig.18 ), so far, all the points to be interpolated have been preliminarily processed.

[0082] In one embodiment, if Figure 5 As shown, S53 includes:

[0083] S531, selecting the number of blocks of the known point block;

[0084] S532, if the selected number of blocks is a first value, establishing a first known point overlapping block for all known point blocks according to a first establishing method;

[0085] S533, if the selected number of blocks is a second value, establishing a second known point overlapping block for all known point blocks according to a second establishing method;

[0086] S534, reducing 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 overlapping block of points to be interpolated;

[0087] S535, reducing the second known point overlapping block by a third predetermined distance toward the center point of the corresponding second known point overlapping block to obtain a second overlapping block of points to be interpolated.

[0088] In this embodiment, overlapping blocks between known point blocks are first established. There are two establishment methods according to the different numbers of blocks in the known point blocks, thereby obtaining two known point overlapping blocks. Then, two to-be-interpolated point overlapping blocks are established based on the two known point overlapping blocks. The establishment method is: respectively reduce the corresponding predetermined distance r to the center point of the corresponding known point overlapping block, that is, the second predetermined distance and the third predetermined distance. These two predetermined distances can be set to be the same or different.

[0089] In a specific implementation, the first value is 2, and the first establishment method is:

[0090] like Fig.14 and Fig.15 As shown, the boundary line of two adjacent known point blocks is used as the symmetry axis, the overlapping block length is set in a direction perpendicular to the symmetry axis, and the overlapping block height is set in a direction parallel to the symmetry axis, and a rectangular overlapping block is established according to the overlapping block length and the overlapping block height, and is used as the first known point overlapping block. Among them, the overlapping block length can be set to 4r (r can be set by itself), and the overlapping block height can be set to the length of the boundary line.

[0091] In a specific implementation, the second value is 4, and the second establishment method is:

[0092] like Fig.16 As shown, with the intersection point of four adjacent known point blocks as the center point, the overlapping block length is set in a direction perpendicular to the symmetry axis, and the overlapping block height is set in a direction parallel to the symmetry axis, and a square overlapping block is established according to the overlapping block length and the overlapping block height, and is used as the second known point overlapping block. Among them, the overlapping block length is equal to the overlapping block height, and both can be set to 4r (r can be set at will).

[0093] S60, using each known point block to interpolate the corresponding point block to be interpolated, to obtain the current stratum interface;

[0094] In this step, each known point block is used to interpolate the corresponding point block to be interpolated one by one, so as to complete the interpolation of all the point blocks to be interpolated, thereby obtaining the entire stratigraphic interface (ie, the current stratigraphic interface). This interpolation method can improve the program running efficiency.

[0095] In one embodiment, if Figure 6 As shown, S60 includes:

[0096] S61, establishing a correspondence relationship for all blocks in the following manner: each known point block corresponds to a block of points to be interpolated, and each known point overlapping block corresponds to an overlapping block of points to be interpolated;

[0097] S62, 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.

[0098] In this embodiment, in order to better interpolate all blocks to be interpolated, it is preferred to establish a corresponding relationship for all blocks, which can be established by numbering, such as known point block 1 corresponds to point block 1 to be interpolated, known point overlapping block 1 corresponds to point overlapping block 1 to be interpolated, and according to the numbering order, the interpolation points of all blocks of points to be interpolated 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 retained first and then waits for subsequent processing.

[0099] In one embodiment, if Figure 7 As shown, S62 includes:

[0100] S621, selecting any of the to-be-interpolated point blocks as the current to-be-interpolated point block, and selecting the first to-be-interpolated point from the current to-be-interpolated point block as the current to-be-interpolated point, and judging whether the current to-be-interpolated point is within the circled range, if not, skipping the current to-be-interpolated point, and if so, judging that the prediction is successful, and performing spatial interpolation on the current to-be-interpolated point through known points within the predetermined range, until all to-be-interpolated points of the current block are predicted and spatial interpolation is performed on the successfully predicted to-be-interpolated points;

[0101] S622: Release the memory of the current block of points to be interpolated, 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.

[0102] In this embodiment, a known point block numbered 1 and a block of points to be interpolated can be selected, and the first point to be interpolated in the block of points to be interpolated can be selected to determine whether the point to be interpolated is within the demarcated range of the project area. Specifically, the ray method can be used to determine whether the point to be interpolated is within the demarcated range. If it is within the demarcated range, the point to be interpolated is skipped. If it is within the demarcated range, the known point is searched in the known point block 1 by the interpolation influence radius r to obtain the range to be interpolated, and the point to be interpolated is spatially interpolated through the known points within the range to be interpolated. Repeat the above process until all the points to be interpolated in the block of points to be interpolated are predicted, and then release the memory of the known points and the points to be interpolated in the block, and then perform spatial interpolation of blocks 2, blocks 3... overlapping blocks 1, overlapping blocks 2... until all blocks are predicted and interpolated.

[0103] S70. Repeat S50-S60 to obtain all stratum interfaces in the engineering area.

[0104] The stratigraphic interfaces are generated in the manner of the aforementioned steps S50-S60 until all stratigraphic interfaces in the project area are generated.

[0105] The embodiment of the present invention interprets geological data and partitions the known points obtained by the interpretation and the interpolation points required for generating stratigraphic interfaces, so that the known point blocks and the interpolation point blocks correspond one to one, converts the global search in the interpolation process into a local search, and then performs spatial interpolation on a partition by partition. When the interpolation of a partition is completed, the memory storing the data points in the partition can be released, thereby improving the program running efficiency and avoiding the situation of over-memory.

[0106] In practical applications, such as Figure 19-21 As shown in the figure, there will be a wind field range in the project area, which will affect the distribution of the survey lines in the project area, making the demarcation range of the project area irregular. Then, according to the above steps, the demarcation range is simplified to the minimum circumscribed rectangle, thereby obtaining a rectangular range, and then the rectangular range is divided to obtain multiple known point blocks, and then the points to be interpolated are divided and various overlapping blocks are established. The division method and establishment method can refer to the above steps and will not be repeated here. Then generate Fig. 22 The stratigraphic interface shown, where Fig. 22 The Z in it represents the depth value of the point to be interpolated.

[0107] The embodiment of the present invention also provides a device for generating a stratum interface using seismic interpretation data, and the device is used to execute any embodiment of the method for generating a stratum interface using seismic interpretation data. Figure 8 , Figure 8 800 is a schematic block diagram of a device for generating a stratum interface using seismic interpretation data provided by an embodiment of the present invention. The device 800 for generating a stratum interface using seismic interpretation data includes:

[0108] The data acquisition unit 810 is used to acquire geological data after geological exploration of the project area through the survey line; wherein the geological exploration includes seismic exploration;

[0109] An interpretation unit 820 is used to perform a comprehensive interpretation on the geological data to obtain a series of data points of the stratigraphic interfaces distributed along the survey line;

[0110] A delineation unit 830 is used 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 including the delineated range along a boundary of the delineated range;

[0111] A grid division unit 840 is used to divide the rectangular range into a plurality of grids, and use the nodes of each grid as the interpolation points of the stratum interface;

[0112] The block division unit 850 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;

[0113] The interpolation unit 860 is used to perform interpolation processing on each block of points to be interpolated and each block of known points to obtain the current stratigraphic interface;

[0114] The repeating unit 870 is used to repeat the block division unit and the interpolation unit to obtain all the stratigraphic interfaces in the project area.

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

[0116] An input unit, used to input multiple vertex coordinates and connect the vertex coordinates to obtain the delineated range of the project area;

[0117] A value obtaining unit, used for obtaining the X values ​​and Y values ​​of all vertex coordinates within the circled range to form an X set and a Y set;

[0118] a combining unit, for respectively calculating the maximum value and the minimum value in the X set and the Y set, and combining 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 used to connect the four coordinate points to obtain a rectangular range.

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

[0121] The first division unit is used to divide the rectangular range along the x direction according to a first predetermined length, and divide the rectangular range along the y direction according to a second predetermined length to obtain multiple grids, and use the nodes of all grids as the interpolation points of the stratum interface.

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

[0123] A second division unit is used to divide the rectangular range from the lower left corner along the x direction according to a third predetermined length, and divide it along the y direction according to a fourth predetermined length to obtain a plurality of known point blocks;

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

[0125] In one embodiment, the block division unit 850 further includes:

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

[0127] A to-be-interpolated point acquisition unit, used to acquire the to-be-interpolated points in all known point blocks that are not covered by the overlapping blocks or the to-be-interpolated point blocks;

[0128] The third division unit is used to divide the to-be-interpolated points into the overlapping blocks or to-be-interpolated point blocks with the closest distance.

[0129] In one embodiment, the overlapping block establishing unit includes:

[0130] A quantity selection unit, used for selecting the block quantity of the known point block;

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

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

[0133] A second reduction unit, configured to reduce the first known point overlapping block by a second predetermined distance toward a center point of the corresponding first known point overlapping block to obtain a first overlapping block of points to be interpolated;

[0134] The third reduction unit is used to reduce the second known point overlapping block by a third predetermined distance toward the center point of the corresponding second known point overlapping block to obtain a second overlapping block of points to be interpolated.

[0135] In one embodiment, the first establishing unit includes:

[0136] The first numerical value is 2, and the boundary line of two adjacent known point blocks is used as the symmetry axis. The overlapping block length is set in a direction perpendicular to the symmetry axis, and the overlapping block height is set in a direction parallel to the symmetry axis. A rectangular overlapping block is established according to the overlapping block length and the overlapping block height, and is used as the first known point overlapping block.

[0137] In one embodiment, the second establishing unit includes:

[0138] The second numerical value is 4, with the intersection point of 4 adjacent known point blocks as the center point, the overlapping block length is set in the direction perpendicular to the symmetry axis, and the overlapping block height is set in the direction parallel to the symmetry axis, and a square overlapping block is established according to the overlapping block length and the overlapping block height, and used as the second known point overlapping block.

[0139] In one embodiment, the interpolation unit 860 includes:

[0140] The third establishing unit is used to establish a corresponding relationship for all blocks in the following manner: each known point block corresponds to a block of points to be interpolated, and each known point overlapping block corresponds to an overlapping block of points to be interpolated;

[0141] The prediction unit is used to predict 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, perform spatial interpolation on all the points to be interpolated.

[0142] In one embodiment, the prediction unit comprises:

[0143] An interpolation subunit is used to select any of the to-be-interpolated point blocks as the current to-be-interpolated point block, select the first to-be-interpolated point from the current to-be-interpolated point block as the current to-be-interpolated point, and determine whether the current to-be-interpolated point is within the circled range. If not, the current to-be-interpolated point is skipped; if so, the prediction is determined to be successful, and the current to-be-interpolated point is spatially interpolated through known points within the predetermined range until all to-be-interpolated points of the current block are predicted and spatial interpolation is performed on the successfully predicted to-be-interpolated points;

[0144] The release unit is used to release the memory of the current block of points to be interpolated, and 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.

[0145] An embodiment of the present invention provides a computer device, including a memory, a processor, and a computer program stored in the memory and executable on the processor. When the processor executes the computer program, the method for generating stratigraphic interfaces using seismic interpretation data as described in the above embodiment is implemented.

[0146] An embodiment of the present invention provides a computer-readable storage medium having a computer program stored thereon. When the computer program is executed by a processor, the method for generating stratigraphic interfaces using seismic interpretation data as described in the above embodiment is implemented.

[0147] The above is only a specific embodiment of the present invention, but the protection scope of the present invention is not limited thereto. Any technician familiar with the technical field can easily think of various equivalent modifications or replacements within the technical scope disclosed by the present invention, and these modifications or replacements should be included in the protection scope of the present invention. Therefore, the protection scope of the present invention shall be based on 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 interpret the geological data to obtain a series of data points of stratigraphic interfaces distributed along the survey line; Step 3: delineate the engineering area to obtain a delineated range of the engineering area, and take a minimum rectangle including 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 be interpolated to obtain the current stratigraphic interface; Step 7: Repeat steps 5-6 to obtain all stratigraphic interfaces in the project area.

2. The method for generating stratigraphic interfaces using seismic interpretation data according to claim 1, characterized in that: The step 3 comprises: Input multiple vertex coordinates and connect them to get the delineated scope of the project area; Taking out the X values ​​and Y values ​​of all vertex coordinates within the circled range to form an X set and a Y set; Calculate the maximum value and the minimum value 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 stratigraphic interfaces using seismic interpretation data according to claim 1, characterized in that: 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 a plurality of grids, and the nodes of all the grids are used as the interpolation points of the stratum interface.

4. The method for generating stratigraphic interfaces using seismic interpretation data according to claim 1, characterized in that: The step 5 comprises: Dividing the rectangular range from the lower left corner along the x direction according to a third predetermined length, and dividing the rectangular range along the y direction according to a fourth predetermined length, to obtain a plurality of known point blocks; Each known point block is reduced by a first predetermined distance toward a center point of the corresponding known point block to obtain a corresponding point block to be interpolated.

5. The method for generating stratigraphic interfaces using seismic interpretation data according to claim 4, characterized in that: The step 5 also includes: Establish overlapping blocks for all blocks of points to be interpolated and blocks of known points respectively; Obtaining the points to be interpolated that are not covered by the overlapping blocks or the points to be interpolated blocks in all known point blocks; The points to be interpolated are divided into the overlapping blocks or the blocks of points to be interpolated that are closest to each other.

6. The method for generating stratigraphic interfaces using seismic interpretation data according to claim 5, characterized in that: The step of respectively establishing overlapping blocks for all the to-be-interpolated point blocks and the known point blocks 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, a second known point overlapping block is established for all known point blocks according to a second establishment method; The first known point overlapping block is reduced 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.

7. The method for generating stratigraphic interfaces using seismic interpretation data according to claim 6, 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.

8. The method for generating stratigraphic interfaces using seismic interpretation data according to claim 7, 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.

9. The method for generating stratigraphic interfaces using seismic interpretation data according to claim 8, characterized in that: The step 6 comprises: Establish a correspondence relationship between all blocks in the following manner: each known point block corresponds to a block of points to be interpolated, and each known point overlapping block 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.

10. The method for generating stratum interfaces using seismic interpretation data according to claim 9, 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 yes, determine that the prediction is successful, and perform spatial interpolation on the current point to be interpolated through known points within the predetermined range until all the points to be interpolated in the current block are predicted and spatial interpolation is performed on the successfully predicted points to be interpolated; The memory of the current block of points to be interpolated is released, and then all the points to be interpolated in the remaining blocks of points to be interpolated are predicted one by one, and spatial interpolation is performed on the successfully predicted points to be interpolated.

11. A device for generating stratigraphic interfaces using seismic interpretation data, used to implement the method for generating stratigraphic interfaces using seismic interpretation data as claimed in any one of claims 1 to 10, characterized in that: include: A data acquisition unit, used to acquire geological data after geological exploration of the engineering area through a survey line; wherein the geological exploration includes seismic exploration; An interpretation unit, used for comprehensively interpreting the geological data to obtain a series of data points of stratigraphic interfaces distributed along the survey line; A delineation unit is used 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 including the delineated range along a boundary of the delineated range; A grid division unit, used for dividing the rectangular range into a plurality of grids, and taking the nodes of each grid as the interpolation points of the stratum interface; A block division unit is used to select any data point of the stratigraphic interface as a known point, and to 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 perform interpolation processing on each block of points to be interpolated and each block of known points to obtain the current stratigraphic interface; Repeating unit is used to repeat the block division unit and interpolation unit to obtain all the stratigraphic interfaces in the project area.

12. 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. 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 10 is implemented.

13. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, the method for generating stratigraphic interfaces using seismic interpretation data according to any one of claims 1 to 10 is implemented.

Images