Static correction method and device for three-dimensional seismic first arrival data

Abstract

The present application provides a static correction method and device for three-dimensional seismic first arrival data. The method includes: acquiring target three-dimensional seismic first arrival data; obtaining the first arrival travel time circle parameter radius and the first arrival travel time error of each of the plurality of receiving points corresponding to the target three-dimensional seismic first arrival data based on a preset roundness error evaluation method; and obtaining the anisotropic static correction amount of each receiving point according to the first arrival travel time circle parameter radius and the first arrival travel time error. The present application can improve the accuracy of static correction, especially the accuracy of anisotropic static correction in complex areas.

Description

Technical Field

[0001] This application relates to the technical field of petroleum seismic exploration, and particularly to a static correction method and device for three-dimensional seismic first arrival data. Background Art

[0002] With the in-depth development of petroleum seismic exploration, the static correction problem has become increasingly prominent, which is an important factor affecting the quality of near-surface anisotropic data.

[0003] At present, there are mainly two types of near-surface anisotropic static corrections: the first is the first arrival travel time medium hypothesis method. This technique is based on the assumption of the anisotropic direction of the underground medium (such as VTI, HTI, TTI, etc. media), deduces the expression of the phase velocity of the medium, conducts anisotropic ray tracing and establishes a multi-parameter tomography equation, completes the velocity and anisotropic parameter distribution of the assumed anisotropic direction, and finally completes the anisotropic static correction. Even by analyzing the tomography sensitivity kernel function of multiple parameters, it is difficult to decouple the velocity and anisotropic parameters and achieve accurate back-projection of the tomography inversion residual, resulting in the anisotropic parameter distribution often changing with the velocity distribution. The other is the first arrival and waveform joint method. This technique realizes the full elastic wave inversion technique, makes no new assumptions about the underground medium, restores the underground elastic parameter distribution through elastic wave forward modeling, converts the velocity and anisotropic parameter distribution, and completes the anisotropic static correction. Due to the great difficulty of accurate elastic wave forward modeling and the large consumption of computing resources in complex areas, this technique is difficult to be widely applied in complex areas. Summary of the Invention

[0004] In view of at least one problem in the prior art, this application proposes a static correction method and device for three-dimensional seismic first arrival data, which can improve the accuracy of static correction, especially the accuracy of anisotropic static correction in complex areas.

[0005] To solve the above technical problems, this application provides the following technical solutions:

[0006] In a first aspect, this application provides a static correction method for three-dimensional seismic first arrival data, including:

[0007] Obtain target three-dimensional seismic first arrival data;

[0008] Based on a preset roundness error evaluation method, obtain the first arrival travel time circle parameter radius and first arrival travel time error of each of the multiple receiving points corresponding to the target three-dimensional seismic first arrival data;

[0009] According to the first arrival travel time circle parameter radius, the first arrival travel time error, and the target three-dimensional seismic first arrival data, obtain the anisotropic static correction amount of each receiving point.

[0010] Further, before obtaining the radius of the first arrival travel time circle parameter and the first arrival travel time error of each of the multiple receiving points corresponding to the target three-dimensional seismic first arrival data based on the preset roundness error evaluation method, the following steps are also included:

[0011] Perform elevation static correction on the first arrival travel times in the target three-dimensional seismic first arrival data, and correct the first arrival travel times to the same reference plane.

[0012] Further, the target three-dimensional seismic first arrival data includes: offset, the receiving point coordinates corresponding to the offset, shot point coordinates, and first arrival travel times;

[0013] Correspondingly, obtaining the radius of the first arrival travel time circle parameter and the first arrival travel time error of each of the multiple receiving points corresponding to the target three-dimensional seismic first arrival data based on the preset roundness error evaluation method includes:

[0014] Divide the offsets of the common shot points into multiple groups;

[0015] Based on the preset roundness error evaluation method, first arrival travel times, shot point coordinates, and receiving point coordinates, obtain the radius of the first arrival travel time circle parameter and the first arrival travel time error corresponding to each group;

[0016] Take the radius of the first arrival travel time circle parameter and the first arrival travel time error corresponding to all groups as the radius of the first arrival travel time circle parameter and the first arrival travel time error of the multiple receiving points corresponding to the target three-dimensional seismic first arrival data.

[0017] Further, obtaining the anisotropic static correction amount of each receiving point according to the radius of the first arrival travel time circle parameter, the first arrival travel time error, and the target three-dimensional seismic first arrival data includes:

[0018] Obtain the first anisotropic velocity correction amount of each receiving point according to the radius of the first arrival travel time circle parameter, the first arrival travel time error, and the offset in the target three-dimensional seismic first arrival data;

[0019] Apply the first anisotropic velocity correction amount to determine the anisotropic static correction amount of each receiving point.

[0020] Further, after obtaining the first anisotropic velocity correction amount of each receiving point according to the radius of the first arrival travel time circle parameter, the first arrival travel time error, and the offset in the target three-dimensional seismic first arrival data, the following steps are also included:

[0021] Based on the preset constraint conditions, optimization algorithm, and the first anisotropic velocity correction amount, obtain the unique second anisotropic velocity correction amount of each receiving point;

[0022] Apply the second anisotropic velocity correction amount to determine the anisotropic static correction amount of each receiving point.

[0023] Further, after obtaining the first arrival travel time circle parameter radius and the first arrival travel time error of each of the multiple receiving points corresponding to the target three-dimensional seismic first arrival data based on the preset roundness error evaluation method, the method further includes:

[0024] Based on the first arrival travel time circle parameter radius, obtaining the isotropic static correction amount of each receiving point.

[0025] Further, the obtaining the isotropic static correction amount of each receiving point based on the first arrival travel time circle parameter radius includes:

[0026] Based on the first arrival travel time circle parameter radius, obtaining the first isotropic velocity of each receiving point;

[0027] Applying the first isotropic velocity to determine the isotropic static correction amount of each receiving point.

[0028] Further, after obtaining the first isotropic velocity of each receiving point based on the first arrival travel time circle parameter radius, the method further includes:

[0029] Based on a preset constraint condition, an optimization algorithm, and the first isotropic velocity, obtaining a unique second isotropic velocity of each receiving point;

[0030] Applying the second isotropic velocity to determine the isotropic static correction amount of each receiving point.

[0031] In a second aspect, the present application provides a static correction device for three-dimensional seismic first arrival data, including:

[0032] An acquisition module, configured to acquire target three-dimensional seismic first arrival data;

