Method and system for acquiring sedimentary channel boundary diagram based on three-dimensional earthquake

Abstract

The invention provides a method and a system for acquiring a sedimentary river boundary diagram based on a three-dimensional earthquake, and belongs to the field of oil and gas field exploration. The acquisition method comprises the following steps: acquiring a plane attribute graph of a sedimentary river channel after frequency division attribute fusion, and performing meshing division on the plane attribute graph to obtain a plurality of single-grid water channel plane attribute graphs; carrying out segmentation according to a preset segmentation step length, determining a water channel superposition type on each segmentation step length, and determining a plane combination type of each water channel in the single grid according to the number of the water channels in the single grid; acquiring boundary points of each color block in the corresponding color code on each segmentation step length, and combining the boundary points of the adjacent segmentation step lengths according to a water channel superposition type and a plane combination type to obtain a boundary diagram of each water channel of the single grid; and combining the boundaries of the water channels in the adjacent planar single grids to obtain a sedimentary river channel boundary diagram. Therefore, the boundary points of the color blocks are accurately picked up by giving each color code water channel superposition type, and then the boundary points are combined, so that the mapping efficiency is improved, and the error probability is reduced.

Description

Technical Field

[0001] The invention relates to the technical field of oil and gas field exploration, and in particular to a method and system for acquiring a sedimentary channel boundary map based on three-dimensional seismic. Background Art

[0002] As the exploration and development of oil and gas fields in the world continues to deepen, the precision of reserve unit characterization is increasing, and the demand for comprehensive analysis of large amounts of three-dimensional seismic data is more urgent. Channel sedimentation, as a microfacies type in river and delta sedimentation, is widely developed in rivers and lobes formed by rivers entering the sea. It not only constitutes the basic skeleton structure of rivers and deltas, but is also an important high-quality reservoir for oil and gas enrichment today. It has always been the focus of research on river-delta sedimentary systems. This type of reservoir is affected by factors such as paleo-structure, paleo-current, sea level rise and fall, and material supply. Multiple water channels overlap each other in space, forming a complex geological structure with vertical superposition and horizontal interlacing. Due to the uneven degree of connectivity between different water channels, the development effect of this type of oil and gas reservoir is seriously affected.

[0003] In this complex waterway system, the sand body thickness of each single waterway is uneven and the width of the waterway varies, making it difficult to form a graph with actual sedimentary waterway significance. The existing characterization methods for studying the sedimentary distribution of complex waterways developed in river-delta sedimentary systems require a lot of manual work to outline point by point and waterway by waterway, which is not only inefficient but also has large errors. Summary of the invention

[0004] The purpose of the embodiments of the present invention is to provide a method and system for acquiring a sedimentary channel boundary map based on three-dimensional seismic data, so as to solve the technical problems existing in the above-mentioned prior art.

[0005] In order to achieve the above object, an embodiment of the present invention provides a method for obtaining a sedimentary channel boundary map based on three-dimensional seismic data, comprising:

[0006] Obtaining a plane attribute map of the sedimentary river channel after frequency division attribute fusion, and gridding the plane attribute map of the sedimentary river channel to obtain a plurality of single-grid waterway plane attribute maps;

[0007] According to a preset segmentation step, each single-grid waterway plane attribute map in the plurality of single-grid waterway plane attribute maps is segmented respectively, and according to a pre-constructed correspondence between the waterway superposition type and the color scale, the waterway superposition type of each waterway at each segmentation step is determined;

[0008] For each single-grid waterway plane attribute map among the plurality of single-grid waterway plane attribute maps, determining the plane combination type of each waterway in the single-grid waterway plane attribute map according to the number of waterways in the single grid and the waterway superposition type of each waterway at each segmentation step in the single grid;

[0009] For each single-grid waterway plane attribute map among the multiple single-grid waterway plane attribute maps, the boundary points of each color block in the color code corresponding to the waterway overlay type of each waterway on each segmentation step are obtained, and according to the waterway overlay type of each waterway on each segmentation step and the plane combination type of each waterway in the single grid, the boundary points of adjacent segmentation steps are combined to obtain the boundary map of each waterway in the single grid; and the boundaries of each waterway in adjacent plane single grids are combined to obtain the boundary map of the sedimentary river channel.

[0010] Optionally, the pre-constructed correspondence between the waterway overlay type and the color code includes:

[0011] If the water channel overlay type is the first overlay type, the corresponding first color code is the first color-second color-first color or the first color-second color-third color; if the water channel overlay type is the second overlay type, the corresponding second color code is the first color-third color / fourth color-first color; if the water channel overlay type is the third overlay type, the corresponding third color code is the first color or the third color; if the water channel overlay type is the fourth overlay type, the corresponding fourth color code is the first color or the third color; if the water channel overlay type is the fifth overlay type, the corresponding fifth color code is the second color, wherein the first color represents thick sand type, the second color represents mud type, and the third color and the fourth color represent medium thick sand type.

[0012] Optionally, before determining the water channel superposition type of each water channel at each segmentation step according to the pre-constructed correspondence between the water channel superposition type and the color code, the acquisition method further includes:

[0013] The water channel superposition type is defined according to the number and distribution characteristics of the water channels at each segmentation step.

[0014] Optionally, the water channel superposition type is defined according to the number and distribution characteristics of the water channels at each segmentation step, including:

[0015] If there are multiple single water channels on each segmentation step and the multiple single water channels are independent of each other, the water channel overlay type is defined as the first overlay type; if there are multiple single water channels on each segmentation step and the multiple single water channels overlap each other, the water channel overlay type is defined as the second overlay type; if there are multiple single water channels on each segmentation step and the multiple single water channels merge with each other, the water channel overlay type is defined as the third overlay type; if there is a single single water channel on each segmentation step, the water channel overlay type is defined as the fourth overlay type; if there is no water channel on each segmentation step, the water channel overlay type is defined as the fifth overlay type.

[0016] Optionally, the determination of the number of waterways in the single grid includes:

[0017] If the number of color blocks in the color mark within the single grid is greater than the first preset threshold, then the number of water channels is multiple; if the number of color blocks in the color mark within the single grid is greater than the second preset threshold and less than or equal to the first preset threshold, and the color of the color blocks is not all the second color representing the mud type, then the number of water channels is single; if the number of color blocks in the color mark within the single grid is equal to the second preset threshold, and the color of the color blocks is the second color, then the number of water channels is zero, wherein the first preset threshold is greater than the second preset threshold.

[0018] Optionally, the determining of the plane combination type of each water channel in the single-grid water channel plane attribute graph according to the number of water channels in the single grid and the water channel superposition type of each water channel at each segmentation step in the single grid includes:

[0019] If there are multiple water channels and the superposition type of two water channels at each segmentation step in the single grid is not the second superposition type, then the plane combination type of the water channels in the single grid is determined to be an independent combination type; if there are multiple water channels and the superposition type of two water channels at each segmentation step in the single grid is not the second superposition type, then the plane combination type of the water channels in the single grid is determined to be a cross combination type; if there is only one water channel, then the plane combination type of the water channels in the single grid is determined to be a single type; if the number of water channels is zero, then the plane combination type of the water channels in the single grid is determined to be a no-water channel type.

[0020] Optionally, combining the adjacent segmentation step boundary points according to the waterway superposition type of each waterway on each segmentation step and the plane combination type of each waterway in the single grid to obtain the boundary map of each waterway in the single grid includes:

[0021] If the plane combination type corresponding to each waterway in the single grid is the independent combination type, and the superposition type of the waterways on each segmentation step is the first superposition type, then the boundary points of each color block in the first color code are connected with the boundary points of the corresponding color blocks on the adjacent segmentation step;

[0022] If the plane combination type corresponding to each waterway in the single grid is the cross combination type, and the superposition type of the waterways on each segmentation step is not less than at least any two of the first superposition type, the second superposition type, and the third superposition type, then the boundary points of the corresponding color blocks in the first color mark, the boundary points of the corresponding color blocks in the second color mark, and the boundary points of the corresponding color blocks in the third color mark corresponding to the adjacent segmentation steps are connected;

[0023] If the plane combination type corresponding to each waterway in the single grid is the single type, and the superposition type of the waterway on each segmentation step is the third superposition type or the fourth superposition type, then the boundary points of each color block in the third color mark or the boundary points of each color block in the fourth color mark are connected with the boundary points of the corresponding color blocks in the third color mark or the boundary points of the corresponding color blocks in the fourth color mark corresponding to the adjacent segmentation step;

[0024] If the plane combination type corresponding to each water channel in the single grid is the type without water channels, and the superposition type of the water channels in each segmentation step is only the fifth superposition type, no connection is performed.

[0025] Optionally, for each single-grid waterway plane attribute map among the plurality of single-grid waterway plane attribute maps, obtaining the boundary points of each color block in the color code corresponding to the waterway superposition type of each waterway at each segmentation step includes:

[0026] According to the water channel overlapping type and the pre-constructed single water channel boundary point acquisition rule corresponding to the water channel overlapping type, the boundary points of each water channel at each segmentation step are determined.

[0027] Optionally, determining the boundary points of each waterway at each segmentation step according to the waterway overlap type and a pre-constructed rule for obtaining a single waterway boundary point corresponding to the waterway overlap type includes:

[0028] If the water channel overlay type is the first overlay type, the boundary points of each color block in the first color mark are obtained; if the water channel overlay type is the second overlay type, the boundary points of each color block in the second color mark are obtained; if the water channel overlay type is the third overlay type or the fourth overlay type, the boundary points of each color block in the third color mark or the fourth color mark are obtained; if the water channel overlay type is the fifth overlay type, the boundary points of the color blocks are not obtained, wherein the numbers of the obtained boundary points of each color block in the first color mark, the boundary points of each color block in the second color mark, and the boundary points of each color block in the third color mark or the fourth color mark are in a decreasing relationship.

[0029] On the other hand, the present invention provides a system for acquiring a sedimentary channel boundary map based on three-dimensional seismic, characterized in that it includes: a memory; and a processor configured to execute the acquisition method described above.

[0030] Through the above technical scheme, the plane attribute map of the sedimentary channel after the frequency division attribute fusion is divided into a single grid, and then the single grid is divided according to the preset segmentation step, and the actual channel superposition type is assigned to each color code corresponding to the waterway on each segmentation step, and the plane combination type of the waterway in the single grid is determined in combination with the number of waterways, and then the boundary points of the color blocks in the color code are accurately picked up, and then the boundary points are combined according to the waterway superposition type and the plane combination type to obtain the sedimentary channel boundary map, which improves the mapping efficiency and reduces the error probability.

