OVT surface element attribute extraction method and device and storage medium

Abstract

The invention provides an OVT surface element attribute extraction method and device and a storage medium, and belongs to the field of seismic exploration and seismic data processing. The method comprises the following steps: acquiring seismic data; carrying out OVT surface element division on the seismic data, and obtaining transverse and longitudinal groups of each OVT surface element; calculating coordinates of the center of each OVT surface element; calculating an offset distance and an azimuth angle of each OVT surface element; and displaying the OVT surface element attribute graph. The OVT piece attribute graph is drawn and output according to the coordinates, the offset distance and the azimuth angle of the surface element center obtained through calculation, the seismic data are converted into visual graphs from abstract numbers, analysis and interpretation are facilitated, a processing result can be displayed to related personnel conveniently, and data communication and sharing are promoted.

Description

Technical Field

[0001] This invention belongs to the field of seismic data processing in seismic exploration, and specifically relates to an OVT surface element attribute extraction method, device, and storage medium. Background Technology

[0002] As oil and gas exploration continues to deepen, the precision requirements for seismic data processing are increasing, leading to a trend towards broadband and full-field processing. Offset Vector Tile (OVT) domain processing has become the preferred method for processing wide-azimuth seismic data. OVT pixel partitioning has become a primary task in OVT domain data processing, and the need for quality control of OVT domain data is also growing.

[0003] Within the cross-shaped arrangement domain, the subsurface cover area is divided into small rectangles at equal intervals according to the shot-receiver distance and receiver line distance; each rectangle is an OVT patch. Chinese public document CN113589361A discloses a rapid method for dividing seismic data shot-receiver distance vector patches. This method calculates the longitudinal and lateral components of the shot-receiver distance and groups the OVTs according to the shot-receiver distance, receiver line distance, shot point distance, and receiver line distance to obtain the OVT patch division results.

[0004] After the face elements are divided, in order to display the overall division of the cross arrangement and attribute information such as offset distance and azimuth angle, it is necessary to calculate the coordinates of the face elements, as well as the offset distance and azimuth angle of each face element. However, there is no corresponding calculation method disclosed in the current technology.

[0005] Based on this, the present invention proposes an OVT (Optical Variable Transformation) element attribute extraction method, which is used to calculate the coordinates of the OVT elements displayed in the cross-shaped arrangement domain, as well as the offset distance and azimuth angle of each element. Summary of the Invention

[0006] The purpose of this invention is to solve the problems existing in the prior art and to provide an OVT surface element attribute extraction method, apparatus and storage medium.

[0007] This invention is achieved through the following technical solution:

[0008] A first aspect of the present invention provides a method for extracting the attributes of OVT surface elements, comprising:

[0009] S100, acquire earthquake data;

[0010] S200: Divide the seismic data into OVT (Outer Threshold) elements and obtain the horizontal and vertical groupings of each OVT element; S300: Calculate the coordinates of the center of each OVT element;

[0011] S400, calculate the offset and azimuth of each OVT element;

[0012] S500 displays the attribute map of the OVT surface element.

[0013] A further improvement of the present invention is that:

[0014] In S200, the horizontal and vertical groupings of each OVT element are obtained, where the horizontal grouping number is denoted as x_num and the vertical grouping number is denoted as y_num;

[0015] The horizontal division size of the face element is x_dis, and the vertical division size of the face element is y_dis.

[0016] A further improvement of the present invention is that:

[0017] In S300, the coordinates of the center of each OVT element are calculated. The specific operations include:

[0018] The center coordinates of each OVT element are calculated using the following formula:

[0019] x-coor: x_coor = (x_num - 1) × x_dis + 1 / 2 × x_dis

[0020] The ordinate is: y_coor = (y_num-1) × y_dis + 1 / 2 × y_dis.

[0021] A further improvement of the present invention is that:

[0022] In S400, the offset and azimuth of each OVT element are calculated. The specific operations include:

[0023] The offset and azimuth of each OVT element are calculated using the following formula:

[0024] Offset distance:

[0025] Azimuth: Azimuth = tan -1 (x_num*x_dis / y_num*y_dis).

[0026] A further improvement of the present invention is that:

