Slope extraction method and device based on geologic horizon guidance

Abstract

The invention provides a slope extraction method and device based on geologic horizon guidance, and belongs to the technical field of geophysical exploration, and the method comprises the steps: obtaining three-dimensional seismic data; and extracting the slope of the nth horizon, the slope of the (n-1) th horizon and the slope between the (n-1) th horizon and the nth horizon by using the three-dimensional seismic data. The device comprises a data acquisition module and a slope extraction module. According to the method, the accuracy of slope extraction is improved, and a reliable data basis is provided for subsequent seismic data processing or seismic data interpretation.

Description

Technical Field

[0001] The present application belongs to the field of geophysical exploration technology, and specifically relates to a slope extraction method and device based on geological layer guidance. Background Art

[0002] In the process of seismic data processing, slope information is crucial. Conventional methods can only obtain the slope of a certain layer, and usually cannot obtain the slope information corresponding to the entire seismic data. Even if the slope information of the entire data body is obtained, the conventional method lacks layer constraints, resulting in low extraction accuracy and cannot meet subsequent accuracy requirements. Therefore, improving the accuracy of slope extraction and providing a reliable data basis for subsequent seismic data processing or seismic data interpretation is a technical problem that needs to be solved urgently. Summary of the invention

[0003] Based on the above technical problems, the present application proposes a slope extraction method and device based on geological layer guidance, which improves the accuracy of slope extraction and provides a reliable data basis for subsequent seismic data processing or seismic data interpretation.

[0004] In the first aspect, the present application proposes a slope extraction method based on geological layer guidance, comprising:

[0005] Acquire 3D seismic data;

[0006] The three-dimensional seismic data is used to extract the slope of the nth layer, the slope of the n-1th layer, and the slope between the n-1th layer and the nth layer of the corresponding point;

[0007] The slope of the nth layer, the slope of the n-1th layer, and the slope between the n-1th layer and the nth layer of all points are constructed as the slope of the three-dimensional seismic data of the entire formation.

[0008] The step of extracting the slope of the nth layer and the slope of the n-1th layer of the corresponding point includes:

[0009] Performing seismic horizon interpretation on the three-dimensional seismic data to obtain the horizon information of the nth horizon and the horizon information of the n-1th horizon;

[0010] The inverse tangent function value is taken for the difference between the layer information of two adjacent points in the same layer, and the slope of the nth layer and the slope of the n-1th layer of the corresponding point are obtained.

[0011] The x-axis direction of the three-dimensional seismic data is data in the direction of the seismic main survey line, the y-axis direction is data in the direction of the seismic contact survey line, and the z-axis direction is data in the direction of the seismic depth.

[0012] The slope between the n-1th layer and the nth layer is extracted by:

[0013] According to the layer information of the nth layer and the layer information of the n-1th layer, a double-parameter inverse distance weighted method is used to calculate the up-sampling coefficient on the non-layer and the down-sampling coefficient on the non-layer;

[0014] According to the up-sampling coefficient, the down-sampling coefficient, the layer information of the n-th layer and the layer information of the n-1-th layer, the slope between the n-1-th layer and the n-th layer of the corresponding point is obtained.

[0015] The inverse tangent function value is taken for the difference between the layer information of two adjacent points in the same layer to obtain the slope of the nth layer and the slope of the n-1th layer of the corresponding point. The calculation formula is as follows:

[0016] S(x,y,z layer(x,y|n) )=arctan[Layer(x,y|n)-Layer(x,y-1|n)]

[0017] Among them, S(x,y,z layer(x,y|n) ) is the slope of point (x, y) in the nth layer, Layer(x, y|n) is the layer information of point (x, y) in the nth layer, Layer(x, y-1|n) is the layer information of point (x, y-1) in the nth layer, and point (x, y) and point (x, y-1) are two adjacent points in the same layer.

[0018] The upsampling coefficient on the non-layer is calculated as follows:

[0019]

[0020] Among them, A up (x,y,z n-1,n ) is the upsampling coefficient on the non-layer, z n-1,n It is the data in the seismic depth direction between the n-1th layer and the nth layer. Layer(x,y|n) is the layer information of the point (x,y) in the nth layer. Layer(x,y|n-1) is the layer information of the point (x,y) in the n-1th layer.

[0021] The downsampling coefficient on the non-layer is calculated as follows:

[0022]

