A data processing method and device based on a mirror coordinate system

Abstract

The specification provides a mirror coordinate system-based data processing method and device. A coordinate group of a target region to be detected in a seismic observation system is obtained. It is determined whether the coordinate system corresponding to the seismic observation system is a mirror coordinate system. In the case where it is determined that the coordinate system corresponding to the seismic observation system is a mirror coordinate system, it is determined whether the geological type of the target region is a tilted transverse anisotropic medium. In the case where the geological type of the target region is a tilted transverse anisotropic medium, the azimuth and the dip angle defined under a preset coordinate system are converted to the mirror coordinate system to obtain a target azimuth and a target dip angle. According to the target azimuth, the target dip angle, a velocity field and an anisotropy field under the mirror coordinate system, a travel time field of the target region is determined, the target region is subjected to imaging processing, and a target imaging result used for representing the geological structure of the target region is obtained. Compared with the method of reconstructing the coordinate system in the prior art, the overall calculation efficiency is improved, and the implementation process is simplified.

Description

Technical Field

[0001] This specification belongs to the field of seismic data technology for petroleum geophysical exploration, and in particular relates to a data processing method and apparatus based on a mirror coordinate system. Background Technology

[0002] Currently, existing seismic detection technologies typically employ the method of reconstructing the coordinate system to redefine the anisotropic field. This involves redefining the left-handed coordinate system according to the right-handed coordinate system for travel time field calculations. However, this approach suffers from cumbersome implementation procedures and low computational efficiency.

[0003] There is currently no effective solution to the above problems. Summary of the Invention

[0004] This specification provides a data processing method and apparatus based on a mirror coordinate system. By converting the coordinate system corresponding to the seismic observation system into a mirror coordinate system, and utilizing the target azimuth, target dip, velocity field, and anisotropic field corresponding to the mirror coordinate system, the travel time field of the target area can be quickly and accurately determined. Based on the travel time field, the target area can be imaged to obtain the target imaging results efficiently and accurately. Compared with methods that require reconstructing the coordinate system to calculate the travel time field, this method improves the overall computational efficiency and simplifies the implementation process.

[0005] This specification provides a data processing method based on a mirror coordinate system, including:

[0006] Obtain the coordinate set of the target area to be detected in the seismic observation system; wherein the coordinate set includes the coordinates of three points, which are composed of the line number, the trace number and the actual geographical coordinates in the seismic observation system, and the coordinates of the three points are not on the same straight line;

[0007] Based on the coordinate set of the target area, determine whether the coordinate system corresponding to the seismic observation system is a mirror coordinate system;

[0008] If the coordinate system corresponding to the seismic observation system is determined to be a mirror coordinate system, it is determined whether the geological type of the target area is a tilted, laterally anisotropic medium.

[0009] When the geological type of the target area is a tilted, laterally anisotropic medium, the azimuth and dip angle defined in the preset coordinate system are transformed to the mirror coordinate system to obtain the target azimuth and target dip angle.

[0010] Based on the target azimuth, target tilt, velocity field, and anisotropic field in the mirror coordinate system, the travel time field of the target region is determined, and the target region is imaged based on the travel time field to obtain target imaging results that characterize the geological structure of the target region.

[0011] In one embodiment, determining whether the coordinate system corresponding to the seismic observation system is a mirror coordinate system based on the coordinate set of the target area includes:

[0012] Based on the coordinate set of the target area, determine the component corresponding to the direction of increase of unit line number, and construct a first planar vector to characterize the direction of increase of line number based on the component corresponding to the direction of increase of unit line number.

[0013] Based on the coordinate set of the target area, determine the component corresponding to the direction of increase of the unit track number, and construct a second plane vector to characterize the direction of increase of the track number based on the component corresponding to the direction of increase of the unit track number.

[0014] Based on the first plane vector and the second plane vector, determine whether the coordinate system corresponding to the seismic observation system is a mirror coordinate system.

[0015] In one embodiment, determining whether the coordinate system corresponding to the seismic observation system is a mirror coordinate system based on the first plane vector and the second plane vector includes:

[0016] The normal vector of the third plane is obtained by taking the cross product of the normal vector of the first plane vector and the normal vector of the second plane vector.

[0017] Based on the third plane normal vector and the preset fixed value, determine whether the coordinate system corresponding to the earthquake observation system is a mirror coordinate system.

[0018] In one embodiment, the preset coordinate system is a geodetic coordinate system. Correspondingly, the step of transforming the azimuth and inclination angles defined in the preset coordinate system to the mirror coordinate system to obtain the target azimuth and target inclination angle includes:

[0019] Based on the first planar vector, determine the target angle between the unit line number increasing direction and the preset axis of the preset coordinate system;

[0020] Based on the target angle, the azimuth and tilt angles defined in the preset coordinate system are transformed to the mirror coordinate system to obtain the target azimuth and the target tilt angle.

[0021] In one embodiment, determining the target angle between the unit line number increment direction and a preset axis of the preset coordinate system based on the first planar vector includes:

[0022] When dyI≥0, the target angle between the direction of increasing unit line number and the preset axis of the preset coordinate system is determined according to the following formula:

[0023]

[0024] When dyI < 0, the target angle between the direction of increasing unit line number and the preset axis of the preset coordinate system is determined according to the following formula:

[0025]

[0026] Where α is the target angle, dxI is the projection of the component corresponding to the direction of increase of the unit line number in the north direction, and dyI is the projection of the component corresponding to the direction of increase of the unit line number in the east direction.

[0027] In one embodiment, the preset coordinate system is a coordinate system constructed based on directional projection angles. Correspondingly, the step of transforming the azimuth and inclination angles defined in the preset coordinate system to the mirror coordinate system to obtain the target azimuth and target inclination angles includes:

[0028] The target azimuth and inclination angles can be obtained using the following formulas:

[0029]

[0030] Where, φ ′ Let θ be the target azimuth angle. ′ Let φ be the target tilt angle, φ be the azimuth angle defined in the preset coordinate system, and θ be the tilt angle defined in the preset coordinate system.

[0031] In one embodiment, the step of imaging the target region based on the travel time field to obtain a target imaging result for characterizing the geological structure of the target region includes:

[0032] Based on the travel time field and the seismic observation data of the target area, the target area is imaged to obtain target imaging results that characterize the geological structure of the target area; wherein, the seismic observation data includes, but is not limited to, seismic signals collected at observation points located in the target area.