[0033] An error evaluation module, configured to obtain the first arrival travel time circle parameter radius and the first arrival travel time error of each of the multiple receiving points corresponding to the target three-dimensional seismic first arrival data based on a preset roundness error evaluation method;

[0034] An anisotropic static correction module, configured to obtain the anisotropic static correction amount of each receiving point according to the first arrival travel time circle parameter radius, the first arrival travel time error, and the target three-dimensional seismic first arrival data.

[0035] Further, the static correction device for three-dimensional seismic first arrival data further includes:

[0036] A correction module, configured to perform elevation static correction on the first arrival travel time in the target three-dimensional seismic first arrival data, and correct the first arrival travel time to the same reference plane.

[0037] Further, the target three-dimensional seismic first arrival data includes: offset, the receiving point coordinates corresponding to the offset, shot point coordinates, and first arrival travel time;

[0038] Correspondingly, the error evaluation module includes:

[0039] A grouping unit, configured to divide the offsets of common shot points into multiple groups;

[0040] An error evaluation unit, configured to obtain the radius of the first arrival travel time circle parameter and the first arrival travel time error corresponding to each group based on the preset roundness error evaluation method, the first arrival travel time, the shot point coordinates, and the receiving point coordinates;

[0041] A determination unit, configured to use the radius of the first arrival travel time circle parameter and the first arrival travel time error corresponding to all groups as the radius of the first arrival travel time circle parameter and the first arrival travel time error of multiple receiving points corresponding to the target three-dimensional seismic first arrival data.

[0042] Further, the anisotropic static correction module includes:

[0043] A velocity correction unit, configured to obtain the first anisotropic velocity correction amount of each receiving point according to the radius of the first arrival travel time circle parameter, the first arrival travel time error, and the offset in the target three-dimensional seismic first arrival data;

[0044] A first anisotropic static correction unit, configured to determine the anisotropic static correction amount of each receiving point by applying the first anisotropic velocity correction amount.

[0045] Further, the anisotropic static correction module further includes:

[0046] A first optimization unit, configured to obtain the unique second anisotropic velocity correction amount of each receiving point based on the preset constraint conditions, optimization algorithm, and the first anisotropic velocity correction amount;

[0047] A second anisotropic static correction unit, configured to determine the anisotropic static correction amount of each receiving point by applying the second anisotropic velocity correction amount.

[0048] Further, the static correction device for three-dimensional seismic first arrival data further includes:

[0049] An isotropic static correction module, configured to obtain the isotropic static correction amount of each receiving point based on the radius of the first arrival travel time circle parameter.

[0050] Further, the isotropic static correction module includes:

[0051] An isotropic velocity determination unit, configured to obtain the first isotropic velocity of each receiving point based on the radius of the first arrival travel time circle parameter;

[0052] A first isotropic static correction unit, configured to determine the isotropic static correction amount of each receiving point by applying the first isotropic velocity.

[0053] Further, the isotropic static correction module further includes:

[0054] A second optimization unit, configured to obtain a unique second isotropic velocity for each receiving point based on a preset constraint condition, an optimization algorithm, and the first isotropic velocity;

[0055] A second isotropic static correction unit, configured to determine the isotropic static correction amount for each receiving point by applying the second isotropic velocity.

[0056] In a third aspect, the present application provides an electronic device, including a memory, a processor, and a computer program stored on the memory and executable on the processor. When the processor executes the program, the static correction method for three-dimensional seismic first arrival data described above is implemented.

[0057] In a fourth aspect, the present application provides a computer-readable storage medium, on which computer instructions are stored. When the instructions are executed, the static correction method for three-dimensional seismic first arrival data described above is implemented.

[0058] As can be seen from the above technical solutions, the present application provides a static correction method and apparatus for three-dimensional seismic first arrival data. Among them, the method includes: obtaining target three-dimensional seismic first arrival data; based on a preset roundness error evaluation method, obtaining the first arrival travel time circle parameter radius and the first arrival travel time error of each of the plurality of receiving points corresponding to the target three-dimensional seismic first arrival data; according to the first arrival travel time circle parameter radius, the first arrival travel time error, and the target three-dimensional seismic first arrival data, obtaining the anisotropic static correction amount for each receiving point, which can improve the accuracy of static correction, especially can improve the accuracy of anisotropic static correction in complex areas. Specifically, it can avoid the assumption of layered distribution of surface layer velocity, especially in the case of severe anisotropy near the surface, can solve the static correction problem in complex areas, can restore the near-surface velocity distribution from seismic first arrival information to achieve high-precision tomographic static correction; fully considers the anisotropy that may exist in all directions, can effectively establish an accurate near-surface velocity model, can complete all-round static correction, can improve the static correction accuracy of three-dimensional seismic data, and thus can be beneficial to the subsequent high-resolution processing of seismic data. It is very practically significant in terms of improving the quality of static correction, shortening the static correction cycle, saving energy, reducing emissions, and reducing costs and increasing efficiency. Description of the Drawings

[0059] In order to more clearly illustrate the technical solutions in the embodiments of the present application or the prior art, the following will briefly introduce the drawings required for use in the description of the embodiments or the prior art. Obviously, the drawings in the following description are only some embodiments of the present application. For those of ordinary skill in the art, without creative efforts, other drawings can be obtained based on these drawings.

[0060] Figure 1 It is a schematic flowchart of the static correction method for three-dimensional seismic first arrival data in an embodiment of the present application;

[0061] Figure 2 It is a schematic flowchart of the static correction method for three-dimensional seismic first arrival data in another embodiment of the present application;

[0062] Figure 3 It is a schematic flowchart of steps 301 to 303 of the static correction method for three-dimensional seismic first arrival data in an embodiment of the present application;

[0063] Figure 4 It is a schematic flowchart of steps 401 and 402 of the static correction method for three-dimensional seismic first arrival data in an embodiment of the present application;

[0064] Figure 5 It is a schematic flowchart of the static correction method for three-dimensional seismic first arrival data in yet another embodiment of the present application;

[0065] Figure 6 It is a schematic flowchart of the static correction method for three-dimensional seismic first arrival data in an application example of the present application;

[0066] Figure 7 It is a schematic diagram of the three-dimensional distribution of first arrivals of a single shot in an example of the present application;

[0067] Figure 8 It is a schematic diagram of the contour distribution of first arrivals in a plane in an example of the present application;

[0068] Figure 9 It is a schematic diagram of the structure of the static correction device for three-dimensional seismic first arrival data in an embodiment of the present application;

[0069] Figure 10 It is a schematic diagram of the structure of the static correction device for three-dimensional seismic first arrival data in another embodiment of the present application;

