A slope extraction method and device based on geological horizon guidance

Abstract

The application provides a slope extraction method and device based on geological horizon guidance, and belongs to the technical field of geophysical exploration. The method comprises the following steps: obtaining three-dimensional seismic data; and extracting the slope of an n-th horizon, the slope of an n-1-th horizon and the slope between the n-1-th horizon and the n-th horizon by using the three-dimensional seismic data. The device comprises a data acquisition module and a slope extraction module. The method provided by the application improves the accuracy of slope extraction and provides reliable data basis for subsequent seismic data processing or seismic data interpretation.

Description

TECHNICAL FIELD

[0001] The application belongs to the technical field of geophysical exploration, and particularly relates to a slope extraction method and device based on geological horizon guidance. BACKGROUND

[0002] In the process of seismic data processing, slope information is crucial, and the conventional method can only obtain the slope of a certain layer, and generally 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 horizon constraint, resulting in low extraction accuracy and failing to meet the subsequent accuracy requirements. Therefore, it is an urgent technical problem to be solved to improve the slope extraction accuracy and provide reliable data basis for subsequent seismic data processing or seismic data interpretation. SUMMARY

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

[0004] In a first aspect, the application provides a slope extraction method based on geological horizon guidance, comprising:

[0005] acquiring three-dimensional seismic data;

[0006] 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 of the corresponding point using the three-dimensional seismic data;

[0007] constructing 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 of all points as the slope of the three-dimensional seismic data of the entire stratum.

[0008] The extraction of the slope of the nth horizon and the slope of the (n-1)th horizon of the corresponding point comprises:

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

[0010] taking the arctangent function value of the difference between the horizon information of the adjacent two points of the same horizon to obtain the slope of the nth horizon and the slope of the (n-1)th horizon of the corresponding point.

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

[0012] The slope between the (n-1)th horizon and the nth horizon comprises:

[0013] According to the layer information of the n-th layer and the layer information of the n-1-th layer, the up-sampling coefficient on the non-layer and the down-sampling coefficient on the non-layer are calculated by using a double-parameter inverse distance weighting method.

[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 difference of the layer information of the two adjacent points on the same layer is taken as an arctangent function value to obtain the slope of the n-th layer and the slope of the n-1-th layer of the corresponding point, and 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] Wherein, S(x,y,z layer(x,y|n) ) is the slope of the n-th layer of the point (x, y), Layer(x,y|n) is the layer information of the n-th layer of the point (x, y), Layer(x,y-1|n) is the layer information of the n-th layer of the point (x, y-1), and the point (x, y) and the point (x, y-1) are two adjacent points on the same layer.

[0018] The up-sampling coefficient on the non-layer is calculated by the following formula:

[0019]

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

[0021] The down-sampling coefficient on the non-layer is calculated by the following formula:

[0022]

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

[0024] The slope between the n-1th layer and the nth layer of the corresponding point is obtained according to the up-sampling coefficient, the down-sampling coefficient, the layer information of the n-1th layer and the layer information of the nth layer, 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] Wherein, 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 up-sampling coefficient on the non-layer, A down (x,y,z n-1,n ) is the down-sampling coefficient on the non-layer, S(x,y,z layer(x,y|n) ) is the slope of the point (x, y) on the nth layer, and S(x,y,z layer(x,y|n-1) ) is the slope of the point (x, y) on the n-1th layer.

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

[0028] A data acquisition module is configured to acquire three-dimensional seismic data.

[0029] A slope extraction module is configured to extract the slope of the n-1th 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] A slope construction module is configured to construct the slope of the n-1th layer, the slope of the n-1th layer and the slope between the n-1th layer and the nth layer of all points as the slope of the three-dimensional seismic data of the entire stratum.

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

[0032] In a fourth aspect, the application provides a computer readable storage medium, which stores executable instructions, and the instructions, when executed, cause the processor to execute the slope extraction method based on the geological layer guidance.

[0033] Advantages:

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

[0035] Figure 1 A slope extraction method based on geological horizon guidance is provided for the embodiment of the application.

[0036] Figure 2 A slope of the nth horizon and a slope of the (n-1)th horizon are provided for the embodiment of the application.

[0037] Figure 3 A slope between the (n-1)th horizon and the nth horizon is provided for the embodiment of the application.