[0033] This specification provides a data processing device based on a mirror coordinate system, including:

[0034] The coordinate acquisition module is used to acquire the coordinate set of the target area to be detected in the seismic observation system; wherein, the coordinate set includes the coordinates of three points, which are composed of the line number, the trace number and the actual geographical coordinates in the seismic observation system, and the coordinates of the three points are not on the same straight line;

[0035] The coordinate system determination module is used to determine whether the coordinate system corresponding to the seismic observation system is a mirror coordinate system based on the coordinate group of the target area.

[0036] The medium type determination module is used to determine whether the geological type of the target area is a tilted transverse anisotropic medium when the coordinate system corresponding to the seismic observation system is determined to be a mirror coordinate system.

[0037] An angle determination module is used to transform the azimuth and dip angle defined in the preset coordinate system to the mirror coordinate system when the geological type of the target area is a tilted transverse anisotropic medium, so as to obtain the target azimuth and target dip angle.

[0038] The imaging data determination module is used to determine the travel time field of the target area based on the target azimuth, target tilt, velocity field and anisotropic field in the mirror coordinate system, and to perform imaging processing on the target area based on the travel time field to obtain target imaging results for characterizing the geological structure of the target area.

[0039] This specification also provides an electronic device, including a processor and a memory for storing processor-executable instructions, wherein the processor, when executing the instructions, implements a data processing method based on a mirror coordinate system.

[0040] This specification also provides a computer-readable storage medium having computer instructions stored thereon, which, when executed, implement a data processing method based on a mirror coordinate system.

[0041] Based on the data processing method based on a mirror coordinate system provided in this specification, the coordinate set of the target area to be detected in the seismic observation system is obtained. The coordinate set includes three coordinates, which are composed of the line number, trace number, and actual geographical location coordinates in the seismic observation system, and the three coordinates are not on the same straight line. Based on the coordinate set of the target area, it is determined whether the coordinate system corresponding to the seismic observation system is a mirror coordinate system. If the coordinate system corresponding to the seismic observation system is determined to be a mirror coordinate system, it is determined whether the geological type of the target area is a tilted laterally anisotropic medium. If the geological type of the target area is a tilted laterally anisotropic medium, the azimuth and dip angle defined under the preset coordinate system are transformed to the mirror coordinate system to obtain the target azimuth and target dip angle. Based on the target azimuth, target dip angle, velocity field, and anisotropic field in the mirror coordinate system, the travel time field of the target area is determined, and based on the travel time field, the target area is imaged to obtain a target imaging result characterizing the geological structure of the target area. In this way, by converting the coordinate system corresponding to the seismic observation system into a mirror coordinate system, and using the target azimuth, target dip, velocity field, and anisotropic field corresponding to the mirror coordinate system, the travel time field of the target area can be quickly and accurately determined. Based on the travel time field, the target area can be imaged to obtain the target imaging results of the target area efficiently and accurately. Compared with the method of calculating the travel time field by reconstructing the coordinate system, this method improves the overall calculation efficiency and simplifies the implementation process. Attached Figure Description

[0042] To more clearly illustrate the embodiments of this specification, the accompanying drawings used in the embodiments will be briefly introduced below. The drawings described below are only some embodiments recorded in this specification. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0043] Figure 1 This is a flowchart illustrating a data processing method based on a mirror coordinate system, provided in one embodiment of this specification.

[0044] Figure 2 This is a schematic diagram of the electronic device structure provided in one embodiment of this specification;

[0045] Figure 3 This is a schematic diagram of the structural composition of a data processing device based on a mirror coordinate system, provided in one embodiment of this specification.

[0046] Figure 4 This is a schematic diagram of a mirrored coordinate system provided in one embodiment of this specification. Detailed Implementation

[0047] To enable those skilled in the art to better understand the technical solutions in this specification, the technical solutions in the embodiments of this specification will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of this specification, and not all embodiments. Based on the embodiments in this specification, all other embodiments obtained by those skilled in the art without creative effort should fall within the scope of protection of this specification.

[0048] Seismic exploration is a method that uses seismic waves to detect underground structures and material distribution. It is widely used in the exploration and development of oil, natural gas, and mineral resources, as well as the prevention and control of geological disasters. During seismic exploration, precise calculations of parameters such as the propagation path and velocity of seismic waves are required to obtain accurate information about the underground structure.

[0049] However, the complexity of the Earth's medium, such as the presence of anisotropic media, complicates seismic wave propagation, increasing the difficulty of seismic detection. Geophysical exploration processing technology utilizes mathematical and physical methods to process and interpret geophysical exploration data. Current seismic detection techniques typically employ a coordinate system reconstruction method to redefine the anisotropic field. This involves redefining the left-handed coordinate system according to a right-handed coordinate system, then calculating the corresponding anisotropic field in the right-handed coordinate system. Furthermore, the observation system grid and anisotropic field grid in the left-handed coordinate system need to be transformed before the offset, followed by the offset in the new coordinate system. After the offset, another transformation is required to return to the initial left-handed coordinate system. This necessitates reordering the three dimensions of the three-dimensional anisotropic field in memory, resulting in low computational efficiency and a relatively cumbersome implementation process.

[0050] To address the root cause of the aforementioned problems, this manual converts the coordinate system corresponding to the seismic observation system into a mirror coordinate system. By utilizing the target azimuth, target dip, velocity field, and anisotropic field corresponding to the mirror coordinate system, the travel time field of the target area can be quickly and accurately determined. Based on the travel time field, the target area can be imaged, resulting in efficient and accurate target imaging. Compared to methods that require reconstructing the coordinate system to calculate the travel time field, this method improves overall computational efficiency and simplifies the implementation process.

[0051] See Figure 1 As shown in the embodiments of this specification, a data processing method based on a mirror coordinate system is provided, wherein the method is specifically applied to the server side. In specific implementation, the method may include the following:

[0052] S101: Obtain the coordinate set of the target area to be detected in the seismic observation system; wherein, the coordinate set includes three point coordinates, the three point coordinates are composed of the line number, the trace number and the actual geographical location coordinates in the seismic observation system, and the three point coordinates are not on the same straight line;

[0053] S102: Based on the coordinate set of the target area, determine whether the coordinate system corresponding to the seismic observation system is a mirror coordinate system;