[0023] Among them, A down (x,y,z n-1,n ) is the downsampling coefficient on the non-layer, z n-1,n It is the data in the seismic depth direction between the n-1th layer and the nth layer. Layer(x,y|n) is the layer information of the point (x,y) in the nth layer. Layer(x,y|n-1) is the layer information of the point (x,y) in the n-1th layer.

[0024] According to the upsampling coefficient, the downsampling coefficient, the layer information of the nth layer and the layer information of the n-1th layer, the slope between the n-1th layer and the nth layer of the corresponding point is obtained, and the calculation formula is as follows:

[0025] S(x,y,z n-1,n )=A up (x,y,z n-1,n )×S(x,y,z layer(x,y|n-1) )+A down (x,y,z n-1,n )×S(x,y,z layer(x,y|n) )

[0026] Among them, S(x,y,z n-1,n ) is the slope between the n-1th layer and the nth layer, A up (x,y,z n-1,n ) is the upsampling coefficient on the non-layer, A down (x,y,z n-1,n ) is the downsampling coefficient on the non-layer, S(x,y,z layer(x,y|n) ) is the slope of the point (x, y) at the nth layer, S(x, y, z layer(x,y|n-1) ) is the slope of the point (x, y) at the n-1th layer.

[0027] In the second aspect, the present application proposes a slope extraction device based on geological layer guidance, comprising:

[0028] A data acquisition module, used for acquiring three-dimensional seismic data;

[0029] A slope extraction module, used to extract the slope of the nth layer, the slope of the n-1th layer, and the slope between the n-1th layer and the nth layer using the three-dimensional seismic data;

[0030] The slope construction module is used to construct the slope of the nth layer, the slope of the n-1th layer, and the slope between the n-1th layer and the nth layer of all points into the slope of the three-dimensional seismic data of the entire formation.

[0031] In a third aspect, the present application proposes an electronic device, comprising: a memory and a processor, wherein the memory stores a computer program, and when the computer program is executed by the processor, the slope extraction method based on geological layer guidance is executed.

[0032] In a fourth aspect, the present application proposes a computer-readable storage medium storing executable instructions, which, when executed, enable a processor to execute the slope extraction method based on geological layer guidance.

[0033] Beneficial effects:

[0034] The present application proposes a slope extraction method and device based on geological horizon guidance. The method of the present application calculates the slope of the data corresponding to the horizon position according to the horizon information. For non-horizontal positions, the slope of each point in the non-horizontal position is obtained by using the horizon information of the nth horizon and the horizon information of the n-1th horizon, and then the slope of the entire data body is obtained by using a double-parameter inverse distance weighted method. The accuracy of slope extraction is improved, providing a reliable data basis for subsequent seismic data processing or seismic data interpretation. BRIEF DESCRIPTION OF THE DRAWINGS

[0035] Figure 1 A flow chart of a slope extraction method based on geological stratum guidance according to an embodiment of the present application;

[0036] Figure 2 This is a flow chart of the slope of the nth layer and the slope of the n-1th layer in an embodiment of the present application;

[0037] Figure 3 This is a flow chart of the slope between the n-1th layer and the nth layer in an embodiment of the present application;

[0038] Figure 4 This is a specific example flow chart of the slope extraction method based on geological layer guidance in this application;

[0039] Figure 5 A schematic diagram of three-dimensional seismic data according to an embodiment of the present application;

[0040] Figure 6 This is a schematic diagram of layer information in an embodiment of the present application;

[0041] Figure 7 A schematic diagram of the slope of three-dimensional seismic data of the entire formation according to an embodiment of the present application;

[0042] Figure 8 This is a principle block diagram of a slope extraction device based on geological layer guidance according to an embodiment of the present application. DETAILED DESCRIPTION

[0043] The present disclosure is further described below in conjunction with the embodiments shown in the accompanying drawings.

[0044] In the process of seismic data processing, slope information is crucial. Conventional methods can only obtain the slope of a certain layer, and usually cannot obtain the slope information corresponding to the entire seismic data. Even if the slope information of the entire data body is obtained, the conventional method lacks stratigraphic constraints, resulting in low extraction accuracy and cannot meet subsequent accuracy requirements.