[0038] Figure 4 A specific example of the slope extraction method based on geological horizon guidance is provided for the application.

[0039] Figure 5 A three-dimensional seismic data schematic diagram is provided for the embodiment of the application.

[0040] Figure 6 A horizon information schematic diagram is provided for the embodiment of the application.

[0041] Figure 7 A slope schematic diagram of three-dimensional seismic data of the entire stratum is provided for the embodiment of the application.

[0042] Figure 8 A principle block diagram of a slope extraction device based on geological horizon guidance is provided for the embodiment of the application. DETAILED DESCRIPTION

[0043] The present disclosure will be further described below in combination with the embodiments shown in the accompanying drawings.

[0044] In the process of seismic data processing, slope information is crucial. The conventional method 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 volume is obtained, the conventional method lacks horizon constraint, resulting in low extraction accuracy and failing to meet the subsequent accuracy requirements.

[0045] The application provides a slope extraction method and device based on geological horizon guidance. First, horizon interpretation is performed on three-dimensional seismic data to obtain interpreted horizon information. On this basis, the slope of data corresponding to horizon positions is calculated according to the horizon information. For non-horizon positions, the slopes calculated by using the upper and lower horizons are inversely distance weighted to obtain the slope of each point of the non-horizon position, and then the horizon slope of the entire data body is obtained, thereby improving the precision of slope extraction and providing reliable data basis for subsequent seismic data processing or seismic data interpretation.

[0046] Embodiment one,

[0047] The embodiment provides a slope extraction method based on geological horizon guidance, as shown in Figure 1 , which comprises the following steps.

[0048] Step S1: obtaining three-dimensional seismic data.

[0049] In the embodiment, the three-dimensional seismic data is obtained in a three-dimensional coordinate formed by taking an x axis as a horizontal axis, a y axis as a vertical axis and a z axis as a vertical axis, wherein the x axis direction of the three-dimensional seismic data is data in a main seismic line direction, the y axis direction is data in a seismic contact line direction, and the z axis direction is data in a seismic depth direction.

[0050] Step S2: extracting the slope of an n-th horizon, the slope of an (n-1)-th horizon and the slope between the (n-1)-th horizon and the n-th horizon by using the three-dimensional seismic data.

[0051] Step S3: constructing the slope of the n-th horizon, the slope of the (n-1)-th horizon and the slope between the (n-1)-th horizon and the n-th horizon of all points as the slope of the three-dimensional seismic data of the entire stratum.

[0052] In a specific implementation, the slope of the n-th horizon and the slope of the (n-1)-th horizon are extracted as shown in Figure 2 , and the extraction process is as follows.

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

[0054] Step S212: taking the arctangent function value of the difference between the horizon information of two adjacent points of the same horizon to obtain the slope of the n-th horizon and the slope of the (n-1)-th horizon of the corresponding points, and 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] wherein S(x,y,z layer(x,y|n)) is the slope of the n-th horizon for the point (x, y), Layer(x, y|n) is the horizon information of the n-th horizon for the point (x, y), Layer(x, y-1|n) is the horizon information of the n-th horizon for the point (x, y-1), and the points (x, y) and (x, y-1) are adjacent points in the same horizon (i.e., the n-th horizon).

[0057] The slope of the n-1-th horizon for the corresponding point can be easily obtained through the slope of the n-th horizon, and 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] wherein S(x, y, z layer(x',y'|n-1) ) is the slope of the n-1-th horizon for the point (x', y'), Layer(x', y'|n-1) is the horizon information of the n-th horizon for the point (x', y'), Layer(x', y'-1|n-1) is the horizon information of the n-1-th horizon for the point (x', y'-1), and the points (x', y') and (x', y'-1) are adjacent points in the same horizon (i.e., the n-1-th horizon).

[0060] In a specific implementation, the slope between the n-1-th horizon and the n-th horizon is as shown in FIG. 4, and the extraction process is as follows: Figure 3

[0061] Step S221: According to the horizon information of the n-th horizon and the horizon information of the n-1-th horizon, the up-sampling coefficient on the non-horizon and the down-sampling coefficient on the non-horizon are calculated by using a two-parameter inverse distance weighting method.

[0062] The up-sampling coefficient on the non-horizon is calculated as follows:

[0063]