[0054] S103: If the coordinate system corresponding to the seismic observation system is determined to be a mirror coordinate system, determine whether the geological type of the target area is a tilted laterally anisotropic medium;

[0055] S104: When the geological type of the target area is a tilted transverse anisotropic medium, the azimuth and dip angle defined under the preset coordinate system are transformed to the mirror coordinate system to obtain the target azimuth and target dip angle.

[0056] S105: Based on the target azimuth, target tilt, velocity field and anisotropic field in the mirror coordinate system, determine the travel time field of the target region, and perform imaging processing on the target region based on the travel time field to obtain target imaging results for characterizing the geological structure of the target region.

[0057] The aforementioned earthquake observation system can refer to a collection of overall equipment layout and working methods used to acquire, record, transmit and analyze seismic wave data. It can form a spatial layout by setting up source points and receiving points (sensor points), and combined with seismic instruments and data processing technology, it can be used to study underground geological structures, detect resources or monitor seismic activity.

[0058] The aforementioned coordinate set includes three coordinate points, which can be used to determine whether the seismic observation system is a left-handed or right-handed system. These three coordinate points are represented by a combination of line and track numbers (cmp) in the seismic observation system and their corresponding geographical coordinates (east coordinate x, north coordinate y). The line number can represent a number along a certain direction in the seismic observation system (usually longitudinal). The track number (cmp) can represent a number along another direction in the seismic observation system (usually transverse).

[0059] For example, the first point (line1, cmp1) has coordinates (x1, y1); the second point (line2, cmp2) has coordinates (x2, y2); and the third point (line3, cmp3) has coordinates (x3, y3). Furthermore, the coordinates of these three points are not on the same straight line; that is, their line numbers and track numbers cannot be located on the same straight line in the coordinate system.

[0060] The aforementioned Tilted Transverse Isotropic Medium (TTI) can be a model of elastic media used to describe the propagation characteristics of seismic waves in complex underground geological conditions. Unlike the Vertically Transverse Isotropy (VTI), which has a tilted axis of symmetry, the TTI medium exhibits isotropy within a plane, but anisotropy in directions perpendicular to that plane. Its axis of symmetry can be defined by its dip angle and azimuth angle.

[0061] The azimuth angle mentioned above refers to the angle between the target point and the north direction on the horizontal plane, usually measured clockwise, ranging from 0° to 360°. The dip angle is the angle between the target point and the observation point in the vertical direction, describing the altitude variation of the target point, ranging from 0° to 90°. The velocity field mentioned above describes the spatial distribution of seismic wave propagation velocity in the target area. The anisotropy field mentioned above can be used to describe the direction-dependent variations of rock physical properties (such as elastic modulus or wave velocity) in the target area.

[0062] The aforementioned travel time field can be obtained by calculating the travel time of seismic waves. It is used to correct for the offset of seismic waves caused by complex geological structures during conventional seismic data processing, reflecting the true spatial location of subsurface reflecting interfaces. Travel time refers to the time required for a seismic wave to propagate from its source (emission point) and receiver point to all locations within the subsurface medium. It reflects the length of time the seismic wave travels through the subsurface medium and is influenced by geological structure, medium type, and wave velocity distribution. The target imaging results are generated based on the travel time field and seismic observation data, used to visually present images of subsurface geological structures, displaying stratigraphic morphology, faults, folds, and other features. These are important data results in seismic exploration, used to guide resource exploration and geological research.

[0063] In some embodiments, when the coordinate system corresponding to the seismic observation system is determined to be a right-handed coordinate system, it is determined whether the geological type of the target area is a tilted laterally anisotropic medium.

[0064] When the geological type of the target area is a tilted, laterally anisotropic medium, the azimuth and dip angles defined in the preset coordinate system are transformed to the right-hand coordinate system to obtain the target azimuth and target dip angles.

[0065] In some embodiments, when the geological type of the target area is a tilted transversely anisotropic medium, the azimuth and dip angle defined in the preset coordinate system are transformed to the mirror coordinate system to obtain the target azimuth and target dip angle.

[0066] Based on the above embodiments, by transforming the azimuth and tilt angles in the preset coordinate system to the target azimuth, target tilt angle, velocity field, and anisotropic field in the mirror coordinate system, it is possible not only to effectively match the symmetry axis direction of the tilted transverse anisotropic medium, but also to simplify the processing flow of azimuth and tilt angles, thereby improving the accuracy during travel.

[0067] In some embodiments, when the geological type of the target area is a laterally isotropic medium, no transformation is required, and the target azimuth and dip angles in the mirror coordinate system are consistent with those in the right-hand coordinate system.

[0068] In some embodiments, the step of imaging the target area based on the travel time field to obtain a target imaging result characterizing the geological structure of the target area may further include:

[0069] Using a preset imaging model, the target area is imaged according to the travel time field to obtain target imaging results that characterize the geological structure of the target area; wherein, the preset imaging model is a model constructed by a fundamental preset machine learning algorithm.

[0070] Based on the above embodiments, by converting the coordinate system corresponding to the seismic observation system into a mirror coordinate system, and utilizing the target azimuth, target dip, velocity field, and anisotropic field corresponding to the mirror coordinate system, the travel time field of the target area can be quickly and accurately determined. Based on the travel time field, the target area can be imaged to obtain the target imaging results of the target area efficiently and accurately. Compared with the method of calculating the travel time field by reconstructing the coordinate system, this method improves the overall calculation efficiency and simplifies the implementation process.

[0071] In some embodiments, the method for determining whether the coordinate system corresponding to the seismic observation system is a mirror coordinate system based on the coordinate group of the target area may further include the following:

[0072] S1: Based on the coordinate set of the target area, determine the component corresponding to the direction of increase of the unit line number, and construct a first planar vector to characterize the direction of increase of the line number based on the component corresponding to the direction of increase of the unit line number.

[0073] S2: Based on the coordinate set of the target area, determine the component corresponding to the direction of increase of the unit track number, and construct a second plane vector to characterize the direction of increase of the track number based on the component corresponding to the direction of increase of the unit track number.

[0074] S3: Based on the first plane vector and the second plane vector, determine whether the coordinate system corresponding to the seismic observation system is a mirror coordinate system.

[0075] The components corresponding to the direction of increase in unit line number can be the components of the direction of increase in line number in the east-west (X-axis) and north-south (Y-axis) directions, used to describe the two-dimensional coordinate changes corresponding to the change in unit line number. For example, assuming that an increment in unit line number corresponds to an eastward coordinate increment of dxI and a northward coordinate increment of dyI, then the components in the direction of increase in unit line number can be represented by the vector (dxI, dyI).