[0045] The present application proposes a slope extraction method and device based on geological horizon guidance. First, the three-dimensional seismic data is interpreted in horizon to obtain the interpreted horizon information. On this basis, the slope of the data corresponding to the horizon position is calculated according to the horizon information. For non-horizontal positions, the slopes calculated by the upper and lower horizons are used for inverse distance weighting to obtain the slope of each point in the non-horizontal position, and then the horizon slope of the entire data body is obtained, which improves the accuracy of slope extraction and provides a reliable data basis for subsequent seismic data processing or seismic data interpretation.

[0046] Embodiment 1

[0047] This embodiment proposes a slope extraction method based on geological layer guidance, such as Figure 1 As shown, including:

[0048] Step S1: Acquire three-dimensional seismic data;

[0049] In this embodiment, the three-dimensional seismic data is three-dimensional seismic data acquired under a three-dimensional coordinate system consisting of an x-axis as the horizontal axis, a y-axis as the vertical axis, and a z-axis as the vertical axis, wherein the x-axis direction of the three-dimensional seismic data is data in the direction of the main seismic survey line, the y-axis direction is data in the direction of the seismic connecting survey line, and the z-axis direction is data in the direction of the seismic depth.

[0050] Step S2: using the three-dimensional seismic data to extract the slope of the nth layer, the slope of the n-1th layer, and the slope between the n-1th layer and the nth layer.

[0051] Step S3: construct the slope of the nth layer, the slope of the n-1th layer, and the slope between the n-1th layer and the nth layer of all points into the slope of the three-dimensional seismic data of the entire formation.

[0052] In a specific implementation, the slope of the nth layer and the slope of the n-1th layer are as follows: Figure 2 As shown, the extraction process is as follows:

[0053] Step S211: performing seismic horizon interpretation on the three-dimensional seismic data to obtain the horizon information of the nth horizon and the horizon information of the n-1th horizon;

[0054] Step S212: Take the inverse tangent function value of the difference between the layer information of two adjacent points in the same layer to obtain the slope of the nth layer and the slope of the n-1th layer of the corresponding point. The calculation formula is as follows:

[0055] S(x,y,z layer(x,y|n) )=arctan[Layer(x,y|n)-Layer(x,y-1|n)]

[0056] Among them, S(x,y,z layer(x,y|n)) is the slope of point (x, y) in the nth layer, Layer(x, y|n) is the layer information of point (x, y) in the nth layer, Layer(x, y-1|n) is the layer information of point (x, y-1) in the nth layer, and point (x, y) and point (x, y-1) are two adjacent points in the same layer (i.e., the nth layer).

[0057] Through the slope of the nth layer, it is easy to get the slope of the n-1th layer of the corresponding point. The calculation formula is as follows:

[0058] S(x,y,z layer(x',y'|n-1) )=arctan[Layer(x',y'|n-1)-Layer(x',y'-1|n-1)]

[0059] Among them, S(x,y,z layer(x',y'|n-1) ) is the slope of the point (x',y') at the n-1th layer, Layer(x',y'|n-1) is the layer information of the point (x',y') at the nth layer, Layer(x',y'-1|n) is the layer information of the point (x',y'-1) at the n-1th layer, and the point (x',y') and the point (x',y'-1) are two adjacent points in the same layer (i.e., the n-1th layer).

[0060] In a specific implementation, the slope between the n-1th layer and the nth layer is as follows: Figure 3 As shown, the extraction process is as follows:

[0061] Step S221: according to the layer information of the nth layer and the layer information of the n-1th layer, a dual-parameter inverse distance weighted method is used to calculate the up-sampling coefficient on the non-layer and the down-sampling coefficient on the non-layer;

[0062] The upsampling coefficient on the non-layer is calculated as follows:

[0063]

[0064] Among them, A up (x,y,z n-1,n ) is the upsampling coefficient on the non-layer, z n-1,n It is the data in the seismic depth direction between the n-1th layer and the nth layer. Layer(x,y|n) is the layer information of the point (x,y) in the nth layer. Layer(x,y|n-1) is the layer information of the point (x,y) in the n-1th layer.

[0065] The downsampling coefficient on the non-layer is calculated as follows:

[0066]

[0067] Among them, Adown (x,y,z n-1,n ) is the downsampling coefficient on the non-layer, z n-1,n It is the data in the seismic depth direction between the n-1th layer and the nth layer. Layer(x,y|n) is the layer information of the point (x,y) in the nth layer. Layer(x,y|n-1) is the layer information of the point (x,y) in the n-1th layer.