[0031] Other features and advantages of the embodiments of the present invention will be described in detail in the subsequent detailed description. BRIEF DESCRIPTION OF THE DRAWINGS

[0032] The accompanying drawings are used to provide a further understanding of the embodiments of the present invention and constitute a part of the specification. Together with the following specific implementations, they are used to explain the embodiments of the present invention, but do not constitute a limitation on the embodiments of the present invention. In the accompanying drawings:

[0033] Figure 1 It is a flowchart of a method for obtaining a sedimentary channel boundary map based on three-dimensional seismic data provided by an embodiment of the present invention;

[0034] FIG. 2( a ) is a schematic diagram of a plane grid coordinate provided by an embodiment of the present invention;

[0035] FIG2( b ) is a plane attribute diagram of multiple single-grid waterways provided by an embodiment of the present invention;

[0036] Figure 3 It is a waterway plane attribute map provided by an embodiment of the present invention after the D23 single grid is segmented according to a preset segmentation step size;

[0037] Figure 4A schematic diagram of the correspondence between the classification of waterway superposition types and color codes provided by an embodiment of the present invention;

[0038] Figure 5 This is a schematic diagram of counting the number of color marks provided by an embodiment of the present invention;

[0039] Figure 6 It is a schematic diagram of a type of plane combination of waterways in a single grid provided by an embodiment of the present invention;

[0040] Figure 7 It is a schematic diagram of a combination type of waterway planes in all grids provided by an embodiment of the present invention;

[0041] Figure 8 It is a schematic diagram of a rule for obtaining a single waterway boundary point corresponding to a pre-constructed waterway overlapping type provided in an embodiment of the present invention;

[0042] Fig. 9 , 10 is a sedimentary river channel boundary map provided by an embodiment of the present invention;

[0043] Fig.11 This is another flow chart for obtaining a sedimentary channel boundary map provided by an embodiment of the present invention. DETAILED DESCRIPTION

[0044] The specific implementation of the embodiment of the present invention is described in detail below in conjunction with the accompanying drawings. It should be understood that the specific implementation described here is only used to illustrate and explain the embodiment of the present invention, and is not used to limit the embodiment of the present invention.

[0045] See also Figure 1 FIG. 1 is a flowchart of a method for obtaining a sedimentary channel boundary map based on three-dimensional seismic data provided by an embodiment of the present invention, which specifically includes the following implementation steps:

[0046] Step 100: obtaining a plane attribute map of the sedimentary river channel after frequency division attribute fusion, and gridding the plane attribute map of the sedimentary river channel to obtain a plurality of single-grid waterway plane attribute maps.

[0047] In some embodiments, before executing step 100, three-dimensional seismic data bodies can be prepared in advance. First, based on the core-logging-seismic calibration results, seismic data bodies with high reservoir thickness sensitivity and good correlation are selected, and then the three-dimensional data body is divided into small layers vertically. Based on a certain selected three-dimensional seismic data body and according to the logging-seismic calibration results, small-layer vertical stratigraphic strata tracking is performed to establish a minimum time scale stratigraphic framework suitable for the work area, and i small layers are divided, where i is the stratigraphic code, i=1, 2,..., n. In this process, the vertical time window scale should be kept as small as possible to ensure that the number of river channels in space during this period can be represented, so as to facilitate the next step of plane boundary identification, and then clarify the basic distribution characteristics of the plane waterway. With the vertical stratigraphic grid division results as constraints, the plane attribute map of the waterway during the deposition period of each small layer is extracted to clarify the basic information such as the distribution direction, shape, and width of the developed waterway. After the data body is optimized, the example work area uses a frequency division attribute fusion body (i.e., RGB attribute fusion) to characterize the waterway distribution characteristics.

[0048] In some embodiments, when executing step 100, the actual working area situation can be analyzed according to the above embodiment based on the core-logging-seismic calibration results and the frequency division fusion attribute map of the working area. The working area range is about 100km 2 The width of a single water channel is 200m-1000m, and can reach 1300m locally. The small layer plane map (i.e., the plane attribute map of the sedimentary channel after the frequency division attributes are fused) is gridded. As shown in Figure 2(a), a rectangular coordinate system is established along the direction of cutting the water channel. According to the actual width and density of the developed water channels, the number of water channels in each grid in this embodiment does not exceed 2, and they are named Dij row by row and column by column, i=1,2,…, j=1,2,…; as shown in Figure 2(b), according to the width and density of the developed water channels in this work area, the plane attribute map after the frequency division attributes are fused is gridded into 16 grids, and they are named respectively. As shown in Figure 2(b), the names of the grids are D11, D12...D44 from left to right and from top to bottom.

[0049] Step 101: According to a preset segmentation step, each single-grid waterway plane attribute map in the plurality of single-grid waterway plane attribute maps is segmented respectively, and according to a pre-constructed correspondence between the waterway superposition type and the color code, the waterway superposition type of each waterway at each segmentation step is determined.

[0050] In some embodiments, the segmentation step length within a single grid is determined according to the actual development of the waterway and the drawing requirements, see Figure 3 As shown, taking D23 as an example, the grid is divided into 12 steps from top to bottom, k = 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12. In practical applications, the grid can be divided according to a certain step size according to the specific application scenario, and there is no limitation here.