[0070] Figure 11 It is a schematic diagram of the structure of the static correction device for three-dimensional seismic first arrival data in yet another embodiment of the present application;

[0071] Figure 12 It is a schematic block diagram of the system composition of the electronic device 9600 in an embodiment of the present application. Detailed implementation manners

[0072] To enable those skilled in the art to better understand the technical solutions in this specification, the technical solutions in the embodiments of this application will be clearly and completely described below in conjunction with the accompanying drawings in the embodiments of this application. Obviously, the described embodiments are only a part of the embodiments of this application, rather than all of the embodiments. Based on the embodiments in this application, all other embodiments obtained by those of ordinary skill in the art without creative efforts shall fall within the protection scope of this application.

[0073] Due to the large near-surface undulation and drastic lateral velocity variation in complex surface areas, static correction processing is often required for the collected seismic data. Conventional static correction methods include: elevation static correction, model-based static correction, refraction static correction, tomography static correction, and static correction methods based on first arrival time, etc.

[0074] To improve the accuracy of static correction, especially the accuracy of anisotropic static correction in complex areas, based on the multiple coverage of seismic data and the characteristic that the first arrival travel time is a comprehensive response of the subsurface physical property distribution, combined with the seismic tomography inversion static correction method in petroleum seismic exploration, this application can first correct the three-dimensional seismic first arrival common shotpoint data to a fixed reference surface through elevation static correction; then group the first arrival travel times of the common shotpoints according to the offset; secondly, complete the estimation of the best isotropic medium and the error distribution for each grouped data, and obtain the distribution of the best isotropic medium parameters and the error distribution at each point within the grouped offset; sequentially complete the estimation of the best isotropic medium and the error distribution for all grouped data within the common shotpoint; combine the estimation of the best isotropic medium and the error distribution of all grouped data within the common shotpoint to complete the estimation of the distribution of the best isotropic medium parameters and the error distribution for all points within the common shotpoint range; finally, combine the estimation of the distribution of the best isotropic medium parameters and the error distribution of all common shotpoints to complete the best estimation within the entire data range. Only the first arrival travel time information is used in the whole process, and no assumption about the anisotropic direction of the medium is required, which can realize the all-round anisotropic velocity modeling within the common shotpoint, and further can complete the all-round velocity modeling within the entire data range, and can achieve the purpose of completing all-round static correction.

[0075] Based on this, to improve the accuracy of static correction, especially the accuracy of anisotropic static correction in complex areas, the embodiments of this application provide a static correction device for three-dimensional seismic first arrival data. The device can be a server or a client device. The client device can include smart phones, tablet electronic devices, network set-top boxes, portable computers, desktop computers, personal digital assistants (PDAs), vehicle-mounted devices, and smart wearable devices, etc. Among them, the smart wearable devices can include smart glasses, smart watches, and smart bracelets, etc.

[0076] In practical applications, the part of performing static correction on three-dimensional seismic first-arrival data can be executed on the server side as described above, or all operations can be completed in the client device. Specifically, it can be selected according to the processing capacity of the client device and the limitations of the user usage scenario, etc. This application does not make any limitations in this regard. If all operations are completed in the client device, the client device may further include a processor.

[0077] The above-mentioned client device may have a communication module (i.e., communication unit), which can communicate with a remote server to realize data transmission with the server. The server may include a server on the task scheduling center side, and in other implementation scenarios, it may also include a server of an intermediate platform, such as a server of a third-party server platform having a communication link with the task scheduling center server. The server may include a single computer device, or may include a server cluster composed of multiple servers, or a server structure of a distributed device.

[0078] Any suitable network protocol can be used for communication between the server and the client device, including network protocols not yet developed on the filing date of this application. The network protocol may, for example, include TCP / IP protocol, UDP / IP protocol, HTTP protocol, HTTPS protocol, etc. Of course, the network protocol may, for example, also include RPC protocol (Remote Procedure Call Protocol) and REST protocol (Representational State Transfer) used on top of the above-mentioned protocols.

[0079] Specifically, it will be described through the following various embodiments.

[0080] In order to improve the accuracy of static correction, especially to improve the accuracy of anisotropic static correction in complex areas, this embodiment provides a method for static correction of three-dimensional seismic first-arrival data, the execution subject of which is a device for static correction of three-dimensional seismic first-arrival data. The device for static correction of three-dimensional seismic first-arrival data includes but is not limited to a server, as Figure 1 shown, the method specifically includes the following content:

[0081] Step 101: Obtain target three-dimensional seismic first-arrival data.

[0082] Specifically, the target three-dimensional seismic first-arrival data may include: first-arrival travel time and information related to the first-arrival travel time. The three-dimensional seismic first-arrival data may be prestack seismic data; prestack seismic data may refer to seismic data that is migrated before stacking. Prestack seismic data has high precision and is more practical for processing seismic data in complex areas.

[0083] Step 102: Based on a preset roundness error evaluation method, obtain the radius of the first arrival travel time circle parameter and the first arrival travel time error for each of the multiple receiving points corresponding to the target 3D seismic first arrival data.

[0084] Specifically, in this embodiment, the roundness error can be used to determine the error between the theoretical first arrival travel time circle and the actual first arrival; the preset roundness error evaluation method can be one of the least squares method, the non-linear least squares method, the approximation method, and the genetic algorithm; the same receiving point can correspond to one or more first arrival travel time circle parameter radii and first arrival travel time errors; the first arrival travel time circle parameter can refer to that based on the assumption of an isotropic medium, the projection of the first arrival travel time with the same offset on the ground should be a circle, and the parameter describing this circle is the first arrival travel time circle parameter, and the radius of the first arrival travel time circle parameter can refer to the radius of this circle; the first arrival travel time error can refer to the error between the actual first arrival travel time and the theoretical first arrival travel time.

[0085] Step 103: Obtain the anisotropic static correction amount for each receiving point according to the radius of the first arrival travel time circle parameter, the first arrival travel time error, and the target 3D seismic first arrival data.

[0086] Specifically, the anisotropic static correction amount can refer to the static correction amount caused by the anisotropy of the velocity.

[0087] To further improve the accuracy of static correction, see Figure 2 , in an embodiment of the present application, before step 102, it further includes:

[0088] Step 201: Perform elevation static correction on the first arrival travel time in the target 3D seismic first arrival data, and correct the first arrival travel time to the same reference plane.

[0089] Specifically, fixed elevation plane static correction can be performed on the first arrival travel time in the entire 3D seismic first arrival data, and all first arrivals are corrected to a unified fixed reference plane.