[0068] Step S222: According to the upsampling coefficient, the downsampling coefficient, the layer information of the nth layer and the layer information of the n-1th layer, the slope between the n-1th layer and the nth layer of the corresponding point is obtained, and the calculation formula is as follows:

[0069] S(x,y,z n-1,n )=A up (x,y,z n-1,n )×S(x,y,z layer(x,y|n-1) )+A down (x,y,z n-1,n )×S(x,y,z layer(x,y|n) )

[0070] Among them, S(x,y,z n-1,n ) is the slope between the n-1th layer and the nth layer, A up (x,y,z n-1,n ) is the upsampling coefficient on the non-layer, A down (x,y,z n-1,n ) is the downsampling coefficient on the non-layer, S(x,y,z layer(x,y|n) ) is the slope of the point (x, y) at the nth layer, S(x, y, z layer(x,y|n-1) ) is the slope of the point (x, y) at the n-1th layer.

[0071] In the specific implementation, it is necessary to calculate the slope S(x, y, z) of the nth layer of all points. layer(x,y|n) ), the slope S(x,y,z layer(x',y'|n-1) ), and the slope S(x,y,z) between the n-1th layer and the nth layer n-1,n ) is constructed as the slope of the 3D seismic data for the entire formation.

[0072] The present embodiment proposes a slope extraction method based on geological layer guidance. First, three-dimensional seismic data is acquired; the slope of the nth layer, the slope of the n-1th layer, and the slope between the n-1th layer and the nth layer are extracted using the three-dimensional seismic data. The slope of the nth layer, the slope of the n-1th layer, and the slope between the n-1th layer and the nth layer of all points are constructed as the slope of the three-dimensional seismic data of the entire formation. The slope extracted by the method of the present embodiment has high accuracy.

[0073] Embodiment 2

[0074] This embodiment is a specific example of the first embodiment. The slope extraction method based on geological layer guidance is described in detail using actual seismic data, and the effectiveness of the method proposed in this embodiment is demonstrated. Figure 4 As shown, including:

[0075] Step S300: Acquire three-dimensional seismic data;

[0076] In this embodiment, the x-axis direction of the three-dimensional seismic data is data in the direction of the seismic main survey line, the y-axis direction is data in the direction of the seismic interconnection survey line, and the z-axis direction is data in the direction of the seismic depth.

[0077] In this example, if Figure 5 As shown, the acquired 3D seismic data is used to test the effectiveness of the method proposed in this example.

[0078] Step S301: performing seismic horizon interpretation on the three-dimensional seismic data to obtain the horizon information of the nth horizon and the horizon information of the n-1th horizon;

[0079] In a specific implementation, the seismic horizon interpretation of three-dimensional seismic data is a technical means commonly used by those skilled in the art, and will not be described in detail in this application.

[0080] like Figure 6 As shown, it is the use of Figure 5 The layer information extracted from the imaging results of the three-dimensional seismic data is obtained through seismic layer interpretation. Figure 6 5 sets of layer information are shown.

[0081] Step S302: Take the inverse tangent function value of the difference between the layer information of two adjacent points in the same layer to obtain the slope of the nth layer and the slope of the n-1th layer of the corresponding point. The calculation formula is as follows:

[0082] S(x,y,z layer(x,y|n) )=arctan[Layer(x,y|n)-Layer(x,y-1|n)]

[0083] Among them, S(x,y,z layer(x,y|n) ) is the slope of point (x, y) in the nth layer, Layer(x, y|n) is the layer information of point (x, y) in the nth layer, Layer(x, y-1|n) is the layer information of point (x, y-1) in the nth layer, and point (x, y) and point (x, y-1) are two adjacent points in the same layer (i.e., the nth layer).

[0084] Through the slope of the nth layer, it is easy to get the slope of the n-1th layer of the corresponding point. The calculation formula is as follows:

[0085] S(x,y,z layer(x',y'|n-1) )=arctan[Layer(x',y'|n-1)-Layer(x',y'-1|n-1)]

[0086] Among them, S(x,y,z layer(x',y'|n-1) ) is the slope of the point (x',y') at the n-1th layer, Layer(x',y'|n-1) is the layer information of the point (x',y') at the nth layer, Layer(x',y'-1|n) is the layer information of the point (x',y'-1) at the n-1th layer, and the point (x',y') and the point (x',y'-1) are two adjacent points in the same layer (i.e., the n-1th layer).