[0027] Converting azimuth radians to degrees involves the following steps:

[0028] First, convert radians to degrees using the following formula:

[0029]

[0030] Then, determine which quadrant the angle belongs to and its specific angle value, including:

[0031] If Azimuth_degree > 0, then when x_num > 0 and y_num > 0, Azimuth_degree belongs to the first quadrant and its value remains unchanged; otherwise, it belongs to the third quadrant, and the azimuth angle value Azimuth_degree = zimuth_degree + 180.

[0032] If Azimuth_degree <= 0, and x_num >= 0 and y_num <= 0, then Azimuth_degree belongs to the second quadrant, and the azimuth angle value Azimuth_degree = Azimuth_degree + 180; otherwise, Azimuth_degree belongs to the fourth quadrant, and the azimuth angle value Azimuth_degree = Azimuth_degree + 360.

[0033] A further improvement of the present invention is that:

[0034] Displaying the OVT element attribute map in S500 includes the following specific operations:

[0035] Draw and output the OVT sheet attribute map based on the center coordinates of the OVT sheet element, as well as the offset distance and azimuth angle of the OVT sheet element;

[0036] If the observation system has a certain angle, the coordinates of the OVT element center are rotated, and the coordinates after rotation are:

[0037]

[0038] Where θ represents the angle of the observation system.

[0039] A second aspect of the present invention provides an OVT (Optical Variable Transformation) surface element attribute extraction device, comprising:

[0040] Acquisition unit, used to acquire seismic data;

[0041] The surface unit is used to divide the seismic data into surface units (OVT) and obtain the horizontal and vertical groupings of each OVT surface unit.

[0042] The first calculation unit is used to calculate the coordinates of the center of each OVT element;

[0043] The second calculation unit is used to calculate the offset and azimuth of each OVT element;

[0044] The display unit is used to display the attribute map of the OVT surface element.

[0045] A further improvement of the present invention is that:

[0046] The first calculation unit performs the following operations:

[0047] The center coordinates of each OVT element are calculated using the following formula:

[0048] x-coor: x_coor = (x_num - 1) × x_dis + 1 / 2 × x_dis

[0049] The ordinate is: y_coor = (y_num-1) × y_dis + 1 / 2 × y_dis.

[0050] A further improvement of the present invention is that:

[0051] The second calculation unit performs the following operations:

[0052] The offset and azimuth of each OVT element are calculated using the following formula:

[0053] Offset distance:

[0054] Azimuth: Azimuth = tan -1 (x_num*x_dis / y_num*y_dis);

[0055] Converting azimuth radians to degrees involves the following steps:

[0056] First, convert radians to degrees using the following formula:

[0057]

[0058] Then, determine which quadrant the angle belongs to and its specific angle value, including:

[0059] If Azimuth_degree > 0, then when x_num > 0 and y_num > 0, Azimuth_degree belongs to the first quadrant and its value remains unchanged; otherwise, it belongs to the third quadrant, and the azimuth angle value Azimuth_degree = zimuth_degree + 180.

[0060] If Azimuth_degree <= 0, and x_num >= 0 and y_num <= 0, then Azimuth_degree belongs to the second quadrant, and the azimuth angle value Azimuth_degree = Azimuth_degree + 180; otherwise, Azimuth_degree belongs to the fourth quadrant, and the azimuth angle value Azimuth_degree = Azimuth_degree + 360.

[0061] A third aspect of the present invention provides a computer-readable storage medium storing at least one computer-executable program, which, when executed by the computer, causes the computer to perform the steps in the OVT surface attribute extraction method.

[0062] Compared with the prior art, the beneficial effects of the present invention are:

[0063] This invention first calculates the coordinates of the center of each OVT surface element based on the horizontal and vertical grouping and the interval of the surface element division. Then, it calculates the azimuth and offset values ​​of each surface element. Finally, it draws the OVT surface element attribute map based on the center coordinates, offset and azimuth of the surface element center. This provides more comprehensive and accurate information for the subsequent processing and analysis of seismic data, which is unprecedented in the prior art.

