Two-dimensional batch chromatography static correction method and device based on three-dimensional chromatography

Abstract

The invention provides a two-dimensional batch chromatography static correction method and device based on three-dimensional chromatography, and the method comprises the steps: picking up the first arrival time of all measuring lines in a target work area through preset picking parameters, so as to obtain the single shot first arrival time of the target work area; wherein all the measuring lines are two-dimensional seismic measuring lines; constructing a surface velocity model of the target work area based on the single shot first arrival time; and determining the tomographic static correction value of the target work area according to the surface velocity model. According to the method, the accurate near-surface velocity model can be obtained, and the closing problem of the measuring line horizon (amplitude energy and time) of the two-dimensional work area is solved. Secondly, the near-surface velocity model is obtained through integral inversion in combination with interpolation and extension, the near-surface velocity precision is effectively improved, and a reliable initial velocity field is provided for subsequent depth domain imaging.

Description

Technical Field

[0001] This application belongs to the technical field of oil and gas exploration, especially the technical field of seismic data processing for oil and gas exploration. Specifically, it relates to a two-dimensional batch tomography static correction method and device based on three-dimensional tomography. Background Art

[0002] It can be understood that the calculation accuracy of static correction directly affects subsequent processing links such as denoising and velocity analysis, restricts the signal-to-noise ratio of the stacked section and the accuracy of the structural form, and is one of the key steps in seismic data processing in areas with complex surface conditions. In the prior art, there are many calculation methods for static correction amounts, such as elevation static correction, refraction static correction, tomography static correction, etc.

[0003] The elevation static correction method is only applicable to work areas where the surface velocity is relatively stable. The refraction static correction requires a relatively stable high-velocity layer interface in the work area where refracted waves can be continuously traced throughout the work area. Tomography static correction can make full use of first arrival information such as direct waves, refracted waves, and diffracted waves, and use the tomography inversion method to invert the near-surface velocity. Compared with the other two static correction methods, in the case of a smaller grid, a more accurate shallow velocity structure can be obtained. Therefore, tomography static correction is a more applicable static correction technology compared to other static correction methods.

[0004] Generally, the calculation of the tomography static correction amount consists of two parts: the weathered layer static correction amount and the elevation correction amount. Among them, the purpose of weathered layer correction is to correct the seismic record from the surface to the defined top boundary of the high-velocity layer to correct the influence of the shallow low-velocity zone on the travel time of depth-reflected waves. This step is based on the shallow velocity obtained by tomography inversion. Since the top interface of the high-velocity layer has large fluctuations and the single-shot time shift is often large after elevation correction, for the need of subsequent data interpretation, the elevation correction method is often used to correct the seismic data from the top interface of the high-velocity layer to a fixed reference surface.

[0005] In the prior art, although the existing two-dimensional tomography static correction technology can better solve the static correction problem of a single survey line. However, due to the existence of multiple survey lines in the work area, affected by differences in acquisition years, surface lithology, first arrival picking, etc., it is difficult for the inversion of a single survey line to solve the unification of the velocity models at the closure points of a large number of survey lines, thus affecting the closure of the horizons of two-dimensional survey lines. In addition, affected by factors such as less micro-logging information in single-line inversion, it is impossible to comprehensively consider micro-logging information on the surface, which will inevitably reduce the accuracy of model inversion. The existence of these factors seriously affects subsequent interpretation work and cannot meet the needs of current fine exploration. Summary of the Invention

[0006] The present invention belongs to the technical field of seismic data processing. One object of the present invention is to solve the technical pain points of non-closure of the multi-line near-surface velocity model and non-closure of seismic results in a 2D work area. Secondly, a relatively reliable 3D near-surface velocity model for the 2D work area is provided through inversion interpolation, so as to provide an accurate near-surface area velocity model for depth-domain imaging.

[0007] Another object of the present invention is to provide a 2D batch tomography static correction device based on 3D tomography. Still another object of the present invention is to provide an electronic device, which includes a memory and a processor. The memory stores a computer program, and when the processor executes the computer program, the steps of the above-mentioned 2D batch tomography static correction method based on 3D tomography are implemented. Still another object of the present invention is to provide a readable medium, on which a computer program is stored, and when the computer program is executed by the processor, the steps of the above-mentioned 2D batch tomography static correction method based on 3D tomography are implemented.

[0008] To solve the technical problems in the background art of the present application, the present invention provides the following technical solutions:

[0009] In a first aspect, the present invention provides a 2D batch tomography static correction method based on 3D tomography, including:

[0010] Picking up the first arrival times of all survey lines in the target work area through preset picking parameters to obtain the single-shot first arrival times of the target work area; wherein, all the survey lines are 2D seismic survey lines;

[0011] Constructing a surface velocity model of the target work area based on the single-shot first arrival times;

[0012] Determining the tomography static correction amount of the target work area according to the surface velocity model.

[0013] In an embodiment of the present invention, before picking up the first arrival times of all survey lines in the target work area through preset picking parameters to obtain the single-shot first arrival times of the target work area, it further includes:

[0014] Preprocessing all the survey lines.

[0015] In an embodiment of the present invention, preprocessing all the survey lines includes:

[0016] Judging whether the elevations of the shot points and the receiving points of all the survey lines are closed;

[0017] If not closed, smoothing the non-closed survey lines.

[0018] In an embodiment of the present invention, constructing a surface velocity model of the target work area based on the single-shot first arrival times includes:

[0019] Based on the first arrival time of a single shot, and according to the shot point coordinates, geophone point coordinates, and surface elevation of all survey lines, the surface velocity model is constructed.

[0020] In an embodiment of the present invention, before determining the tomographic static correction amount of the target work area according to the surface velocity model, it further includes:

[0021] Correcting the surface velocity model according to the micro-log data in the target work area.

[0022] In an embodiment of the present invention, determining the tomographic static correction amount of the target work area according to the surface velocity model includes:

[0023] Determining the three-dimensional high-velocity layer top interface of the target work area according to the surface velocity model;

[0024] Determining the tomographic static correction amount according to the three-dimensional high-velocity layer top interface.

[0025] In an embodiment of the present invention, before determining the tomographic static correction amount according to the three-dimensional high-velocity layer top interface, it further includes:

[0026] Smoothing the three-dimensional high-velocity layer top interface.

[0027] In a second aspect, the present invention provides a two-dimensional batch tomographic static correction device based on three-dimensional tomography, and the device includes:

[0028] The first arrival time picking module is used to pick the first arrival time of all survey lines in the target work area through preset picking parameters to obtain the first arrival time of a single shot in the target work area; wherein, all the survey lines are two-dimensional seismic survey lines;

[0029] The surface velocity model construction module is used to construct the surface velocity model of the target work area based on the first arrival time of a single shot;