[0076] The components corresponding to the direction of increase in the unit track number can be the components of the direction of increase in the track number in the east-west (X-axis) and north-south (Y-axis) directions, used to describe the two-dimensional coordinate changes corresponding to the change in the unit track number. For example, if the increment of one unit track number corresponds to an increment of dxX in the east and an increment of dyX in the north, then the components of the direction of increase in the unit track number can be represented by the vector (dxX, dyX).

[0077] Based on the above embodiments, according to the coordinate group of the target area, the components of the unit line number increase direction and the unit trace number increase direction are determined in sequence, and the first plane vector and the second plane vector are constructed. Then, it is determined whether the coordinate system of the seismic observation system is a mirror coordinate system. This can effectively solve the problem of inconsistent or misjudged coordinate system definitions in the seismic observation system and improve the accuracy of seismic data processing.

[0078] In some embodiments, the method for determining whether the coordinate system corresponding to the seismic observation system is a mirror coordinate system based on the first plane vector and the second plane vector may further include the following:

[0079] S1: Take the cross product of the normal vector of the first plane vector and the normal vector of the second plane vector to obtain the normal vector of the third plane.

[0080] S2: Based on the third plane normal vector and the preset fixed value, determine whether the coordinate system corresponding to the seismic observation system is a mirror coordinate system.

[0081] Specifically, for example, suppose the target area has the following three points: Point A corresponds to line number and track number (1,1), and its coordinates are (x0, y0). Point B corresponds to line number and track number (2,1), and its coordinates are (x1, y1). Point C corresponds to line number and track number (1,2), and its coordinates are (x2, y2).

[0082] The direction of increase for a unit line number refers to the direction in which the line number increases from 1 to 2, i.e., from point A to point B; its east-west component is dxI = x1 - x0, and its north-south component is dyI = y1 - y0. Therefore, the vector corresponding to the direction of increase for a unit line number is (dxI, dyI). The direction of increase for a unit track number refers to the direction in which the track number increases from 1 to 2, i.e., from point A to point C; its east-west component is dxX = x2 - x0, and its north-south component is dyX = y2 - y0. Therefore, the vector corresponding to the direction of increase for a unit track number is (dxX, dyX).

[0083] Construct a first plane vector (dxI, dyI, 0) and a second plane vector (dxX, dyX, 0) based on the direction of increase in unit line number and the direction of increase in unit track number, respectively. Calculate the cross product of the two vectors to obtain the third plane normal vector v3 = (0, 0, dxI·dyX - dxX·dyI). Determine whether the coordinate system is right-handed (positive normal vector) or left-handed (negative normal vector) based on the sign of the third plane normal vector. If it is left-handed, it can be considered as a mirror image of the left-handed coordinate system about the z = 0 plane.

[0084] In some embodiments, determining whether the coordinate system corresponding to the seismic observation system is a mirror coordinate system based on the third plane normal vector and a preset fixed value may further include:

[0085] Taking the aforementioned preset fixed value of 0 as an example, assuming the direction of line number increase in the target area is (1,0,0) and the direction of track number increase is (0,1,0), the plane normal vector formed by the two is calculated as (0,0,1). When compared with the preset fixed value of 0, the third component of the normal vector is 1, which is greater than 0, so the coordinate system is determined to be a conventional coordinate system; conversely, if the direction of line number increase is (1,0,0) and the direction of track number increase is (0,-1,0), its normal vector is (0,0,-1), and the third component is less than 0, then the coordinate system is a mirrored coordinate system.

[0086] Based on the above embodiments, by determining whether the coordinate system corresponding to the seismic observation system is a mirror coordinate system according to the relationship between the normal vector of the third plane and a preset fixed value, the efficiency and accuracy of coordinate system type identification can be greatly improved. This method simplifies complex geometric analysis into simple symbol comparison, facilitating automated judgment.

[0087] In some embodiments, the preset coordinate system is a geodetic coordinate system. Accordingly, the method of transforming the azimuth and inclination angles defined in the preset coordinate system to the mirror coordinate system to obtain the target azimuth and target inclination angle may further include the following in a specific implementation:

[0088] S1: Determine the target angle between the unit line number increasing direction and the preset axis of the preset coordinate system based on the first plane vector;

[0089] S2: Based on the target angle, the azimuth and tilt angle defined in the preset coordinate system are transformed to the mirror coordinate system to obtain the target azimuth and the target tilt angle.

[0090] In some embodiments, the method for determining the target angle between the unit line number increment direction and the preset axis of the preset coordinate system based on the first planar vector may further include the following:

[0091] When dyI≥0, the target angle between the direction of increasing unit line number and the preset axis of the preset coordinate system is determined according to the following formula:

[0092]

[0093] When dyI < 0, the target angle between the direction of increasing unit line number and the preset axis of the preset coordinate system is determined according to the following formula:

[0094]

[0095] Where α is the target angle, dxI is the projection of the component corresponding to the direction of increase of the unit line number in the north direction, and dyI is the projection of the component corresponding to the direction of increase of the unit line number in the east direction.

[0096] Specifically, the above coordinate system uses north as the positive x-axis, east as the positive y-axis, and depth as the positive z-axis. This coordinate system forms a right-handed coordinate system. The azimuth angle is defined as the angle of the projection direction of the xy plane in this coordinate system, that is, the angle relative to the x-axis. The angle range is 0°-360°. The angle between the axis of symmetry and the z-axis in the coordinate system is defined as the inclination angle θ. The angle range is 0°-90°. Assume that the angle of the direction of increasing line number relative to the x-axis is α.

[0097] Furthermore, the target azimuth and inclination angles in the mirror coordinate system can be obtained using the following transformation formulas:

[0098] φ ′ =360°+α-φ

[0099] θ ′ =θ

[0100] Where, φ ′ Let θ be the target azimuth angle. ′ The target tilt angle is given.

[0101] Furthermore, the target azimuth and inclination angles in the right-hand coordinate system can be obtained using the following conversion formulas:

[0102] φ ′ =φ-α

[0103] θ ′ =θ

[0104] Where, φ ′ Let θ be the target azimuth angle. ′ The target tilt angle is given.