[0064] This invention calculates the coordinates, offset, and azimuth of the surface element center, and then plots and outputs an OVT (Operational Vibration Tunneling) image attribute map based on these coordinates, offset, and azimuth. This helps processing personnel to more intuitively understand the spatial distribution characteristics of seismic data. The graphical display of offset and azimuth in this invention transforms seismic data from abstract numbers into intuitive graphics, facilitating analysis and interpretation, and making it easier to present processing results to relevant personnel, promoting data exchange and sharing. It can also evaluate the performance and layout of the observation system, providing optimization suggestions for subsequent data acquisition and processing.

[0065] This invention also provides a method for displaying and processing the attributes of an observation system when it has a certain angle, which is suitable for the needs of complex observation systems. Attached Figure Description

[0066] Figure 1 A flowchart of an OVT surface element attribute extraction method provided by the present invention;

[0067] Figure 2 This is the result of OVT surface element division in Embodiment 6 of the present invention;

[0068] Figure 3 This is a property diagram of the offset distance of the OVT surface element in Embodiment 6 of the present invention;

[0069] Figure 4 This is a property diagram of the azimuth angle of the OVT surface element in Embodiment 6 of the present invention;

[0070] Figure 5 This is a property diagram of the offset distance of the OVT surface element when the observation system has a certain angle in Embodiment 6 of the present invention;

[0071] Figure 6 This is a property diagram of the azimuth angle of the OVT surface element when the observation system has a certain angle in Embodiment 6 of the present invention. Detailed Implementation

[0072] The present invention will now be described in further detail with reference to the accompanying drawings:

[0073] Within the cross-shaped arrangement domain, the underground coverage area is divided into small rectangles at equal distances according to the shot-receiver distance and receiver line distance. Each rectangle is an OVT sheet, or OVT element. After the OVT elements are divided, in order to display the overall division of the cross-shaped arrangement and attribute information such as offset distance and azimuth angle, it is necessary to calculate the coordinates of the displayed elements, as well as the offset distance and azimuth angle of each element. However, there is currently no corresponding calculation method disclosed in the existing technology.

[0074] Based on this, the present invention proposes an OVT (Optical Variable Transformation) element attribute extraction method, which is used to calculate the coordinates of the OVT elements displayed in the cross-shaped arrangement domain, as well as the offset distance and azimuth angle of each element.

[0075]

Example 1

[0076] This invention proposes a method for extracting the attributes of OVT surface elements, such as... Figure 1 As shown, the specific steps include:

[0077] S100, acquire earthquake data;

[0078] S200 divides the seismic data into OVT (Outer Threshold) elements and obtains the horizontal and vertical groupings of each OVT element;

[0079] S300, calculate the coordinates of the center of each OVT element;

[0080] S400, calculate the offset and azimuth of each OVT element;

[0081] S500 displays the attribute map of the OVT surface element.

[0082] This invention first calculates the coordinates of the center of each OVT surface element based on the horizontal and vertical grouping and the interval of the surface element division. Then, it calculates the azimuth and offset values ​​of the center of each surface element. Finally, it draws the OVT surface element attribute map based on the center coordinates, offset and azimuth of the surface element center. This provides more comprehensive and accurate information for the subsequent processing and analysis of seismic data, which is unprecedented in the prior art.

[0083]

Example 2

[0084] In S200, seismic data is divided into OVT (Outer Threshold) elements, and the horizontal and vertical groups of each OVT element are obtained. Specific operations include:

[0085] Existing techniques are used to divide seismic data into OVT (Outer Virtual Unit) cells. For example, the method disclosed in Chinese Patent Application No. 202010361237.8 can be used to divide the data into cells and obtain the horizontal and vertical groupings of each OVT cell.

[0086] The horizontal group number is denoted as x_num, and the vertical group number is denoted as y_num;

[0087] The horizontal division size of the face element is x_dis, and the vertical division size of the face element is y_dis.

[0088]

Example 3

[0089] In S300, the coordinates of the center of each OVT element are calculated. The specific operations include:

[0090] The information of each OVT element is represented by the information of its center point. The center coordinates of each OVT element are calculated using the following formula:

[0091] x-coor: x_coor = (x_num - 1) × x_dis + 1 / 2 × x_dis