[0090] To further improve the accuracy of static correction, in an embodiment of the present application, the target 3D seismic first arrival data includes: offset, the receiving point coordinates corresponding to the offset, the shot point coordinates, and the first arrival travel time; see Figure 3 , correspondingly, step 102 includes:

[0091] Step 301: Divide the offsets of common shot points into multiple groups.

[0092] Specifically, within common shot points, the first arrival travel times of the 3D seismic first arrival data can be sorted according to the offset from small to large, and the offsets within common shot points are grouped as needed. The more groups, the higher the accuracy of tomographic static correction and the greater the computational amount; the minimum grouping interval is the field acquisition offset interval, and the maximum grouping interval is the maximum offset.

[0093] Step 302: Based on the preset roundness error evaluation method, first arrival travel time, shot point coordinates, and receiving point coordinates, obtain the corresponding first arrival travel time circle parameter radius and first arrival travel time error for each group.

[0094] Specifically, based on the preset roundness error evaluation method, the first arrival travel time corresponding to each group, shot point coordinates, and receiving point coordinates, obtain the optimal first arrival travel time circle parameter radius and first arrival travel time error of the receiving point corresponding to each offset in each group.

[0095] Step 303: Use the first arrival travel time circle parameter radius and first arrival travel time error corresponding to all groups as the first arrival travel time circle parameter radius and first arrival travel time error of multiple receiving points corresponding to the target three-dimensional seismic first arrival data.

[0096] Specifically, all the first arrival travel time circle parameter radius and first arrival travel time error may include: the optimal first arrival travel time circle parameter radius and first arrival travel time error of the receiving points corresponding to multiple shot points respectively; different shot points may correspond to the same receiving point, so the same receiving point may correspond to multiple first arrival travel time circle parameter radius and first arrival travel time error.

[0097] In order to achieve omnidirectional velocity modeling on the basis of avoiding assumptions about the anisotropic direction of the medium, and then achieve anisotropic static correction to improve the reliability of static correction in complex areas, refer to Figure 4 , in an embodiment of the present application, step 103 includes:

[0098] Step 401: According to the first arrival travel time circle parameter radius, first arrival travel time error, and offset in the target three-dimensional seismic first arrival data, obtain the first anisotropic velocity correction amount for each receiving point.

[0099] Specifically, according to the first arrival travel time circle parameter radius, first arrival travel time error, and the distance from the receiving point to the shot point, obtain the isotropic velocity and the first anisotropic velocity correction amount for each receiving point. The isotropic velocity and the first anisotropic velocity correction amount for each receiving point can be equivalent to the velocity between the shot point and the receiving point; since the same receiving point may correspond to multiple first arrival travel time circle parameter radius and first arrival travel time error, each receiving point may correspond to multiple first anisotropic velocity correction amounts.

[0100] Step 402: Apply the first anisotropic velocity correction amount to determine the anisotropic static correction amount for each receiving point.

[0101] Specifically, in combination with conventional static correction calculation methods, the first anisotropic velocity correction amount can be applied to determine the anisotropic static correction amount for each receiving point. The conventional static correction calculation methods can be tomographic static correction methods. Tomographic static correction methods can utilize the first arrival wave to invert the surface velocity that varies longitudinally and transversely. The anisotropic velocity correction amount can refer to the velocity correction amount caused by anisotropy. The first anisotropic velocity correction amount can represent the anisotropic velocity correction amount before constrained solution, and there may be multiple first anisotropic velocity correction amounts at the same position.

[0102] Due to the overlap in the acquisition range, there are different values for the velocity distribution at some positions, and there may also be different values for the anisotropic velocity correction amount at the same position. To improve the accuracy of anisotropic static correction, in an embodiment of the present application, after step 401, the following is further included:

[0103] Step 403: Based on the preset constraint conditions, optimization algorithm, and the first anisotropic velocity correction amount, obtain the unique second anisotropic velocity correction amount for each receiving point.

[0104] Step 404: Apply the second anisotropic velocity correction amount to determine the anisotropic static correction amount for each receiving point.

[0105] Specifically, the preset constraint conditions can be at least one of distance constraint, micro-logging constraint, and structural constraint. The optimization algorithm can be one of the comprehensive averaging method, Kalman filtering method, Bayesian estimation method, statistical decision method, fuzzy logic method, and neural network method. After performing constrained solution on the multiple anisotropic velocity correction amounts for each receiving point, the unique second anisotropic velocity correction amount for each receiving point can be obtained. In combination with conventional static correction calculation methods, the second anisotropic velocity correction amount can be applied to determine the anisotropic static correction amount for each receiving point.

[0106] To improve the accuracy and efficiency of isotropic static correction on the basis of improving the accuracy and efficiency of anisotropic static correction, see Figure 5 , in an embodiment of the present application, after step 102, the following is further included:

[0107] Step 501: Based on the first arrival travel time circle parameter radius, obtain the isotropic static correction amount for each receiving point.

[0108] Specifically, the isotropic static correction amount can refer to the static correction amount corresponding to the isotropic velocity. To improve the accuracy of obtaining the static correction amount, the isotropic static correction amount can be obtained first. This part corresponds to the main part, and the anisotropic static correction amount can make appropriate supplements and modifications to the isotropic static correction amount.

[0109] In order to achieve isotropic static correction and improve the reliability of static correction in complex areas without making assumptions about the direction of medium anisotropy, in an embodiment of the present application, step 501 includes:

[0110] Step 601: Obtain the first isotropic velocity of each receiving point based on the first arrival travel time circle parameter radius.

[0111] Step 602: Apply the first isotropic velocity to determine the isotropic static correction amount of each receiving point.

[0112] Specifically, the isotropic velocity can be the background velocity, the first isotropic velocity can be the isotropic velocity before constraint solution, and there may be multiple first isotropic velocities at the same position.

[0113] Since there is an overlap in the acquisition range, different values may exist for the velocity distribution at some of the same positions. To improve the accuracy of isotropic static correction, in an embodiment of the present application, after step 601, it further includes:

[0114] Step 603: Based on a preset constraint condition, optimization algorithm, and the first isotropic velocity, obtain a unique second isotropic velocity for each receiving point.

[0115] Step 604: Apply the second isotropic velocity to determine the isotropic static correction amount of each receiving point.

[0116] Specifically, the second isotropic velocity can be the isotropic velocity after constraint solution, and there is a unique second isotropic velocity at the same position.

[0117] Specifically, the isotropic static correction amount of each receiving point can be determined by combining a conventional static correction calculation method and applying the second isotropic velocity; the conventional static correction calculation method can be the tomographic static correction method.