[0105] Based on the above embodiments, by transforming the azimuth and dip angles defined in the preset coordinate system to the target azimuth and dip angles in the mirror coordinate system, the impact of coordinate system differences on seismic data processing can be effectively eliminated. This method ensures the consistency of angle information in different coordinate systems, improving the accuracy and reliability of seismic imaging.

[0106] In some embodiments, the preset coordinate system is a coordinate system constructed based on directional projection angles. Accordingly, the method of transforming the azimuth and inclination angles defined in the preset coordinate system to the mirror coordinate system to obtain the target azimuth and target inclination angles may further include the following in a specific implementation:

[0107] The target azimuth and inclination angles can be obtained using the following formulas:

[0108]

[0109] Where, φ ′ Let θ be the target azimuth angle. ′ Let φ be the target tilt angle, φ be the azimuth angle defined in the preset coordinate system, and θ be the tilt angle defined in the preset coordinate system.

[0110] The above conversion formulas all correspond to the symmetry axis of an inclined transversely anisotropic medium pointing towards the positive depth direction, i.e., when the inclination angle is 0°, the symmetry axis direction is the positive z-axis direction. When the symmetry axis direction points in the opposite direction to the depth direction, i.e., when the inclination angle is 0°, the symmetry axis is the negative z-axis direction, the azimuth angle is converted into the target azimuth angle φ pointing towards the positive depth direction using the following formula. ′ :

[0111] φ ′ =φ+180°

[0112] Then, the range of the target tilt angle is adjusted, and the angle outside 0°-360° is adjusted to 0°-360°, finally obtaining the target tilt angle and target azimuth angle in the mirror right-hand coordinate system corresponding to the left-hand observation system.

[0113] In some embodiments, the method of imaging the target area based on the travel time field to obtain target imaging results for characterizing the geological structure of the target area may further include the following:

[0114] Based on the travel time field and the seismic observation data of the target area, the target area is imaged to obtain target imaging results that characterize the geological structure of the target area; wherein, the seismic observation data includes, but is not limited to, seismic signals collected at observation points located in the target area.

[0115] In some embodiments, determining the travel time field corresponding to the target region based on the target azimuth, target tilt, velocity field, and anisotropic field in the mirror coordinate system may further include:

[0116] Using a preset wavefield propagation model, combined with the velocity distribution and path relationship of seismic wave propagation, and based on a preset ray tracing method, the propagation time of the seismic wave from the source to the target area is calculated, and the propagation time is used as the travel time field corresponding to the target area.

[0117] For example, given that the target azimuth of a point in the target area is 60° (representing the horizontal angle between the epicenter and the point) and the dip angle is 30° (representing the angle between the point and the horizontal plane), and combining the velocity distribution and path relationship of seismic wave propagation, the travel time from the epicenter to the point is calculated to be 3.5 seconds through ray tracing.

[0118] In some embodiments, the offset imaging result of the target region can also be determined based on the travel time field. Specifically, this may include:

[0119] Using a preset ray migration algorithm, the recorded seismic data is processed to correct the migration phenomenon of the seismic waves caused by the complexity of the underground structure, determine the true location of the underground interface, and generate migration imaging results to characterize the geological structure of the target area.

[0120] For example, based on the travel time, a preset ray-based migration algorithm is used to reposition the recorded seismic signal to the propagation path of the line connecting the source and the point, correcting the unrealistic reflections of seismic waves caused by tilted interfaces or velocity changes, and generating migration imaging results to characterize the geological structure of the target area.

[0121] Based on the above embodiments, by accurately calculating the travel time of the target area according to the target azimuth and dip angles in the mirror coordinate system, and determining the migration imaging results based on this, the accuracy of seismic data processing can be effectively improved. This method can correct deviations in the propagation paths of seismic waves in complex geological structures, restore the true location of subsurface interfaces, and significantly enhance the resolution and accuracy of seismic imaging. This not only helps to accurately describe the geological structure of the target area, but also provides more scientific data support for oil and gas resource exploration and subsurface structure analysis, thereby improving exploration efficiency and reducing exploration risks.

[0122] As can be seen from the above, the data processing method based on a mirror coordinate system provided in this specification obtains the coordinate set of the target area to be detected in the seismic observation system. The coordinate set includes three coordinates, which are composed of the line number, trace number, and actual geographical location coordinates in the seismic observation system, and the three coordinates are not on the same straight line. Based on the coordinate set of the target area, it is determined whether the coordinate system corresponding to the seismic observation system is a mirror coordinate system. If the coordinate system corresponding to the seismic observation system is determined to be a mirror coordinate system, it is determined whether the geological type of the target area is a tilted laterally anisotropic medium. If the geological type of the target area is a tilted laterally anisotropic medium, the azimuth and dip angle defined under the preset coordinate system are transformed to the mirror coordinate system to obtain the target azimuth and target dip angle. Based on the target azimuth, target dip angle, velocity field, and anisotropic field in the mirror coordinate system, the travel time field of the target area is determined, and based on the travel time field, the target area is imaged to obtain a target imaging result characterizing the geological structure of the target area. In this way, by converting the coordinate system corresponding to the seismic observation system into a mirror coordinate system, and using the target azimuth, target dip, velocity field, and anisotropic field corresponding to the mirror coordinate system, the travel time field of the target area can be quickly and accurately determined. Based on the travel time field, the target area can be imaged to obtain the target imaging results of the target area efficiently and accurately. Compared with the method of calculating the travel time field by reconstructing the coordinate system, this method improves the overall calculation efficiency and simplifies the implementation process.

[0123] See Figure 2 As shown in the embodiments of this specification, a specific electronic device is also provided, wherein the electronic device includes a network communication port 201, a processor 202 and a memory 203, and the above structures are connected by internal cables so that the various structures can perform specific data interaction.

[0124] Specifically, the network communication port 201 can be used to acquire the coordinate set of the target area to be detected in the seismic observation system; wherein the coordinate set includes three-point coordinates, which are composed of the line number, the track number and the actual geographical location coordinates in the seismic observation system, and the three-point coordinates are not on the same straight line.