[0064] wherein A up (x, y, z n-1,n ) is the up-sampling coefficient on the non-horizon, z n-1,n is the data in the seismic depth direction between the n-1-th horizon and the n-th horizon, Layer(x, y|n) is the horizon information of the n-th horizon for the point (x, y), and Layer(x, y|n-1) is the horizon information of the n-1-th horizon for the point (x, y).

[0065] The down-sampling coefficient on the non-horizon is calculated as follows:

[0066]

[0067] wherein A​down (x,y,z n-1,n ) is a down-sampling coefficient on non-layer position, z n-1,n is the data in the seismic depth direction between the n-1th layer position and the nth layer position, Layer(x,y|n) is the layer position information of point (x, y) in the nth layer position, and Layer(x,y|n-1) is the layer position information of point (x, y) in the n-1th layer position.

[0068] Step S222: According to the up-sampling coefficient, the down-sampling coefficient, the layer position information of the nth layer position, and the layer position information of the n-1th layer position, the slope between the n-1th layer position and the nth layer position 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] Wherein, S(x,y,z n-1,n ) is the slope between the n-1th layer position and the nth layer position, A up (x,y,z n-1,n ) is the up-sampling coefficient on non-layer position, A down (x,y,z n-1,n ) is the down-sampling coefficient on non-layer position, S(x,y,z layer(x,y|n) ) is the slope of point (x, y) in the nth layer position, and S(x,y,z layer(x,y|n-1) ) is the slope of point (x, y) in the n-1th layer position.

[0071] In a specific implementation, the slope S(x,y,z layer(x,y|n) ) of the nth layer position, the slope S(x,y,z layer(x',y'|n-1) ) of the n-1th layer position, and the slope S(x,y,z n-1,n ) between the n-1th layer position and the nth layer position of all points are constructed as the slope of the three-dimensional seismic data of the entire stratum.

[0072] The slope extraction method based on geological layer position guidance provided in the embodiment first acquires three-dimensional seismic data; extracts the slope of the nth layer position, the slope of the n-1th layer position, and the slope between the n-1th layer position and the nth layer position by using the three-dimensional seismic data, constructs the slope of the nth layer position, the slope of the n-1th layer position, and the slope between the n-1th layer position and the nth layer position of all points as the slope of the three-dimensional seismic data of the entire stratum, and the slope extracted by the method of the embodiment is high in accuracy.

[0073] Embodiment Two,

[0074] This embodiment is a specific example of Embodiment One, which uses actual seismic data to describe in detail the slope extraction method based on geological horizon guidance, and proves the effectiveness of the method proposed in this embodiment, as shown in Figure 4 , including:

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

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

[0077] In this example, as shown in Figure 5 , the obtained three-dimensional seismic data is used to test the effectiveness of the method proposed in this example.

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

[0079] In specific implementation, performing seismic horizon interpretation on three-dimensional seismic data is a common technical means for those skilled in the art, and will not be described herein.

[0080] As shown in Figure 6 , the horizon information is extracted using the imaging result of the three-dimensional seismic data Figure 5 , and through seismic horizon interpretation, five sets of horizon information as shown in Figure 6 are obtained.

[0081] Step S302: Take the arctangent value of the difference between the horizon information of the adjacent two points of the same horizon to obtain the slope of the nth horizon and the slope of the n-1th horizon of the corresponding point, and 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] Where S(x,y,z layer(x,y|n) ) is the slope of the nth horizon of point (x, y), Layer(x,y|n) is the horizon information of point (x, y) in the nth horizon, Layer(x,y-1|n) is the horizon information of point (x, y-1) in the nth horizon, and point (x, y) and point (x, y-1) are adjacent two points in the same horizon (i.e. the nth horizon).

[0084] The slope of the n-1th horizon of the corresponding point is easily obtained through the slope of the nth horizon, and 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] Wherein, S(x,y,z layer(x',y'|n-1) ) is the slope of the point (x',y') in the n-1 layer, Layer(x',y'|n-1) is the layer information of the point (x',y') in the n layer, Layer(x',y'-1|n-1) is the layer information of the point (x',y'-1) in the n-1 layer, the point (x',y') and the point (x',y'-1) are adjacent points in the same layer (i.e. the n-1 layer).

