An OVT spiral track set automatic sequencing analysis method, device, equipment and medium
By constructing an OVT surface mesh and coordinate system, and performing spiral numbering and encoding, the problem of low generation efficiency of OVT spiral gathers is solved, and efficient spiral gather sorting and coherence improvement are achieved.
Patent Information
- Authority / Receiving Office
- CN · China
- Patent Type
- Applications(China)
- Current Assignee / Owner
- CHINA PETROLEUM & CHEMICAL CORP
- Filing Date
- 2024-12-27
- Publication Date
- 2026-06-30
AI Technical Summary
Existing OVT spiral gather generation technologies are inefficient, manual coding is complex and prone to errors, while automatic coding disrupts the continuity and coherence of spiral coding.
By constructing an OVT surface grid and surface coordinate system, spiral numbering and encoding are performed, the relative coordinates of seismic traces are calculated, spiral trace numbers are matched and written into the trace head field, and sorted according to spiral trace numbers.
It improves the efficiency of encoding and sorting, and enhances the continuity of in-phase axes and the coherence of adjacent channels within the imaging gather.
Smart Images
Figure CN122307654A_ABST
Abstract
Description
Technical Field
[0001] This disclosure relates to the field of seismic exploration OVT data processing technology, and in particular to an automatic sorting analysis method, apparatus, equipment and medium for OVT spiral gathers. Background Technology
[0002] Before processing OVT data in seismic exploration, the seismic data needs to be divided into OVT elements, sorted into OVT gathers, and then migrated for imaging. Finally, the OVG imaging gathers are sorted into spiral gathers according to the spiral order of the OVT element distribution.
[0003] Existing spiral gather generation technologies include a manual coding and sorting method, in which the user manually calculates, codes, and adjusts the order of each OVT slice. However, the coding is complex, slow, and prone to errors. Another method is an automatic coding and sorting method, but it strictly follows the incremental coding of azimuth and offset, resulting in some traces "skipping" and disrupting the continuity of spiral coding. This reduces the coherence of adjacent traces within the gather, leading to a lower spiral coding speed and thus lower efficiency in spiral gather sorting. Summary of the Invention
[0004] This disclosure provides an automatic sorting and analysis method, apparatus, equipment, and medium for OVT spiral gathers to solve the problem of low efficiency in spiral gather sorting.
[0005] Firstly, this disclosure provides an automatic sorting and analysis method for OVT spiral gathers, including:
[0006] An OVT (Optical Transformer Dimension) grid is constructed based on the cross arrangement information and OVT grid division parameters in the pre-acquired seismic data, and an OVT grid coordinate system is constructed based on the OVT grid.
[0007] In the OVT surface element coordinate system, all data points in the OVT surface element grid are spirally numbered, and a spiral code table is generated based on the spirally numbered data points.
[0008] Calculate the relative coordinates of each seismic trace in the OVT surface coordinate system of the pre-acquired seismic trace set;
[0009] The spiral code table is used to match the spiral trace number corresponding to the relative coordinates of each seismic trace, and the spiral trace number corresponding to each seismic trace is written into the predefined trace header field.
[0010] The seismic gathers are sorted according to the spiral gather numbers written in the trace header field to obtain spiral gathers.
[0011] In some embodiments, constructing an OVT (Optical Virtual Transformer) mesh based on the cross-shaped arrangement information and OVT mesh division parameters in pre-acquired seismic data includes:
[0012] The shot line length and detector length are determined based on the cross arrangement information.
[0013] The shot spacing and detector spacing are determined based on the OVT element division parameters.
[0014] The number of surface elements in the direction of the shot line is calculated using the shot line length and the shot line spacing.
[0015] The number of detector directional elements is calculated using the detector length and the detector spacing.
[0016] The number of surface elements in the shot direction and the number of surface elements in the detector direction are combined to form an OVT surface element grid.
[0017] In some embodiments, the step of spirally numbering all data points in the OVT element grid in the OVT element coordinate system includes:
[0018] In the OVT surface coordinate system, determine the coordinate points corresponding to all data points in the OVT surface grid, and determine the spiral path corresponding to all data points based on the coordinate points;
[0019] When the number of data points is less than a preset threshold, the spiral tracks corresponding to all data points are numbered in turn to obtain the spiral number corresponding to the spiral track.
[0020] When the number of data points is greater than or equal to a preset threshold, the spiral path corresponding to all data points is encoded using preset encoding logic to obtain the spiral code corresponding to the spiral path.
[0021] In some embodiments, the step of sequentially numbering the spiral tracks corresponding to all data points to obtain the spiral number corresponding to the spiral track includes:
[0022] Construct a square divergence sequence according to the order of the first quadrant, fourth quadrant, third quadrant, and second quadrant in the OVT surface element coordinate system;
[0023] In the OVT surface element coordinate system, the target circle diverging outward from the origin of the coordinate system is determined according to the square divergence order, and it is detected whether there are any missing OVT surface elements in the target circle;
[0024] When a missing OVT element exists, the missing OVT element is skipped, and the data points in the target circle are sequentially numbered from the origin of the coordinate system outward in a square divergence order.
[0025] When the missing OVT element does not exist, the data points in the target circle are sequentially numbered from the origin of the coordinate system outward in a square divergence order.
[0026] When the target circle after outward divergence ends, extend outward to the starting position of the next circle corresponding to the target circle, based on the ending position of the target circle;
[0027] The origin of the coordinate system is updated to the starting position of the next circle, and the process returns to the step of determining the target circle diverging outward from the origin of the coordinate system in the OVT surface element coordinate system according to the square divergence order, until all data points are encoded and the spiral number corresponding to the spiral path is obtained.
[0028] In some embodiments, encoding the spiral paths corresponding to all data points using preset encoding logic to obtain the spiral code corresponding to the spiral path includes:
[0029] Extract the relative coordinates corresponding to the data points according to the encoding logic;
[0030] Calculate the maximum relative coordinate among the relative coordinates, and calculate the number of coding cycles based on the maximum relative coordinate;
[0031] When the number of coding cycles is a preset number of cycles threshold, the spiral code corresponding to the data point is output;
[0032] When the number of coding cycles is not a preset cycle threshold, the total coding value is calculated based on the number of coding cycles;
[0033] Match the relative coordinates with the coordinate rules corresponding to the spiral path in the OVT surface element coordinate system;
[0034] Once a match is successful, the spiral path coding logic formula is determined according to the coordinate rules, and the spiral code corresponding to the spiral path is calculated according to the spiral path coding logic formula and the total coding value.
[0035] In some embodiments, calculating the relative coordinates of each seismic trace in the pre-acquired seismic trace set within the OVT element coordinate system includes:
[0036] Calculate the offset and azimuth of each seismic trace in the seismic trace set;
[0037] The rotation angle of each seismic trace is calculated based on the azimuth angle and the predefined work area azimuth angle;
[0038] Calculate the gun line projection data of the offset distance in the gun line direction, and calculate the detector line projection data of the offset distance in the detector line direction;
[0039] The relative coordinates of each seismic trace in the OVT element coordinate system are calculated based on the shot line projection data, the receiver projection data, and the predefined OVT element size. The formula for calculating the relative coordinates is as follows:
[0040] X OVTi =Offset XL / Cell Detector +0.5
[0041] Y OVTi =Offset IL / Cell Source +0.5
[0042] Among them, X OVTi Y is the x-coordinate of the i-th seismic trace in the OVT element coordinate system. OVTi Let be the ordinate of the i-th seismic trace in the OVT element coordinate system, and Offset be the ordinate of the i-th seismic trace. XL For the projection data of the detector line, Offset IL For the projected data of the gun line, Cell Detector Cell is the detector line spacing within the OVT element size. Source The shot spacing in the OVT element size.
[0043] In some embodiments, sorting the seismic gathers according to the spiral gather number written in the trace header field to obtain spiral gathers includes:
[0044] Extract the horizontal survey line number, vertical survey line number, and spiral trace number corresponding to each seismic trace in the seismic trace set;
[0045] Generate key fields from the horizontal survey line direction number, the vertical survey line direction number, and the spiral track number;
[0046] Each seismic trace in the seismic trace set is sorted according to the trace number encoding order of the spiral trace number in the key field to obtain the spiral trace set.
[0047] Secondly, this disclosure provides an OVT spiral gather automatic sorting and analysis device, comprising:
[0048] The OVT surface coordinate system construction module is used to construct an OVT surface grid based on the cross arrangement information and OVT surface division parameters in the pre-acquired seismic data, and to construct an OVT surface coordinate system based on the OVT surface grid.
[0049] The spiral numbering module is used to spiral number all data points in the OVT surface element grid in the OVT surface element coordinate system, and generate a spiral code table based on the spiral numbered data points.
[0050] The relative coordinate calculation module is used to calculate the relative coordinates of each seismic trace in the pre-acquired seismic trace set in the OVT surface element coordinate system;
[0051] The spiral trace number matching module is used to match the spiral trace number corresponding to the relative coordinates of each seismic trace through the spiral encoding table, and write the spiral trace number corresponding to each seismic trace into a predefined trace header field;
[0052] The spiral gather generation module is used to sort the seismic gathers according to the spiral gather number written in the trace header field to obtain spiral gathers.
[0053] Thirdly, this disclosure provides a computer device including a memory, a processor, and a computer program stored in the memory, wherein the processor executes the computer program to implement the steps of the OVT spiral gather automatic sorting analysis method described in the above aspects.
[0054] Fourthly, this disclosure provides a computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the steps of the OVT spiral gather automatic sorting analysis method described above.
[0055] This disclosure provides an automatic sorting and analysis method, apparatus, device, and medium for OVT spiral gathers. It constructs an OVT grid based on the cross-shaped arrangement information and OVT grid division parameters from pre-acquired seismic data, and then constructs an OVT coordinate system based on the OVT grid. All data points in the OVT grid are spirally numbered within the OVT coordinate system, and a spiral coding table is generated based on these spirally numbered data points. The relative coordinates of each seismic trace in the pre-acquired seismic gather are calculated within the OVT coordinate system. The spiral coding table is used to match the spiral trace number corresponding to the relative coordinates of each seismic trace, and the spiral trace number corresponding to each seismic trace is written into a predefined trace header field. The seismic gathers are sorted according to the spiral trace numbers written into the trace header field to obtain spiral gathers. This avoids manual coding, improves the efficiency of coding and sorting, and enhances the continuity of in-phase axes and the coherence of adjacent traces within the imaging gather.
[0056] 1. The technical feature of constructing the OVT surface element coordinate system solves the problem of inaccurate data point positioning.
[0057] 2. The technical feature generates spiral numbering and spiral coding tables. The technical effect is to solve the problem of partial "track skipping" phenomenon, which disrupts the continuity of spiral coding.
[0058] 3. The technical feature calculation of the relative coordinates of seismic traces solves the problems of unclear seismic trace identification and difficulty in rapid positioning. Attached Figure Description
[0059] The present disclosure will be described in more detail below based on embodiments and with reference to the accompanying drawings:
[0060] Figure 1 A flowchart illustrating the automatic sorting and analysis method for OVT spiral gathers provided in this embodiment of the disclosure;
[0061] Figure 2 This is a schematic diagram of the OVT surface element coordinate system and spiral encoding provided in the embodiments of this disclosure;
[0062] Figure 3 This is a schematic diagram of the cross arrangement and OVT element division provided in the embodiments of this disclosure;
[0063] Figure 4 A schematic diagram of the pseudocode for calculating the sequence number of a point in the OVT surface element coordinate system provided in this embodiment of the disclosure;
[0064] Figure 5 A schematic diagram of the automatic spiral gather sorting process provided in the embodiments of this disclosure;
[0065] Figure 6 A schematic diagram of the automatic encoding and sorting parameters for the OVT spiral gather provided in the embodiments of this disclosure;
[0066] Figure 7 A schematic diagram comparing the gathers before and after spiral gather sorting, provided in an embodiment of this disclosure;
[0067] Figure 8 A schematic diagram of the azimuth angle and offset distance of the spiral gather provided in the embodiments of this disclosure;
[0068] Figure 9 A functional block diagram of the OVT spiral gather automatic sorting and analysis device provided in the embodiments of this disclosure.
[0069] In the accompanying drawings, the same parts are referred to by the same reference numerals, and the drawings are not drawn to scale. Detailed Implementation
[0070] To enable those skilled in the art to better understand the technical solutions of this disclosure, and to fully understand and implement the process of how this disclosure applies technical means to solve technical problems and achieve corresponding technical effects, the technical solutions in the embodiments of this disclosure will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of this disclosure, not all embodiments. The embodiments of this disclosure and the various features within them can be combined with each other without conflict, and the resulting technical solutions are all within the protection scope of this disclosure. All other embodiments obtained by those skilled in the art based on the embodiments of this disclosure without creative effort should fall within the protection scope of this disclosure.
[0071] It should be noted that the terms "first," "second," etc., in the specification, claims, and accompanying drawings of this disclosure are used to distinguish similar objects and are not necessarily used to describe a specific order or sequence. It should be understood that such data can be interchanged where appropriate so that the embodiments of this disclosure described herein can be implemented in orders other than those illustrated or described herein. Furthermore, the terms "comprising" and "having," and any variations thereof, are intended to cover non-exclusive inclusion; for example, a process, method, system, product, or apparatus that comprises a series of steps or units is not necessarily limited to those steps or units explicitly listed, but may include other steps or units not explicitly listed or inherent to such processes, methods, products, or apparatus.
[0072] It should be noted that the steps shown in the flowchart in the accompanying drawings can be executed in a computer system such as a set of computer-executable instructions, and although a logical order is shown in the flowchart, in some cases the steps shown or described may be executed in a different order than that shown here.
[0073] Example 1
[0074] Figure 1 This is a flowchart illustrating an automatic sorting and analysis method for OVT spiral gathers provided in an embodiment of this disclosure. Figure 1 As shown, an automatic sorting and analysis method for OVT spiral gathers includes:
[0075] S1. Construct an OVT (Optical Variable Transformer) grid based on the cross arrangement information and OVT grid division parameters in the pre-acquired seismic data, and construct an OVT grid coordinate system based on the OVT grid.
[0076] In this embodiment of the invention, the cross-shaped arrangement information includes the shot line length and the geophone length. The shot line length refers to the straight line length along which the shot points are arranged at the seismic data acquisition site, determining the total span of data coverage in the shot line direction. The geophone length is the straight line length along which the geophones are arranged, defining the coverage range in the geophone line direction. The OVT (Offset Vector Tile) surface division parameters include the shot line spacing and the geophone line spacing. The shot line spacing refers to the distance between two adjacent shot lines, determining the density of the grid in the shot line direction. The geophone line spacing refers to the distance between two adjacent geophone lines, determining the density of the grid in the geophone line direction. The OVT surface grid is a spatial division method used to organize and analyze seismic data, dividing the seismic data acquisition area into many small rectangular units, each with its specific location and range.
[0077] In this embodiment of the invention, the step of constructing an OVT (Optical Transformer Node) mesh based on the cross-shaped arrangement information and OVT mesh division parameters in the pre-acquired seismic data includes:
[0078] The shot line length and detector length are determined based on the cross arrangement information.
[0079] The shot spacing and detector spacing are determined based on the OVT element division parameters.
[0080] The number of surface elements in the direction of the shot line is calculated using the shot line length and the shot line spacing.
[0081] The number of detector directional elements is calculated using the detector length and the detector spacing.
[0082] The number of surface elements in the shot direction and the number of surface elements in the detector direction are combined to form an OVT surface element grid.
[0083] In detail, during seismic data acquisition, a cross-shaped arrangement is typically used. The shot line is a straight line through which shot points (locations that generate seismic waves, such as the explosion point of explosives or the location of artificial seismic sources) are placed, and the shot line length is the actual length of this straight line. The geophone length is the straight line through which geophones (used to receive seismic waves) are placed. This length information determines the spatial coverage of the seismic data and forms the basis for subsequent construction of the surface grid. The shot line spacing is the distance between two adjacent shot lines, and the geophone spacing is the distance between two adjacent geophones. These spacing parameters are used to define the size of each surface grid.
[0084] Specifically, the number of elements in the shot direction = (shot length / shot spacing) + 1, where 1 is added to include the boundary conditions at both ends. The number of elements in the detector direction = (detector length / detector spacing) + 1. After determining the number of elements in the shot and detector directions, the number of detector elements is determined in the X direction, and the number of shot elements is determined in the Y direction. The number of elements in the X and Y directions is intersected to divide the grid into small rectangles, i.e., the OVT element grid. Each element has a unique position in this grid, which can be identified by its index in the shot and detector directions.
[0085] Furthermore, in order to analyze the spiral numbers corresponding to all data points in the OVT surface grid, it is necessary to establish an OVT surface coordinate system based on the OVT surface grid, and analyze the spiral numbers of the data points in the coordinate system, so as to facilitate the location and processing of seismic data within this surface.
[0086] In this embodiment of the invention, the OVT element coordinate system is a coordinate system that describes the spatial position of the OVT element, providing a unique coordinate identifier for each OVT element, enabling accurate positioning and manipulation of specific elements when processing seismic data.
[0087] In this embodiment of the invention, constructing the OVT surface element coordinate system based on the OVT surface element mesh includes:
[0088] Identify the directional dimension of the OVT surface mesh;
[0089] The center of the OVT surface element mesh is taken as the origin of the coordinate system;
[0090] An OVT surface coordinate system is constructed on the OVT surface grid based on the coordinate origin and the direction dimension.
[0091] In detail, determine the dimensions NX of the mesh in both directions. OVT NY OVT Using the detector direction as the X dimension NX OVT Using the direction of the shot line as the Y dimension NY OVT An OVT element coordinate system is established with the center of the OVT element mesh as the origin (0, 0). In the shot direction, the x-coordinate of the element to the left of the origin is negative, and the x-coordinate to the right is positive. In the receiver direction, the y-coordinate of the element below the origin is negative, and the y-coordinate above is positive. Once the origin and direction dimensions are determined, the OVT element coordinate system can be constructed with the shot direction as the X-axis and the receiver direction as the Y-axis. The positive direction of the coordinate axes is defined, and then coordinates are assigned to each element according to its relative position to the origin in the shot and receiver directions.
[0092] For example, suppose a surface element is 3 surface element spacings away from the origin in the shot direction (and in the positive direction), and 2 surface element spacings away from the origin in the detector direction (also in the positive direction). If the surface element spacings in the shot direction and the detector direction are 100 meters and 80 meters respectively, then the coordinates of this surface element in the OVT surface element coordinate system may be (300, 160). In this way, a coordinate system in the OVT surface element coordinate system is constructed for each surface element, so that all surface elements can be conveniently located, analyzed and processed in this coordinate system.
[0093] Furthermore, in the traditional automatic encoding and sorting method, encoding is strictly performed according to the increasing azimuth angle and offset distance. This results in some traces "skipping traces", which disrupts the continuity of spiral encoding and reduces the coherence of adjacent traces within the trace set. Therefore, spiral encoding of the OVT surface element mesh should be performed in the OVT surface element coordinate system.
[0094] S2. In the OVT surface element coordinate system, all data points in the OVT surface element grid are spirally numbered, and a spiral code table is generated based on the spirally numbered data points.
[0095] In this embodiment of the invention, spiral numbering refers to assigning a unique identifier or code to each data point in the OVT surface grid. The purpose is to facilitate the differentiation, management, and subsequent data processing and analysis of numerous data points.
[0096] In this embodiment of the invention, the step of spirally numbering all data points in the OVT surface element grid in the OVT surface element coordinate system includes:
[0097] In the OVT surface coordinate system, determine the coordinate points corresponding to all data points in the OVT surface grid, and determine the spiral path corresponding to all data points based on the coordinate points;
[0098] When the number of data points is less than a preset threshold, the spiral tracks corresponding to all data points are numbered in turn to obtain the spiral number corresponding to the spiral track.
[0099] When the number of data points is greater than or equal to a preset threshold, the spiral path corresponding to all data points is encoded using preset encoding logic to obtain the spiral code corresponding to the spiral path.
[0100] In detail, in the previously constructed OVT surface element coordinate system, each data point can determine its corresponding coordinate point based on its location, and then determine the corresponding spiral track based on the coordinate point. For example, let the origin of the OVT surface element coordinate system (0, 0) be the first spiral track, and (0, 1) be the second spiral track.
[0101] Specifically, when the number of data points is less than a preset threshold, the spiral tracks can be numbered sequentially. If there are 5 spiral tracks, the first spiral track is numbered 1, the second spiral track is numbered 2, and so on, until the 5th spiral track is numbered 5. This way, each spiral track receives a spiral number, allowing data points to be distinguished and processed based on their spiral track number. However, when the number of data points is large, sequentially numbering the spiral tracks is inefficient. Therefore, it is necessary to utilize specific functional logic relationships to generate corresponding spiral codes for numerous spiral tracks, such as... Figure 2 The diagram shows the OVT surface element coordinate system and spiral coding. With coordinate (0, 0) as the origin, the spirals are numbered sequentially in the order of quadrants 1, 4, 3, 2, 1, 4, 3, 2, 1... Then, the spiral number corresponding to (0, 0) is 1, the spiral number corresponding to (0, 1) is 2, the spiral number corresponding to (1, 1) is 3, the spiral number corresponding to (1, 0) is 4, and so on, until all data points are encoded, thus generating a spiral coding table.
[0102] Furthermore, after generating corresponding spiral numbers for all data points in the OVT surface mesh, an index table is formed between the relative coordinates of the OVT mesh and the spiral sequence number, such as ((0,0),1), ((0,1),2), ((1,1),3), ((1,0),4), etc.
[0103] Furthermore, in order to calculate the spiral trace number corresponding to the actual seismic trace, it is necessary to divide each seismic trace in the seismic trace set into OVT surface elements to calculate the relative coordinates corresponding to each seismic trace, and then calculate the spiral trace number corresponding to each seismic trace based on the relative coordinates.
[0104] S3. Calculate the relative coordinates of each seismic trace in the OVT surface coordinate system of the pre-acquired seismic trace set.
[0105] In this embodiment of the invention, the relative coordinates refer to the coordinate values of each OVT surface element after the seismic trace is divided in the OVT surface element coordinate system.
[0106] In this embodiment of the invention, calculating the relative coordinates of each seismic trace in the pre-acquired seismic trace set within the OVT surface coordinate system includes:
[0107] Calculate the offset and azimuth of each seismic trace in the seismic trace set;
[0108] The rotation angle of each seismic trace is calculated based on the azimuth angle and the predefined work area azimuth angle;
[0109] Calculate the gun line projection data of the offset distance in the gun line direction, and calculate the detector line projection data of the offset distance in the detector line direction;
[0110] The relative coordinates of each seismic trace in the OVT element coordinate system are calculated based on the shot line projection data, the receiver projection data, and the predefined OVT element size. The formula for calculating the relative coordinates is as follows:
[0111] X OVTi =Offset XL / Cell Detector +0.5
[0112] Y OVTi =Offset IL / Cell Source +0.5
[0113] Among them, X OVTi Y is the x-coordinate of the i-th seismic trace in the OVT element coordinate system. OVTi Let be the ordinate of the i-th seismic trace in the OVT element coordinate system, and Offset be the ordinate of the i-th seismic trace. XL For the projection data of the detector line, Offset IL For the projected data of the gun line, Cell Detector Cell is the detector line spacing within the OVT element size. Source The shot spacing in the OVT element size.
[0114] In detail, the offset of each seismic trace in the seismic trace set is calculated. Trace and azimuth angle Azi Trace According to the work area azimuth Azi defined by the grid Grid Calculate the rotation angle for each seismic trace, i.e., the rotation angle is Azi. Rotate =Azi Trace -Azi Grid And calculate the projection of the offset of each seismic trace onto the cross-shaped shot line and receiver line, then the projection in the shot line direction is Offset. IL The projection of the detector line direction is Offset XL Then Offset IL =Offset Trace *cos(Azi Rotate ), Offset XL =Offset Trace *sin(Azi Rotate Then, based on the defined OVT element size, the relative coordinates (X, Y, F) of each seismic trace in the OVT element coordinate system are calculated. OVTi Y OVTi ).
[0115] Specifically, such as Figure 3 The diagram shown is a cross-shaped arrangement and a schematic diagram of the OVT (Outer Threshold Variable) element division. Figure 3 The OVT (Out-of-Touch) element represents the area between the shot line and the receiver. Multiple dashed lines connect the shot and receiver points, and these dashed lines represent the path of seismic waves propagating from the shot to the receiver. Each dashed line represents a seismic trace, that is, the signal path recorded when a seismic wave generated by a shot is received by a receiver. The OVT point is connected to some shot and receiver points through dashed lines, indicating the seismic wave propagation path involved in this OVT and the corresponding shot and receiver points, that is, the OVT element divided by each seismic trace.
[0116] Furthermore, based on the coordinates (X... OVTi Y OVTi Search for and match the "spiral code table", find the corresponding spiral channel number, and write it into the channel header word OVT_NUM.
[0117] S4. Match the spiral trace number corresponding to the relative coordinates of each seismic trace using the spiral coding table, and write the spiral trace number corresponding to each seismic trace into the predefined trace head field.
[0118] In this embodiment of the invention, for each seismic trace, the relative coordinates (X, X) in the determined OVT surface element coordinate system are used. OVTi Y OVTi Find the spiral trace number corresponding to this relative coordinate in the spiral coding table. For example, if the relative coordinate is (1, 0) and the spiral trace number corresponding to (1, 0) in the spiral coding table is 4, then the spiral number corresponding to this seismic trace is 4. Thus, the spiral number corresponding to each seismic trace is matched. Once the spiral trace number corresponding to each seismic trace is matched, write this spiral trace number into the predefined trace head field.
[0119] In detail, in seismic data formats, the trace header contains a lot of information describing the attributes of seismic traces, such as sampling rate, trace number, line number, etc. The trace header field is a predefined location specifically for storing spiral trace numbers. By writing the spiral trace number into it, the seismic trace can be easily identified and manipulated according to this spiral trace number when processing seismic data. For example, seismic traces can be sorted and grouped according to the spiral trace number, which is helpful for processing such as azimuth anisotropy analysis.
[0120] Furthermore, in order to make the seismic gathers usable for subsequent azimuth anisotropy processing, all seismic gathers need to be reordered after they are numbered in spiral order and written into the gather header to obtain the required spiral gathers.
[0121] S5. Sort the seismic gathers according to the spiral gather numbers written in the trace header field to obtain spiral gathers.
[0122] In this embodiment of the invention, the spiral gather refers to a set of seismic traces organized in a specific order (spiral order). In seismic data processing, by sorting the seismic traces according to the spiral trace number, the seismic data can exhibit certain regularities in subsequent processing (such as azimuth anisotropy analysis).
[0123] In this embodiment of the invention, sorting the seismic gathers according to the spiral gather numbers written in the trace header field to obtain spiral gathers includes:
[0124] Extract the horizontal survey line number, vertical survey line number, and spiral trace number corresponding to each seismic trace in the seismic trace set;
[0125] Generate key fields from the horizontal survey line direction number, the vertical survey line direction number, and the spiral track number;
[0126] Each seismic trace in the seismic trace set is sorted according to the trace number encoding order of the spiral trace number in the key field to obtain the spiral trace set.
[0127] In detail, during seismic data acquisition, the horizontal survey line number represents the line number along the horizontal direction (generally parallel to the acquisition arrangement direction), so the horizontal survey line number is Inline. The vertical survey line number represents the line number perpendicular to the horizontal survey line direction, so the vertical survey line number is Xline. Then, the horizontal survey line number, the vertical survey line number, and the spiral trace number are used to generate key fields, so the key fields are (Inline, Xline, OVT_NUM).
[0128] Specifically, each seismic trace in the seismic trace set is sorted according to the trace number encoding order of the spiral trace number in the key field. That is, the seismic traces are rearranged according to the encoding order of the spiral trace number. For example, if the spiral trace number corresponding to seismic trace A is 1, then seismic trace A is placed at the first position in the spiral trace set. In this way, all seismic traces are sorted according to the spiral trace number order in the seismic trace set, and the sorted seismic trace set is used as the spiral trace set.
[0129] Example 2
[0130] Based on the above embodiments, the step of sequentially numbering the spiral tracks corresponding to all data points to obtain the spiral number corresponding to the spiral track includes:
[0131] Construct a square divergence sequence according to the order of the first quadrant, fourth quadrant, third quadrant, and second quadrant in the OVT surface element coordinate system;
[0132] In the OVT surface element coordinate system, the target circle diverging outward from the origin of the coordinate system is determined according to the square divergence order, and it is detected whether there are any missing OVT surface elements in the target circle;
[0133] When a missing OVT element exists, the missing OVT element is skipped, and the data points in the target circle are sequentially numbered from the origin of the coordinate system outward in a square divergence order.
[0134] When the missing OVT element does not exist, the data points in the target circle are sequentially numbered from the origin of the coordinate system outward in a square divergence order.
[0135] When the target circle after outward divergence ends, extend outward to the starting position of the next circle corresponding to the target circle, based on the ending position of the target circle;
[0136] The origin of the coordinate system is updated to the starting position of the next circle, and the process returns to the step of determining the target circle diverging outward from the origin of the coordinate system in the OVT surface element coordinate system according to the square divergence order, until all data points are encoded and the spiral number corresponding to the spiral path is obtained.
[0137] In detail, during spiral numbering, the data points need to be numbered sequentially in the order of the first quadrant, fourth quadrant, third quadrant, and second quadrant, and then gradually diverge outwards from the origin in a square pattern, i.e., a square divergence sequence. Within the OVT surface element coordinate system, the target circles diverging outwards from the origin are determined according to this square divergence sequence. Centered on the origin, the outward divergence forms concentric circles, each of which can be considered a target circle. For example, the data points in the circle closest to the origin constitute the first target circle, the data points in the next outermost layer belong to the second target circle, and so on. Determining these target circles is for the purpose of numbering the data points within each circle according to a set order and rules.
[0138] Specifically, due to various reasons (such as omissions during data acquisition, filtering during data processing, etc.), there may be OVT facets that should exist but are actually missing in a certain target circle. That is, when a circle is not completed, the quadrants are still numbered in the order of 2→1→4→3→2…, but the missing OVT facets are skipped, and the sequence number is uninterrupted. When the data points in a target circle are numbered according to the rules, that is, when the target circle after outward divergence ends, it is necessary to naturally extend outward to the starting position of the next circle (that is, the next larger circle corresponding to the target circle) based on the ending position of this target circle. For example, when each circle ends (-Ni, Ni), it extends outward to the starting position of the next circle (-Ni, Ni+1). Figure 2Extending outward from (-1, 1) to (-1, 2), a new loop begins. Starting from the beginning of the new loop, the data points in the new loop are sequentially numbered until all data points of the OVT surface mesh have been traversed, thus obtaining the spiral numbering corresponding to all data points of the OVT surface mesh.
[0139] Example 3
[0140] Based on the above embodiments, the step of encoding the spiral path corresponding to all data points using preset encoding logic to obtain the spiral code corresponding to the spiral path includes:
[0141] Extract the relative coordinates corresponding to the data points according to the encoding logic;
[0142] Calculate the maximum relative coordinate among the relative coordinates, and calculate the number of coding cycles based on the maximum relative coordinate;
[0143] When the number of coding cycles is a preset number of cycles threshold, the spiral code corresponding to the data point is output;
[0144] When the number of coding cycles is not a preset cycle threshold, the total coding value is calculated based on the number of coding cycles;
[0145] Match the relative coordinates with the coordinate rules corresponding to the spiral path in the OVT surface element coordinate system;
[0146] Once a match is successful, the spiral path coding logic formula is determined according to the coordinate rules, and the spiral code corresponding to the spiral path is calculated according to the spiral path coding logic formula and the total coding value.
[0147] In detail, when (Mx, My) is a coordinate point in the OVT surface element coordinate system, the sequence number of the point is calculated through pseudocode logic, that is, the relative coordinate is calculated based on (Mx, My). The relative X coordinate of the data point is determined by ABSX = Abs(Mx), the relative Y coordinate of the data point is determined by ABSY = Abs(My), the maximum relative coordinate in the relative coordinates is calculated by M = Max(ABSX, ABSY), and the current lap number is calculated by N = M + 1, which represents the Nth lap.
[0148] Specifically, based on the decision branch, it checks whether the current coding circle number N is equal to 1. If (N == 1), the spiral coding is OVT_NUM = 1; otherwise, it calculates the total coding value based on the current coding circle number, i.e., total = (2M-1)*(2M-1). Then, according to the coordinate rules, it compares the coordinates of the data points with the maximum relative coordinates, first comparing them with the coordinate values corresponding to the first track of each circle. That is, if (-Mx == M-1 && My == M), the spiral coding is OVT. NUM = total + 1; If the relative coordinate value does not conform to the coordinate rule of the first lane of each lap, it is compared with the coordinate value corresponding to the last lane of each lap, i.e., if (-Mx == M && My == M), the spiral code is OVT_NUM = (2N-1)*(2N-1); If the relative coordinate value also does not conform to the coordinate rule of the last lane of each lap, it is compared with the coordinate value corresponding to the top, i.e., if (My == M), the distance to the first lane of the current lap is calculated, i.e., dis = Mx + (M-1), and the spiral code is determined to be OVT NUM = total + 1 + dis; If the relative coordinate value also does not conform to the coordinate rule of the top, it is compared with the coordinate value corresponding to the left, i.e., if (Mx == -M), the distance to the last lane of the current lap is calculated, i.e., dis = My - M, and the spiral code is determined to be OVT. NUM = N*N + dis; If the relative coordinate value does not conform to the coordinate rules of the left part, then compare it with the corresponding coordinate value of the right part. That is, if (Mx == M), calculate the distance to the upper right inflection point, i.e., dis = M - My, and determine the spiral code as OVT_NUM = total + 1 + (2 * M - 1) + dis; If the relative coordinate value does not conform to the coordinate rules of the right part, then compare it with the corresponding coordinate value of the bottom part. That is, if (My == - M), calculate the distance to the lower right inflection point, i.e., dis = M - Mx, and determine the spiral code as total + 1 + (4 * M - 1) + dis, such as Figure 4 The diagram shows a pseudocode for calculating the sequence number of a point in the OVT surface coordinate system. The diagram calculates the spiral sequence number corresponding to a point in the OVT surface coordinate system based on different judgment conditions.
[0149] Example 4
[0150] Based on the above embodiments, this embodiment provides an application example.
[0151] The automatic encoding and sorting function of spiral gathers was tested using OVT imaging gather data from a certain work area. Figure 5 The diagram shows the flowchart of automatic spiral gather sorting. Based on the cross arrangement information and OVT element partitioning parameters, a spiral gather coordinate system and spiral code table are established using spiral gather coding. One trace is taken, and the relative coordinates (X, Y, X) of the current trace in the OVT element coordinate system are calculated.OVTi Y OVTi ), query the "spiral coding table" to obtain the current trace number, and write it into the trace head OVT_NUM. After processing all seismic traces, reorder the internal structures of the imaging gathers according to the new trace head OVT_NUM to obtain the spiral gathers; for example Figure 6 The diagram shows the automatic encoding and sorting parameters for the OVT spiral gather. It configures the spacing between detector lines and source lines, outputs cross-expansion attributes, and sets the maximum offset (Max Offset) to 10400 and 8000, and the tile size to 520 and 400. Figure 7 As shown, this is a comparative diagram of the gathers before and after spiral gather sorting. The continuity of in-phase axes and the coherence of adjacent gathers are significantly enhanced after spiral gather sorting, which is beneficial for subsequent anisotropy analysis and processing; as shown... Figure 8 The diagram shown illustrates the azimuth and offset of the spiral gather. Figure 8 It can be seen that the azimuth of the spiral gather increases from 0 to 360 degrees, while the offset generally shows an increasing trend, but "oscillates" within one azimuth period, which is consistent with the arrangement pattern of the spiral gather.
[0152] Example 5
[0153] like Figure 9 As shown in the figure, this embodiment also provides a functional block diagram of an OVT spiral gather automatic sorting and analysis device.
[0154] The OVT spiral gather automatic sorting analysis and generation device 100 described in this embodiment can be installed in an electronic device. Depending on the functions implemented, the OVT spiral gather automatic sorting analysis device 100 may include an OVT surface element coordinate system construction module 101, a spiral numbering module 102, a relative coordinate calculation module 103, a spiral gather number matching module 104, and a spiral gather generation module 105. The module described in this invention can also be called a unit, referring to a series of computer program segments that can be executed by the processor of an electronic device and can perform a fixed function, stored in the memory of the electronic device.
[0155] In this embodiment, the functions of each module / unit are as follows:
[0156] The OVT surface coordinate system construction module 101 is used to construct an OVT surface grid based on the cross arrangement information and OVT surface division parameters in the pre-acquired seismic data, and to construct an OVT surface coordinate system based on the OVT surface grid.
[0157] The spiral numbering module 102 is used to spiral number all data points in the OVT surface element grid in the OVT surface element coordinate system, and generate a spiral code table based on the spiral numbered data points.
[0158] The relative coordinate calculation module 103 is used to calculate the relative coordinates of each seismic trace in the pre-acquired seismic trace set in the OVT surface element coordinate system;
[0159] The spiral trace number matching module 104 is used to match the spiral trace number corresponding to the relative coordinates of each seismic trace through the spiral encoding table, and write the spiral trace number corresponding to each seismic trace into a predefined trace header field;
[0160] The spiral gather generation module 105 is used to sort the seismic gathers according to the spiral gather number written in the trace header field to obtain spiral gathers.
[0161] In detail, each module in the OVT spiral gather automatic sorting and analysis device 100 described in the embodiments of the present invention adopts the same technical means as the OVT spiral gather automatic sorting and analysis method described in Embodiments 1 to 4, and can produce the same technical effect, which will not be repeated here.
[0162] Example 6
[0163] Based on the above embodiments, this embodiment provides a computer device, including a memory, a processor, and a computer program stored in the memory, wherein the processor executes the computer program to implement the steps of the method described in the above embodiments.
[0164] In some embodiments of this example, a computer-readable storage medium is provided, on which a computer program is stored, characterized in that the computer program, when executed by a processor, implements the steps of the method described in the above embodiments.
[0165] In some embodiments of this example, a computer program product is provided, including a computer program / instructions, characterized in that the computer program, when executed by a processor, implements the steps of the method described in the above embodiments.
[0166] The processor may include, but is not limited to, one or more processors or microprocessors. Each processor may be implemented as an Application Specific Integrated Circuit (ASIC), Digital Signal Processor (DSP), Digital Signal Processing Device (DSPD), Programmable Logic Device (PLD), Field Programmable Gate Array (FPGA), controller, microcontroller, microprocessor, or other electronic component, for executing the methods in the above embodiments.
[0167] Computer-readable storage media can be implemented by any type of volatile or non-volatile storage device or a combination thereof. Computer-readable storage media may include, but are not limited to, random access memory (RAM), read-only memory (ROM), flash memory, EPROM memory, EEPROM memory, registers, and computer storage media (e.g., hard disks, floppy disks, solid-state drives, removable disks, CD-ROMs, DVD-ROMs, Blu-ray discs, etc.).
[0168] Computer-readable storage media may also store at least one computer-executable program / instruction, such as computer-readable instructions. Computer-readable storage media include, but are not limited to, volatile memory and / or non-volatile memory. Volatile memory may include, for example, random access memory (RAM) and / or cache memory. Computer-readable storage media may include, for example, read-only memory (ROM), hard disk, flash memory, etc. For example, a non-transitory computer-readable storage medium may be connected to a computing device such as a computer, and then, when the computing device executes the computer-readable instructions stored on the computer-readable storage medium, the various methods described above can be performed.
[0169] In addition, the computer device may include (but is not limited to) a data bus, an input / output (I / O) bus, a display, and input / output devices (e.g., keyboard, mouse, speakers, etc.).
[0170] The processor can communicate with external devices via the I / O bus through wired or wireless networks.
[0171] In one embodiment, the at least one computer-executable instruction may also be compiled into or comprise a software product / computer program product, wherein one or more computer-executable instructions are executed by a processor to perform the steps of the various functions and / or methods in the embodiments described herein.
[0172] In the embodiments provided in this disclosure, it should be understood that the disclosed apparatus and methods can also be implemented in other ways. The apparatus embodiments described above are merely illustrative; for example, the flowcharts and block diagrams in the accompanying drawings illustrate the architecture, functionality, and operation of possible implementations of apparatus, methods, and computer program products according to various embodiments of this disclosure. In this regard, each block in a flowchart or block diagram may represent a module, segment, or portion of code containing one or more executable instructions for implementing a specified logical function. It should also be noted that in some alternative implementations, the functions marked in the blocks may occur in a different order than those marked in the drawings. For example, two consecutive blocks may actually be executed substantially in parallel, and they may sometimes be executed in reverse order, depending on the functions involved. It should also be noted that each block in a block diagram and / or flowchart, and combinations of blocks in block diagrams and / or flowcharts, can be implemented using a dedicated hardware-based system that performs the specified function or action, or using a combination of dedicated hardware and computer instructions.
[0173] It should be noted that, in this disclosure, the terms "comprising," "including," or any other variations thereof are intended to cover non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements includes not only those elements but also other elements not expressly listed, or elements inherent to such a process, method, article, or apparatus. Without further limitation, an element limited by the phrase "comprising one..." does not exclude the presence of other identical elements in the process, method, article, or apparatus that includes that element.
[0174] While the embodiments disclosed herein are as described above, the foregoing content is merely for the purpose of facilitating understanding of this disclosure and is not intended to limit this disclosure. Any person skilled in the art to which this disclosure pertains may make any modifications and changes in form and detail of the implementation without departing from the spirit and scope of this disclosure; however, the scope of patent protection of this disclosure shall still be determined by the scope defined in the appended claims.
Claims
1. An automatic sorting and analysis method for OVT spiral gathers, characterized in that, include: An OVT (Optical Transformer Dimension) grid is constructed based on the cross arrangement information and OVT grid division parameters in the pre-acquired seismic data, and an OVT grid coordinate system is constructed based on the OVT grid. In the OVT surface element coordinate system, all data points in the OVT surface element grid are spirally numbered, and a spiral code table is generated based on the spirally numbered data points. Calculate the relative coordinates of each seismic trace in the OVT surface coordinate system of the pre-acquired seismic trace set; The spiral code table is used to match the spiral trace number corresponding to the relative coordinates of each seismic trace, and the spiral trace number corresponding to each seismic trace is written into the predefined trace header field. The seismic gathers are sorted according to the spiral gather numbers written in the trace header field to obtain spiral gathers.
2. The automatic sorting and analysis method for OVT spiral gathers according to claim 1, characterized in that, The step of constructing an OVT (Optical Virtual Transformer) mesh based on the cross-shaped arrangement information and OVT mesh division parameters in the pre-acquired seismic data includes: The shot line length and detector length are determined based on the cross arrangement information. The shot spacing and detector spacing are determined based on the OVT element division parameters. The number of surface elements in the direction of the shot line is calculated using the shot line length and the shot line spacing. The number of detector directional elements is calculated using the detector length and the detector spacing. The number of surface elements in the shot direction and the number of surface elements in the detector direction are combined to form an OVT surface element grid.
3. The automatic sorting and analysis method for OVT spiral gathers according to claim 1, characterized in that, The step of spirally numbering all data points in the OVT element grid in the OVT element coordinate system includes: In the OVT surface coordinate system, determine the coordinate points corresponding to all data points in the OVT surface grid, and determine the spiral path corresponding to all data points based on the coordinate points; When the number of data points is less than a preset threshold, the spiral tracks corresponding to all data points are numbered in turn to obtain the spiral number corresponding to the spiral track. When the number of data points is greater than or equal to a preset threshold, the spiral path corresponding to all data points is encoded using preset encoding logic to obtain the spiral code corresponding to the spiral path.
4. The automatic sorting and analysis method for OVT spiral gathers according to claim 3, characterized in that, The step of sequentially numbering the spiral paths corresponding to all data points to obtain the spiral number corresponding to the spiral path includes: Construct a square divergence sequence according to the order of the first quadrant, fourth quadrant, third quadrant, and second quadrant in the OVT surface element coordinate system; In the OVT surface element coordinate system, the target circle diverging outward from the origin of the coordinate system is determined according to the square divergence order, and it is detected whether there are any missing OVT surface elements in the target circle; When a missing OVT element exists, the missing OVT element is skipped, and the data points in the target circle are sequentially numbered from the origin of the coordinate system outward in a square divergence order. When the missing OVT element does not exist, the data points in the target circle are sequentially numbered from the origin of the coordinate system outward in a square divergence order. When the target circle after outward divergence ends, extend outward to the starting position of the next circle corresponding to the target circle, based on the ending position of the target circle; The origin of the coordinate system is updated to the starting position of the next circle, and the process returns to the step of determining the target circle diverging outward from the origin of the coordinate system in the OVT surface element coordinate system according to the square divergence order, until all data points are encoded and the spiral number corresponding to the spiral path is obtained.
5. The automatic sorting and analysis method for OVT spiral gathers according to claim 3, characterized in that, The step of encoding the spiral path corresponding to all data points through a preset encoding logic to obtain the spiral code corresponding to the spiral path includes: Extract the relative coordinates corresponding to the data points according to the encoding logic; Calculate the maximum relative coordinate among the relative coordinates, and calculate the number of coding cycles based on the maximum relative coordinate; When the number of coding cycles is a preset number of cycles threshold, the spiral code corresponding to the data point is output; When the number of coding cycles is not a preset cycle threshold, the total coding value is calculated based on the number of coding cycles; Match the relative coordinates with the coordinate rules corresponding to the spiral path in the OVT surface element coordinate system; Once a match is successful, the spiral path coding logic formula is determined according to the coordinate rules, and the spiral code corresponding to the spiral path is calculated according to the spiral path coding logic formula and the total coding value.
6. The automatic sorting and analysis method for OVT spiral gathers according to claim 1, characterized in that, The calculation of the relative coordinates of each seismic trace in the OVT surface coordinate system of the pre-acquired seismic trace set includes: Calculate the offset and azimuth of each seismic trace in the seismic trace set; The rotation angle of each seismic trace is calculated based on the azimuth angle and the predefined azimuth angle of the work area; Calculate the gun line projection data of the offset distance in the gun line direction, and calculate the detector line projection data of the offset distance in the detector line direction; The relative coordinates of each seismic trace in the OVT element coordinate system are calculated based on the shot line projection data, the receiver projection data, and the predefined OVT element size. The formula for calculating the relative coordinates is as follows: X OVTi =Offset XL / Cell Detector +0.5 Y OVTi =Offset IL / Cell Source +0.5 Among them, X OVTi Y is the x-coordinate of the i-th seismic trace in the OVT element coordinate system. OVTi Let be the ordinate of the i-th seismic trace in the OVT element coordinate system, and Offset be the ordinate of the i-th seismic trace. XL For the projection data of the detector line, Offset IL For the projected data of the shot line, Cell Detector Cell is the detector line spacing within the OVT element size. Source The shot spacing in the OVT element size.
7. The automatic sorting and analysis method for OVT spiral gathers according to claim 1, characterized in that, The process of sorting the seismic gathers according to the spiral gather number written in the trace header field to obtain spiral gathers includes: Extract the horizontal survey line number, vertical survey line number, and spiral trace number corresponding to each seismic trace in the seismic trace set; Generate key fields from the horizontal survey line direction number, the vertical survey line direction number, and the spiral track number; Each seismic trace in the seismic trace set is sorted according to the trace number encoding order of the spiral trace number in the key field to obtain the spiral trace set.
8. An OVT spiral gather automatic sorting and analysis device, characterized in that, include: O The VT surface coordinate system construction module is used to construct an OVT surface grid based on the cross arrangement information and OVT surface division parameters in the pre-acquired seismic data, and to construct an OVT surface coordinate system based on the OVT surface grid. The spiral numbering module is used to spiral number all data points in the OVT surface element grid in the OVT surface element coordinate system, and generate a spiral code table based on the spiral numbered data points. The relative coordinate calculation module is used to calculate the relative coordinates of each seismic trace in the pre-acquired seismic trace set in the OVT surface element coordinate system; The spiral trace number matching module is used to match the spiral trace number corresponding to the relative coordinates of each seismic trace through the spiral encoding table, and write the spiral trace number corresponding to each seismic trace into a predefined trace header field; The spiral gather generation module is used to sort the seismic gathers according to the spiral gather number written in the trace header field to obtain spiral gathers.
9. A computer device, comprising a memory, a processor, and a computer program stored in the memory, characterized in that, The processor executes the computer program to implement the steps of the OVT spiral gather automatic sorting analysis method according to any one of claims 1 to 7.
10. A computer-readable storage medium having a computer program stored thereon, characterized in that, When executed by a processor, the computer program implements the steps of the OVT spiral gather automatic sorting analysis method according to any one of claims 1 to 7.