[0118] To further illustrate the present solution, the present application also provides an application example of a static correction method for three-dimensional seismic first arrival data. See Figure 6 , and the specific description is as follows:

[0119] S100: Fixed surface static correction; that is, perform elevation static correction on the first arrival travel times of the entire three-dimensional seismic first arrival data, and correct all first arrival travel times to a unified fixed reference surface H0.

[0120] S200: Common shot-offset sorting; that is, arrange the first arrival travel times of the three-dimensional seismic first arrival data according to common shot points, and sort the first arrival travel times within each common shot point from small to large according to offset: sort the first arrival travel times of the three-dimensional seismic first arrival data from small to large according to common shot-offset.

[0121] S300: Common shot offset grouping; the offsets within a common shot are automatically divided into N groups according to a preset offset interval ΔD (the minimum of ΔD is the field acquisition offset interval, and the maximum is the maximum offset); the function implemented by the common shot in this application example can be equivalent to the function implemented by a single shot. The stereo distribution of the first arrivals of a single shot can be as shown in Figure 7 as follows, Figure 7 where the abscissa and ordinate in it represent the geodetic coordinates, with the unit of meter; the vertical coordinate represents the first arrival travel time, with the unit of millisecond.

[0122] S400: Evaluation of the roundness error within the offset group to obtain the isotropic velocity within the offset group and the anisotropic velocity correction amount within the offset group; that is, select all the first arrival travel times within a common shot offset grouping. According to the shot point coordinates (sx, sy, sz) and the receiver point coordinates (rx, ry, rz) of the first arrival travel times, perform a roundness error evaluation within an offset grouping (the evaluation methods include but are not limited to the least squares method, the nonlinear least squares method, the approximation method, and the genetic algorithm), and obtain the radius T0 of the best first arrival travel time circle parameter for each point in the current offset grouping i,j and D0 i,j and the first arrival travel time error Tε for each point i,j,k , where i is the number of the first arrival travel time of a common shot, j is the number of the current offset grouping, and k is the point number within the offset grouping. For example, the equal value distribution of the first arrival plane can be as shown in Figure 8 as follows, Figure 8 where the abscissa and ordinate represent the geodetic coordinates, with the unit of meter; 250, 500, 750, 1000, and 1250 represent the first arrival travel times, with the unit of millisecond. The lines where 250, 500, 750, 1000, and 1250 are located respectively represent the actual first arrival equal value lines, and the lines where a1 to a4 are located respectively represent the best circle estimates corresponding to the radius of the best first arrival travel time circle parameter.

[0123] Combined with the distance D from the receiver point to the shot point i,j,k and the offset grouping j, calculate:

[0124]

[0125]

[0126] and T0 i,-1 =0, D0 i,-1 =0; the velocity between the shot point and the receiver point is the isotropic velocity V0 i,j and the anisotropic velocity correction amount Vε i,j,k , and the isotropic velocity V0 of each point within the same offset grouping is i,j the same, and the anisotropic velocity corrections of each point are not exactly the same.

[0127] S500: Determine whether all common shot points are completed. If so, execute step S600; otherwise, return to execute step S400, that is, repeat step S400 until all offset groups within all common shot points are completed, obtaining the velocity distribution V0 of the entire space x,y,z and the anisotropic velocity correction amount Vε x,y,z , due to the overlap of the acquisition range, some of V0 x,y,z has different values at the same position. Similarly, Vε at the same position x,y,z also has different values.

[0128] S600: Constraint solution to obtain the global isotropic velocity and the global anisotropic velocity correction amount; introduce constraint conditions (distance constraint, micro-logging constraint, structural constraint, etc.), and use optimization algorithms (comprehensive averaging method, Kalman filtering method, Bayesian estimation method, statistical decision method, fuzzy logic method, neural network method, etc.) to perform constraint solution on the entire V0 x,y,z and Vε x,y,z to obtain the unique isotropic velocity V0 at each position x,y,z and the anisotropic velocity correction amount Vε x,y,z , where the anisotropic velocity correction amount can be equivalent to the omnidirectional velocity correction amount.

[0129] S700: Obtain the global isotropic static correction amount.

[0130] S800: Obtain the global anisotropic static correction amount.

[0131] Specifically, the isotropic static correction amount corr_0 and the anisotropic static correction amount corr_w can be obtained according to the conventional method for obtaining static correction amounts.

[0132] S900: Comprehensive application of static correction amounts; comprehensively apply the isotropic static correction amount corr_0 and the anisotropic static correction amount corr_w, where the anisotropic static correction amount can be equivalent to the omnidirectional static correction amount; the anisotropic static correction amounts of all receiving points can form the global omnidirectional static correction amount.

[0133] As can be seen from the above description, the static correction method for three-dimensional seismic first arrival data provided by this application example decouples the first arrival travel time in the data domain on a unified reference surface, performs azimuthal anisotropic near-surface velocity modeling on the three-dimensional pre-stack seismic first arrival travel time, can fully consider the possible anisotropy in all directions, can effectively establish an accurate near-surface velocity model, complete omnidirectional static correction, can improve the static correction accuracy of three-dimensional seismic data, and is beneficial to the high-resolution processing of subsequent seismic data.

[0134] At the software level, in order to improve the accuracy of static correction, especially the accuracy of anisotropic static correction in complex areas, this application provides an embodiment of a static correction device for three-dimensional seismic first arrival data, which implements all or part of the content of the static correction method for the three-dimensional seismic first arrival data. See Figure 9 The static correction device for the three-dimensional seismic first arrival data specifically includes the following:

[0135] An acquisition module 01, configured to acquire target three-dimensional seismic first arrival data.

[0136] An error evaluation module 02, configured to obtain the first arrival travel time circle parameter radius and the first arrival travel time error of each of the plurality of receiving points corresponding to the target three-dimensional seismic first arrival data based on a preset roundness error evaluation method.

[0137] An anisotropic static correction module 03, configured to obtain the anisotropic static correction amount of each receiving point according to the first arrival travel time circle parameter radius, the first arrival travel time error, and the target three-dimensional seismic first arrival data.

[0138] See Figure 10 In an embodiment of this application, the static correction device for the three-dimensional seismic first arrival data further includes:

[0139] A correction module 04, configured to perform elevation static correction on the first arrival travel time in the target three-dimensional seismic first arrival data and correct the first arrival travel time to the same reference plane.

[0140] In an embodiment of this application, the target three-dimensional seismic first arrival data includes: offset, the receiving point coordinates corresponding to the offset, the shot point coordinates, and the first arrival travel time; correspondingly, the error evaluation module includes:

[0141] A grouping unit, configured to divide the offsets of common shot points into multiple groups;

[0142] An error evaluation unit, configured to obtain the first arrival travel time circle parameter radius and the first arrival travel time error corresponding to each group based on the preset roundness error evaluation method, the first arrival travel time, the shot point coordinates, and the receiving point coordinates;

[0143] A determination unit, configured to use the first arrival travel time circle parameter radius and the first arrival travel time error corresponding to all groups as the first arrival travel time circle parameter radius and the first arrival travel time error of the plurality of receiving points corresponding to the target three-dimensional seismic first arrival data.

[0144] In an embodiment of this application, the anisotropic static correction module includes:

[0145] A velocity correction unit, configured to obtain the first anisotropic velocity correction amount of each receiving point according to the first arrival travel time circle parameter radius, the first arrival travel time error, and the offset in the target three-dimensional seismic first arrival data;

[0146] The first anisotropic static correction unit is configured to determine the anisotropic static correction amount of each receiving point by applying the first anisotropic velocity correction amount.

[0147] In an embodiment of the present application, the anisotropic static correction module further includes:

[0148] The first optimization unit is configured to obtain a unique second anisotropic velocity correction amount for each receiving point based on a preset constraint condition, an optimization algorithm, and the first anisotropic velocity correction amount;

[0149] The second anisotropic static correction unit is configured to determine the anisotropic static correction amount of each receiving point by applying the second anisotropic velocity correction amount.

[0150] See Figure 11 , in an embodiment of the present application, the static correction device for three-dimensional seismic first arrival data further includes:

[0151] The isotropic static correction module 05 is configured to obtain the isotropic static correction amount of each receiving point based on the radius of the first arrival travel time circle parameter.

[0152] In an embodiment of the present application, the isotropic static correction module includes:

[0153] The isotropic velocity determination unit is configured to obtain the first isotropic velocity of each receiving point based on the radius of the first arrival travel time circle parameter.

[0154] The first isotropic static correction unit is configured to determine the isotropic static correction amount of each receiving point by applying the first isotropic velocity.

[0155] In an embodiment of the present application, the isotropic static correction module further includes:

[0156] The second optimization unit is configured to obtain a unique second isotropic velocity for each receiving point based on a preset constraint condition, an optimization algorithm, and the first isotropic velocity.

[0157] The second isotropic static correction unit is configured to determine the isotropic static correction amount of each receiving point by applying the second isotropic velocity.

[0158] The embodiments of the static correction device for three-dimensional seismic first arrival data provided in this specification can specifically be used to execute the processing flow of the embodiments of the above-mentioned static correction method for three-dimensional seismic first arrival data, and its functions will not be elaborated here. For details, reference can be made to the detailed description of the embodiments of the static correction method for three-dimensional seismic first arrival data.

[0159] As can be seen from the above description, the static correction method and device for three-dimensional seismic first arrival data provided by the present application can improve the accuracy of static correction, especially the accuracy of anisotropic static correction in complex areas. Specifically, it fully considers the possible anisotropy in all directions, can effectively establish an accurate near-surface velocity model, can complete all-round static correction, can improve the static correction accuracy of three-dimensional seismic data, and thus is conducive to the subsequent high-resolution processing of seismic data. It is of great practical significance in improving the quality of static correction, shortening the static correction period, saving energy, reducing emissions, and reducing costs and increasing efficiency.

[0160] From the hardware level, in order to improve the accuracy of static correction, especially the accuracy of anisotropic static correction in complex areas, the present application provides an embodiment of an electronic device for implementing all or part of the content in the static correction method for the three-dimensional seismic first arrival data. The electronic device specifically includes the following:

[0161] A processor, a memory, a communication interface, and a bus; wherein, the processor, the memory, and the communication interface complete communication with each other through the bus; the communication interface is used to implement information transmission between the static correction device for the three-dimensional seismic first arrival data and related devices such as user terminals. The electronic device can be a desktop computer, a tablet computer, a mobile terminal, etc., and this embodiment is not limited thereto. In this embodiment, the electronic device can be implemented with reference to the embodiments for implementing the static correction method for the three-dimensional seismic first arrival data and the embodiments for implementing the static correction device for the three-dimensional seismic first arrival data, and the content is incorporated herein, and the repeated parts will not be elaborated.

[0162] Figure 12 It is a schematic block diagram of the system composition of the electronic device 9600 according to an embodiment of the present application. As Figure 12 shown, the electronic device 9600 may include a central processing unit 9100 and a memory 9140; the memory 9140 is coupled to the central processing unit 9100. It should be noted that this Figure 12 is exemplary; other types of structures can also be used to supplement or replace this structure to achieve telecommunication functions or other functions.

[0163] In one or more embodiments of the present application, the static correction function for the three-dimensional seismic first arrival data can be integrated into the central processing unit 9100. Among them, the central processing unit 9100 can be configured to perform the following controls:

[0164] Step 101: Obtain target three-dimensional seismic first arrival data.

[0165] Step 102: Based on a preset roundness error evaluation method, obtain the radius of the first arrival travel time circle parameter and the first arrival travel time error of each of the multiple receiving points corresponding to the target 3D seismic first arrival data.

[0166] Step 103: According to the radius of the first arrival travel time circle parameter, the first arrival travel time error, and the target 3D seismic first arrival data, obtain the anisotropic static correction amount of each receiving point.

[0167] As can be seen from the above description, the electronic device provided by the embodiments of the present application can improve the accuracy of static correction, especially the accuracy of anisotropic static correction in complex areas.

[0168] In another embodiment, the static correction device for 3D seismic first arrival data can be separately configured from the central processing unit 9100. For example, the static correction device for 3D seismic first arrival data can be configured as a chip connected to the central processing unit 9100, and the static correction function of the 3D seismic first arrival data is realized through the control of the central processing unit.

[0169] As Figure 12 shown, the electronic device 9600 may further include: a communication module 9110, an input unit 9120, an audio processor 9130, a display 9160, and a power supply 9170. It should be noted that the electronic device 9600 does not necessarily have to include all the components shown in Figure 12 ; in addition, the electronic device 9600 may further include components not shown in Figure 12 , and reference can be made to the prior art.

[0170] As Figure 12 shown, the central processing unit 9100 is sometimes also referred to as a controller or an operation control, and may include a microprocessor or other processor devices and / or logic devices. The central processing unit 9100 receives inputs and controls the operations of the various components of the electronic device 9600.