[0087] By traversing the entire layer data, the slope of the layer corresponding to all layer positions can be calculated. Then the slope of the layer corresponding to the non-layer position in the data is calculated. This example uses the dual-parameter inverse distance weighted method to obtain it.

[0088] Step S303: according to the layer information of the nth layer and the layer information of the n-1th layer, a dual-parameter inverse distance weighted method is used to calculate the up-sampling coefficient on the non-layer and the down-sampling coefficient on the non-layer;

[0089] The upsampling coefficient on the non-layer is calculated as follows:

[0090]

[0091] Among them, A up (x,y,z n-1,n ) is the upsampling coefficient on the non-layer, z n-1,n It is the data in the seismic depth direction between the n-1th layer and the nth layer. Layer(x,y|n) is the layer information of the point (x,y) in the nth layer. Layer(x,y|n-1) is the layer information of the point (x,y) in the n-1th layer.

[0092] The downsampling coefficient on the non-layer is calculated as follows:

[0093]

[0094] Among them, A down (x,y,z n-1,n ) is the downsampling coefficient on the non-layer, z n-1,n It is the data in the seismic depth direction between the n-1th layer and the nth layer. Layer(x,y|n) is the layer information of the point (x,y) in the nth layer. Layer(x,y|n-1) is the layer information of the point (x,y) in the n-1th layer.

[0095] Step S304: According to the upsampling coefficient, the downsampling coefficient, the layer information of the nth layer and the layer information of the n-1th layer, the slope between the n-1th layer and the nth layer of the corresponding point is obtained, and the calculation formula is as follows:

[0096] S(x,y,z n-1,n )=A up (x,y,z n-1,n )×S(x,y,z layer(x,y|n-1) )+A down (x,y,z n-1,n )×S(x,y,z layer(x,y|n) )

[0097] Among them, S(x,y,z n-1,n ) is the slope between the n-1th layer and the nth layer, A up (x,y,z n-1,n ) is the upsampling coefficient on the non-layer, A down (x,y,z n-1,n ) is the downsampling coefficient on the non-layer, S(x,y,z layer(x,y|n) ) is the slope of the point (x, y) at the nth layer, S(x, y, z layer(x,y|n-1) ) is the slope of the point (x, y) at the n-1th layer.

[0098] Step S305: constructing the slope of the nth layer, the slope of the n-1th layer, and the slope between the n-1th layer and the nth layer of all points into the slope of the three-dimensional seismic data of the entire formation;

[0099] like Figure 7 As shown, it is the use of Figure 6 The slope of the 3D seismic data of the entire formation is calculated based on the 5 sets of stratigraphic information.

[0100] This embodiment proposes a slope extraction method based on geological layer guidance. First, three-dimensional seismic data is obtained; the slope of the nth layer, the slope of the n-1th layer, and the slope between the n-1th layer and the nth layer are extracted using the three-dimensional seismic data, and the slope of the nth layer, the slope of the n-1th layer, and the slope between the n-1th layer and the nth layer of all points are constructed as the slope of the three-dimensional seismic data of the entire formation. Among them, the slope of the nth layer and the slope of the n-1th layer are obtained by taking the inverse tangent function value of the difference between the layer information of two adjacent points of the same layer, and the two-parameter inverse distance weighted method is abandoned to calculate the upsampling coefficient on the non-layer and the downsampling coefficient on the non-layer. According to the upsampling coefficient, the downsampling coefficient, the layer information of the nth layer and the layer information of the n-1th layer, the slope between the n-1th layer and the nth layer is obtained. The slope extraction method of this embodiment has high accuracy and provides a reliable data basis for subsequent seismic data processing or seismic data interpretation.

[0101] Embodiment 3

[0102] This embodiment proposes a slope extraction device based on geological layer guidance, such as Figure 8 As shown, it includes: a data acquisition module, a slope extraction module and a slope construction module, wherein the data acquisition module is connected to the slope extraction module; the slope extraction module is connected to the slope construction module;

[0103] A data acquisition module, used for acquiring three-dimensional seismic data;

[0104] A slope extraction module, used to extract the slope of the nth layer, the slope of the n-1th layer, and the slope between the n-1th layer and the nth layer using the three-dimensional seismic data;

[0105] The slope construction module is used to construct the slope of the nth layer, the slope of the n-1th layer, and the slope between the n-1th layer and the nth layer of all points into the slope of the three-dimensional seismic data of the entire formation.