[0051] In some embodiments, the pre-constructed correspondence between the waterway overlay type and the color code includes: if the waterway overlay type is the first overlay type, the corresponding first color code is the first color-second color-first color or the first color-second color-third color; if the waterway overlay type is the second overlay type, the corresponding second color code is the first color-third color / fourth color-first color; if the waterway overlay type is the third overlay type, the corresponding third color code is the first color or the third color; if the waterway overlay type is the fourth overlay type, the corresponding fourth color code is the first color or the third color; if the waterway overlay type is the fifth overlay type, the corresponding fifth color code is the second color, wherein the first color represents the thick sand type, the second color represents the mud type, the third color and the fourth color represent the medium-thick sand type. In this embodiment, the first color is red, the second color is black, the third color is green, and the fourth color is blue. In actual applications, the colors can be set according to the specific application scenarios, which are not limited here.

[0052] In some implementations, before executing step 101 to determine the channel superposition type of each waterway at each segmentation step according to the pre-constructed correspondence between the channel superposition type and the color code, the following steps may also be performed:

[0053] The water channel superposition type is defined according to the number and distribution characteristics of the water channels at each segmentation step.

[0054] In some embodiments, the water channel overlay type is defined according to the number and distribution characteristics of the water channels on each segmentation step, including: if there are multiple single water channels on each segmentation step, and the multiple single water channels are independent of each other, the water channel overlay type is defined as the first overlay type; if there are multiple single water channels on each segmentation step, and the multiple single water channels overlap each other, the water channel overlay type is defined as the second overlay type; if there are multiple single water channels on each segmentation step, and the multiple single water channels merge with each other, the water channel overlay type is defined as the third overlay type; if there is a single single water channel on each segmentation step, the water channel overlay type is defined as the fourth overlay type; if there is no water channel on each segmentation step, the water channel overlay type is defined as the fifth overlay type.

[0055] In some embodiments, the classification of the waterway overlapping type is described by taking two waterways C1 and C2 as an example. Figure 4As shown, if there are two single water channels (C1, C2), there are three superposition modes on the section. The first one is that C1 and C2 are independent of each other, and the type is defined as J1, which can be subdivided into J1a, J1b, and J1c according to the water channel thickness and mutual interval; the second one is that C1 and C2 overlap each other, and the type is defined as J2, which can be subdivided into J2a and J2b according to the water channel thickness and mutual interval; the third one is that one of C1 and C2 contains the other, and the type is defined as J3, which can be subdivided into J3a, J3b, and J3c according to the water channel thickness and relative distance; if there is only one single water channel (C1), then C1 exists independently, and the type is defined as J4, which can be subdivided into J4a and J4b according to the water channel thickness; if there is no single water channel, the number is 0, and the type is defined as J5.

[0056] In some embodiments, the correspondence between the channel overlay type and the color code combination on a single step length is as follows: Figure 4 As shown in the figure, due to the complexity of the actual sedimentation process, thin sand bodies at the edge of the river channel are often missing, so the main body of the thick and medium-thick channels (red and green) are mainly determined; the channel superposition type is J1, corresponding to the color code red-black-red or red-black-green (thick sand-mud-thick sand, thick sand-mud-medium-thick sand); the channel superposition type is J2, corresponding to the color code red-green / blue-red (thick sand-medium-thick sand-thick sand); the channel superposition type is J3, corresponding to the color code red or green (thick sand, medium-thick sand); the channel superposition type is J4, corresponding to the color code red or green (thick sand, medium-thick sand); the channel superposition type is J5, corresponding to the color code black.

[0057] In some embodiments, the number of waterways in the single grid can be determined according to the following conditions:

[0058] Condition 1: If the number of color blocks in the color mark in the single grid is greater than a first preset threshold, then the number of the waterways is multiple;

[0059] Condition 2: If the number of color blocks in the color code in the single grid is greater than the second preset threshold and less than or equal to the first preset threshold, and the color of the color blocks is not all the second color representing the mud type, then the number of the waterway is single;

[0060] Condition 3: If the number of color blocks in the color mark within the single grid is equal to the second preset threshold, and the color of the color blocks is the second color, then the number of waterways is zero, wherein the first preset threshold is greater than the second preset threshold.

[0061] In some embodiments, see Figure 5 As shown, if the number of color blocks is greater than 3, the number of water channels is equal to 2; when the number of color blocks is greater than 1 and less than or equal to 3, and the colors are not all black, the number of water channels is equal to 1; when the mean number of color blocks contains only 1 color block, and the color is black, the number of water channels is equal to 0.

[0062] Step 102: For each single-grid waterway plane attribute map among the multiple single-grid waterway plane attribute maps, determine the plane combination type of each waterway in the single-grid waterway plane attribute map according to the number of waterways in the single grid and the waterway superposition type of each waterway at each segmentation step in the single grid.

[0063] In some implementations, when executing step 102 to determine the plane combination type of each waterway in the single-grid waterway plane attribute map, the determination may be made according to the following conditions:

[0064] Condition 1: if there are multiple waterways, and at each segmentation step in the single grid, the superposition type of two waterways is not the second superposition type, then the plane combination type of the waterways in the single grid is determined to be an independent combination type;

[0065] Condition 2: if there are multiple waterways, and at each segmentation step in the single grid, the superposition type of two waterways is not the second superposition type, then the plane combination type of the waterways in the single grid is determined to be a cross combination type;

[0066] Condition 3: If the number of the water channel is single, then the plane combination type of the water channel in the single grid is determined to be a single type;

[0067] Condition 4: If the number of the water channels is zero, the plane combination type of the water channels in the single grid is determined to be a type without water channels.