[0125] The processor 202 can be specifically used to determine whether the coordinate system corresponding to the seismic observation system is a mirror coordinate system based on the coordinate group of the target area; if the coordinate system corresponding to the seismic observation system is determined to be a mirror coordinate system, it can determine whether the geological type of the target area is a tilted laterally anisotropic medium; if the geological type of the target area is a tilted laterally anisotropic medium, it can transform the azimuth and dip angle defined under the preset coordinate system to the mirror coordinate system to obtain the target azimuth and target dip angle; based on the target azimuth, target dip angle, velocity field and anisotropic field under the mirror coordinate system, it can determine the travel time field of the target area, and perform imaging processing on the target area based on the travel time field to obtain target imaging results used to characterize the geological structure of the target area.

[0126] The memory 203 can be used to store the corresponding instruction program.

[0127] Based on the above method, the relevant structural performance of electronic devices can be effectively utilized to improve the data processing speed of electronic devices and efficiently realize the data processing method based on the mirror coordinate system.

[0128] In this embodiment, the network communication port 201 can be a virtual port bound to different communication protocols, thereby enabling the sending or receiving of different data. For example, the network communication port can be a port responsible for web data communication, a port responsible for FTP data communication, or a port responsible for email data communication. Furthermore, the network communication port can also be a physical communication interface or communication chip. For example, it can be a wireless mobile network communication chip, such as GSM or CDMA; it can also be a Wi-Fi chip; or it can be a Bluetooth chip.

[0129] In this embodiment, the processor 202 can be implemented in any suitable manner. For example, the processor can take the form of a microprocessor or processor and a computer-readable medium storing computer-readable program code (e.g., software or firmware) executable by the (micro)processor, logic gates, switches, application-specific integrated circuits (ASICs), programmable logic controllers, and embedded microcontrollers, etc. This specification is not limiting.

[0130] In this embodiment, the memory 203 may include multiple layers. In a digital system, anything that can store binary data can be a memory. In an integrated circuit, a circuit with storage function but no physical form is also called a memory, such as RAM, FIFO, etc. In a system, a storage device with a physical form is also called a memory, such as a memory stick, TF card, etc.

[0131] This specification also provides a computer-readable storage medium based on the above-described data processing method using a mirror coordinate system. The method acquires a coordinate set of a target area to be detected within a seismic observation system. The coordinate set includes three coordinates, each composed of a line number, a trace number, and the actual geographical location coordinates within the seismic observation system, and these three coordinates are not on the same straight line. Based on the coordinate set of the target area, it is determined whether the coordinate system corresponding to the seismic observation system is a mirror coordinate system. If the coordinate system is determined to be a mirror coordinate system, it is determined whether the geological type of the target area is a tilted laterally anisotropic medium. If the geological type of the target area is a tilted laterally anisotropic medium, the azimuth and dip angles defined under a preset coordinate system are transformed to the mirror coordinate system to obtain the target azimuth and target dip angle. Based on the target azimuth, target dip angle, velocity field, and anisotropic field in the mirror coordinate system, the travel time field of the target area is determined. Based on the travel time field, the target area is imaged to obtain a target imaging result characterizing the geological structure of the target area.

[0132] In this embodiment, the storage medium includes, but is not limited to, Random Access Memory (RAM), Read-Only Memory (ROM), cache, hard disk drive (HDD), or memory card. The memory can be used to store computer program instructions. The network communication unit can be an interface configured according to standards specified in the communication protocol for network connection communication.

[0133] In this embodiment, the specific functions and effects implemented by the program instructions stored in the computer-readable storage medium can be explained in comparison with other embodiments, and will not be repeated here.

[0134] See Figure 3 At the software level, embodiments of this specification also provide a data processing device based on a mirror coordinate system, which may specifically include the following structural modules:

[0135] The coordinate acquisition module 301 is used to acquire the coordinate set of the target area to be detected in the seismic observation system; wherein, the coordinate set includes three point coordinates, the three point coordinates are composed of the line number, the trace number and the actual geographical location coordinates in the seismic observation system, and the three point coordinates are not on the same straight line;

[0136] The coordinate system determination module 302 is used to determine whether the coordinate system corresponding to the seismic observation system is a mirror coordinate system based on the coordinate group of the target area.

[0137] The medium type determination module 303 is used to determine whether the geological type of the target area is a tilted transverse anisotropic medium when the coordinate system corresponding to the seismic observation system is determined to be a mirror coordinate system.

[0138] Angle determination module 304 is used to transform the azimuth and dip angle defined in the preset coordinate system to the mirror coordinate system when the geological type of the target area is a tilted transverse anisotropic medium, so as to obtain the target azimuth and target dip angle.

[0139] The imaging data determination module 305 is used to determine the travel time field of the target area based on the target azimuth, target tilt, velocity field and anisotropic field in the mirror coordinate system, and to perform imaging processing on the target area based on the travel time field to obtain target imaging results for characterizing the geological structure of the target area.

[0140] In some embodiments, the coordinate system determination module 302, in specific implementation, determines the component corresponding to the increasing direction of the unit line number based on the coordinate group of the target area, and constructs a first plane vector representing the increasing direction of the line number based on the component corresponding to the increasing direction of the unit line number; determines the component corresponding to the increasing direction of the unit trace number based on the coordinate group of the target area, and constructs a second plane vector representing the increasing direction of the trace number based on the component corresponding to the increasing direction of the unit trace number; the coordinate system judgment module is used to determine whether the coordinate system corresponding to the seismic observation system is a mirror coordinate system based on the first plane vector and the second plane vector.

[0141] In some embodiments, the coordinate system determination module, in its specific implementation, performs a cross product of the normal vector of the first plane vector and the normal vector of the second plane vector to obtain the normal vector of the third plane; and determines whether the coordinate system corresponding to the seismic observation system is a mirror coordinate system based on the normal vector of the third plane and a preset fixed value.

[0142] In some embodiments, the preset coordinate system is a geodetic coordinate system. Accordingly, the angle determination module 304 is specifically implemented by determining, based on the first plane vector, the target angle between the unit line number increasing direction and the preset axis of the preset coordinate system; and based on the target angle, converting the azimuth and inclination angles defined in the preset coordinate system to the mirror coordinate system to obtain the target azimuth and the target inclination angle.

[0143] In some embodiments, the angle determination module, in specific implementation, determines the target angle between the unit line number increasing direction and the preset axis of the preset coordinate system according to the following formula when dyI≥0:

[0144]

[0145] When dyI < 0, the target angle between the direction of increasing unit line number and the preset axis of the preset coordinate system is determined according to the following formula:

[0146]

[0147] Where α is the target angle, dxI is the projection of the component corresponding to the direction of increase of the unit line number in the north direction, and dyI is the projection of the component corresponding to the direction of increase of the unit line number in the east direction.