[0171] Among them, the memory 9140 can be, for example, one or more of a buffer, a flash memory, a hard drive, a removable medium, a volatile memory, a non-volatile memory, or other suitable devices. The above information related to failures can be stored, and in addition, programs for executing relevant information can also be stored. And the central processing unit 9100 can execute the programs stored in the memory 9140 to implement information storage or processing, etc.

[0172] The input unit 9120 provides inputs to the central processing unit 9100. The input unit 9120 is, for example, a key or a touch input device. The power supply 9170 is used to supply power to the electronic device 9600. The display 9160 is used to display display objects such as images and texts. The display can be, for example, an LCD display, but is not limited thereto.

[0173] The memory 9140 may be a solid-state memory, for example, a read-only memory (ROM), a random access memory (RAM), a SIM card, etc. It may also be a memory that stores information even when power is off, can be selectively erased and has more data. An example of this memory is sometimes referred to as an EPROM, etc. The memory 9140 may also be some other type of device. The memory 9140 includes a buffer memory 9141 (sometimes referred to as a buffer). The memory 9140 may include an application / function storage unit 9142, which is used to store application programs and function programs or the processes for operating the electronic device 9600 through the central processor 9100.

[0174] The memory 9140 may also include a data storage unit 9143, which is used to store data, such as contacts, digital data, pictures, sounds, and / or any other data used by the electronic device. The driver storage unit 9144 of the memory 9140 may include various drivers of the electronic device for communication functions and / or for performing other functions of the electronic device (such as a messaging application, an address book application, etc.).

[0175] The communication module 9110 is a transmitter / receiver 9110 that transmits and receives signals via the antenna 9111. The communication module (transmitter / receiver) 9110 is coupled to the central processor 9100 to provide input signals and receive output signals, which may be the same as in the case of a conventional mobile communication terminal.

[0176] Based on different communication technologies, multiple communication modules 9110 may be provided in the same electronic device, such as a cellular network module, a Bluetooth module, and / or a wireless local area network module, etc. The communication module (transmitter / receiver) 9110 is also coupled to the speaker 9131 and the microphone 9132 via the audio processor 9130 to provide an audio output via the speaker 9131 and receive an audio input from the microphone 9132, so as to implement the usual telecommunication functions. The audio processor 9130 may include any suitable buffer, decoder, amplifier, etc. In addition, the audio processor 9130 is also coupled to the central processor 9100, so that recording can be performed on the local machine through the microphone 9132, and the sound stored on the local machine can be played through the speaker 9131.

[0177] As can be seen from the above description, the electronic device provided by the embodiments of the present application can improve the accuracy of static correction, especially the accuracy of anisotropic static correction in complex areas.

[0178] Embodiments of the present application also provide a computer-readable storage medium capable of implementing all steps in the static correction method for three-dimensional seismic first arrival data in the above embodiments. A computer program is stored on the computer-readable storage medium, and when the computer program is executed by a processor, all steps of the static correction method for three-dimensional seismic first arrival data in the above embodiments are implemented. For example, when the processor executes the computer program, the following steps are implemented:

[0179] Step 101: Obtain target three-dimensional seismic first arrival data.

[0180] Step 102: Based on a preset roundness error evaluation method, obtain the first arrival travel time circle parameter radius and the first arrival travel time error of each of the multiple receiving points corresponding to the target three-dimensional seismic first arrival data.

[0181] Step 103: According to the first arrival travel time circle parameter radius, the first arrival travel time error, and the target three-dimensional seismic first arrival data, obtain the anisotropic static correction amount of each receiving point.

[0182] As can be seen from the above description, the computer-readable storage medium provided by the embodiments of the present application can improve the accuracy of static correction, especially the accuracy of anisotropic static correction in complex areas.

[0183] In the present application, each embodiment of the above method is described in a progressive manner. The same or similar parts between each embodiment can be referred to each other, and the key points of each embodiment are the differences from other embodiments. For the relevant parts, refer to the partial description of the method embodiment.

[0184] Those skilled in the art should understand that the embodiments of the present application can be provided as a method, a system, or a computer program product. Therefore, the present application can take the form of a complete hardware embodiment, a complete software embodiment, or an embodiment combining software and hardware aspects. Moreover, the present application can take the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) containing computer-usable program code.

[0185] The present application is described with reference to the flowcharts and / or block diagrams of methods, devices (systems), and computer program products according to the embodiments of the present application. It should be understood that each process and / or block in the flowchart and / or block diagram can be implemented by computer program instructions, and the combination of processes and / or blocks in the flowchart and / or block diagram can also be implemented by computer program instructions. These computer program instructions can be provided to the processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing devices to generate a machine, so that the instructions executed by the processor of the computer or other programmable data processing devices generate for implementing in the process Figure 1 one process or multiple processes and / or blocksFigure 1 means for the functions specified in one or more boxes.

[0186] These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to work in a particular manner, such that the instructions stored in the computer-readable memory produce a manufacture including an instruction means that implements the functions specified in one Figure 1 process or processes and / or boxes Figure 1 or more boxes.

[0187] These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus, such that a series of operational steps are performed on the computer or other programmable apparatus to produce a computer-implemented process, whereby the instructions executed on the computer or other programmable apparatus provide steps for implementing the functions specified in one Figure 1 process or processes and / or boxes Figure 1 or more boxes.

[0188] Specific embodiments are used in this application to elaborate on the principles and implementation manners of this application. The description of the above embodiments is only used to help understand the method and its core idea of this application; at the same time, for those of ordinary skill in the art, based on the idea of this application, there will be changes in the specific implementation manners and application scopes. In summary, the content of this specification should not be construed as a limitation to this application.

Claims

1. A static correction method for three-dimensional seismic first arrival data, characterized in that, Including: Obtaining target three-dimensional seismic first arrival data; Based on a preset roundness error evaluation method, obtaining the first arrival travel time circle parameter radius and the first arrival travel time error of each of the multiple receiving points corresponding to the target three-dimensional seismic first arrival data; According to the first arrival travel time circle parameter radius, the first arrival travel time error, and the target three-dimensional seismic first arrival data, obtaining the anisotropic static correction amount of each receiving point; The target three-dimensional seismic first arrival data includes: first arrival travel time, receiving point coordinates, shot point coordinates, and offset; correspondingly, the obtaining the first arrival travel time circle parameter radius and the first arrival travel time error of each of the multiple receiving points corresponding to the target three-dimensional seismic first arrival data based on a preset roundness error evaluation method includes: Dividing the offsets of common shot points into multiple groups; Based on the preset roundness error evaluation method, the first arrival travel time, the shot point coordinates, and the receiving point coordinates, obtaining the first arrival travel time circle parameter radius and the first arrival travel time error corresponding to each group; Taking the first arrival travel time circle parameter radius and the first arrival travel time error corresponding to all groups as the first arrival travel time circle parameter radius and the first arrival travel time error of the multiple receiving points corresponding to the target three-dimensional seismic first arrival data; The obtaining the anisotropic static correction amount of each receiving point according to the first arrival travel time circle parameter radius, the first arrival travel time error, and the target three-dimensional seismic first arrival data includes: According to the first arrival travel time circle parameter radius, the first arrival travel time error, and the offset in the target three-dimensional seismic first arrival data, obtaining the first anisotropic velocity correction amount of each receiving point; Applying the first anisotropic velocity correction amount to determine the anisotropic static correction amount of each receiving point.