[0030] The tomographic static correction amount determination module is used to determine the tomographic static correction amount of the target work area according to the surface velocity model.

[0031] In an embodiment of the present invention, a two-dimensional batch tomographic static correction device based on three-dimensional tomography further includes:

[0032] The survey line preprocessing module is used to preprocess all the survey lines.

[0033] In an embodiment of the present invention, the survey line preprocessing module includes:

[0034] The closure judgment unit is used to judge whether the shot point elevation and geophone point elevation of all the survey lines are closed;

[0035] A survey line smoothing unit, configured to smooth the unclosed survey lines if they are not closed.

[0036] In an embodiment of the present invention, the surface velocity model construction module includes:

[0037] A surface velocity model construction unit, configured to construct the surface velocity model based on the first arrival time of a single shot, and according to the shot point coordinates, geophone point coordinates, and surface elevation of all the survey lines.

[0038] In an embodiment of the present invention, a two-dimensional batch tomographic static correction device based on three-dimensional tomography further includes:

[0039] A velocity model correction module, configured to correct the surface velocity model according to the micro-log data in the target work area.

[0040] In an embodiment of the present invention, the tomographic static correction amount determination module includes:

[0041] A high-velocity layer top interface determination unit, configured to determine the three-dimensional high-velocity layer top interface of the target work area according to the surface velocity model;

[0042] A tomographic static correction amount determination unit, configured to determine the tomographic static correction amount according to the three-dimensional high-velocity layer top interface.

[0043] In an embodiment of the present invention, a two-dimensional batch tomographic static correction device based on three-dimensional tomography further includes:

[0044] A high-velocity layer top interface smoothing module, configured to smooth the three-dimensional high-velocity layer top interface.

[0045] In a third aspect, the present invention provides a computer program product, including a computer program / instructions, which when executed by a processor, implement the steps of a two-dimensional batch tomographic static correction method based on three-dimensional tomography.

[0046] In a fourth aspect, the present invention 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, it implements the steps of a two-dimensional batch tomographic static correction method based on three-dimensional tomography.

[0047] In a fifth aspect, the present invention provides a computer-readable storage medium, on which a computer program is stored. When the computer program is executed by a processor, it implements the steps of a two-dimensional batch tomographic static correction method based on three-dimensional tomography.

[0048] As can be seen from the above description, the embodiments of the present invention provide a two-dimensional batch tomography static correction method and device based on three-dimensional tomography. The corresponding two-dimensional batch tomography static correction method based on three-dimensional tomography includes: First, pick up the first arrival times of all survey lines in the target work area through preset picking parameters to obtain the single-shot first arrival times of the target work area; where all survey lines are two-dimensional seismic survey lines; then, construct the surface velocity model of the target work area based on the single-shot first arrival times; finally, determine the tomography static correction amount of the target work area according to the surface velocity model.

[0049] The corresponding two-dimensional batch tomography static correction device based on three-dimensional tomography includes: a first arrival time picking module, configured to pick up the first arrival times of all survey lines in the target work area through preset picking parameters to obtain the single-shot first arrival times of the target work area; where all survey lines are two-dimensional seismic survey lines; a surface velocity model construction module, configured to construct the surface velocity model of the target work area based on the single-shot first arrival times; a tomography static correction amount determination module, configured to determine the tomography static correction amount of the target work area according to the surface velocity model.

[0050] The two-dimensional batch tomography static correction method and device based on three-dimensional tomography provided by the embodiments of the present invention are used for near-surface velocity modeling in complex two-dimensional work areas, can obtain an accurate near-surface velocity model, and solve the problem of line layer position (amplitude energy, time) closure in two-dimensional work areas, especially the closure problem in the depth domain.

[0051] Secondly, obtain the near-surface velocity model through overall inversion combined with interpolation and extrapolation, effectively improving the accuracy of the near-surface velocity and providing a reliable initial velocity field for subsequent depth domain imaging. BRIEF DESCRIPTION OF THE DRAWINGS

[0052] In order to more clearly illustrate the technical solutions in the embodiments of the present invention 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 some embodiments of the present invention. For those of ordinary skill in the art, other drawings can be obtained based on these drawings without creative efforts.

[0053] Figure 1 It is a schematic flowchart of a two-dimensional batch tomography static correction method based on three-dimensional tomography in an embodiment of the present invention;

[0054] Figure 2 It is another schematic flowchart of a two-dimensional batch tomography static correction method based on three-dimensional tomography in an embodiment of the present invention;

[0055] Figure 3 It is a schematic flowchart of step 900 of a two-dimensional batch tomography static correction method based on three-dimensional tomography in an embodiment of the present invention;

[0056] Figure 4 Schematic diagram of the third process of the two-dimensional batch tomography static correction method based on three-dimensional tomography in the embodiments of the present invention;

[0057] Figure 5 Schematic diagram of one process of step 200 of the two-dimensional batch tomography static correction method based on three-dimensional tomography in the embodiments of the present invention;

[0058] Figure 6 Schematic diagram of another process of step 200 of the two-dimensional batch tomography static correction method based on three-dimensional tomography in the embodiments of the present invention;

[0059] Figure 7 Schematic diagram of the process of the two-dimensional batch tomography static correction method based on three-dimensional tomography in the specific embodiments of the present invention;

[0060] Figure 8 Schematic diagram of the near-surface velocity model of single-line inversion in the specific embodiments of the present invention;

[0061] Figure 9 Schematic diagram of the initial near-surface velocity model constructed after fusing the first arrival time files of all survey lines in the work area in the specific embodiments of the present invention;

[0062] Figure 10 Schematic diagram of the near-surface velocity model obtained by unified inversion of all survey lines in the work area in the specific embodiments of the present invention;

[0063] Figure 11 Schematic diagram of the comparison of the stacked profiles at the intersection of single-line tomography inversion in the specific embodiments of the present invention (survey line A (left), survey line B (right));

[0064] Figure 12 Schematic diagram of the comparison of the stacked profiles at the intersection of overall three-dimensional tomography inversion in the specific embodiments of the present invention (survey line A (left), survey line B (right));

[0065] Figure 13 Schematic diagram of the overall three-dimensional velocity model of the two-dimensional work area in the specific embodiments of the present invention (yellow and above is the near-surface velocity model);

[0066] Figure 14 Block diagram of the two-dimensional batch tomography static correction device based on three-dimensional tomography in the specific embodiments of the present invention;

[0067] Figure 15 Schematic diagram of the structure of the electronic device in the embodiments of the present invention. Specific embodiments