[0148] In some embodiments, the preset coordinate system is a coordinate system constructed based on the directional projection angle. Accordingly, the angle determination module 304, in its specific implementation, obtains the target azimuth and target tilt angle according to the following formula:

[0149]

[0150] Where, φ ′ Let θ be the target azimuth angle. ′ Let φ be the target tilt angle, φ be the azimuth angle defined in the preset coordinate system, and θ be the tilt angle defined in the preset coordinate system.

[0151] In some embodiments, the imaging data determination module 305, in specific implementation, performs imaging processing on the target area based on the travel time field and the seismic observation data of the target area to obtain target imaging results for characterizing the geological structure of the target area; wherein, the seismic observation data includes, but is not limited to, seismic signals collected at the observation points deployed in the target area.

[0152] It should be noted that the units, devices, or modules described in the above embodiments can be implemented by computer chips or physical entities, or by products with certain functions. For ease of description, the above devices are described by dividing them into various modules according to their functions. Of course, in implementing this specification, the functions of each module can be implemented in one or more software and / or hardware, or the module that implements the same function can be implemented by a combination of multiple sub-modules or sub-units, etc. The device embodiments described above are merely illustrative. For example, the division of units is only a logical functional division, and there may be other division methods in actual implementation. For example, multiple units or components can be combined or integrated into another system, or some features can be ignored or not executed. Furthermore, the coupling or direct coupling or communication connection between the devices or units shown or discussed can be through some interfaces, and the indirect coupling or communication connection between devices or units can be electrical, mechanical, or other forms.

[0153] As can be seen from the above, the data processing device based on the mirror coordinate system provided in the embodiments of this specification can quickly and accurately determine the travel time field of the target area by converting the coordinate system corresponding to the seismic observation system into the mirror coordinate system and using the target azimuth, target dip, velocity field and anisotropic field corresponding to the mirror coordinate system. Based on the travel time field, the target area can be imaged to obtain the target imaging result of the target area efficiently and accurately. Compared with the method of calculating the travel time field by reconstructing the coordinate system, this method improves the overall calculation efficiency and simplifies the implementation process.

[0154] In a specific scenario example, the data processing method and apparatus based on a mirror coordinate system provided in this specification can be applied. By converting the coordinate system corresponding to the seismic observation system into a mirror coordinate system, and utilizing the target azimuth, target dip, velocity field, and anisotropic field corresponding to the mirror coordinate system, the travel time field of the target area can be quickly and accurately determined. Based on the travel time field, imaging processing of the target area can be performed to obtain the target imaging results efficiently and accurately. Compared with methods that require reconstructing the coordinate system to calculate the travel time field, this method improves the overall computational efficiency and simplifies the implementation process. The specific implementation process may include the following:

[0155] S1: Acquire seismic data, migration velocity field, coordinates of three points on the seismic observation system, and migration parameters;

[0156] S2: Based on the coordinates of three points in the seismic observation system, A1(line1, cmp1), (x1, y1), A2(line2, cmp2), (x2, y2), and A3(line3, cmp4), (x3, y3), determine whether the seismic observation system is left-handed or right-handed. (linei, cmpi) represent the line and trace number of one point, and (xi, yi) represent the east and north coordinates of the line and trace number positions, respectively. The three points are not on the same straight line. Assume that the coordinates of the line and trace number (1,1) are (x0, y0), the direction of increase of the unit line number is (dxI, dyI), and the direction of increase of the unit trace number is (dxX, dyX). By solving the system of equations, the values ​​of the six unknowns (x0, y0, dxI, dyI, dxX, dyX) can be calculated:

[0157]

[0158] Solving for:

[0159]

[0160] x0=[x1-(line1-1)*dxI-(cmp1-1)*dxX]*0.5

[0161]

[0162] y0=[y1-(line1-1)*dyI-(cmp1-1)*dyX]*0.5

[0163] Based on the increasing direction of the unit line number and unit track number, two three-dimensional vectors can be constructed: the first plane vector v1 = (dxI, dyI, 0), and the second plane vector v2 = (dxX, dyX, 0), where v1 and v2 represent the increasing direction of the line number and track number, respectively. Then, the normal vectors of the v1 and v2 planes are calculated: the third plane vector v3 = v1 × v2 = (0, 0, dxI*dyX - dxX*dyI), where × represents the cross product. When the third component of v3 is greater than zero, it indicates that the third plane vector v3 is aligned with the positive z-direction. In this case, v1, v2, and the positive z-axis form a right-handed coordinate system. When the third component of v3 is less than zero, it indicates that vector v3 is aligned with the negative z-direction. In this case, v1, v2, and the positive z-axis form a left-handed coordinate system. When v1, v2, and the positive z-axis form a left-handed coordinate system, this coordinate system can be considered as a mirror image of the original left-handed coordinate system about the z=0 plane. (See [reference needed]). Figure 4 As shown, in the mirror coordinate system, the coordinates on the v1, v2, z axes and their order in memory remain unchanged, and the travel time calculated in the mirror coordinate system is the same as that in the original left-handed coordinate system.

[0164] While this specification provides the steps of operation for the methods described in the embodiments or flowcharts, more or fewer steps may be included based on conventional or non-inventive means. The order of steps listed in the embodiments is merely one possible order of execution among many steps and does not represent the only possible order. In actual device or client product execution, the methods shown in the embodiments or drawings may be executed sequentially or in parallel (e.g., in a parallel processor or multi-threaded processing environment, or even a distributed data processing environment). The terms "comprising," "including," or any other variations thereof are intended to cover a non-exclusive inclusion, such that a process, method, product, 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, product, or apparatus. Without further limitations, the presence of other identical or equivalent elements in a process, method, product, or apparatus that includes said elements is not excluded. The terms "first," "second," etc., are used to denote names and do not indicate any particular order.

[0165] Those skilled in the art will also know that, besides implementing the controller using purely computer-readable program code, the same functions can be achieved by logically programming the method steps, making the controller function as logic gates, switches, application-specific integrated circuits (ASICs), programmable logic controllers (PLCs), and embedded microcontrollers. Therefore, such a controller can be considered a hardware component, and the devices within it used to implement various functions can also be considered structures within that hardware component. Alternatively, the devices used to implement various functions can be considered as both software modules implementing the method and structures within a hardware component.