[0087] The slope of the layer corresponding to all layer positions can be calculated by traversing the entire layer data. Then the slope of the layer corresponding to the non-layer position in the data is calculated by using the double parameter inverse distance weighting method.

[0088] Step S303: according to the layer information of the n layer and the layer information of the n-1 layer, the up-sampling coefficient on the non-layer and the down-sampling coefficient on the non-layer are calculated by using the double parameter inverse distance weighting method;

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

[0090]

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

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

[0093]

[0094] Wherein, A down (x,y,z n-1,n ) is the down-sampling coefficient on the non-layer, z n-1,n is the data in the seismic depth direction between the n-1 layer and the n layer, Layer(x,y|n) is the layer information of the point (x,y) in the n layer, Layer(x,y|n-1) is the layer information of the point (x,y) in the n-1 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 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 3D 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 using the 5 sets of horizon information.

[0100] The embodiment provides a slope extraction method based on geological horizon guidance, first, three-dimensional seismic data is acquired; the slope of the nth horizon, the slope of the n-1th horizon and the slope between the n-1th horizon and the nth horizon are extracted from the three-dimensional seismic data, and the slope of the nth horizon, the slope of the n-1th horizon and the slope between the n-1th horizon and the nth horizon of all points are constructed as the slope of the three-dimensional seismic data of the entire stratum. Wherein, the slope of the nth horizon and the slope of the n-1th horizon are obtained by taking the arctangent function value of the difference of horizon information of two adjacent points of the same horizon, and the up-sampling coefficient on the non-horizon and the down-sampling coefficient on the non-horizon are obtained by abandoning the double-parameter inverse distance weighting method, the slope between the n-1th horizon and the nth horizon is obtained according to the up-sampling coefficient, the down-sampling coefficient, the horizon information of the nth horizon and the horizon information of the n-1th horizon, and the slope obtained by the slope extraction method of the embodiment has high precision, and provides reliable data basis for subsequent seismic data processing or seismic data interpretation.

[0101] Embodiment three,

[0102] The embodiment provides a slope extraction device based on geological horizon guidance, as shown in the figure, comprising: a data acquisition module, a slope extraction module and a slope construction module, the data acquisition module is connected with the slope extraction module; the slope extraction module is connected with the slope construction module; Figure 8

[0103] The data acquisition module is used for acquiring three-dimensional seismic data;

[0104] The slope extraction module is used for extracting the slope of the nth horizon, the slope of the n-1th horizon and the slope between the n-1th horizon and the nth horizon from the three-dimensional seismic data;

[0105] The slope construction module is used for constructing the slope of the nth horizon, the slope of the n-1th horizon and the slope between the n-1th horizon and the nth horizon of all points as the slope of the three-dimensional seismic data of the entire stratum.

[0106] The slope extraction module comprises a first slope extraction unit and a second slope extraction unit; the first slope extraction unit is connected with the second slope extraction unit;

[0107] The first slope extraction unit is used for calculating the slope of the nth horizon and the slope of the n-1th horizon of the corresponding point, and comprises:

[0108] The three-dimensional seismic data is interpreted by seismic horizon, and the horizon information of the nth horizon and the horizon information of the n-1th horizon are obtained;

[0109] ​The difference value of the layer information of two adjacent points in the same layer is taken as an arctangent function value to obtain the slope of the n-th layer of the corresponding point and the slope of the n-1-th layer, and 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] Wherein, S(x,y,z layer(x,y|n) ) is the slope of the n-th layer of the point (x, y), Layer(x,y|n) is the layer information of the n-th layer of the point (x, y), Layer(x,y-1|n) is the layer information of the n-th layer of the point (x, y-1), and the point (x, y) and the point (x, y-1) are two adjacent points in the same layer (i.e. the n-th layer).

[0112] The slope of the n-1-th layer of the corresponding point is easily obtained through the slope of the n-th layer, and 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] Wherein, S(x,y,z layer(x',y'|n-1) ) is the slope of the n-1-th layer of the point (x', y'), Layer(x',y'|n-1) is the layer information of the n-th layer of the point (x', y'), Layer(x',y'-1|n-1) is the layer information of the n-1-th layer of the point (x', y'-1), and the point (x', y') and the point (x', y'-1) are two adjacent points in the same layer (i.e. the n-1-th layer).

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