[0106] The slope extraction module includes: a first slope extraction unit and a second slope extraction unit; the first slope extraction unit is connected to the second slope extraction unit;

[0107] The first slope extraction unit is used to calculate the slope of the nth layer and the slope of the n-1th layer of the corresponding point, including:

[0108] Performing seismic horizon interpretation on the three-dimensional seismic data to obtain the horizon information of the nth horizon and the horizon information of the n-1th horizon;

[0109] Take the inverse tangent function value of the difference between the layer information of two adjacent points in the same layer to obtain the slope of the nth layer and the slope of the n-1th layer of the corresponding point. The calculation formula is as follows:

[0110] S(x,y,z layer(x,y|n) )=arctan[Layer(x,y|n)-Layer(x,y-1|n)]

[0111] Among them, S(x,y,z layer(x,y|n) ) is the slope of point (x, y) in the nth layer, Layer(x, y|n) is the layer information of point (x, y) in the nth layer, Layer(x, y-1|n) is the layer information of point (x, y-1) in the nth layer, and point (x, y) and point (x, y-1) are two adjacent points in the same layer (i.e., the nth layer).

[0112] Through the slope of the nth layer, it is easy to get the slope of the n-1th layer of the corresponding point. The calculation formula is as follows:

[0113] S(x,y,z layer(x',y'|n-1) )=arctan[Layer(x',y'|n-1)-Layer(x',y'-1|n-1)]

[0114] Among them, S(x,y,z layer(x',y'|n-1) ) is the slope of the point (x',y') at the n-1th layer, Layer(x',y'|n-1) is the layer information of the point (x',y') at the nth layer, Layer(x',y'-1|n) is the layer information of the point (x',y'-1) at the n-1th layer, and the point (x',y') and the point (x',y'-1) are two adjacent points in the same layer (i.e., the n-1th layer).

[0115] The second slope extraction unit is used to calculate the slope between the n-1th layer and the nth layer of the corresponding point, including:

[0116] According to the layer information of the nth layer and the layer information of the n-1th layer, a double-parameter inverse distance weighted method is used to calculate the up-sampling coefficient on the non-layer and the down-sampling coefficient on the non-layer;

[0117] The upsampling coefficient on the non-layer is calculated as follows:

[0118]

[0119] Among them, A up (x,y,z n-1,n ) is the upsampling coefficient on the non-layer, z n-1,nIt is the data in the seismic depth direction between the n-1th layer and the nth layer. Layer(x,y|n) is the layer information of the point (x,y) in the nth layer. Layer(x,y|n-1) is the layer information of the point (x,y) in the n-1th layer.

[0120] The downsampling coefficient on the non-layer is calculated as follows:

[0121]

[0122] Among them, A down (x,y,z n-1,n ) is the downsampling coefficient on the non-layer, z n-1,n It is the data in the seismic depth direction between the n-1th layer and the nth layer. Layer(x,y|n) is the layer information of the point (x,y) in the nth layer. Layer(x,y|n-1) is the layer information of the point (x,y) in the n-1th layer.

[0123] According to the upsampling coefficient, the downsampling coefficient, the layer information of the nth layer and the layer information of the n-1th layer, the slope between the n-1th layer and the nth layer of the corresponding point is obtained, and the calculation formula is as follows:

[0124] S(x,y,z n-1,n )=A up (x,y,z n-1,n )×S(x,y,z layer(x,y|n-1) )+A down (x,y,z n-1,n )×S(x,y,z layer(x,y|n) )

[0125] Among them, S(x,y,z n-1,n ) is the slope between the n-1th layer and the nth layer, A up (x,y,z n-1,n ) is the upsampling coefficient on the non-layer, A down (x,y,z n-1,n ) is the downsampling coefficient on the non-layer, S(x,y,z layer(x,y|n) ) is the slope of the point (x, y) at the nth layer, S(x, y, z layer(x,y|n-1) ) is the slope of the point (x, y) at the n-1th layer.

[0126] The slope construction unit is used to construct the slope of the nth layer, the slope of the n-1th layer, and the slope between the n-1th layer and the nth layer of all points into the slope of the three-dimensional seismic data of the entire formation.