[0166] As can be seen from the above description of the embodiments, those skilled in the art can clearly understand that this specification can be implemented by means of software plus necessary general-purpose hardware platforms. Based on this understanding, the technical solutions of this specification can essentially be embodied in the form of a software product. This computer software product can be stored in a storage medium, such as ROM / RAM, magnetic disk, optical disk, etc., and includes several instructions to cause a computer device (which may be a personal computer, mobile terminal, server, or network device, etc.) to execute the methods described in the various embodiments or some parts of the embodiments of this specification.

[0167] Although this specification has been described by way of examples, those skilled in the art will recognize that many variations and modifications are possible without departing from the spirit of this specification, and it is intended that the appended claims cover such variations and modifications without departing from the spirit of this specification.

Claims

1. A data processing method based on a mirror coordinate system, characterized in that, include: Obtain the coordinate set of the target area to be detected in the seismic observation system; wherein the coordinate set includes the coordinates of three points, which are composed of the line number, the trace number and the actual geographical coordinates in the seismic observation system, and the coordinates of the three points are not on the same straight line; Based on the coordinate set of the target area, determine whether the coordinate system corresponding to the seismic observation system is a mirror coordinate system; If the coordinate system corresponding to the seismic observation system is determined to be a mirror coordinate system, it is determined whether the geological type of the target area is a tilted, laterally anisotropic medium. When the geological type of the target area is a tilted, laterally anisotropic medium, the azimuth and dip angle defined in the preset coordinate system are transformed to the mirror coordinate system to obtain the target azimuth and target dip angle. Based on the target azimuth, target tilt, velocity field, and anisotropic field in the mirror coordinate system, the travel time field of the target region is determined, and the target region is imaged based on the travel time field to obtain target imaging results that characterize the geological structure of the target region.

2. The method according to claim 1, characterized in that, The step of determining whether the coordinate system corresponding to the seismic observation system is a mirror coordinate system based on the coordinate set of the target area includes: Based on the coordinate set of the target area, determine the component corresponding to the direction of increase of unit line number, and construct a first planar vector to characterize the direction of increase of line number based on the component corresponding to the direction of increase of unit line number. Based on the coordinate set of the target area, determine the component corresponding to the direction of increase of the unit track number, and construct a second plane vector to characterize the direction of increase of the track number based on the component corresponding to the direction of increase of the unit track number. Based on the first plane vector and the second plane vector, determine whether the coordinate system corresponding to the seismic observation system is a mirror coordinate system.

3. The method according to claim 2, characterized in that, The step of determining whether the coordinate system corresponding to the seismic observation system is a mirror coordinate system based on the first plane vector and the second plane vector includes: The normal vector of the third plane is obtained by performing a cross product of the normal vector of the first plane vector and the normal vector of the second plane vector. Based on the third plane normal vector and the preset fixed value, determine whether the coordinate system corresponding to the earthquake observation system is a mirror coordinate system.

4. The method according to claim 3, characterized in that, The preset coordinate system is a geodetic coordinate system. Correspondingly, the step of transforming the azimuth and dip angles defined in the preset coordinate system to the mirror coordinate system to obtain the target azimuth and dip angles includes: Based on the first planar vector, determine the target angle between the unit line number increasing direction and the preset axis of the preset coordinate system; Based on the target angle, the azimuth and tilt angles defined in the preset coordinate system are transformed to the mirror coordinate system to obtain the target azimuth and the target tilt angle.

5. The method according to claim 4, characterized in that, Determining the target angle between the unit line number increment direction and the preset axis of the preset coordinate system based on the first planar vector includes: when At that time, the target angle between the direction of increase of the unit line number and the preset axis of the preset coordinate system is determined according to the following formula: when At that time, the target angle between the direction of increase of the unit line number and the preset axis of the preset coordinate system is determined according to the following formula: in, Let dxI be the target angle, dyI be the projection of the component corresponding to the direction of the unit line number in the north direction, and dyI be the projection of the component corresponding to the direction of the unit line number in the east direction.

6. The method according to claim 1, characterized in that, The preset coordinate system is a coordinate system constructed based on directional projection angles. Correspondingly, the step of transforming the azimuth and inclination angles defined in the preset coordinate system to the mirror coordinate system to obtain the target azimuth and target inclination angles includes: The target azimuth and inclination angles can be obtained using the following formulas: in, The target azimuth angle is... The target tilt angle, The azimuth angle defined in the preset coordinate system. The tilt angle is defined under the preset coordinate system.

7. The method according to claim 1, characterized in that, The step of imaging the target region based on the travel time field to obtain target imaging results for characterizing the geological structure of the target region includes: Based on the travel time field and the seismic observation data of the target area, the target area is imaged to obtain target imaging results that characterize the geological structure of the target area; wherein, the seismic observation data includes, but is not limited to, seismic signals collected at observation points located in the target area.

8. A data processing device based on a mirror coordinate system, characterized in that, include: The coordinate acquisition module is used to acquire the coordinate set of the target area to be detected in the seismic observation system; wherein, the coordinate set includes the coordinates of three points, which are composed of the line number, the trace number and the actual geographical coordinates in the seismic observation system, and the coordinates of the three points are not on the same straight line; The coordinate system determination module is used to determine whether the coordinate system corresponding to the seismic observation system is a mirror coordinate system based on the coordinate group of the target area. The medium type determination module is used to determine whether the geological type of the target area is a tilted transverse anisotropic medium when the coordinate system corresponding to the seismic observation system is determined to be a mirror coordinate system. An angle determination module is used to transform the azimuth and dip angle defined in the preset coordinate system to the mirror coordinate system when the geological type of the target area is a tilted transverse anisotropic medium, so as to obtain the target azimuth and target dip angle. The imaging data determination module is used to determine the travel time field of the target area based on the target azimuth, target tilt, velocity field and anisotropic field in the mirror coordinate system, and to perform imaging processing on the target area based on the travel time field to obtain target imaging results for characterizing the geological structure of the target area.

9. An electronic device, characterized in that, It includes a processor and a memory for storing processor-executable instructions, wherein the processor, when executing the instructions, implements the steps of the data processing method based on a mirror coordinate system as described in any one of claims 1 to 7.

10. A computer-readable storage medium, characterized in that, It stores computer instructions that, when executed by a processor, implement the steps of the data processing method based on a mirror coordinate system as described in any one of claims 1 to 7.

Images