[0116] According to the layer information of the n-th layer and the layer information of the n-1-th layer, the up-sampling coefficient on the non-layer and the down-sampling coefficient on the non-layer are calculated by using a double-parameter inverse distance weighting method;

[0117] The up-sampling coefficient on the non-layer is calculated by the following formula:

[0118]

[0119] Wherein, A up (x,y,z n-1,n ) is the up-sampling coefficient on the non-layer, and z n-1,nThe data in the depth direction between the n-1th horizon and the nth horizon, Layer(x, y|n) is the horizon information of the point (x, y) in the nth horizon, and Layer(x, y|n-1) is the horizon information of the point (x, y) in the n-1th horizon.

[0120] The down-sampling coefficient on the non-horizon, and the calculation formula is as follows:

[0121]

[0122] Wherein, A down (x, y, z n-1,n ) is the up-sampling coefficient on the non-horizon, z n-1,n is the data in the depth direction between the n-1th horizon and the nth horizon, Layer(x, y|n) is the horizon information of the point (x, y) in the nth horizon, and Layer(x, y|n-1) is the horizon information of the point (x, y) in the n-1th horizon.

[0123] According to the up-sampling coefficient, the down-sampling coefficient, the horizon information of the n-1th horizon and the horizon information of the n-1th horizon, the slope between the n-1th horizon and the nth horizon 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] Wherein, S(x, y, z n-1,n ) is the slope between the n-1th horizon and the nth horizon, A up (x, y, z n-1,n ) is the up-sampling coefficient on the non-horizon, A down (x, y, z n-1,n ) is the down-sampling coefficient on the non-horizon, S(x, y, z layer(x,y|n) ) is the slope of the point (x, y) in the nth horizon, and S(x, y, z layer(x,y|n-1) ) is the slope of the point (x, y) in the n-1th horizon.

[0126] The slope construction unit is configured to construct the slope of the n-1th horizon, the slope of the nth horizon, and the slope between the n-1th horizon and the nth horizon of all points as the slope of the three-dimensional seismic data of the entire stratum.

[0127] This embodiment proposes a slope extraction device based on geological horizon guidance. First, a data acquisition module is used to acquire three-dimensional seismic data. Second, the three-dimensional seismic data is used to extract the slope of the nth horizon, the slope of the n-1th horizon, and the slope between the n-1th and nth horizons. Finally, a slope construction unit is used to construct the slope of the nth horizon, the slope of the n-1th horizon, and the slope between the n-1th and nth horizons at all points into the slope of the three-dimensional seismic data for the entire formation. The accuracy of the slope extracted by the device of this embodiment is improved, providing a reliable data foundation for subsequent seismic data processing or interpretation.

[0128] Example 4:

[0129] This embodiment provides an electronic device, including: 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 performed.

[0130] The electronic device can be a mobile phone, computer, or tablet computer, and includes a memory and a processor. The memory stores a computer program that, when executed by the processor, implements the slope extraction method based on geological horizon guidance as described in the embodiments. It is understood that the electronic device may also include an input / output (I / O) interface and a communication component.

[0131] The processor is configured to execute all or part of the steps of the slope extraction method based on geological horizon guidance as described in the above embodiment. The memory is configured to store various types of data, such as instructions for 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] Example 5

[0134] The embodiment provides a computer readable storage medium storing executable instructions, the instructions, when executed, causing a processor to perform the slope extraction method guided by a geological horizon.

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

[0136] Based on such understanding, the technical solution of the present application or the part of the technical solution which is essential or contributes to the prior art can be embodied in the form of a software product, which is stored in a storage medium and includes instructions for causing a computer device (which can be a personal computer, a server, or a network device, etc.) to execute all or part of the slope extraction method guided by a geological horizon described in various embodiments of the present application.

[0137] The aforementioned storage medium includes a flash memory, a hard disk, a multimedia card, a card-type memory (for example, an SD (Secure Digital Memory Card) or a DX (Memory Data Register, MDR) memory, etc.), a random access memory (RAM), a static random access memory (SRAM), a read-only memory (ROM), an electrically erasable programmable read-only memory (EEPROM), a programmable read-only memory (PROM), a magnetic memory, a magnetic disk, an optical disk, a server, an APP (Application) application store, and various media capable of storing program check codes, on which a computer program is stored, the computer program being executable by a processor to implement each step of the slope extraction method guided by a geological horizon.