[0127] The present embodiment proposes a slope extraction device based on geological layer guidance. First, a data acquisition module is used to acquire three-dimensional seismic data; secondly, the three-dimensional seismic data is used to extract the slope of the nth layer, the slope of the n-1th layer, and the slope between the n-1th layer and the nth layer; finally, a slope construction unit is used to construct the slope of the nth layer, the slope of the n-1th layer, and the slope between the n-1th layer and the nth layer of all points into the slope of the three-dimensional seismic data of the entire stratum. The accuracy of the slope extracted by the device of the present embodiment is improved, providing a reliable data basis for subsequent seismic data processing or seismic data interpretation.

[0128] Embodiment 4:

[0129] This embodiment provides an electronic device, including: a memory and a processor, wherein a computer program is stored in the memory, and when the computer program is executed by the processor, the slope extraction method based on geological layer guidance is executed.

[0130] The electronic device may be a mobile phone, a computer or a tablet computer, etc., including a memory and a processor, wherein a computer program is stored in the memory, and when the computer program is executed by the processor, the slope extraction method based on geological layer guidance as described in the embodiment is implemented. It is understood that the electronic device may also include an input / output (I / O) interface and a communication component.

[0131] The processor is used to execute all or part of the steps in the slope extraction method based on geological layer guidance as described in the above embodiment. The memory is used to store various types of data, which may include instructions of any application or method in the electronic device, as well as data related to the application.

[0132] The processor can be an application specific integrated circuit (ASIC), a digital signal processor (DSP), a programmable logic device (PLD), a field programmable gate array (FPGA), a controller, a microcontroller, a microprocessor or other electronic components, and is used to execute the slope extraction method based on geological layer guidance described in the above embodiment.

[0133] Embodiment 5

[0134] This embodiment provides a computer-readable storage medium storing executable instructions, which, when executed, enable a processor to execute the slope extraction method based on geological layer guidance.

[0135] If implemented in the form of a software functional unit and sold or used as an independent product, it can be stored in a computer-readable storage medium.

[0136] Based on this understanding, the technical solution of the present application, or the part that contributes to the prior art, or the part of the technical solution, can be embodied in the form of a software product. The computer software product is stored in a storage medium and includes a number of instructions for enabling a computer device (which can be a personal computer, a server, or a network device, etc.) to execute all or part of the steps of the slope extraction method based on geological layer guidance described in each embodiment of the present application.

[0137] The aforementioned storage media include: flash memory, hard disk, multimedia card, card-type memory (for example, SD (Secure Digital Memory Card) or DX (Memory Data Register, MDR abbreviation, memory data register) memory, etc.), random access memory (RAM), static random access memory (SRAM), read-only memory (ROM), electrically erasable programmable read-only memory (EEPROM), programmable read-only memory (PROM), magnetic memory, disk, CD, server, APP (Application, abbreviation of application software) application store and other media that can store program verification codes, on which computer programs are stored. When the computer program is executed by the processor, the various steps of the above-mentioned slope extraction method based on geological layer guidance can be implemented.

[0138] The various embodiments in the present disclosure are described in a progressive manner, and the same or similar parts between the various embodiments can be referenced to each other, and each embodiment focuses on the differences from other embodiments.

[0139] The protection scope of the present disclosure is not limited to the above-mentioned embodiments. Obviously, those skilled in the art can make various changes and modifications to the present disclosure without departing from the scope and spirit of the present disclosure. If these changes and modifications fall within the scope of the claims of the present disclosure and their equivalents, the intention of the present disclosure also includes these changes and modifications.

Claims

1. A slope extraction method based on geological layer guidance, characterized in that: include: Acquire 3D seismic data; The three-dimensional seismic data is used to extract the slope of the nth layer, the slope of the n-1th layer, and the slope between the n-1th layer and the nth layer of the corresponding point; The slope of the nth layer, the slope of the n-1th layer, and the slope between the n-1th layer and the nth layer of all points are constructed as the slope of the three-dimensional seismic data of the entire formation.

2. The slope extraction method based on geological layer guidance according to claim 1 is characterized in that: The step of extracting the slope of the nth layer and the slope of the n-1th layer of the corresponding point includes: Performing seismic horizon interpretation on the three-dimensional seismic data to obtain the horizon information of the nth horizon and the horizon information of the n-1th horizon; The inverse tangent function value is taken for the difference between the layer information of two adjacent points in the same layer, and the slope of the nth layer and the slope of the n-1th layer of the corresponding point are obtained.