[0068] To make the objectives, technical solutions, and advantages of the embodiments of the present invention clearer, the following will clearly and completely describe the technical solutions in the embodiments of the present invention in conjunction with the accompanying drawings in the embodiments of the present invention. Obviously, the described embodiments are part of the embodiments of the present invention, rather than all of them. All other embodiments obtained by those of ordinary skill in the art based on the embodiments of the present invention without creative efforts belong to the scope of protection of the present invention.

[0069] Those skilled in the art should understand that the embodiments of the present invention can be provided as methods, systems, or computer program products. Therefore, the present invention can take the form of a completely hardware embodiment, a completely software embodiment, or an embodiment combining software and hardware aspects. Moreover, the present invention 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.

[0070] It should be noted that the terms "including" and "having" in the specification and claims of this application and any variations thereof in the above-mentioned drawings are intended to cover non-exclusive inclusion. For example, a process, method, system, product, or device that includes a series of steps or units does not necessarily have to be limited to those steps or units clearly listed, but may include other steps or units not clearly listed or inherent to these processes, methods, products, or devices. Without conflict, the embodiments in this application and the features in the embodiments can be combined with each other. The following will detail this application with reference to the drawings and in conjunction with the embodiments.

[0071] In the technical solutions of this application, the acquisition, storage, use, processing, etc. of data all comply with the relevant regulations of laws and regulations.

[0072] Embodiment 1:

[0073] Although the existing two-dimensional tomography static correction technology can better solve the static correction problem of a single survey line. However, due to the existence of multiple survey lines in the work area, affected by differences in acquisition years, surface lithology, first arrival picking, etc., it is difficult to unify the velocity models at the closed points of a large number of survey lines by the inversion of a single survey line, thus affecting the closure of the horizons of the two-dimensional survey lines and unable to meet the needs of current fine exploration. Based on this, the embodiments of the present invention provide a specific implementation manner of a two-dimensional batch tomography static correction method based on three-dimensional tomography. See Figure 1 , which specifically includes the following content:

[0074] Step 100: Pick the first arrival times of all survey lines in the target work area through preset picking parameters to obtain the single-shot first arrival times of the target work area; wherein, all the survey lines are two-dimensional seismic survey lines;

[0075] Step 200: Construct a surface velocity model of the target work area based on the single-shot first arrival time;

[0076] Step 300: Determine the tomographic static correction amount of the target work area according to the surface velocity model.

[0077] As can be seen from the above description, the embodiment of the present invention provides a two-dimensional batch tomographic static correction method based on three-dimensional tomography, including: first, pick up the first arrival time of all survey lines in the target work area through preset picking parameters to obtain the single-shot first arrival time of the target work area; wherein, all survey lines are two-dimensional seismic survey lines; then, construct a surface velocity model of the target work area based on the single-shot first arrival time; finally, determine the tomographic static correction amount of the target work area according to the surface velocity model.

[0078] The two-dimensional batch tomographic static correction method based on three-dimensional tomography provided by the embodiment of the present invention is used for near-surface velocity modeling in complex two-dimensional work areas. This method first solves the problems of non-closure of the multi-line near-surface velocity model and non-closure of seismic results in two-dimensional work areas. Secondly, it provides a relatively reliable three-dimensional near-surface velocity model for the two-dimensional work area through inversion interpolation, providing an accurate near-surface regional velocity model for depth-domain imaging, thereby solving the problem of closure of the survey line horizons (amplitude energy, time) in two-dimensional work areas, especially the closure problem in the depth domain.

[0079] Secondly, the near-surface velocity model is obtained through overall inversion combined with interpolation and extrapolation, effectively improving the near-surface velocity accuracy and providing a reliable initial velocity field for subsequent depth-domain imaging.

[0080] Embodiment 2:

[0081] In seismic exploration, when the seismic wavefront reaches a certain observation point, the moment when the geophone first detects the particle vibration at the observation point is called the first arrival time of the seismic wave, simply referred to as the first arrival (First arrival, First break). In the daily communication among industry insiders, the first arrival is also often used as an abbreviation for "first arrival wave".

[0082] When implementing Step 100, pick up the first arrival time corresponding to all survey lines with unified picking parameters. Specifically, set a time window (length), divide the time window into two equal parts before and after, calculate the ratio of the energy sum of the two parts before and after the time window (after / before); determine the time window where the first arrival time is located according to the change of the ratio, usually take the maximum value; pick up the first arrival according to relevant criteria within this time window, preferably take the midpoint of the time window.

[0083] The single-shot first arrival time is the time when the seismic waves first reach the seismic receiver after being generated by a single seismic source (or shot). Therefore, the single-shot first arrival time in Step 100 is a substitute value of the theoretical single-shot first arrival time.

[0084] By measuring the first arrival time of single-shot seismic waves, the distance and velocity of seismic wave propagation can be deduced, and then the subsurface geological structure can be inferred, including information such as the depth and thickness of rock layers and the dip of strata.

[0085] Two-dimensional seismic lines collect seismic data by deploying a series of seismic sources (such as vibrating trucks or air guns) and seismic receivers (such as seismographs or seismic probes) on the earth's surface. These seismic sources and receivers are deployed along a straight line (i.e., the seismic line).

[0086] In a typical two-dimensional seismic line exploration process, first, seismic waves are generated on the earth's surface by seismic sources. These seismic waves propagate downward and are reflected when they encounter subsurface rock layers. The reflected seismic waves are received and recorded by seismic receivers. By analyzing the recorded seismic data (including the first arrival time, amplitude, frequency, etc. of seismic waves), the subsurface geological structure can be inferred, including information such as the depth and thickness of rock layers and the dip of strata.

[0087] The result of this method is a two-dimensional seismic profile, which shows the subsurface geological structure under the seismic line. However, this two-dimensional seismic profile can only provide geological information under one seismic line. If more comprehensive subsurface geological information is to be obtained, three-dimensional seismic survey is required. And three-dimensional seismic survey requires higher costs and is not suitable for work areas where two-dimensional seismic survey data already exist.

[0088] For step 200, the surface velocity model describes the velocity distribution of seismic waves at different depths within and below the earth's crust. The propagation velocity of seismic waves is affected by geological media (such as rock type, temperature, pressure, etc.), so the velocity model can reflect the subsurface geological conditions and geological structure.

[0089] Specifically, the surface velocity model is established by analyzing seismic waveform data. The propagation velocity and path of seismic waves can be calculated through the travel time (the time from the seismic source to the seismic station) and the first arrival time of seismic waves. By analyzing a large amount of seismic waveform data, the velocity distribution of seismic waves at different depths can be constructed, that is, the surface velocity model.