[0068] In some embodiments, taking two waterways as an example, the plane combination type of the waterways in a single grid is shown in Figure 6 As shown in the figure, in a single grid, there are several types of combinations of two or less waterways:

[0069] There are two single water channels (C1, C2), and there are two combination modes on the plane. The first one is that C1 and C2 do not intersect, and the type is defined as M1. According to the difference in water channel position and width, it can be subdivided into M1a, M1b, and M1c; the second one is that C1 and C2 cross-cut and overlap, and the type is defined as M2, which can be subdivided into M2a and M2b according to the crossing method; if there is only one single water channel (C1), C1 exists independently, and the type is defined as M3, which can be subdivided into M3a, M3b, and M3c according to the difference in water channel position and width; if there is no single water channel, the number is 0, and the type is defined as M4.

[0070] In some embodiments, Figure 3 As shown, the segmentation step k = 1, 2...8, the number of color blocks is 5; the segmentation step k = 9, 10...12, the number of color blocks is 3, the color blocks are greater than 3, so it is determined that the number of D23 waterways is equal to 2; and according to D23, when k = 8, it is J2 type, so it is determined to be type M2 cross-type waterway combination.

[0071] In some embodiments, the plane combination type determination result of each waterway in the plane attribute map of all single-grid waterways in the work area is as follows: Figure 7 shown.

[0072] Step 103: For each single-grid waterway plane attribute map among the multiple single-grid waterway plane attribute maps, obtain the boundary points of each color block in the color code corresponding to the waterway overlay type of each waterway on each segmentation step, and combine the adjacent segmentation step boundary points according to the waterway overlay type of each waterway on each segmentation step and the plane combination type of each waterway in the single grid to obtain the boundary map of each waterway in the single grid.

[0073] In some implementations, when executing step 103, a single-grid waterway boundary map may be obtained according to the following steps:

[0074] S1031: If the plane combination type corresponding to each waterway in the single grid is the independent combination type, and the superposition type of the waterway on each segmentation step is the first superposition type, then the boundary points of each color block in the first color code are connected with the boundary points of the corresponding color blocks on the adjacent segmentation step.

[0075] In some embodiments, if the plane combination type corresponding to each water channel in the single grid is the independent combination type M1, and the superposition types of all water channels on the segmentation step k are all of type J1, then the boundary points (Dij, k-1, x1) and (Dij, k, x1), (Dij, k-1, x2) and (Dij, k, x2), (Dij, k-1, x3) and (Dij, k, x3), and (Dij, k, x4) and (Dij, k, x4) are connected respectively.

[0076] S1032: If the plane combination type corresponding to each waterway in the single grid is the cross combination type, and the superposition type of the waterway on each segmentation step is not less than at least any two of the first superposition type, the second superposition type and the third superposition type, then the boundary points of the corresponding color blocks in the first color mark, the boundary points of the corresponding color blocks in the second color mark and the boundary points of the corresponding color blocks in the third color mark corresponding to the adjacent segmentation steps are connected.

[0077] In some embodiments, if the plane combination type corresponding to each waterway in the single grid is the cross combination type M2, the segmentation step k includes at most three waterway overlapping types J1, J2 and J3, and no less than two types, such as Figure 6 As shown, there are two combination methods. The proportion of cross-replacement type and cross-separation type can be set according to the development of waterways in the work area, and then one type is randomly selected to continue the following steps;

[0078] If it is M2a type, the combination is as follows:

[0079] When the water channel superposition type of segmentation step k-1 is determined to be J1 and the water channel superposition type of segmentation step k is determined to be J2, when the two are combined, J1(Dij, k-1, x1) is connected with J2(Dij, k, x1), J1(Dij, k-1, x2) is connected with J2(Dij, k, x2), J1(Dij, k-1, x3) is connected with J2(Dij, k, x3), and J1(Dij, k, x4) is connected with J2(Dij, k, x4);

[0080] When the water channel superposition type of segmentation step k-1 is determined to be J2 and the water channel superposition type of segmentation step k is determined to be J3, when the two are combined, J2(Dij, k-1, x1) is connected with J3(Dij, k, x1), J2(Dij, k-1, x2) is connected with J3(Dij, k, x2), J2(Dij, k-1, x3) is connected with J3(Dij, k, x1), and J2(Dij, k, x4) is connected with J3(Dij, k, x2);

[0081] When the water channel superposition type of segmentation step k-1 is determined to be J3 and the water channel superposition type of segmentation step k is determined to be J2, when the two are combined, J3(Dij, k-1, x1) is connected with J2(Dij, k, x1), J3(Dij, k-1, x2) is connected with J2(Dij, k, x2), J3(Dij, k-1, x1) is connected with J2(Dij, k, x3), and J3(Dij, k, x2) is connected with J2(Dij, k, x4) respectively;

[0082] If it is M2b type, the combination is as follows:

[0083] When the water channel superposition type of segmentation step k-1 is determined to be J1 and the water channel superposition type of segmentation step k is determined to be J2, when the two are combined, J1(Dij, k-1, x1) is connected with J2(Dij, k, x1), J1(Dij, k-1, x2) is connected with J2(Dij, k, x2), J1(Dij, k-1, x3) is connected with J2(Dij, k, x3), and J1(Dij, k, x4) is connected with J2(Dij, k, x4);