[0092] The ordinate is: y_coor = (y_num-1) × y_dis + 1 / 2 × y_dis.

[0093]

Example 4

[0094] In S400, the offset and azimuth of each OVT element are calculated. The specific operations include:

[0095] The offset and azimuth of each OVT element are calculated using the following formula:

[0096] Offset distance:

[0097] Azimuth: Azimuth = tan -1 (x_num*x_dis / y_num*y_dis).

[0098] The azimuth angle value obtained from the above azimuth angle calculation formula is limited to a range of 0 to 360 degrees to more clearly and intuitively display the azimuth angle attribute. Therefore, the azimuth angle is further processed as follows:

[0099] First, convert radians to degrees:

[0100] Then, determine which quadrant the angle belongs to and its specific angle value. Specific operations include:

[0101] If Azimuth_degree > 0, then when x_num > 0 and y_num > 0, Azimuth_degree belongs to the first quadrant and its value remains unchanged; otherwise (that is, when x_num > 0 and y_num ≯ 0, or x_num ≯ 0 and y_num > 0, or x_num ≯ 0 and y_num ≯ 0), it belongs to the third quadrant, and the azimuth angle value Azimuth_degree = zimuth_degree + 180;

[0102] If Azimuth_degree <= 0, and x_num >= 0, y_num <= 0, then Azimuth_degree belongs to the second quadrant, and the azimuth angle value Azimuth_degree = Azimuth_degree + 180; otherwise (that is, when x_num < 0, y_num <= 0 or x_num >= 0, y_num > 0 or x_num < 0, y_num > 0), Azimuth_degree belongs to the fourth quadrant, and the azimuth angle value Azimuth_degree = Azimuth_degree + 360.

[0103]

Example 5

[0104] Displaying the OVT element attribute map in S500 includes the following specific operations:

[0105] Draw and output the OVT sheet attribute map based on the center coordinates of the OVT sheet element, as well as the offset distance and azimuth angle of the OVT sheet element;

[0106] If the observation system has a certain angle, the coordinates of the OVT element center will be rotated, and the coordinates after rotation will be:

[0107]

[0108] Where θ represents the angle of the observation system.

[0109] The method described in this invention calculates the coordinates, offset, and azimuth of the surface element center, and then plots and outputs an OVT (Operational Vibration Tunneling) image attribute map based on these coordinates, offset, and azimuth. This helps processing personnel to more intuitively understand the spatial distribution characteristics of seismic data. The graphical display of offset and azimuth transforms seismic data from abstract numbers into intuitive graphics, facilitating analysis and interpretation, and enabling the presentation of processing results to relevant personnel. This promotes data exchange and sharing, allows for the evaluation of the observation system's performance and layout, and provides optimization suggestions for subsequent data acquisition and processing.

[0110] This invention also provides a method for displaying and processing the attributes of an observation system when it has a certain angle, which is suitable for the needs of complex observation systems.

[0111]

Example 6

[0112] In this embodiment of the invention, for seismic data in a cross-shaped arrangement domain, the OVT (Optical Variable Transform) grouping number of each data track is calculated according to Chinese patent application number 202010361237.8.

[0113] like Figure 2 As shown in Figure 2.2, (2.2) represents a face cell with horizontal group number x_num = 2 and vertical group number y_num = 2. The horizontal division size of the face cell is defined as x_dis = 500m and y_dis = 300.

[0114] The information of each surface element is represented by the information of its center point in the OVT surface element. The coordinates of the center of each OVT surface element are obtained. Taking surface element (2,2) as an example, the coordinates displayed by its attributes are calculated as follows:

[0115] x_coor=(2-1)×500+1 / 2×500=750,

[0116] y_coor=(2-1)×300+1 / 2×300=650.

[0117] Calculate the offset and azimuth of this element:

[0118] Offset distance:

[0119] Azimuth: Azimuth = tan -1 (5 / 3)≈0.9273.

[0120] The azimuth angle value obtained from the above azimuth angle calculation formula is limited to a range of 0 to 360 degrees to more clearly and intuitively display the azimuth angle attribute. Therefore, the azimuth angle is further processed as follows:

[0121] Convert radians to degrees: Azimuth_degree is greater than 0, and x_num>0, y_num>0, which belongs to the first quadrant. After converting to 0~360 degrees, Azimuth_degree is still 53.13 degrees.

[0122] After calculating all OVT elements according to the above steps, draw the OVT patch attribute map based on the element center coordinates, offset distance, and azimuth angle, as shown below. Figure 3 and Figure 4 As shown.

[0123] When the observation system has a certain angle, the calculated coordinates are rotated, and the result after rotation is as follows. Figure 5 and Figure 6 As shown.

[0124]

Example 7

[0125] This invention provides an OVT (Optical Variable Transformation) surface element attribute extraction device, comprising:

[0126] Acquisition unit, used to acquire seismic data;

[0127] The surface unit is used to divide the seismic data into surface units (OVT) and obtain the horizontal and vertical groupings of each OVT surface unit. Specifically, the following operations are performed:

[0128] Existing techniques are used to divide seismic data into OVT (Outer Virtual Unit) cells. For example, the method disclosed in Chinese Patent Application No. 202010361237.8 can be used to divide the data into cells and obtain the horizontal and vertical groupings of each OVT cell.

[0129] The horizontal group number is denoted as x_num, and the vertical group number is denoted as y_num;

[0130] The horizontal division size of the face element is x_dis, and the vertical division size of the face element is y_dis.

[0131] The first calculation unit is used to calculate the coordinates of the center of each OVT element, and specifically performs the following operations:

[0132] The information of each OVT element is represented by the information of its center point. The center coordinates of each OVT element are calculated using the following formula:

[0133] x-coor: x_coor = (x_num - 1) × x_dis + 1 / 2 × x_dis

[0134] The ordinate is: y_coor = (y_num-1) × y_dis + 1 / 2 × y_dis.

[0135] The second calculation unit is used to calculate the offset and azimuth of each OVT element, and specifically performs the following operations:

[0136] The offset and azimuth of each OVT element are calculated using the following formula:

[0137] Offset distance:

[0138] Azimuth: Azimuth = tan -1 (x_num*x_dis / y_num*y_dis).

[0139] The azimuth angle value obtained from the above azimuth angle calculation formula is limited to a range of 0 to 360 degrees to more clearly and intuitively display the azimuth angle attribute. Therefore, the azimuth angle is further processed as follows:

[0140] First, convert radians to degrees:

[0141] Then, determine which quadrant the angle belongs to and its specific angle value. Specific operations include:

[0142] If Azimuth_degree > 0, then when x_num > 0 and y_num > 0, Azimuth_degree belongs to the first quadrant and its value remains unchanged; otherwise, it belongs to the third quadrant, and the azimuth angle value Azimuth_degree = zimuth_degree + 180.

[0143] If Azimuth_degree <= 0, and x_num >= 0 and y_num <= 0, then Azimuth_degree belongs to the second quadrant, and the azimuth angle value Azimuth_degree = Azimuth_degree + 180; otherwise, Azimuth_degree belongs to the fourth quadrant, and the azimuth angle value Azimuth_degree = Azimuth_degree + 360.

[0144] The display unit is used to display the OVT element attribute map and performs the following operations:

[0145] Draw and output the OVT sheet attribute map based on the center coordinates of the OVT sheet element, as well as the offset distance and azimuth angle of the OVT sheet element;

[0146] If the observation system has a certain angle, the coordinates of the OVT element center will be rotated, and the coordinates after rotation will be:

[0147]

[0148] Where θ represents the angle of the observation system.

[0149]

Example 8

[0150] This invention provides a computer-readable storage medium storing at least one computer-executable program, which, when executed by the computer, causes the computer to perform the following steps:

[0151] S100, acquire earthquake data;

[0152] S200 divides the seismic data into OVT (Outer Threshold) elements and obtains the horizontal and vertical groupings of each OVT element;

[0153] in,

[0154] The horizontal group number is denoted as x_num, and the vertical group number is denoted as y_num;

[0155] The horizontal division size of the face element is x_dis, and the vertical division size of the face element is y_dis.

[0156] S300, calculates the coordinates of the center of each OVT element. Specific operations include:

[0157] The information of each OVT element is represented by the information of its center point. The center coordinates of each OVT element are calculated using the following formula:

[0158] x-coor: x_coor = (x_num - 1) × x_dis + 1 / 2 × x_dis

[0159] The ordinate is: y_coor = (y_num-1) × y_dis + 1 / 2 × y_dis.

[0160] S400 calculates the offset and azimuth of each OVT element. Specific operations include:

[0161] The offset and azimuth of each OVT element are calculated using the following formula:

[0162] Offset distance:

[0163] Azimuth: Azimuth = tan -1 (x_num*x_dis / y_num*y_dis).

[0164] The azimuth angle value obtained from the above azimuth angle calculation formula is limited to a range of 0 to 360 degrees to more clearly and intuitively display the azimuth angle attribute. Therefore, the azimuth angle is further processed as follows:

[0165] First, convert radians to degrees:

[0166] Then, determine which quadrant the angle belongs to and its specific angle value. Specific operations include:

[0167] If Azimuth_degree > 0, then when x_num > 0 and y_num > 0, Azimuth_degree belongs to the first quadrant and its value remains unchanged; otherwise, it belongs to the third quadrant, and the azimuth angle value Azimuth_degree = Azimuth_degree + 180.

[0168] If Azimuth_degree <= 0, and x_num >= 0 and y_num <= 0, then Azimuth_degree belongs to the second quadrant, and the azimuth angle value Azimuth_degree = Azimuth_degree + 180; otherwise, Azimuth_degree belongs to the fourth quadrant, and the azimuth angle value Azimuth_degree = Azimuth_degree + 360.

[0169] S500 displays the OVT element attribute map. Specific operations include:

[0170] Draw and output the OVT sheet attribute map based on the center coordinates of the OVT sheet element, as well as the offset distance and azimuth angle of the OVT sheet element;

[0171] If the observation system has a certain angle, the coordinates of the OVT element center will be rotated, and the coordinates after rotation will be:

[0172]

[0173] Where θ represents the angle of the observation system.

[0174] Those skilled in the art will understand that all or part of the processes in the methods of the above embodiments can be implemented by a computer program instructing related hardware. The computer program can be stored in a non-volatile computer-readable storage medium, and when executed, it can include the processes of the embodiments of the above methods. Any references to memory, storage, databases, or other media used in the embodiments provided in this application can include non-volatile and / or volatile memory. Non-volatile memory can include read-only memory (ROM), programmable ROM (PROM), electrically programmable ROM (EPROM), electrically erasable programmable ROM (EEPROM), or flash memory. Volatile memory can include random access memory (RAM) or external cache memory. By way of illustration and not limitation, RAM is available in various forms, such as static RAM (SRAM), dynamic RAM (DRAM), synchronous DRAM (SDRAM), dual data rate SDRAM (DDRSDRAM), enhanced SDRAM (ESDRAM), synchronous link DRAM (SLDRAM), Rambus direct RAM (RDRAM), direct memory bus dynamic RAM (DRDRAM), and memory bus dynamic RAM (RDRAM), etc.

[0175] The technical features of the above embodiments can be combined in any way. For the sake of brevity, not all possible combinations of the technical features in the above embodiments are described. However, as long as there is no contradiction in the combination of these technical features, they should be considered to be within the scope of this specification.

[0176] The above technical solution is only one embodiment of the present invention. For those skilled in the art, based on the principles disclosed in the present invention, it is easy to make various types of improvements or modifications, and not limited to the technical solutions described in the specific embodiments of the present invention. Therefore, the foregoing description is only a preferred option and is not restrictive.

Claims

1. A method for extracting attributes of OVT surface elements, characterized in that, include: S100, acquire earthquake data; S200 divides the seismic data into OVT (Outer Threshold) elements and obtains the horizontal and vertical groupings of each OVT element; S300, calculate the coordinates of the center of each OVT element; S400, calculate the offset and azimuth of each OVT element; S500 displays the attribute map of the OVT surface element.

2. The OVT surface element attribute extraction method according to claim 1, characterized in that, In step S200, the horizontal and vertical groups of each OVT element are obtained, where the horizontal group number is denoted as x_num and the vertical group number is denoted as y_num; The horizontal division size of the face element is x_dis, and the vertical division size of the face element is y_dis.