[0090] Constructing a surface velocity model includes the following steps: Based on seismic data and experience, construct an initial velocity model. Use tomography algorithms to invert the seismic data with the initial velocity model to obtain a preliminary image. Evaluate and adjust the velocity model according to the preliminary image and single-shot first-arrival time data. According to the results of evaluation and adjustment, perform iterative inversion to continuously optimize the velocity model. The accuracy and reliability of the velocity model can be judged by comparing the fitting degree between the inversion result and the observed data. Finally, compare and verify the inverted velocity model with the actual geological situation. Other geological data and logging data can be used to verify the accuracy of the velocity model.

[0091] In some embodiments of the present invention, referring to Figure 2 , before step 100, it further includes:

[0092] Step 90: Preprocess all the survey lines.

[0093] In some embodiments of the present invention, referring to Figure 3 , step 90 includes:

[0094] Step 901: Determine whether the elevation of the shot points and the elevation of the geophone points of all the survey lines are closed;

[0095] The elevation of the shot point refers to the altitude or ground elevation of the position where the shot point (seismic source) is located in seismic exploration. It is the vertical reference height used when specifying the position of the shot point.

[0096] The elevation of the shot point is usually measured using surveying instruments such as total stations to measure the height of the shot point position on the surface. The accurate measurement of the elevation of the shot point is very important for seismic data processing and interpretation, and can help determine the propagation path and velocity of seismic waves underground.

[0097] The elevation of the geophone point refers to the altitude or ground elevation of the position where the receiver (geophone) is arranged in seismic exploration. It is the vertical reference height of the ground position for recording the underground signals reflected and refracted during the propagation of seismic waves.

[0098] The measurement method of the elevation of the geophone point is similar to that of the elevation of the shot point, and a total station or other surveying instruments can be used to measure the height of the geophone point position on the surface. The accurate measurement of the elevation of the geophone point can help correct the time shift in seismic data and improve the interpretation and imaging quality of seismic data.

[0099] The specific implementation method of step 901 is as follows:

[0100] For a complete survey line, if the shot point positions at the beginning and end are on the same horizontal plane, that is, the elevation or ground elevation is equal, then the shot point elevation can be considered closed. This means that from the starting point to the end point of the survey line, there is no obvious change in the elevation of the shot points along the survey line direction.

[0101] For a complete survey line, if the geophone point positions at the beginning and end are on the same horizontal plane, that is, the elevation or ground elevation is equal, then the geophone point elevation can be considered closed. This means that from the starting point to the end point of the survey line, there is no obvious change in the elevation of the geophone points along the survey line direction.

[0102] Closed shot point elevation and geophone point elevation can provide more accurate seismic data, help reduce the errors in time migration and depth conversion, and improve the reliability of seismic imaging and interpretation. Therefore, in actual survey line design and data processing, survey lines with closed elevation are usually selected as much as possible.

[0103] Step 902: If it is not closed, smooth the unclosed survey line.

[0104] Specifically, check whether the elevations of the shot and geophone points of each survey line in the work area are closed. If there are small non-closures, smoothing treatment should be done to ensure the accuracy of the measurement results at the same location in different periods.

[0105] When implementing Step 902, it specifically includes: First, for the shot point elevation or geophone point elevation data on the survey line, perform differential analysis to calculate the elevation difference between adjacent points. By analyzing the distribution and trend of the elevation differences, the areas with elevation changes can be determined. For the areas with large elevation changes, smoothing filtering techniques can be used to process the data. Smoothing filtering methods include moving average filtering, Gaussian filtering, etc. These filtering methods can remove some high-frequency noise and local mutations, making the data smoother.

[0106] For the areas with large elevation changes, interpolation methods can be used to fill the gaps between the data. Commonly used interpolation methods include linear interpolation, spline interpolation, etc. Interpolation methods can infer the elevation of unknown points through known data points, thereby achieving data smoothing.

[0107] If there is an obvious trend or pattern in the overall survey line data, global optimization methods can be considered for data smoothing. Global optimization methods can perform overall fitting and adjustment of the data through least squares fitting or other mathematical optimization algorithms, making the data smoother.

[0108] It should be noted that when smoothing the unclosed survey line data, the influence of data boundaries and endpoints needs to be carefully handled to avoid introducing unnecessary errors. In addition, during the data processing, parameter selection and result verification should be carried out according to the actual situation to ensure the effect and reliability of the smoothing process.

[0109] In some embodiments of the present invention, referring to Figure 4 , step 200 includes:

[0110] Step 201: Based on the first arrival time of a single shot, and according to the shot point coordinates, geophone coordinates, and surface elevation of all survey lines, construct the surface velocity model.

[0111] In some embodiments of the present invention, referring to Figure 4 , before step 300, the two-dimensional batch tomographic static correction method based on three-dimensional tomography further includes:

[0112] Step 290: Correct the surface velocity model according to the micro-logging data in the target work area.

[0113] Micro-logging is a logging method used to obtain the formation properties and wellbore characteristics around the wellbore. Usually, a logging tool with a smaller diameter is used, and high-resolution, local area measurements can be carried out in the wellbore. By monitoring the propagation and reflection of seismic waves, physical parameters such as the velocity, density, and wave impedance of the underground formation can be inferred.

[0114] During seismic micro-logging, first drill a hole of a certain depth from the ground underground. At certain intervals and quantities in this hole, then place an excitation (detonator or explosive) at preset positions, or place geophones for receiving the signals generated by the excitation at these preset positions. Also, place a certain number of geophones at a certain distance from the wellhead on the ground to receive, or excitation (if excitation is placed in the hole, then geophones are placed on the ground; if geophones are placed in the hole, then excitation is placed on the ground). Only the method of placing excitation in the hole and geophones on the ground will be described below. By detonating the excitation points at each depth from deep to shallow, the first arrival information of seismic waves from each excitation to each geophone (the geophones can be multiple or one. To improve efficiency and accuracy, usually multiple geophones are set simultaneously) is obtained. After vertical correction processing of the received first arrival information of seismic waves, draw the vertical time-depth curve of this seismic micro-logging. Using the fitting and interpretation of this time-depth curve, obtain the structural change amounts of the thickness and velocity of each low-velocity and velocity-reducing zone in the shallow surface layer. Finally, correct the surface velocity model according to the structural change amounts of the thickness and velocity of each low-velocity and velocity-reducing zone in the shallow surface layer.

[0115] It is understandable that single-line micro-logging is rare or non-existent, while there are relatively more in the entire work area. This realizes the change of micro-logging data application from point (line) to surface, thereby improving the accuracy of the initial velocity model.

[0116] In some embodiments of the present invention, referring to Figure 5 , step 200 includes:

[0117] Step 201: Determine the top interface of the three-dimensional high-velocity layer in the target work area according to the surface velocity model;