3. The slope extraction method based on geological layer guidance according to claim 1 or 2, characterized in that: The x-axis direction of the three-dimensional seismic data is data in the direction of the seismic main survey line, the y-axis direction is data in the direction of the seismic contact survey line, and the z-axis direction is data in the direction of the seismic depth.

4. The slope extraction method based on geological layer guidance according to claim 2 is characterized in that: The slope between the n-1th layer and the nth layer is extracted by: According to the layer information of the nth layer and the layer information of the n-1th layer, a double-parameter inverse distance weighted method is used to calculate the up-sampling coefficient on the non-layer and the down-sampling coefficient on the non-layer; According to the up-sampling coefficient, the down-sampling coefficient, the layer information of the n-th layer and the layer information of the n-1-th layer, the slope between the n-1-th layer and the n-th layer of the corresponding point is obtained.

5. The slope extraction method based on geological layer guidance according to claim 2 is characterized in that: The inverse tangent function value is taken for the difference between the layer information of two adjacent points in the same layer to obtain the slope of the nth layer and the slope of the n-1th layer of the corresponding point. The calculation formula is as follows: S(x,y,z layer(x,y|n) )=arctan[Layer(x,y|n)-Layer(x,y-1|n)] Among them, S(x,y,z layer(x,y|n) ) is the slope of point (x, y) in the nth layer, Layer(x, y|n) is the layer information of point (x, y) in the nth layer, Layer(x, y-1|n) is the layer information of point (x, y-1) in the nth layer, and point (x, y) and point (x, y-1) are two adjacent points in the same layer.

6. The slope extraction method based on geological layer guidance according to claim 4 is characterized in that: The upsampling coefficient on the non-layer is calculated as follows: Among them, A up (x,y,z n-1,n ) is the upsampling coefficient on the non-layer, z n-1,n It is the data in the seismic depth direction between the n-1th layer and the nth layer. Layer(x,y|n) is the layer information of the point (x,y) in the nth layer. Layer(x,y|n-1) is the layer information of the point (x,y) in the n-1th layer.

7. The slope extraction method based on geological layer guidance according to claim 4 is characterized in that: The downsampling coefficient on the non-layer is calculated as follows: Among them, A down (x,y,z n-1,n ) is the downsampling coefficient on the non-layer, z n-1,n It is the data in the seismic depth direction between the n-1th layer and the nth layer. Layer(x,y|n) is the layer information of the point (x,y) in the nth layer. Layer(x,y|n-1) is the layer information of the point (x,y) in the n-1th layer.

8. The slope extraction method based on geological layer guidance according to claim 4 is characterized in that: According to the upsampling coefficient, the downsampling coefficient, the layer information of the nth layer and the layer information of the n-1th layer, the slope between the n-1th layer and the nth layer of the corresponding point is obtained, and the calculation formula is as follows: S(x,y,z n-1,n )=A up (x,y,z n-1,n )×S(x,y,z layer(x,y|n-1) )+A down (x,y,z n-1,n )×S(x,y,z layer(x,y|n) ) Among them, S(x,y,z n-1,n ) is the slope between the n-1th layer and the nth layer, A up (x,y,z n-1,n ) is the upsampling coefficient on the non-layer, A down (x,y,z n-1,n ) is the downsampling coefficient on the non-layer, S(x,y,z layer(x,y|n) ) is the slope of the point (x, y) at the nth layer, S(x, y, z layer(x,y|n-1) ) is the slope of the point (x, y) at the n-1th layer.

9. A slope extraction device based on geological layer guidance, characterized in that: include: A data acquisition module, used for acquiring three-dimensional seismic data; A slope extraction module, used to extract the slope of the nth layer, the slope of the n-1th layer, and the slope between the n-1th layer and the nth layer using the three-dimensional seismic data; The slope construction module is used to construct the slope of the nth layer, the slope of the n-1th layer, and the slope between the n-1th layer and the nth layer of all points into the slope of the three-dimensional seismic data of the entire formation.

10. An electronic device, characterized in that: include: A memory and a processor, wherein a computer program is stored in the memory, and when the computer program is executed by the processor, the slope extraction method based on geological layer guidance according to any one of claims 1 to 8 is executed.

11. A computer-readable storage medium, characterized in that: It stores executable instructions, which, when executed, enable the processor to execute the slope extraction method based on geological layer guidance as described in any one of claims 1-8.

Images