[0138] Each of the embodiments in the present disclosure is described in a progressive manner, and the same or similar parts between the embodiments can be referred to each other, and each embodiment mainly explains the difference from other embodiments.

[0139] The protection scope of the present disclosure is not limited to the above-described embodiments, and it is obvious that those skilled in the art can make various modifications and changes to the present disclosure without departing from the scope and spirit of the present disclosure. If these modifications and changes belong to the scope of the claims of the present disclosure and equivalent technologies, the present disclosure also intends to include these modifications and changes.

Claims

1. A slope extraction method based on geological horizon 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; wherein: The extracting of the slope of the nth layer and the slope of the n-1th layer of the corresponding point comprises: performing seismic layer interpretation on the three-dimensional seismic data to obtain layer information of the nth layer and layer information of the n-1th layer; taking 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 slope between the n-1th layer and the nth layer is extracted by: using a dual-parameter inverse distance weighted method to calculate an upsampling coefficient on a non-layer and a downsampling coefficient on a non-layer according to the layer information of the nth layer and the layer information of the n-1th layer; obtaining the slope between the n-1th layer and the nth layer of the corresponding point 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 upsampling coefficient on the non-layer is calculated as follows: in, is the upsampling coefficient on the non-layer, is the data of seismic depth direction between the n-1th layer and the nth layer, is the layer information of point (x, y) at the nth layer, is the layer information of point (x, y) at the n-1th layer; The downsampling coefficient on the non-layer is calculated as follows: in, is the downsampling coefficient on the non-layer, is the data of seismic depth direction between the n-1th layer and the nth layer, is the layer information of point (x, y) at the nth layer, is the layer information of point (x, y) at the n-1th 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 of the corresponding point is obtained, and the calculation formula is as follows: in, is the slope between the n-1th layer and the nth layer, is the upsampling coefficient on the non-layer, is the downsampling coefficient on the non-layer, is the slope of point (x, y) at the nth layer, is the slope of the point (x, y) at the n-1th layer; 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, 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 interconnection survey line, and the z-axis direction is data in the direction of the seismic depth.

3. The slope extraction method based on geological layer guidance according to claim 1, 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: in, is the slope of point (x, y) at the nth layer, is the layer information of point (x, y) at the nth layer, is the layer information of point (x, y-1) in the nth layer. Point (x, y) and point (x, y-1) are two adjacent points in the same layer.

4. A slope extraction device based on geological layer guidance, characterized in that: include: A data acquisition module, used for acquiring three-dimensional seismic data; The slope extraction module 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 using the three-dimensional seismic data; wherein: The extracting of the slope of the nth layer and the slope of the n-1th layer comprises: performing seismic layer interpretation on the three-dimensional seismic data to obtain layer information of the nth layer and layer information of the n-1th layer; taking 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 points; The slope between the n-1th layer and the nth layer is extracted by: using a dual-parameter inverse distance weighted method to calculate an upsampling coefficient on a non-layer and a downsampling coefficient on a non-layer according to the layer information of the nth layer and the layer information of the n-1th layer; obtaining the slope between the n-1th layer and the nth layer of the corresponding point 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 upsampling coefficient on the non-layer is calculated as follows: in, is the upsampling coefficient on the non-layer, is the data of seismic depth direction between the n-1th layer and the nth layer, is the layer information of point (x, y) at the nth layer, is the layer information of point (x, y) at the n-1th layer; The downsampling coefficient on the non-layer is calculated as follows: in, is the downsampling coefficient on the non-layer, is the data of seismic depth direction between the n-1th layer and the nth layer, is the layer information of point (x, y) at the nth layer, is the layer information of point (x, y) at the n-1th 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 of the corresponding point is obtained, and the calculation formula is as follows: in, is the slope between the n-1th layer and the nth layer, is the upsampling coefficient on the non-layer, is the downsampling coefficient on the non-layer, is the slope of point (x, y) at the nth layer, is the slope of the point (x, y) at the n-1th layer; 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.

5. 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 3 is executed.

6. A computer-readable storage medium, characterized in that The device stores executable instructions, which, when executed, enable a processor to execute the slope extraction method based on geological layer guidance according to any one of claims 1 to 3.

Images