[0118] Different from the prior art, step 201 no longer constructs the top interface of the high-velocity layer in a single-line manner, but adopts an overall three-dimensional top interface of the high-velocity layer in the whole area (moderate smoothing processing is required according to actual conditions). On this basis, the tomographic static correction amount of the whole work area is calculated.

[0119] When implementing step 201, first perform tomographic imaging based on the surface velocity model and seismic data to obtain the structural information (velocity change and interface position) of the underground formation, namely velocity change and interface position. Then, according to the tomographic imaging results, perform inversion to infer the position and shape of the top interface of the high-velocity layer. Inversion methods can use methods such as least squares fitting and regularization. Verify and adjust the constructed top interface model of the high-velocity layer. Other geological data (such as borehole data, geological profiles, etc.) can be used for comparison and verification to ensure the accuracy and reliability of the model. Finally, according to the information of the top interface of the high-velocity layer obtained by inversion, combine it with the velocity model to construct a three-dimensional model of the top interface of the high-velocity layer.

[0120] Step 202: Determine the tomographic static correction amount according to the three-dimensional top interface of the high-velocity layer.

[0121] In some embodiments of the present invention, referring to Figure 6 , the two-dimensional batch tomographic static correction method based on three-dimensional tomography further includes before step 202:

[0122] Step 203: Smooth the three-dimensional top interface of the high-velocity layer.

[0123] When implementing step 203, it can be achieved through the following several ways:

[0124] Adopt a window with a fixed size, average the data within the window, and then use the average value as the new value of the center point of the window. This can reduce the influence of local noise and make the interface smoother.

[0125] Adopt a window with a fixed size, sort the data within the window, and take the median value as the new value of the center point of the window. This method has a good noise suppression effect and can retain boundary and detail information.

[0126] By performing weighted averaging on the data, the new value at the center point of the window is affected by neighboring points, while points farther away have less influence. This can smooth the interface and retain details.

[0127] Performing edge-preserving filtering on the top interface of the high-velocity layer can retain boundary information while smoothing. This processing method usually uses a gradient-based approach to perform different degrees of smoothing on edge points and non-edge points.

[0128] As can be seen from the above description, the embodiment of the present invention provides a two-dimensional batch tomography static correction method based on three-dimensional tomography, including: first, picking the first arrival times of all survey lines in the target work area through preset picking parameters to obtain the single-shot first arrival times of the target work area; where all survey lines are two-dimensional seismic survey lines; then, constructing the surface velocity model of the target work area based on the single-shot first arrival times; and finally determining the tomography static correction amount of the target work area according to the surface velocity model.

[0129] The present invention solves the closure problem of the survey line horizons (amplitude energy, time) in the two-dimensional work area, especially the closure problem in the depth domain. By merging the first arrival files of multiple lines in the whole area and inversely calculating the near-surface velocity model within the entire work area in a three-dimensional manner, the tomography static correction amount for the entire work area can be obtained, and the corresponding static correction amount is extracted for each survey line through the positions of shot and receiver points. The near-surface velocity model and amount obtained in this way are closed at any intersection position. At the same time, it provides a relatively accurate near-surface velocity model for the entire work area for subsequent three-dimensional construction, provides certain reference for near-surface investigation; and also provides a more accurate overall near-surface velocity for depth-domain velocity modeling.

[0130] Embodiment Three:

[0131] In a specific implementation manner, the present invention also provides a specific implementation manner of a two-dimensional batch tomography static correction method based on three-dimensional tomography, see Figure 7 , which specifically includes the following steps.

[0132] S1: Check whether the elevations of shot and receiver points of each survey line in the work area are closed;

[0133] If there are small non-closures, smoothing should be performed to ensure the accuracy of measurement results at the same position in different periods.

[0134] S2: Merge the data of all survey lines in the work area and pick the initial value time using unified parameters;

[0135] In this way, the single-shot first arrival times of the entire fused work area can be obtained. This file contains information such as the coordinates and surface elevations of all shot points and receiver points in the work area, as well as the first arrival time information of each trace.

[0136] S3: Perform a global coordinate transformation on the entire first arrival time file after multi-line fusion.

[0137] Affected by the complexity of near-surface excitation and reception conditions, there are often some anomalies in the first arrival time. Therefore, it is necessary to uniformly edit the first arrival time file. Specifically, step S3 includes the following aspects:

[0138] Shot point coordinate transformation: The shot point coordinate refers to the position coordinate where the seismic source is placed in seismic exploration. In some cases, it is necessary to convert the shot point coordinate from one coordinate system to another, such as from the geographic coordinate system (latitude and longitude) to the engineering coordinate system (such as the UTM coordinate system).

[0139] Geophone point coordinate transformation: The geophone point coordinate refers to the position coordinate where the seismic wave signal is received in seismic exploration. Similar to the shot point coordinate, the geophone point coordinate may also need to be converted between coordinate systems.

[0140] Rotation and translation of the coordinate system: In some cases, in order to register or integrate with other data, it may be necessary to rotate and translate the coordinate system of the seismic data. This usually requires adjusting the coordinates of the seismic data based on known control points or registration points.

[0141] Conversion of projection coordinate systems: Seismic exploration data is usually processed and analyzed using projection coordinate systems, such as the Universal Transverse Mercator projection (UTM) or the equal area projection (such as the Lambert projection). When converting between different projection coordinate systems, corresponding coordinate transformations are required.

[0142] S4: Construct an initial velocity model.

[0143] Figure 8 It is the near-surface velocity model for single-line inversion. There is no micro-logging information at this line, and it relies only on the inversion of the first arrival time. Figure 9 It is the initial near-surface velocity model constructed after fusing the first arrival time files of all survey lines in the work area. This velocity model can be regarded as the three-dimensional near-surface initial velocity model of the entire work area.

[0144] S5: Calibrate the surface velocity model according to the micro-logging data in the target work area.

[0145] There are very few or no single-line micro-logging wells, while there are relatively more in the entire work area. This realizes the change of applying micro-logging data from point (line) to surface, and improves the accuracy of the initial velocity model. After uniformly constructing the initial velocity model, tomography inversion work based on micro-logging constraints is adopted. A high-precision near-surface velocity model of the entire work area is obtained.

[0146] Figure 10 For Figure 9Based on this, a near-surface velocity model obtained by post-inversion with micro-logging constraints is adopted. This model can be used for the calculation of tomographic static correction for the entire area and the fusion velocity of the near-surface in the depth domain.

[0147] S6: Determine the tomographic static correction amount for the target work area based on the corrected surface velocity model, and perform tomographic static correction.

[0148] In step S6, instead of constructing the top interface of the high-velocity layer in a single-line manner, a holistic three-dimensional top interface of the high-velocity layer is adopted for the entire area (moderate smoothing processing is required according to the actual situation). Based on this, the tomographic static correction amount for the entire work area is calculated.