[0084] When the water channel superposition type of segmentation step k-1 is determined to be J2 and the water channel superposition type of segmentation step k is determined to be J3, when the two are combined, J2(Dij, k-1, x1) is connected with J3(Dij, k, x1), J2(Dij, k-1, x2) is connected with J3(Dij, k, x2), J2(Dij, k-1, x3) is connected with J3(Dij, k, x1), and J2(Dij, k, x4) is connected with J3(Dij, k, x2);

[0085] When the water channel overlay type with segmentation step k-1 is determined to be J3 and the water channel overlay type with segmentation step k is determined to be J2, when the two are combined, J3(Dij, k-1, x1) and J2(Dij, k, x3), J3(Dij, k-1, x2) and J2(Dij, k, x4), J3(Dij, k-1, x1) and J2(Dij, k, x3), J3(Dij, k, x2) and J2(Dij, k, x4) are connected respectively.

[0086] S1033: If the plane combination type corresponding to each waterway in the single grid is the single type, and the overlay type of the waterway on each segmentation step is the third overlay type or the fourth overlay type, then the boundary points of each color block in the third color mark or the boundary points of each color block in the fourth color mark are connected with the boundary points of the corresponding color blocks in the third color mark or the boundary points of the corresponding color blocks in the fourth color mark on the adjacent segmentation step.

[0087] In some embodiments, if the plane combination type corresponding to each water channel in the single grid is a single type M3: the segmentation step k includes one of the water channel superposition types J3 or J4, the boundary points (Dij, k-1, x1) and (Dij, k, x1), (Dij, k-1, x2) and (Dij, k, x2) are connected respectively.

[0088] S1034: If the plane combination type corresponding to each water channel in the single grid is the type without water channels, and the superposition type of the water channels in each segmentation step is only the fifth superposition type, no connection is performed.

[0089] In some embodiments, if the plane combination type corresponding to each water channel in the single grid is the non-water channel type M4, and the superposition type of each water channel at the segmentation step k is a type J5, there is no connection line.

[0090] In some implementations, before executing step 103, for each single-grid waterway plane attribute map in the plurality of single-grid waterway plane attribute maps, obtaining the boundary points of each color block in the color code corresponding to the waterway superposition type of each waterway at each segmentation step includes:

[0091] According to the water channel overlapping type and the pre-constructed single water channel boundary point acquisition rule corresponding to the water channel overlapping type, the boundary points of each water channel at each segmentation step are determined.

[0092] In some implementations, the boundary points of each waterway at each segmentation step may be determined by the following acquisition rules, including:

[0093] Rule 1: If the waterway overlay type is the first overlay type, then obtain the boundary points of each color block in the first color code;

[0094] Rule 2: If the waterway overlay type is the second overlay type, then obtain the boundary points of each color block in the second color code;

[0095] Rule 3: If the waterway superposition type is the third superposition type or the fourth superposition type, then obtain the boundary points of each color block in the third color code or the fourth color code;

[0096] Rule 4: If the waterway overlay type is the fifth overlay type, the boundary points of the color blocks are not obtained, wherein the numbers of the boundary points of each color block in the first color code, the boundary points of each color block in the second color code, and the boundary points of each color block in the third color code or the fourth color code are in a decreasing relationship.

[0097] In some embodiments, see Figure 8 As shown, when the water channel overlay type is J1, 3-4 boundary points of the three color blocks are picked up in sequence, which are defined as (Dij, k, x1), (Dij, k, x2), (Dij, k, x3), and (Dij, k, x4), respectively, where x2 can be equal to x3; when the water channel overlay type is J2, 3-4 boundary points of the three color blocks are picked up in sequence, which are defined as (Dij, k, x1), (Dij, k, x3), (Dij, k, x2), and (Dij, k, x4), respectively; when the water channel overlay type is J3 and J4, 2 boundary points of the two color blocks are picked up in sequence, which are defined as (Dij, k, x1) and (Dij, k, x2), respectively; when the water channel overlay type is J5, 0 boundary points of the color blocks are picked up.

[0098] Step 104: Combine the boundaries of each waterway in adjacent plane single grids to obtain a sedimentary river channel boundary map.

[0099] In some embodiments, the boundaries of each waterway in the adjacent plane single grids in the y-axis direction can be combined first, and then the boundaries of each waterway in the adjacent plane single grids in the x-axis direction can be combined to obtain a sedimentary river boundary map, such as Fig. 9 , Fig.10 shown.

[0100] In some embodiments, see Fig.11As shown, it is another flow chart for obtaining a sedimentary channel boundary map provided by an embodiment of the present invention. First, a three-dimensional seismic data body is optimized and prepared, see step 100, and then plane grid division (Dnm) is performed on it, and the segmentation step (K) is determined. For a single step, the picked color code boundary points, color code combination patterns and the number of color blocks are output one by one, and the channel superposition type (Jn) is determined according to the color code boundary points and the color code combination patterns, see step 101, and a single channel boundary point is defined, and the channel plane combination type (Mn) is determined in combination with the channel superposition type (Jn) and the number of color blocks. According to the single channel boundary point and the channel plane combination type (Mn), see step 102, adjacent step points are combined one by one by determining the boundary point plane combination method, see step 103, and then the grid mosaic combination is performed to finally obtain a complex channel boundary map, see step 104, wherein each Dnm is independently run from determining the segmentation step (K) to combining adjacent step points one by one, see steps 101-103.

