Quasi-three-dimensional loess tableland complex earth surface multi-line two-dimensional seismic data static correction and velocity analysis method and device, electronic equipment and medium
By using a pseudo-3D processing method, a 3D mesh was established for seismic data processing, which solved the accuracy problems of static correction and velocity analysis in 2D seismic data of the complex surface of the Loess Plateau, and improved the imaging accuracy and interpretation reliability of the seismic data.
Patent Information
- Authority / Receiving Office
- CN · China
- Patent Type
- Applications(China)
- Current Assignee / Owner
- Filing Date
- 2024-09-23
- Publication Date
- 2026-03-24
AI Technical Summary
In the static correction and velocity analysis of the complex surface of the Loess Plateau, 2D seismic data suffers from problems such as difficulty in closing near-surface velocities and inconsistent velocity points, leading to misinterpretations and affecting the accuracy and precision of tectonic interpretation.
A quasi-3D processing method was adopted to establish a 3D grid covering all 2D survey lines in the entire area. Seismic data was loaded, and first arrival picking and noise attenuation were performed. Stacking velocity analysis and residual static correction were carried out line by line to ensure consistent velocity at the intersection of survey lines. Imaging accuracy was improved by pre-stack time migration and migration velocity analysis.
It improves the imaging accuracy and interpretability of seismic data, ensures velocity consistency at the intersection of survey lines, reduces the risk of misjudgment, and provides more reliable exploration data.
Smart Images
Figure CN121721728A_ABST
Abstract
Description
TECHNICAL FIELD
[0001] The present application relates to the field of oil and gas geophysical exploration, and more particularly, to a quasi-three-dimensional loess tableland complex surface multi-line two-dimensional seismic data static correction and velocity analysis method, device, electronic equipment and medium. BACKGROUND
[0002] The two-dimensional survey line single line is closed by the static correction amount, and on this basis, the static correction and velocity analysis are performed, and the near-surface and deep velocity are often difficult to close, which easily leads to difficulties in structural interpretation closure; and the non-uniform velocity points easily lead to misjudgment of the interpreter and cause difficulties for variable velocity mapping.
[0003] The loess tableland surface is complex, the near-surface velocity changes greatly, and the signal-to-noise ratio of seismic data is low, and in view of these problems, a quasi-three-dimensional loess tableland complex surface multi-line two-dimensional seismic data static correction and velocity analysis method is formed through research, which leads to spatial discontinuity of seismic data caused by two-dimensional seismic data, and the spatial discontinuity leads to the problems of difficult surveying and seismic interpretation layer closure and difficult fault closure in multi-line seismic data interpretation.
[0004] At present, the two-dimensional seismic data processing generally performs the closure of the static correction amount in the static correction work by the artificial method, and the artificial comparison method is used for velocity analysis to partially solve the closure problem, but this method often faces difficulties in implementation, the static correction amount reflects the static correction stripping amount of the shot point and the geophone position, and cannot reflect the common depth point and the reflection point static correction at the intersection of the two-dimensional survey line, in this case, the separation of the high-frequency static correction amount and the low-frequency static correction amount is generally performed, the static correction amount separated out at the common depth point is matched, the near-surface velocity model has not been truly closed, the accuracy is low, and many difficulties are brought to the later velocity analysis.
[0005] Therefore, it is necessary to develop a quasi-three-dimensional loess tableland complex surface multi-line two-dimensional seismic data static correction and velocity analysis method, device, electronic equipment and medium.
[0006] The information disclosed in the background section of this application is only intended to deepen the understanding of the general background of the application, and should not be regarded as acknowledging or implying in any form that the information constitutes prior art known to those skilled in the art. SUMMARY
[0007] The present application proposes a quasi-three-dimensional loess tableland complex surface multi-line two-dimensional seismic data static correction and velocity analysis method, device, electronic equipment and medium, which can improve the imaging accuracy of seismic data, increase the interpretability of data, and provide more reliable data for exploration.
[0008] In a first aspect, the present application provides a quasi-three-dimensional loess tableland complex surface multi-line two-dimensional seismic data static correction and velocity analysis method, comprising:
[0009] establishing a three-dimensional grid containing all two-dimensional survey line shot points in the whole area and loading seismic data;
[0010] picking up first arrivals of seismic data in the whole area and calculating tomographic static correction amounts;
[0011] inverting a near-surface velocity model according to the first arrivals based on the three-dimensional grid after loading seismic data;
[0012] preprocessing noise attenuation data, performing stacking velocity analysis line by line, and obtaining a velocity spectrum;
[0013] performing velocity analysis and residual static correction on the velocity spectrum, performing pre-stack time migration and migration velocity analysis line by line;
[0014] checking survey line intersections and outputting migration results.
[0015] As a specific implementation manner of the embodiment of the present disclosure, the method further comprises:
[0016] After loading seismic data, it is checked whether the seismic data exceeds the grid range.
[0017] As a specific implementation manner of the embodiment of the present disclosure, the noise attenuation data comprises surface waves, linear interference, abnormal amplitude and 50Hz industrial electric noise.
[0018] As a specific implementation manner of the embodiment of the present disclosure, the preprocessing comprises attenuation processing, surface consistency amplitude compensation and surface consistency deconvolution.
[0019] As a specific implementation manner of the embodiment of the present disclosure, the checking of the survey line intersections comprises:
[0020] checking whether there is a time difference of the same set of strata at the intersections of different survey lines.
[0021] In the second aspect, the embodiment of the present disclosure also provides a quasi-three-dimensional loess tableland complex surface multi-line two-dimensional seismic data static correction and velocity analysis device, comprising:
[0022] a grid establishing module, which establishes a three-dimensional grid containing all two-dimensional survey line shot points in the whole area and loads seismic data;
[0023] a first arrival picking module, which picks up first arrivals of seismic data in the whole area and calculates tomographic static correction amounts;
[0024] an inversion module, which inverts a near-surface velocity model according to the first arrivals based on the three-dimensional grid after loading seismic data;
[0025] a preprocessing module, which preprocesses noise attenuation data, performs stacking velocity analysis line by line, and obtains a velocity spectrum;
[0026] The processing module performs velocity analysis and residual static correction on the velocity spectrum, and performs pre-stack time migration, migration velocity analysis line by line;
[0027] The output module checks the intersection of the survey lines and outputs the migration result.
[0028] As a specific implementation manner of the embodiment of the present disclosure, the method further comprises:
[0029] After loading the seismic data, it is checked whether the seismic data exceeds the grid range.
[0030] As a specific implementation manner of the embodiment of the present disclosure, the noise attenuation data comprises surface waves, linear interference, abnormal amplitude and 50Hz industrial electric noise.
[0031] As a specific implementation manner of the embodiment of the present disclosure, the preprocessing comprises attenuation processing, surface consistency amplitude compensation and surface consistency deconvolution.
[0032] As a specific implementation manner of the embodiment of the present disclosure, the checking of the intersection of the survey lines comprises:
[0033] Checking whether there is a time difference in the same set of strata at the intersection of different survey lines.
[0034] In a third aspect, the embodiment of the present disclosure further provides an electronic device, which comprises:
[0035] A memory storing executable instructions;
[0036] A processor running the executable instructions in the memory to implement the quasi-three-dimensional loess table complex surface multi-line two-dimensional seismic data static correction and velocity analysis method.
[0037] In a fourth aspect, the embodiment of the present disclosure further provides a computer readable storage medium storing a computer program, which is executed by a processor to implement the quasi-three-dimensional loess table complex surface multi-line two-dimensional seismic data static correction and velocity analysis method.
[0038] The method and device of the present application have other characteristics and advantages, which will be apparent or will be described in detail in the accompanying drawings and subsequent specific embodiments incorporated herein, which together serve to explain the specific principles of the present application. BRIEF DESCRIPTION OF DRAWINGS
[0039] The above and other objects, features and advantages of the present application will become more apparent from the following detailed description when taken in conjunction with the accompanying drawings in which like reference characters refer to like parts throughout the figures, and in which:
[0040] Figure 1 A flow chart showing the steps of a pseudo-three-dimensional loess tableland complex surface multi-line two-dimensional seismic data static correction and velocity analysis method according to an embodiment of the present application is shown.
[0041] Figure 2 A schematic diagram showing a grid design according to an embodiment of the present application is shown.
[0042] Figure 3 A schematic diagram showing first arrival picking according to an embodiment of the present application is shown.
[0043] Figure 4a , Figure 4b A comparative schematic diagram showing a near-surface inversion model and a near-surface data formation model according to an embodiment of the present application is shown.
[0044] Figure 5a , Figure 5b A schematic diagram showing single-shot effects before and after static correction according to an embodiment of the present application is shown.
[0045] Figure 6a , Figure 6b A schematic diagram showing profile comparison of elevation static correction and tomographic static correction according to an embodiment of the present application is shown.
[0046] Figure 7 A schematic diagram showing stack velocity analysis according to an embodiment of the present application is shown.
[0047] Figure 8 A migration velocity analysis point distribution diagram according to an embodiment of the present application is shown.
[0048] Figure 9a , Figure 9b A schematic diagram showing migration velocity closure analysis diagrams of east-west lines and north-south lines according to an embodiment of the present application is shown.
[0049] Figure 10 A schematic diagram showing line closure analysis according to an embodiment of the present application is shown.
[0050] Figure 11 A block diagram of a pseudo-three-dimensional loess tableland complex surface multi-line two-dimensional seismic data static correction and velocity analysis device according to an embodiment of the present application is shown.
[0051] BRIEF DESCRIPTION OF THE DRAWINGS
[0052] 201, a grid establishing module; 202, a first arrival picking module; 203, an inversion module; 204, a preprocessing module; 205, a processing module; 206, an output module. DETAILED DESCRIPTION
[0053] The preferred embodiments of the present application will be described in more detail below. Although the preferred embodiments of the present application are described below, it is understood that the present application can be implemented in various forms and should not be limited by the embodiments set forth herein.
[0054] To make the scheme and effects of the embodiments of the present application more comprehensible, six specific application examples are given below. It should be understood by those skilled in the art that the examples are only for the convenience of understanding the present application, and any specific details thereof are not intended to limit the present application in any way.
[0055] Example 1
[0056] Figure 1 A flowchart showing the steps of a pseudo-three-dimensional loess tableland complex surface multi-line two-dimensional seismic data static correction and velocity analysis method according to an embodiment of the present application is shown.
[0057] As shown in Figure 1 The pseudo-three-dimensional loess tableland complex surface multi-line two-dimensional seismic data static correction and velocity analysis method includes:
[0058] Step 101, a three-dimensional grid containing all two-dimensional line shot points in the whole area is established and seismic data is loaded;
[0059] Step 102, first arrival picking is performed on seismic data in the whole area, and tomographic static correction is calculated;
[0060] Step 103, based on the three-dimensional grid after loading the seismic data, a near-surface velocity model is inverted according to the first arrival wave;
[0061] Step 104, noise attenuation data is preprocessed, and stacking velocity analysis is performed line by line to obtain a velocity spectrum;
[0062] Step 105, velocity analysis and residual static correction are performed on the velocity spectrum, and pre-stack time migration and migration velocity analysis are performed line by line;
[0063] Step 106, the intersection points of the lines are checked, and the migration results are output.
[0064] In one example, it further includes:
[0065] After loading the seismic data, it is checked whether the seismic data exceeds the grid range.
[0066] In one example, the noise attenuation data includes: surface wave, linear interference, abnormal amplitude, 50Hz industrial electric noise.
[0067] In one example, the pre-processing includes attenuation processing, surface consistency amplitude compensation, surface consistency deconvolution.
[0068] In one example, the inspection for the intersection of the survey lines includes:
[0069] The inspection is performed on whether the same set of strata at the intersection of different survey lines exists time difference.
[0070] Specifically, after collecting data such as an observation system file, data, and shallow seismic data, and loading the observation system, the degree of dispersion of shot-receiver points is analyzed.
[0071] According to the observation system, a three-dimensional processing grid is established. In seismic data processing, in the face of two-dimensional acquisition seismic data, the common method is to establish a processing grid for each two-dimensional survey line separately, and to perform a series of subsequent processing on a single line. However, there is one or more intersections between the east-west and north-south two-dimensional survey lines in a work area. The static correction amount is calculated by establishing a grid for a single line, and the static correction amount of different two-dimensional survey lines at the intersection may be different, which causes time difference problems of the same set of strata on different survey lines and affects subsequent interpretation.
[0072] The application adopts pseudo-three-dimensional processing, that is, a large grid covering all two-dimensional survey lines in the work area is established, and all shot-receiver point information in the work area is used to calculate a unified static correction amount for the whole work area, so as to avoid the possible time difference problem at the intersection of survey lines. Specifically, the shot-receiver point coordinates and common midpoint are distributed in a plane, and a processing grid is designed. The bin grid size should be consistent with the design bin of the two-dimensional survey line.
[0073] Observation system and grid loading: the data is loaded with the processing grid, and whether the seismic data exceeds the grid range after the quality control of the loaded grid is checked.
[0074] The first arrival is picked up for the whole area data, and is used to calculate the tomographic static correction amount. The near-surface model is inverted according to the first arrival. The near-surface velocity model is quality controlled to check whether there is an abnormal value and whether it is consistent with the geological rule. The static correction amount is quality controlled to analyze whether the static correction amount is reasonable and applied to the seismic data.
[0075] The noise such as surface wave, linear interference, abnormal amplitude, and 50Hz industrial electric noise is attenuated to improve the signal-to-noise ratio of the data, and the surface consistency amplitude compensation, surface consistency deconvolution and other data preprocessing work are completed.
[0076] After the velocity line is selected according to the distribution of the two-dimensional survey lines, the above pseudo-three-dimensional static correction processing and the preprocessing work of improving the signal-to-noise ratio are completed, and the stacking velocity analysis is performed line by line.
[0077] Through multiple iterations of velocity analysis and residual static correction for surface consistency, the accuracy of velocity was improved. Simultaneously, during line-by-line velocity analysis, emphasis was placed on the velocity at intersections between survey lines to ensure velocity consistency at these points. Pre-stack migration velocity analysis was extracted from 2D lines, and pre-stack time migration and migration velocity analysis were performed line-by-line. The presence of time differences within the same stratigraphic set at intersections of different survey lines was checked to confirm the effectiveness of the pseudo-3D processing. The migration results were output for comprehensive interpretation.
[0078] Example 2
[0079] This invention also provides a device for static correction and velocity analysis of multi-line two-dimensional seismic data of complex loess plateau surfaces (pseudo-three-dimensional), comprising:
[0080] The grid creation module creates a 3D grid containing all 2D seismic line shot points in the entire area and loads seismic data.
[0081] The first arrival picking module performs first arrival picking on the seismic data for the entire region and calculates the tomographic static correction.
[0082] The inversion module, based on a 3D grid after loading seismic data, inverts the near-surface velocity model according to the first arrival wave;
[0083] The preprocessing module preprocesses the noise attenuation data and performs line-by-line velocity analysis to obtain the velocity spectrum.
[0084] The processing module performs velocity analysis and residual static correction on the velocity spectrum, and performs pre-stack time migration and migration velocity analysis line by line.
[0085] The output module checks the intersection points of the survey lines and outputs the offset results.
[0086] In one example, it also includes:
[0087] After loading the seismic data, check whether the quality control seismic data exceeds the grid range.
[0088] In one example, noise attenuation data includes: surface waves, linear interference, anomalous amplitude, and 50Hz industrial electrical noise.
[0089] In one example, preprocessing includes attenuation processing, surface consistency amplitude compensation, and surface consistency deconvolution.
[0090] In one example, checking the intersection of survey lines includes:
[0091] Check whether there is a time difference at the intersection of different survey lines for the same set of strata.
[0092] Specifically, data such as observation system files, data, and shallow seismic data are collected, and after the observation system is loaded, the dispersion of shot-receiver points is analyzed.
[0093] According to the observation system, a three-dimensional processing grid is established. In seismic data processing, in the face of two-dimensional acquisition seismic data, the common method is to establish a processing grid for each two-dimensional line separately, and then a series of subsequent processing is carried out for a single line. However, there is one or more intersection points between the east-west and north-south two-dimensional lines in a work area. If a single line is used to establish a grid and calculate a static correction amount, the static correction amounts of different two-dimensional lines at the intersection points may be different, which causes the existence of time difference of the same set of strata on different lines, and affects the subsequent interpretation.
[0094] The application adopts pseudo-three-dimensional processing, that is, a large grid covering all two-dimensional lines in the work area is established, and all shot point information in the work area is used to calculate a unified static correction amount for the whole work area, so as to avoid the possible time difference problem at the intersection points of the lines. Specifically, the shot point coordinates and common midpoint are distributed in a plane, a processing grid is designed, and the bin grid size should be consistent with the design of the two-dimensional line bin.
[0095] Selecting an observation system and loading a grid: load the processing grid to the data, and check whether the seismic data loaded after the quality control of the grid exceeds the grid range.
[0096] First arrival picking is performed on the whole area data to calculate the tomographic static correction amount. The near-surface model is inverted according to the first arrival wave. The near-surface velocity model is quality controlled to check whether there are abnormal values and whether it conforms to the geological law. The static correction amount is quality controlled to analyze whether the static correction amount is reasonable and applied to the seismic data.
[0097] Attenuate the noises such as surface wave, linear interference, abnormal amplitude, 50Hz industrial electricity and the like to improve the signal-to-noise ratio of the data, and complete the surface consistency amplitude compensation, surface consistency deconvolution and other data preprocessing work.
[0098] After selecting the velocity line according to the distribution of the two-dimensional line and completing the above-mentioned pseudo-three-dimensional static correction processing and improving the signal-to-noise ratio and other preprocessing work, the stacking velocity analysis is performed line by line.
[0099] Through the velocity analysis and the multiple iterations of the surface consistency residual static correction, the accuracy of the velocity is improved. At the same time, attention is paid to the velocity at the intersection of the lines during the line-by-line velocity analysis to ensure the consistency of the velocity at the intersection. The two-dimensional line prestack migration velocity analysis is extracted, and the prestack time migration and migration velocity analysis are performed line by line. The existence of time difference of the same set of strata at the intersection points of different lines is checked to confirm the effect of the pseudo-three-dimensional processing. The migration result is output for comprehensive interpretation.
[0100] Example 3
[0101] Collect the observation system file, data, shallow seismic data and the like, load the observation system, and analyze the dispersion degree of the shot point.
[0102] Figure 2 A schematic diagram of a grid design according to an embodiment of the present invention is shown, with red indicating the design grid range.
[0103] A three-dimensional processing grid is established based on the observation system. In a certain work area, the observation system is 3590-10-(20)-10-3590, designed as a variable surface element, requiring processing with 20m surface elements. Therefore, the design grid for the entire area is 20m x 20m. The center point of the grid should be consistent with the center point of the existing design surface element as much as possible. All shot and receiver coordinates should be included in the grid, such as... Figure 2 As shown.
[0104] After loading the observation system, the data is loaded into a grid. After loading, the line and point numbers of each survey line are checked and marked. At the same time, it is checked whether there is any data leakage. If leakage is found, the grid needs to be adjusted.
[0105] Figure 3 A schematic diagram of initial arrival pickup according to an embodiment of the present invention is shown, with red indicating the picked-up initial arrival.
[0106] like Figure 3 As shown, initial arrival data is picked up, the picking range is defined based on the initial arrival velocity, and an energy threshold is determined. The time that meets the conditions within the detection threshold is used to obtain the initial arrival data. Correlation analysis is performed on adjacent tracks, an adjacent track time difference threshold is set, non-compliant items are removed, and initial arrival information that meets the requirements is obtained.
[0107] Near-surface models are inverted through grid tomography; all survey lines are jointly inverted through grid tomography to obtain a three-dimensional near-surface model, thus avoiding errors caused by inconsistent inversion models between survey lines due to separate inversion of two-dimensional survey lines.
[0108] Figure 4a , Figure 4b A comparative schematic diagram of a near-surface inversion model and a near-surface data formation model according to an embodiment of the present invention is shown.
[0109] like Figure 4a , Figure 4b As shown, the rationality of the near-surface velocity model inverted is analyzed by comparing it with near-surface data. The near-surface model inverted from the survey line is compared and analyzed with near-surface data (small refraction, micrologging) to determine the rationality of the model (if the models are basically the same, but the details are better, it indicates that the model inversion is correct).
[0110] Figure 5a , Figure 5b Schematic diagrams showing the effects of a single shot before and after static correction according to an embodiment of the present invention are shown respectively.
[0111] Figure 6a , Figure 6b respectively show a profile comparison of elevation statics and tomographic statics according to an embodiment of the present application.
[0112] As shown in Figure 5a , Figure 5b , Figure 6a , Figure 6b , the statics loading is performed, the statics rationality is analyzed, the statics loading is performed, the single shot before and after the statics is compared, and the statics effect is analyzed according to the stack profile comparison. Generally, the single shot first arrival is leveled, the hyperbolic characteristics are obvious, the stack profile continuity is good, the structure understanding is consistent, and the statics can be applied.
[0113] According to the original data, the noise attenuation and deconvolution are performed. According to the two-dimensional line distribution, the velocity-line direction and the trace direction are selected, and the line direction and the trace direction data corresponding to the two-dimensional line are extracted.
[0114] The velocity spectrum of the selected line is formed, and the residual statics is completed. According to the conventional method, the velocity spectrum is formed, the velocity is picked up in combination with the line direction and the trace direction velocity trend, and the residual boundary positive iteration is performed.
[0115] The residual statics and the velocity analysis iteration processing are performed line by line. In the processing process, the velocity is an important parameter of the seismic data processing, and the accuracy of the velocity directly affects the effect of the seismic data stacking and migration. In order to ensure the accuracy of the seismic data processing, the accurate velocity must be obtained.
[0116] Figure 7 A schematic diagram of stack velocity analysis according to an embodiment of the present application is shown.
[0117] As shown in Figure 7 , in the velocity analysis, the velocity analysis and the surface consistency residual statics are iterated for multiple times to improve the accuracy of the velocity. At the same time, the velocity at the intersection of the lines is paid attention to in the line-by-line velocity analysis process to ensure the consistency of the velocity at the intersection.
[0118] The two-dimensional line pre-stack migration velocity analysis is extracted. The pseudo-three-dimensional pre-stack migration is performed. The target line is selected according to the two-dimensional line, and the migration velocity analysis is performed. The line direction and the trace direction velocity (two-dimensional line position) are analyzed, the pre-stack migration velocity analysis and iteration are completed, which ensures the closure of the line direction and the trace direction velocity, and provides a guarantee for the reliability of the velocity.
[0119] Figure 8 A migration velocity analysis point distribution diagram according to an embodiment of the present application is shown.
[0120] As shown in Figure 8The pre-stack migration velocity analysis is shown. The process of the pre-stack time migration is to image each common offset profile separately, and then stack all the results to form a migration profile. In the case of not severe lateral velocity variation, the pre-stack time migration transfers the reflection energy existing in each trace to its real underground position.
[0121] Figure 9a 、 Figure 9b The schematic diagrams of the east-west line and the north-south line of the migration velocity closure analysis according to one embodiment of the present application are shown respectively.
[0122] As Figure 9a 、 Figure 9b shown, the preprocessing in the time domain can provide high-quality prestack gather data for the pre-stack time migration, and the establishment of the pre-stack time migration velocity field needs to be analyzed, and the analysis method is similar to the stacking velocity analysis. In the process of the line-by-line migration velocity analysis, attention should also be paid to the velocity at the intersection of the lines, and the velocity at the intersection should be ensured to be consistent.
[0123] Figure 10 The schematic diagram of the line closure analysis according to one embodiment of the present application is shown.
[0124] The pseudo-three-dimensional pre-stack migration is performed according to the line position, the line closure is analyzed, and as Figure 10 shown, the lines 1, 2 and 3 are all closed well, and the expected effect is achieved, and the migration output result is obtained.
[0125] Example 4
[0126] Figure 11 The block diagram of a pseudo-three-dimensional loess table complex surface multi-line two-dimensional seismic data static correction and velocity analysis device according to one embodiment of the present application is shown.
[0127] As Figure 11 shown, the pseudo-three-dimensional loess table complex surface multi-line two-dimensional seismic data static correction and velocity analysis device comprises:
[0128] The grid establishment module 201 establishes a three-dimensional grid containing all two-dimensional line shot points in the whole area and loads seismic data;
[0129] The first arrival picking module 202 picks the first arrival of the seismic data in the whole area, and calculates the tomographic static correction amount;
[0130] The inversion module 203 inverts the near-surface velocity model based on the three-dimensional grid after the seismic data is loaded according to the first arrival;
[0131] The preprocessing module 204 pre-processes the noise attenuation data, performs stacking velocity analysis line by line, and obtains a velocity spectrum;
[0132] The processing module 205 performs velocity analysis and residual static correction on the velocity spectrum, performs pre-stack time migration line by line, and performs migration velocity analysis.
[0133] The output module 206 checks the intersection of the survey lines and outputs the migration result.
[0134] In one example, the method further comprises:
[0135] After loading the seismic data, it is checked whether the seismic data exceeds the grid range.
[0136] In one example, the noise attenuation data includes surface waves, linear interference, abnormal amplitude, and 50Hz industrial noise.
[0137] In one example, the preprocessing includes attenuation processing, surface consistency amplitude compensation, and surface consistency deconvolution.
[0138] In one example, the checking of the intersection of the survey lines includes:
[0139] Checking whether there is a time difference in the same set of strata at the intersection of different survey lines.
[0140] Example 5
[0141] The present disclosure provides an electronic device, comprising: a memory storing executable instructions; and a processor running the executable instructions in the memory to implement the above-mentioned pseudo-three-dimensional loess tableland complex surface multi-line two-dimensional seismic data static correction and velocity analysis method.
[0142] The electronic device according to the embodiments of the present disclosure comprises a memory and a processor.
[0143] The memory is configured to store non-transitory computer-readable instructions. Specifically, the memory can include one or more computer program products, which can include various forms of computer-readable storage media, such as volatile memory and / or non-volatile memory. The volatile memory may, for example, include random access memory (RAM), cache memory, and / or the like. The non-volatile memory may, for example, include read-only memory (ROM), hard disk, flash memory, and / or the like.
[0144] The processor can be a central processing unit (CPU) or other forms of processing units having data processing and / or instruction execution capabilities, and can control other components in the electronic device to perform desired functions. In one embodiment of the present disclosure, the processor is configured to run the computer-readable instructions stored in the memory.
[0145] Those skilled in the art shall understand that, in order to solve the technical problem of how to obtain a good user experience effect, the embodiment can also include well-known structures such as a communication bus, an interface, and the like, which shall also be included in the protection scope of the present disclosure.
[0146] The detailed description of the embodiment can refer to the corresponding description in the foregoing embodiments, and will not be repeated here.
[0147] Example 6
[0148] The embodiment of the present disclosure provides a computer readable storage medium, which stores a computer program. The computer program is executed by a processor to implement the pseudo-three-dimensional loess table complex surface multi-line two-dimensional seismic data static correction and velocity analysis method.
[0149] The computer readable storage medium according to the embodiment of the present disclosure has non-transitory computer readable instructions stored thereon. When the non-transitory computer readable instructions are run by a processor, all or part of the steps of the method of the embodiment of the present disclosure are executed.
[0150] The computer readable storage medium includes but is not limited to optical storage media (for example, CD-ROM and DVD), magneto-optical storage media (for example, MO), magnetic storage media (for example, magnetic tape or mobile hard disk), media with built-in rewritable non-volatile memory (for example, memory card), and media with built-in ROM (for example, ROM cartridge).
[0151] Those skilled in the art shall understand that the purpose of the above description of the embodiments of the present application is only to exemplarily illustrate the beneficial effects of the embodiments of the present application, and is not intended to limit the embodiments of the present application to any examples given.
[0152] The above has described the embodiments of the present application, and the above description is exemplary, is not exhaustive, and is not limited to the disclosed embodiments. Many modifications and changes are obvious to those skilled in the art without departing from the scope and spirit of the described embodiments.
Claims
1. A method for static correction and velocity analysis of multi-line two-dimensional seismic data on complex surfaces of the Loess Plateau (quasi-three-dimensional), characterized in that, include: A three-dimensional grid containing all two-dimensional seismic line shot points in the entire area was established and seismic data was loaded. First arrival picking is performed on the seismic data for the entire region, and tomographic static correction is calculated. Based on the three-dimensional grid after loading seismic data, a near-surface velocity model is retrieved according to the first arrival wave. The noise attenuation data is preprocessed, and velocity analysis is performed line by line to obtain the velocity spectrum. Velocity analysis and residual static correction were performed on the velocity spectrum, and pre-stack time migration and migration velocity analysis were performed line by line. The intersection points of the survey lines are checked, and the offset results are output.
2. The method for static correction and velocity analysis of multi-line two-dimensional seismic data for complex loess plateau surfaces according to claim 1, wherein, Also includes: After loading the seismic data, check whether the quality control seismic data exceeds the grid range.
3. The method for static correction and velocity analysis of multi-line two-dimensional seismic data for complex loess plateau surfaces according to claim 1, wherein, Noise attenuation data includes: surface waves, linear interference, abnormal amplitude, and 50Hz industrial electrical noise.
4. The method for static correction and velocity analysis of multi-line two-dimensional seismic data for complex loess plateau surfaces according to claim 1, wherein, Preprocessing includes attenuation processing, surface uniformity amplitude compensation, and surface uniformity deconvolution.
5. The method for static correction and velocity analysis of multi-line two-dimensional seismic data for complex loess plateau surfaces according to claim 1, wherein, The inspection of the intersection of the survey lines includes: Check whether there is a time difference at the intersection of different survey lines for the same set of strata.
6. A device for static correction and velocity analysis of multi-line two-dimensional seismic data of complex loess plateau surfaces (pseudo-three-dimensional), characterized in that, include: The grid creation module creates a 3D grid containing all 2D seismic line shot points in the entire area and loads seismic data. The first arrival picking module performs first arrival picking on the seismic data for the entire region and calculates the tomographic static correction. The inversion module, based on a 3D grid after loading seismic data, inverts the near-surface velocity model according to the first arrival wave; The preprocessing module preprocesses the noise attenuation data and performs line-by-line velocity analysis to obtain the velocity spectrum. The processing module performs velocity analysis and residual static correction on the velocity spectrum, and performs pre-stack time migration and migration velocity analysis line by line. The output module checks the intersection points of the survey lines and outputs the offset results.
7. The device for static correction and velocity analysis of multi-line two-dimensional seismic data of complex loess plateau surfaces according to claim 6, wherein, Also includes: After loading the seismic data, check whether the quality control seismic data exceeds the grid range.
8. The device for static correction and velocity analysis of multi-line two-dimensional seismic data of complex loess plateau surfaces according to claim 6, wherein, The inspection of the intersection of the survey lines includes: Check whether there is a time difference at the intersection of different survey lines for the same set of strata.
9. An electronic device, characterized in that, The electronic device includes: Memory, which stores executable instructions; A processor that executes the executable instructions in the memory to implement the method for static correction and velocity analysis of multi-line two-dimensional seismic data of complex loess plateau surfaces as described in any one of claims 1-5.
10. A computer-readable storage medium, characterized in that, The computer-readable storage medium stores a computer program that, when executed by a processor, implements the method for static correction and velocity analysis of multi-line two-dimensional seismic data for complex loess plateau surfaces as described in any one of claims 1-5.