[0149] The tomographic static correction amount in step S6 includes the shot and receiver point information for the entire work area and can be applied to the entire work area at once or line by line.

[0150] Figure 11 For the closure situation of line A and line B at the added points after applying the layer static correction amount calculated by single-line inversion, it can be seen that there are obvious closure problems (the black frame area is the intersection position and both sides).

[0151] Figure 12 For the closure situation of line A and line B at the intersection after applying the correction amount obtained by unified inversion calculation of all two-dimensional lines in the work area using three-dimensional tomography technology, it can be seen that this method significantly eliminates the closure problems of two-dimensional lines (the black frame area is the intersection position and both sides).

[0152] Figure 13 For the overall three-dimensional near-surface model of the two-dimensional work area (the dotted line and above is the tomographic inversion model, which is smoothly fused through the top interface of the high-velocity layer), the near-surface velocity model is obtained through overall inversion combined with interpolation and extrapolation, effectively improving the accuracy of the near-surface velocity and providing a reliable initial velocity field for subsequent depth-domain imaging.

[0153] The following conclusions can be drawn from the above description:

[0154] Figure 8 For the near-surface velocity model obtained by single-line inversion, there is no micro-logging information at this line, and it relies only on the inversion of the first arrival time. The initial near-surface velocity model constructed by fusing the first arrival time files of all lines in the work area can be regarded as the three-dimensional near-surface initial velocity model of the entire work area, and the micro-logging data of the entire work area can be used, which has obvious advantages compared with single-line inversion. The near-surface velocity model obtained by post-inversion with micro-logging constraints can be used for the calculation of tomographic static correction for the entire area and the fusion velocity of the near-surface in the depth domain.

[0155] From Figure 11To test the closure of Line A and Line B at the added points after applying the layer static correction amount calculated by one - line inversion, it can be seen that there are obvious closure problems (the black - framed area is the intersection position and both sides). After applying the correction amount calculated by unified inversion of all 2D survey lines in the work area using 3D tomography technology, the closure of Line A and Line B at the intersection can be seen, and it can be seen that this method significantly eliminates the closure problems of 2D survey lines (the black - framed area is the intersection position and both sides).

[0156] The overall 3D near - surface model of the 2D work area (the yellow line and above is the tomography inversion model, smoothed and fused through the top interface of the high - velocity layer). By combining overall inversion with interpolation and extrapolation, a near - surface velocity model is obtained, effectively improving the accuracy of the near - surface velocity and providing a reliable initial velocity field for subsequent depth - domain imaging.

[0157] As can be seen from the above description, a specific application example of the present invention provides a 2D batch tomography static correction method based on 3D tomography, including: first, picking the first - arrival times of all survey lines in the target work area through preset picking parameters to obtain the single - shot first - arrival times of the target work area; where all survey lines are 2D seismic survey lines; then, constructing the surface velocity model of the target work area based on the single - shot first - arrival times; and finally, determining the tomography static correction amount of the target work area according to the surface velocity model.

[0158] The present invention realizes the unified inversion of the near - surface velocity model and the calculation of the tomography static correction amount for 2D seismic data processing at one time. It effectively solves the closure problems of 2D seismic data (time domain) in terms of the in - phase axis time and amplitude energy; provides a relatively accurate near - surface velocity model for the entire work area for subsequent 3D construction, provides certain reference for near - surface investigation; and provides a more accurate overall near - surface velocity for depth - domain velocity modeling.

[0159] A 2D batch tomography static correction method based on 3D tomography provided by a specific application example of the present invention is used for near - surface velocity modeling in a complex 2D work area, thereby obtaining an accurate near - surface velocity model, solving the closure problems of the survey line horizons (amplitude energy, time) in the 2D work area, especially the closure problems in the depth domain. By merging the first - arrival files of multiple lines in the whole area, editing them uniformly in a 3D manner, and using appropriate offset segments to invert the near - surface velocity model within the entire work area, the tomography static correction amount for the entire work area is obtained, and the corresponding static correction amount is extracted for each survey line through the positions of shot - build points. The near - surface velocity model and amount obtained in this way are closed at any intersection position.

[0160] On the other hand, the present invention fills the gap in two-dimensional multi-line processing based on three-dimensional tomography. The issues of closure and static correction accuracy in two-dimensional seismic work areas have always been a very important part of two-dimensional seismic exploration. This technical solution should be able to solve this problem well from the source; it provides a relatively accurate near-surface velocity model for the entire work area for subsequent three-dimensional construction, provides certain reference for near-surface surveys; provides a more accurate overall near-surface velocity for subsequent depth-domain velocity modeling; and has a wide range of application prospects.

[0161] Embodiment 4:

[0162] Based on the same inventive concept, an embodiment of the present application also provides a two-dimensional batch tomography static correction device based on three-dimensional tomography, which can be used to implement the method described in the above embodiment, as in the following embodiment. Since the principle of the two-dimensional batch tomography static correction device based on three-dimensional tomography to solve problems is similar to that of the two-dimensional batch tomography static correction method based on three-dimensional tomography, the implementation of the two-dimensional batch tomography static correction device based on three-dimensional tomography can refer to the implementation of the two-dimensional batch tomography static correction method based on three-dimensional tomography, and the repeated parts will not be elaborated. As used below, the term "unit" or "module" can be a combination of software and / or hardware that can achieve a predetermined function. Although the system described in the following embodiments is preferably implemented in software, implementation in hardware, or a combination of software and hardware is also possible and contemplated.

[0163] An embodiment of the present invention provides a specific implementation manner of a two-dimensional batch tomography static correction device based on three-dimensional tomography that can implement the two-dimensional batch tomography static correction method based on three-dimensional tomography. Refer to Figure 14 , a two-dimensional batch tomography static correction device based on three-dimensional tomography includes:

[0164] The first arrival time picking module 10 is used to pick the first arrival time of all survey lines in the target work area through preset picking parameters to obtain the single-shot first arrival time of the target work area; wherein, all the survey lines are two-dimensional seismic survey lines;

[0165] The surface velocity model construction module 20 is used to construct the surface velocity model of the target work area based on the single-shot first arrival time;

[0166] The tomography static correction amount determination module 30 is used to determine the tomography static correction amount of the target work area according to the surface velocity model.

[0167] In an embodiment of the present invention, a two-dimensional batch tomography static correction device based on three-dimensional tomography further includes:

[0168] The survey line preprocessing module is used to preprocess all the survey lines.

[0169] In an embodiment of the present invention, the survey line preprocessing module includes:

[0170] A closed judgment unit, configured to judge whether the shot elevations and geophone elevations of all survey lines are closed;