[0101] On the other hand, the present invention provides a system for acquiring a sedimentary channel boundary map based on three-dimensional seismic, characterized in that it includes: a memory; and a processor, configured to execute the acquisition method described in any of the above embodiments.

[0102] Those skilled in the art will appreciate that the embodiments of the present application may be provided as methods, systems, or computer program products. Therefore, the present application may adopt the form of a complete hardware embodiment, a complete software embodiment, or an embodiment in combination with software and hardware. Moreover, the present application may adopt the form of a computer program product implemented in one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) that include computer-usable program code.

[0103] The present application is described with reference to the flowcharts and / or block diagrams of the methods, devices (systems), and computer program products according to the embodiments of the present application. It should be understood that each process and / or box in the flowchart and / or block diagram, as well as the combination of the processes and / or boxes in the flowchart and / or block diagram, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing device to generate a machine, so that the instructions executed by the processor of the computer or other programmable data processing device generate instructions for implementing the processes in the flowchart and / or block diagram. Figure 1 A process or multiple processes and / or boxes Figure 1 A device that provides the functions specified in a block or multiple blocks.

[0104] These computer program instructions may also be stored in a computer-readable memory capable of directing a computer or other programmable data processing device to operate in a specific manner, so that the instructions stored in the computer-readable memory produce an article of manufacture comprising an instruction device, which implements the process Figure 1 A process or multiple processes and / or boxes Figure 1 A function specified in one or more boxes.

[0105] These computer program instructions can also be loaded onto a computer or other programmable data processing device so that a series of operating steps are executed on the computer or other programmable device to produce a computer-implemented process, thereby providing instructions for implementing the process. Figure 1 A process or multiple processes and / or boxes Figure 1 A step that specifies a function in one or more boxes.

[0106] In a typical configuration, a computing device includes one or more processors (CPU), input / output interfaces, network interfaces, and memory.

[0107] The memory may include non-permanent memory in a computer-readable medium, random access memory (RAM) and / or non-volatile memory in the form of read-only memory (ROM) or flash RAM. The memory is an example of a computer-readable medium.

[0108] Computer readable media include permanent and non-permanent, removable and non-removable media that can be implemented by any method or technology to store information. Information can be computer readable instructions, data structures, program modules or other data. Examples of computer storage media include, but are not limited to, phase change memory (PRAM), static random access memory (SRAM), dynamic random access memory (DRAM), other types of random access memory (RAM), read-only memory (ROM), electrically erasable programmable read-only memory (EEPROM), flash memory or other memory technology, compact disk read-only memory (CD-ROM), digital versatile disk (DVD) or other optical storage, magnetic cassettes, magnetic tape magnetic disk storage or other magnetic storage devices or any other non-transmission media that can be used to store information that can be accessed by a computing device. As defined herein, computer readable media does not include temporary computer readable media (transitory media), such as modulated data signals and carrier waves.

[0109] It should also be noted that the terms "include", "comprises" or any other variations thereof are intended to cover non-exclusive inclusion, so that a process, method, commodity or device including a series of elements includes not only those elements, but also other elements not explicitly listed, or also includes elements inherent to such process, method, commodity or device. In the absence of more restrictions, the elements defined by the sentence "comprises a ..." do not exclude the existence of other identical elements in the process, method, commodity or device including the elements.

[0110] The above are only embodiments of the present application and are not intended to limit the present application. For those skilled in the art, the present application may have various changes and variations. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present application should be included within the scope of the claims of the present application.

Claims

1. A method for obtaining a sedimentary channel boundary map based on three-dimensional seismic, characterized in that: include: Obtaining a plane attribute map of the sedimentary river channel after frequency division attribute fusion, and gridding the plane attribute map of the sedimentary river channel to obtain a plurality of single-grid waterway plane attribute maps; According to a preset segmentation step, each single-grid waterway plane attribute map in the plurality of single-grid waterway plane attribute maps is segmented respectively, and according to a pre-constructed correspondence between the waterway superposition type and the color scale, the waterway superposition type of each waterway at each segmentation step is determined; For each single-grid waterway plane attribute map among the plurality of single-grid waterway plane attribute maps, determining the plane combination type of each waterway in the single-grid waterway plane attribute map according to the number of waterways in the single grid and the waterway superposition type of each waterway at each segmentation step in the single grid; For each single-grid waterway plane attribute map among the multiple single-grid waterway plane attribute maps, the boundary points of each color block in the color code corresponding to the waterway superposition type of each waterway on each segmentation step are obtained, and according to the waterway superposition type of each waterway on each segmentation step and the plane combination type of each waterway in the single grid, the boundary points of adjacent segmentation steps are combined to obtain the boundary map of each waterway in the single grid; and The boundaries of each channel in adjacent plane single grids are combined to obtain the sedimentary channel boundary map.

2. The acquisition method according to claim 1, characterized in that: The correspondence between the pre-constructed waterway overlay types and color codes includes: If the waterway superposition type is the first superposition type, the corresponding first color code is the first color-the second color-the first color or the first color-the second color-the third color; If the waterway superposition type is the second superposition type, the corresponding second color code is the first color-third color / fourth color-first color; If the waterway superposition type is the third superposition type, the corresponding third color code is the first color or the third color; If the waterway superposition type is the fourth superposition type, the corresponding fourth color code is the first color or the third color; If the waterway superposition type is the fifth superposition type, the corresponding fifth color code is the second color, wherein the first color represents the thick sand type, the second color represents the mud type, and the third color and the fourth color represent the medium-thick sand type.