2. The static correction method for three-dimensional seismic first arrival data according to claim 1, characterized in that, Before the obtaining the first arrival travel time circle parameter radius and the first arrival travel time error of each of the multiple receiving points corresponding to the target three-dimensional seismic first arrival data based on a preset roundness error evaluation method, it further includes: Performing elevation static correction on the first arrival travel time in the target three-dimensional seismic first arrival data, and correcting the first arrival travel time to the same reference plane.

3. The static correction method for three-dimensional seismic first arrival data according to claim 1, characterized in that, After the obtaining the first anisotropic velocity correction amount of each receiving point according to the first arrival travel time circle parameter radius, the first arrival travel time error, and the offset in the target three-dimensional seismic first arrival data, it further includes: Based on a preset constraint condition, an optimization algorithm, and the first anisotropic velocity correction amount, obtaining a unique second anisotropic velocity correction amount of each receiving point; Applying the second anisotropic velocity correction amount to determine the anisotropic static correction amount of each receiving point.

4. The static correction method for three-dimensional seismic first arrival data according to claim 1, wherein, After the obtaining the first arrival travel time circle parameter radius and the first arrival travel time error of each of the multiple receiving points corresponding to the target three-dimensional seismic first arrival data based on a preset roundness error evaluation method, it further includes: Based on the first arrival travel time circle parameter radius, obtaining the isotropic static correction amount of each receiving point.

5. The static correction method for three-dimensional seismic first arrival data according to claim 4, characterized in that, The obtaining the isotropic static correction amount of each receiving point based on the first arrival travel time circle parameter radius includes: Based on the first arrival travel time circle parameter radius, obtaining the first isotropic velocity of each receiving point; Applying the first isotropic velocity to determine the isotropic static correction amount of each receiving point.

6. The static correction method for three-dimensional seismic first arrival data according to claim 5, characterized in that, After the obtaining the first isotropic velocity of each receiving point based on the first arrival travel time circle parameter radius, it further includes: Based on the preset constraint conditions, optimization algorithm, and the first isotropic velocity, obtain the unique second isotropic velocity for each receiving point; Apply the second isotropic velocity to determine the isotropic static correction amount for each receiving point.

7. A static correction device for three-dimensional seismic first arrival data, characterized in that It includes: An acquisition module for acquiring the target three-dimensional seismic first arrival data; An error evaluation module for obtaining the first arrival travel time circle parameter radius and first arrival travel time error of each of the multiple receiving points corresponding to the target three-dimensional seismic first arrival data based on a preset roundness error evaluation method; An anisotropic static correction module for obtaining the anisotropic static correction amount for each receiving point according to the first arrival travel time circle parameter radius, first arrival travel time error, and the target three-dimensional seismic first arrival data; The target three-dimensional seismic first arrival data includes: offset, the receiving point coordinates corresponding to the offset, shot point coordinates, and first arrival travel time; correspondingly, the error evaluation module includes: A grouping unit for dividing the offsets of common shot points into multiple groups; An error evaluation unit for obtaining the first arrival travel time circle parameter radius and first arrival travel time error corresponding to each group based on the preset roundness error evaluation method, first arrival travel time, shot point coordinates, and receiving point coordinates; A determination unit for using the first arrival travel time circle parameter radius and first arrival travel time error corresponding to all groups as the first arrival travel time circle parameter radius and first arrival travel time error of the multiple receiving points corresponding to the target three-dimensional seismic first arrival data; The anisotropic static correction module includes: A velocity correction unit for obtaining the first anisotropic velocity correction amount for each receiving point according to the first arrival travel time circle parameter radius, first arrival travel time error, and the offset in the target three-dimensional seismic first arrival data; A first anisotropic static correction unit for applying the first anisotropic velocity correction amount to determine the anisotropic static correction amount for each receiving point.

8. The static correction device for three-dimensional seismic first arrival data according to claim 7, characterized in that, It further includes: A correction module for performing elevation static correction on the first arrival travel time in the target three-dimensional seismic first arrival data and correcting the first arrival travel time to the same reference plane.

9. The static correction device for three-dimensional seismic first arrival data according to claim 7, characterized in that, The anisotropic static correction module further includes: A first optimization unit for obtaining the unique second anisotropic velocity correction amount for each receiving point based on the preset constraint conditions, optimization algorithm, and the first anisotropic velocity correction amount; A second anisotropic static correction unit for applying the second anisotropic velocity correction amount to determine the anisotropic static correction amount for each receiving point.

10. The static correction device for three-dimensional seismic first arrival data according to claim 7, characterized in that, It further includes: An isotropic static correction module for obtaining the isotropic static correction amount for each receiving point based on the first arrival travel time circle parameter radius.

11. The static correction device for three-dimensional seismic first arrival data according to claim 10, characterized in that, The isotropic static correction module includes: An isotropic velocity determination unit for obtaining the first isotropic velocity for each receiving point based on the first arrival travel time circle parameter radius; A first isotropic static correction unit for applying the first isotropic velocity to determine the isotropic static correction amount for each receiving point.

12. The static correction device for three-dimensional seismic first arrival data according to claim 11, wherein The isotropic static correction module further includes: A second optimization unit for obtaining the unique second isotropic velocity for each receiving point based on the preset constraint conditions, optimization algorithm, and the first isotropic velocity; A second isotropic static correction unit, configured to determine the isotropic static correction amount for each receiving point by applying the second isotropic velocity.

13. An electronic device, comprising a memory, a processor, and a computer program stored on the memory and executable on the processor, characterized in that, When the processor executes the program, it implements the static correction method for three-dimensional seismic first arrival data according to any one of claims 1 to 6.

14. A computer-readable storage medium having computer instructions stored thereon, characterized in that, When the instructions are executed, they implement the static correction method for three-dimensional seismic first arrival data according to any one of claims 1 to 6.

Images