3. The OVT surface element attribute extraction method according to claim 2, characterized in that, Step S300 calculates the coordinates of the center of each OVT element. The specific operations include: The center coordinates of each OVT element are calculated using the following formula: x-coor: x_coor = (x_num - 1) × x_dis + 1 / 2 × x_dis The ordinate is: y_coor = (y_num-1) × y_dis + 1 / 2 × y_dis.

4. The OVT surface element attribute extraction method according to claim 2, characterized in that, Step S400 calculates the offset distance and azimuth angle of each OVT element. Specific operations include: The offset and azimuth of each OVT element are calculated using the following formula: Offset distance: New:Azimuth=tan -1 (x_num*x_dis / y_num*y_dis)。 5. The OVT surface element attribute extraction method according to claim 4, characterized in that, Converting azimuth radians to degrees involves the following steps: First, convert radians to degrees using the following formula: Then, determine which quadrant the angle belongs to and its specific angle value, specifically including: If Azimuth_degree > 0, when x_num > 0 and y_num > 0, then Azimuth_degree belongs to the first quadrant and its value remains unchanged; otherwise, Azimuth_degree belongs to the third quadrant, and the azimuth angle value Azimuth_degree = zimuth_degree + 180. If Azimuth_degree <= 0, and x_num >= 0 and y_num <= 0, then Azimuth_degree belongs to the second quadrant, and the azimuth angle value Azimuth_degree = Azimuth_degree + 180; otherwise, Azimuth_degree belongs to the fourth quadrant, and the azimuth angle value Azimuth_degree = Azimuth_degree + 360.

6. The OVT surface element attribute extraction method according to claim 3, characterized in that, Step S500 displays the OVT element attribute map, and the specific operations include: Draw and output the OVT sheet attribute map based on the center coordinates of the OVT sheet element, as well as the offset distance and azimuth angle of the OVT sheet element; If the observation system has a certain angle, the coordinates of the OVT element center are rotated, and the coordinates after rotation are: Where θ represents the angle of the observation system.

7. An OVT (Optical Variable Transformation) surface element attribute extraction device, characterized in that, include: Acquisition unit, used to acquire seismic data; The surface unit is used to divide the seismic data into surface units (OVT) and obtain the horizontal and vertical groupings of each OVT surface unit. The first calculation unit is used to calculate the coordinates of the center of each OVT element; The second calculation unit is used to calculate the offset and azimuth of each OVT element; The display unit is used to display the attribute map of the OVT surface element.

8. The OVT surface element attribute extraction device according to claim 7, characterized in that, The first calculation unit performs the following operations: The center coordinates of each OVT element are calculated using the following formula: x-coor: x_coor = (x_num - 1) × x_dis + 1 / 2 × x_dis The ordinate is: y_coor = (y_num-1) × y_dis + 1 / 2 × y_dis.

9. The OVT surface element attribute extraction device according to claim 7, characterized in that, The second calculation unit performs the following operations: The offset and azimuth of each OVT element are calculated using the following formula: Offset distance: Azimuth angle: Azimuth=tan-1(x_num*x_dis / y_num*y_dis); Converting azimuth radians to degrees involves the following steps: First, convert radians to degrees using the following formula: Then, determine which quadrant the angle belongs to and its specific angle value, including: If Azimuth_degree > 0, and x_num > 0 and y_num > 0, then Azimuth_degree belongs to the first quadrant and its value remains unchanged; otherwise, it belongs to the third quadrant, and the azimuth angle value Azimuth_degree = zimuth_degree + 180. If Azimuth_degree <= 0, and x_num >= 0 and y_num <= 0, then Azimuth_degree belongs to the second quadrant, and the azimuth angle value Azimuth_degree = Azimuth_degree + 180; otherwise, Azimuth_degree belongs to the fourth quadrant, and the azimuth angle value Azimuth_degree = Azimuth_degree + 360.

10. A computer-readable storage medium storing at least one computer-executable program, which, when executed by the computer, causes the computer to perform the steps of the OVT surface element attribute extraction method as claimed in any one of claims 1-6.

Images