3. The acquisition method according to claim 1, characterized in that: Before determining the water channel superposition type of each water channel at each segmentation step according to the pre-constructed correspondence between the water channel superposition type and the color code, the acquisition method further includes: The water channel superposition type is defined according to the number and distribution characteristics of the water channels at each segmentation step.

4. The acquisition method according to claim 3, characterized in that: Defining the waterway superposition type according to the number and distribution characteristics of the waterways at each segmentation step includes: If there are multiple single water channels on each segmentation step, and the multiple single water channels are independent of each other, then the water channel superposition type is defined as the first superposition type; If there are multiple single water channels on each segmentation step, and the multiple single water channels overlap each other, then the water channel overlap type is defined as the second overlap type; If there are multiple single water channels on each segmentation step, and the multiple single water channels merge with each other, the water channel superposition type is defined as the third superposition type; If there is a single water channel on each segmentation step, the water channel superposition type is defined as the fourth superposition type; If no water channel exists on each segmentation step, the water channel superposition type is defined as the fifth superposition type.

5. The acquisition method according to claim 1, characterized in that: Determining the number of waterways in a single grid includes: If the number of color blocks in the color mark in the single grid is greater than a first preset threshold, then the number of waterways is multiple; If the number of color blocks in the color code in the single grid is greater than the second preset threshold and less than or equal to the first preset threshold, and the color of the color blocks is not all the second color representing the mud type, then the number of the waterway is single; If the number of color blocks in the color mark within the single grid is equal to a second preset threshold, and the color of the color blocks is the second color, then the number of waterways is zero, wherein the first preset threshold is greater than the second preset threshold.

6. The acquisition method according to claim 2, characterized in that: Determining the plane combination type of each water channel in the single grid water channel plane attribute graph according to the number of water channels in the single grid and the water channel superposition type of each water channel at each segmentation step in the single grid includes: If there are multiple waterways, and the superposition type of two waterways in each segmentation step in the single grid is not the second superposition type, then determining the plane combination type of the waterways in the single grid to be an independent combination type; If there are multiple waterways, and the superposition type of two waterways in each segmentation step in the single grid is not the second superposition type, then determining the plane combination type of the waterways in the single grid to be a cross combination type; If the number of the water channel is single, determining the plane combination type of the water channel in the single grid to be a single type; If the number of the water channels is zero, it is determined that the plane combination type of the water channels in the single grid is a type without water channels.

7. The acquisition method according to claim 6, characterized in that: The method combines the adjacent segmentation step boundary points according to the waterway superposition type of each waterway on each segmentation step and the plane combination type of each waterway in the single grid to obtain the boundary map of each waterway in the single grid, including: If the plane combination type corresponding to each waterway in the single grid is the independent combination type, and the superposition type of the waterways on each segmentation step is the first superposition type, then the boundary points of each color block in the first color code are connected with the boundary points of the corresponding color blocks on the adjacent segmentation step; If the plane combination type corresponding to each waterway in the single grid is the cross combination type, and the superposition type of the waterways on each segmentation step is not less than at least any two of the first superposition type, the second superposition type, and the third superposition type, then the boundary points of the corresponding color blocks in the first color mark, the boundary points of the corresponding color blocks in the second color mark, and the boundary points of the corresponding color blocks in the third color mark corresponding to the adjacent segmentation steps are connected; If the plane combination type corresponding to each waterway in the single grid is the single type, and the superposition type of the waterway on each segmentation step is the third superposition type or the fourth superposition type, then the boundary points of each color block in the third color mark or the boundary points of each color block in the fourth color mark are connected with the boundary points of the corresponding color blocks in the third color mark or the boundary points of the corresponding color blocks in the fourth color mark corresponding to the adjacent segmentation step; If the plane combination type corresponding to each water channel in the single grid is the type without water channels, and the superposition type of the water channels in each segmentation step is only the fifth superposition type, no connection is performed.

8. The acquisition method according to claim 2, characterized in that: The step of obtaining, for each single-grid waterway plane attribute map in the plurality of single-grid waterway plane attribute maps, the boundary points of each color block in the color scale corresponding to the waterway superposition type of each waterway at each segmentation step includes: According to the water channel overlapping type and the pre-constructed single water channel boundary point acquisition rule corresponding to the water channel overlapping type, the boundary points of each water channel at each segmentation step are determined.

9. The acquisition method according to claim 8, characterized in that: The step of determining the boundary points of each waterway at each segmentation step according to the waterway overlap type and the pre-constructed single waterway boundary point acquisition rule corresponding to the waterway overlap type comprises: If the waterway superposition type is the first superposition type, obtaining the boundary points of each color block in the first color code; If the waterway overlay type is the second overlay type, obtaining the boundary points of each color block in the second color code; If the waterway superposition type is the third superposition type or the fourth superposition type, obtaining the boundary points of each color block in the third color code or the fourth color code; If the waterway overlay type is the fifth overlay type, the boundary points of the color blocks are not obtained, wherein the numbers of the boundary points of each color block in the first color code, the boundary points of each color block in the second color code, and the boundary points of each color block in the third color code or the fourth color code are in a decreasing relationship.

10. A system for acquiring sedimentary channel boundary maps based on three-dimensional seismic, characterized in that: include: Memory; as well as A processor, configured to execute the acquisition method as described in any one of claims 1-9.

Images