[0171] A survey line smoothing unit, configured to perform smoothing processing on the unclosed survey lines if they are not closed.

[0172] In an embodiment of the present invention, the surface velocity model construction module includes:

[0173] A surface velocity model construction unit, configured to construct the surface velocity model based on the first arrival time of a single shot and according to the shot coordinates, geophone coordinates, and surface elevation of all survey lines.

[0174] In an embodiment of the present invention, a two-dimensional batch tomography static correction device based on three-dimensional tomography further includes:

[0175] A velocity model correction module, configured to correct the surface velocity model according to the micro-logging data in the target work area.

[0176] In an embodiment of the present invention, the tomography static correction amount determination module includes:

[0177] A high-velocity layer top interface determination unit, configured to determine the three-dimensional high-velocity layer top interface of the target work area according to the surface velocity model;

[0178] A tomography static correction amount determination unit, configured to determine the tomography static correction amount according to the three-dimensional high-velocity layer top interface.

[0179] In an embodiment of the present invention, a two-dimensional batch tomography static correction device based on three-dimensional tomography further includes:

[0180] A high-velocity layer top interface smoothing module, configured to perform smoothing processing on the three-dimensional high-velocity layer top interface.

[0181] As can be seen from the above description, an embodiment of the present invention provides a two-dimensional batch tomography static correction device based on three-dimensional tomography, including: a three-parameter determination module, configured to respectively determine the energy intensity similarity, energy distribution contrast, and structural similarity between a target imaging profile and at least one reference imaging profile; wherein, the target imaging profile and the reference imaging profile are generated from seismic data measured by different observation systems; a target imaging profile evaluation module, configured to evaluate the target imaging profile according to the energy intensity similarity, energy distribution contrast, structural similarity, and at least one reference imaging profile.

[0182] The two-dimensional batch tomography static correction device based on three-dimensional tomography provided by the embodiments of the present invention provides a reliable quantitative index to analyze the change trend of the imaging results of different observation systems, and solves the problem that the subtle differences in the imaging results cannot be visually distinguished.

[0183] Embodiment 5:

[0184] The embodiments of the present application also provide a specific implementation manner of an electronic device that can implement all the steps in the two-dimensional batch tomography static correction method based on three-dimensional tomography in the above embodiments. Refer to Figure 15 , and the electronic device specifically includes the following content:

[0185] A processor 1201, a memory 1202, a communication interface 1203, and a bus 1204;

[0186] Among them, the processor 1201, the memory 1202, and the communication interface 1203 complete mutual communication through the bus 1204; the communication interface 1203 is used to implement information transmission between related devices such as server-side devices and client-side devices;

[0187] The processor 1201 is used to call the computer program in the memory 1202. When the processor executes the computer program, all the steps in the two-dimensional batch tomography static correction method in the above embodiments are implemented. For example, when the processor executes the computer program, the following steps are implemented:

[0188] Pick up the first arrival times of all survey lines in the target work area through preset picking parameters to obtain the single-shot first arrival times of the target work area; wherein, all the survey lines are two-dimensional seismic survey lines;

[0189] Construct the surface velocity model of the target work area based on the single-shot first arrival times;

[0190] Determine the tomography static correction amount of the target work area according to the surface velocity model.

[0191] In one embodiment, before picking up the first arrival times of all survey lines in the target work area through preset picking parameters to obtain the single-shot first arrival times of the target work area, it further includes:

[0192] Preprocess all the survey lines.

[0193] In one embodiment, preprocessing all the survey lines includes:

[0194] Judge whether the elevation of the shot points and the elevation of the geophone points of all the survey lines are closed;

[0195] If it is not closed, smooth the unclosed survey lines.

[0196] In one embodiment, constructing the surface velocity model of the target work area based on the first arrival time of a single shot includes:

[0197] Based on the first arrival time of a single shot, and according to the shot point coordinates, geophone point coordinates, and surface elevation of all the survey lines, construct the surface velocity model.

[0198] In one embodiment, before determining the tomographic static correction amount of the target work area according to the surface velocity model, it further includes:

[0199] Correct the surface velocity model according to the micro-log data in the target work area.

[0200] In one embodiment, determining the tomographic static correction amount of the target work area according to the surface velocity model includes:

[0201] Determine the three-dimensional high-velocity layer top interface of the target work area according to the surface velocity model;

[0202] Determine the tomographic static correction amount according to the three-dimensional high-velocity layer top interface.

[0203] In one embodiment, before determining the tomographic static correction amount according to the three-dimensional high-velocity layer top interface, it further includes:

[0204] Smooth the three-dimensional high-velocity layer top interface.

[0205] Embodiment Six:

[0206] An embodiment of the present application further provides a computer-readable storage medium capable of implementing all steps in the above-mentioned two-dimensional batch tomographic static correction method based on three-dimensional tomography. A computer program is stored on the computer-readable storage medium. When the computer program is executed by a processor, all steps of the two-dimensional batch tomographic static correction method based on three-dimensional tomography in the above-mentioned embodiment are implemented. For example, when the processor executes the computer program, the following steps are implemented:

[0207] Pick up the first arrival time of all survey lines in the target work area through preset picking parameters to obtain the first arrival time of a single shot in the target work area; wherein, all the survey lines are two-dimensional seismic survey lines;

[0208] Construct the surface velocity model of the target work area based on the first arrival time of a single shot;

[0209] Determine the tomographic static correction amount of the target work area according to the surface velocity model.

[0210] In one embodiment, before picking up the first arrival times of all survey lines in the target work area through preset picking parameters to obtain the single-shot first arrival times of the target work area, the following steps are further included:

[0211] Preprocess all the survey lines.

[0212] In one embodiment, preprocessing all the survey lines includes:

[0213] Judge whether the elevations of the shot points and the geophone points of all the survey lines are closed;

[0214] If not closed, perform smoothing processing on the unclosed survey lines.

[0215] In one embodiment, constructing the surface velocity model of the target work area based on the single-shot first arrival times includes:

[0216] Based on the single-shot first arrival times, and according to the shot point coordinates, geophone point coordinates, and surface elevation of all the survey lines, construct the surface velocity model.

[0217] In one embodiment, before determining the tomographic static correction amount of the target work area according to the surface velocity model, the following steps are further included:

[0218] Correct the surface velocity model according to the micro-log data in the target work area.

[0219] In one embodiment, determining the tomographic static correction amount of the target work area according to the surface velocity model includes:

[0220] Determine the three-dimensional high-velocity layer top interface of the target work area according to the surface velocity model;

[0221] Determine the tomographic static correction amount according to the three-dimensional high-velocity layer top interface.

[0222] In one embodiment, before determining the tomographic static correction amount according to the three-dimensional high-velocity layer top interface, the following steps are further included:

[0223] Perform smoothing processing on the three-dimensional high-velocity layer top interface.

[0224] Each embodiment in this specification is described in a progressive manner. For the same or similar parts among the embodiments, reference can be made to each other. Each embodiment focuses on the differences from other embodiments. In particular, for the hardware + program type embodiments, since they are basically similar to the method embodiments, the description is relatively simple, and the relevant parts can refer to the partial description of the method embodiments.

[0225] The above describes specific embodiments of this specification. Other embodiments are within the scope of the appended claims. In some cases, the acts or steps recited in the claims may be performed in a different order than in the embodiments and still achieve the desired result. Additionally, the processes depicted in the figures do not necessarily require the particular order or sequential order shown to achieve the desired result. In certain embodiments, multitasking and parallel processing are also possible or may be advantageous.

[0226] Although this application provides method operation steps such as in the embodiments or flowcharts, based on routine or non-creative labor, there may be more or fewer operation steps. The order of steps listed in the embodiments is only one way among many execution orders of the steps and does not represent the only execution order. When the actual device or client product is executed, it can be executed in the order shown in the embodiments or the figures or in parallel (e.g., in an environment with parallel processors or multi-threaded processing).

[0227] For convenience of description, when describing the above device, it is divided into various modules according to functions for separate description. Of course, when implementing the embodiments of this specification, the functions of each module can be implemented in the same or multiple software and / or hardware, or the modules implementing the same function can be realized by a combination of multiple sub-modules or sub-units, etc. The device embodiments described above are merely illustrative. For example, the division of the units is only a logical function 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. Another point is that the displayed or discussed coupling or direct coupling or communication connection to each other can be through some interfaces, and the indirect coupling or communication connection of the device or unit can be in electrical, mechanical or other forms.

[0228] Those skilled in the art also know that in addition to implementing the controller in the form of pure computer-readable program code, the method steps can be logically programmed to enable the controller to be implemented in the form of logic gates, switches, application-specific integrated circuits, programmable logic controllers, embedded microcontrollers, etc. to achieve the same function. Therefore, such a controller can be considered a hardware component, and the devices included therein for implementing various functions can also be regarded as the structures within the hardware component. Or even, the devices for implementing various functions can be regarded as both software modules for implementing the method and structures within the hardware component.

[0229] In a typical configuration, a computing device includes one or more processors (CPUs), an input / output interface, a network interface, and memory.

[0230] The memory may include non-permanent memory in the form of computer-readable media, such as random access memory (RAM) and / or non-volatile memory, such as read-only memory (ROM) or flash RAM. The memory is an example of computer-readable media.

[0231] Each embodiment in this specification is described in a progressive manner. For the same or similar parts among the embodiments, reference can be made to each other. Each embodiment focuses on the differences from other embodiments. In particular, for the system embodiment, since it is basically similar to the method embodiment, the description is relatively simple. For the relevant parts, reference can be made to the description of the method embodiment. In the description of this specification, the description with reference to the terms "one embodiment", "some embodiments", "example", "specific example", or "some examples", etc. means that the specific features, structures, materials, or characteristics described in connection with the embodiment or example are included in at least one embodiment or example of the embodiments of this specification. In this specification, the schematic expressions of the above terms do not necessarily refer to the same embodiment or example. Moreover, the specific features, structures, materials, or characteristics described can be combined in a suitable manner in any one or more embodiments or examples. In addition, without contradiction, those skilled in the art can combine and combine the different embodiments or examples described in this specification and the features of different embodiments or examples.

[0232] The above is only for the embodiments of this specification and is not used to limit the embodiments of this specification. For those skilled in the art, various changes and modifications can be made to the embodiments of this specification. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the embodiments of this specification shall be included within the scope of the claims of the embodiments of this specification.

Claims

1. A two-dimensional batch tomography static correction method based on three-dimensional tomography, characterized in that, it includes: Pick up the first arrival time of all survey lines in the target work area through preset picking parameters to obtain the single-shot first arrival time of the target work area; wherein, all the survey lines are two-dimensional seismic survey lines; Construct the surface velocity model of the target work area based on the single-shot first arrival time; Determine the tomography static correction amount of the target work area according to the surface velocity model.

2. The two-dimensional batch tomography static correction method according to claim 1, characterized in that, Before picking up the first arrival time of all survey lines in the target work area through preset picking parameters to obtain the single-shot first arrival time of the target work area, it further includes: Preprocess all the survey lines.

3. The two-dimensional batch tomography static correction method according to claim 2, characterized in that, Preprocessing all the survey lines includes: Judge whether the elevation of the shot points and the elevation of the geophone points of all the survey lines are closed; If not closed, smooth the unclosed survey lines.

4. The two-dimensional batch tomography static correction method according to claim 1, characterized in that, Constructing the surface velocity model of the target work area based on the single-shot first arrival time includes: Construct the surface velocity model based on the single-shot first arrival time and according to the shot point coordinates, geophone point coordinates and surface elevation of all the survey lines.

5. The two-dimensional batch tomography static correction method according to claim 1, characterized in that, Before determining the tomography static correction amount of the target work area according to the surface velocity model, it further includes: Correct the surface velocity model according to the micro-log data in the target work area.

6. The two-dimensional batch tomography static correction method according to any one of claims 1 to 5, characterized in that, Determining the tomography static correction amount of the target work area according to the surface velocity model includes: Determine the top interface of the three-dimensional high-velocity layer of the target work area according to the surface velocity model; Determine the tomography static correction amount according to the top interface of the three-dimensional high-velocity layer.

7. The two-dimensional batch tomography static correction method according to claim 6, characterized in that, Before determining the tomography static correction amount according to the top interface of the three-dimensional high-velocity layer, it further includes: Smooth the top interface of the three-dimensional high-velocity layer.

8. A two-dimensional batch tomography static correction device based on three-dimensional tomography, characterized in that, it includes: The first arrival time picking module is used to pick up the first arrival time of all survey lines in the target work area through preset picking parameters to obtain the single-shot first arrival time of the target work area; wherein, all the survey lines are two-dimensional seismic survey lines; The surface velocity model construction module is used to construct the surface velocity model of the target work area based on the single-shot first arrival time; The tomography static correction amount determination module is used to determine the tomography static correction amount of the target work area according to the surface velocity model.

9. An electronic device, including 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 steps of a two-dimensional batch tomography static correction method based on three-dimensional tomography according to any one of claims 1 to 7.

10. A computer-readable storage medium, on which a computer program is stored, characterized in that, when the computer program is executed by a processor, it implements the steps of a two-dimensional batch tomography static correction method based on three-dimensional tomography according to any one of claims 1 to 7.

Images