An efficient quality control optimization method and device for RMS velocity picking in the time domain

Abstract

The present disclosure relates to an efficient quality control optimization method and device for time-domain RMS velocity picking, including: determining the corresponding migration profile according to the obtained velocity picking file; performing horizon interpretation on the migration profile to determine the reflection standard horizons on the migration profile; according to the velocity picking file, on each reflection standard horizon, statistically calculating the corresponding velocity picking values; respectively calculating the differences between every two adjacent velocity picking values after sorting; according to the differences, respectively establishing a velocity difference plan view for each reflection standard horizon, and finding the difference abnormal points on the velocity difference plan view; correcting the velocity picking values corresponding to the difference abnormal points, and after all the difference abnormal points are corrected, obtaining the velocity picking file after quality control optimization. To solve the problems that there will be differences when manually picking velocity values in the past, resulting in low accuracy of the velocity model, and the method for velocity quality control of the picked velocity values in the past is relatively cumbersome and cannot meet the timeliness requirements of data processing.

Description

Technical Field

[0001] The present disclosure relates to the technical field of seismic data processing, and particularly to an efficient quality control optimization method and device for RMS velocity picking in the time domain. Background Art

[0002] In recent years, with the continuous development and technological progress in the exploration field, the exploration target has gradually shifted from macroscopic to microscopic. The imaging quality of seismic data is mainly affected by the quality of CMP gathers and the accuracy of the velocity model. Therefore, higher requirements are placed on the accuracy of the velocity model in the time domain.

[0003] However, at present, velocity picking mainly relies on manual picking and the accumulation of work experience. The picked velocities may vary among different staff members, and there may also be certain differences in the velocities picked by the same staff member twice. Velocity quality control is a relatively cumbersome process, and there is a lack of a simple and efficient velocity quality control means. Summary of the Invention

[0004] The present disclosure proposes an efficient quality control optimization method and device for RMS velocity picking in the time domain to solve the problems that there are differences when picking velocity values manually in the past, resulting in low accuracy of the velocity model, and the method for velocity quality control of the picked velocity values in the past is relatively cumbersome and cannot meet the timeliness requirements of data processing.

[0005] According to one aspect of the present disclosure, there is provided an efficient quality control optimization method for RMS velocity picking in the time domain, including:

[0006] Obtain a velocity picking file for the work area, and determine the corresponding migration profile according to the velocity picking file;

[0007] Perform horizon interpretation on the migration profile to determine the reflection standard horizons on the migration profile;

[0008] According to the velocity picking file, on each of the reflection standard horizons, count the corresponding velocity pick values;

[0009] Sort the velocity pick values counted for each reflection standard horizon respectively, and calculate the difference between every two adjacent velocity pick values after sorting;

[0010] According to the differences, establish a velocity difference planar graph for each reflection standard horizon respectively, and find the difference abnormal points on the velocity difference planar graph;

[0011] Correct the velocity pick values corresponding to the difference abnormal points. After all the difference abnormal points are corrected, obtain a velocity picking file after quality control optimization.

[0012] Preferably, before obtaining the velocity picking file for the work area, it further includes: creating the velocity picking file, and the method includes:

[0013] Obtain the CMP gather of the work area;

[0014] Perform maximum amplitude coherent energy scanning on the CMP gather to obtain a velocity spectrum;

[0015] Perform velocity picking on the energy cluster center of the velocity spectrum to form a velocity picking file.

[0016] Preferably, the method for determining the corresponding migration profile according to the velocity picking file includes:

[0017] Perform interpolation smoothing processing on the velocity picking file to form an initial RMS velocity;

[0018] Perform pre-stack time migration operation using the RMS velocity and the CMP gather to obtain a CRP gather;

[0019] Perform stacking processing on the CRP gather to obtain a migration profile.

[0020] Preferably, the method for statistically calculating the corresponding velocity picking values on each reflection standard layer according to the velocity picking file includes:

[0021] According to the velocity picking file, perform anti-normal moveout on the CRP gather to obtain an anti-normal moveout gather;

[0022] Perform maximum amplitude coherent energy scanning on the anti-normal moveout gather to obtain a velocity spectrum;

[0023] On the velocity spectrum, use the reflection standard layer as a constraint condition to statistically calculate the velocity picking values.

[0024] Preferably, the method for respectively establishing a velocity difference planar graph for each reflection standard layer according to the difference includes:

[0025] Use the trace number of the velocity picking value as the vertical coordinate and the line number as the horizontal coordinate to establish a velocity difference planar graph for each reflection standard layer.

[0026] Preferably, the method for correcting the velocity picking value corresponding to the difference anomaly point includes:

[0027] Refer to the overall structural change trend and velocity change trend of the reflection standard layer corresponding to the difference anomaly point, and adjust the velocity picking value corresponding to the difference anomaly point to obtain an adjusted velocity picking value;

[0028] Judge whether the adjusted velocity picking value satisfies the flattening of the common reflection point gather. If so, determine the difference corresponding to the adjusted velocity picking value;

[0029] Using the velocity error calculation formula, determine the error of the difference corresponding to the adjusted velocity pick value;

[0030] Determine whether the magnitude of the error is within a predetermined percentage of the average velocity of its corresponding reflection standard layer. If so, the adjusted velocity pick value is the corrected velocity pick value.

[0031] Preferably, the velocity error calculation formula includes:

[0032] Wherein, In the formula, Ea is the velocity error amount, Δv ti is the velocity difference at the i-th point within the time window t, and v av is the average value of all picked velocities within the time window t.

[0033] According to one aspect of the present disclosure, there is provided a high-efficiency quality control optimization device for time-domain RMS velocity picking, including:

[0034] An acquisition unit for acquiring a velocity picking file of the work area and determining the corresponding migration profile according to the velocity picking file;

[0035] An interpreted horizon determination unit for performing horizon interpretation on the migration profile and determining the reflection standard layers on the migration profile;

[0036] A velocity pick value determination unit for statistically calculating the corresponding velocity pick values on each of the reflection standard layers according to the velocity picking file;

[0037] A difference determination unit for sorting the velocity pick values statistically calculated for each reflection standard layer and calculating the difference between every two adjacent velocity pick values after sorting;

[0038] A difference plane map determination unit for respectively establishing a velocity difference plane map for each reflection standard layer according to the differences and finding the difference anomaly points on the velocity difference plane map;

[0039] An outlier correction unit for correcting the velocity pick values corresponding to the difference anomaly points. After all the difference anomaly points are corrected, a quality control optimized velocity picking file is obtained.

[0040] The present invention has at least the following beneficial effects:

[0041] The present disclosure provides an efficient quality control optimization method and device for time-domain RMS velocity picking. By using a reflection standard layer for constraint, the velocity picking values at the same interpreted horizon are statistically analyzed. By calculating the differences between adjacent points and establishing a difference plan view, abnormal points are found and corrected, so that the efficiency and accuracy of velocity value picking and quality control are higher, and the accuracy of each velocity picking file can be quality-controlled without performing final migration imaging processing. BRIEF DESCRIPTION OF THE DRAWINGS

[0042] The accompanying drawings are incorporated herein and form a part of this specification. These drawings illustrate embodiments consistent with the present disclosure and, together with the specification, are used to explain the technical solutions of the present disclosure.

[0043] Figure 1 A flowchart showing an efficient quality control optimization method for time-domain RMS velocity picking according to an embodiment of the present disclosure;

[0044] Figure 2 Showing a CMP gather and the corresponding velocity spectrum according to an embodiment of the present disclosure;

[0045] Figure 3 Showing a comparison graph of the abnormal velocity change at a certain velocity picking point and the velocity values of its adjacent points according to an embodiment of the present disclosure;

[0046] Figure 4 Showing a velocity slice corresponding to the abnormal velocity picking point before correction according to an embodiment of the present disclosure;

[0047] Figure 5 Showing a structural diagram of the velocity slice corresponding to the abnormal velocity picking point after correction and the reflection standard horizon corresponding to the abnormal velocity picking point according to an embodiment of the present disclosure. DETAILED DESCRIPTION OF THE EMBODIMENTS

[0048] Various exemplary embodiments, features, and aspects of the present disclosure will be described in detail below with reference to the accompanying drawings. The same reference numerals in the drawings denote elements having the same or similar functions. Although various aspects of the embodiments are shown in the drawings, the drawings are not necessarily drawn to scale unless otherwise specified.

[0049] The term "exemplary" used herein means "serving as an example, embodiment, or illustration". Any embodiment described herein as "exemplary" is not necessarily to be construed as superior to or better than other embodiments.

[0050] As used herein, the term "and / or" is merely a description of the associated relationship between associated objects, indicating that there can be three relationships. For example, A and / or B can represent three situations: A exists alone, A and B exist simultaneously, and B exists alone. In addition, the term "at least one" as used herein means any one of a plurality or any combination of at least two of a plurality. For example, including at least one of A, B, and C can mean including any one or more elements selected from the set composed of A, B, and C.

[0051] In addition, for a better illustration of the present disclosure, numerous specific details are given in the following detailed implementation manners. Those skilled in the art should understand that the present disclosure can also be implemented without some specific details. In some instances, methods, means, elements, and circuits well-known to those skilled in the art are not described in detail so as to highlight the gist of the present disclosure.

[0052] Figure 1 A flowchart showing a method for efficiently optimizing quality control of time-domain RMS velocity picking according to an embodiment of the present disclosure; Figure 2 Showing a CMP gather and its corresponding velocity spectrum according to an embodiment of the present disclosure; Figure 3 Showing a comparison graph of abnormal velocity changes at a certain velocity picking point and the velocity values of adjacent points according to an embodiment of the present disclosure; Figure 4 Showing a velocity slice corresponding to an abnormal velocity picking point before correction according to an embodiment of the present disclosure; Figure 5 Showing a structural diagram of a velocity slice corresponding to an abnormal velocity picking point after correction and a reflection standard horizon corresponding to the abnormal velocity picking point according to an embodiment of the present disclosure. As Figures 1-5 As shown, a method for efficiently optimizing quality control of time-domain RMS velocity picking includes: Step S01: Obtain a velocity picking file for a work area, and determine the corresponding migration profile according to the velocity picking file; Step S02: Interpret horizons for the migration profile to determine the reflection standard horizons on the migration profile; Step S03: According to the velocity picking file, on each of the reflection standard horizons, count the corresponding velocity picking values; Step S04: Sort the velocity picking values counted for each reflection standard horizon respectively, and calculate the difference between every two adjacent velocity picking values after sorting; Step S05: According to the differences, establish a velocity difference planar graph for each reflection standard horizon respectively, and find the difference abnormal points on the velocity difference planar graph; Step S06: Correct the velocity picking values corresponding to the difference abnormal points. After all the difference abnormal points are corrected, a velocity picking file with optimized quality control is obtained.

[0053] The method for efficiently optimizing quality control of time-domain RMS velocity picking provided by the embodiments of the present invention specifically includes the following steps:

[0054] Step S01: Obtain the velocity picking file of the work area, and determine the corresponding migration profile according to the velocity picking file.

[0055] In the present disclosure, before obtaining the velocity picking file of the work area, it further includes: creating the velocity picking file, and the method includes: obtaining the CMP gather of the work area; performing maximum amplitude coherent energy scanning on the CMP gather to obtain a velocity spectrum; picking velocities at the centers of energy clusters of the velocity spectrum to form a velocity picking file.

[0056] In the present disclosure, the method for determining the corresponding migration profile according to the velocity picking file includes: performing interpolation and smoothing processing on the velocity picking file to form an initial RMS velocity; performing pre-stack time migration operation using the RMS velocity and the CMP gather to obtain a CRP gather; performing stacking processing on the CRP gather to obtain a migration profile.

[0057] In an embodiment of the present disclosure, the CMP gather (common midpoint gather) is obtained by preprocessing the original seismic data. Perform maximum amplitude coherent energy scanning on the obtained CMP gather to obtain a velocity spectrum, and pick velocity values on the velocity spectrum to form a velocity picking file.

[0058] Perform interpolation and smoothing processing on the velocity picking file to form an RMS (root mean square) velocity. Use the RMS velocity and the CMP gather to perform pre-stack time migration operation to obtain a CRP gather (common reflection point gather); perform stacking processing on the CRP gather to obtain a migration profile.

[0059] Figure 2 Show the CMP gather (right figure) and the corresponding velocity spectrum (left figure) in an embodiment of the present disclosure. In Figure 2 the velocity spectrum, the black dots are the centers of energy clusters. When picking velocities, make the picked points coincide with the centers of energy clusters, and at the same time ensure that the CMP gather is flattened.

[0060] Step S02: Perform horizon interpretation on the migration profile to determine the reflection standard horizons on the migration profile.

[0061] In an embodiment of the present disclosure, the method for determining the reflection standard horizons on the migration profile includes: obtaining the well logging data of the work area, and determining all the positions of the reflection standard horizons on the geological horizons at the well points in the well logging data; on the migration profile, find the positions corresponding to the reflection standard horizons on the geological horizons at the well points, and from shallow to deep, perform horizon interpretation and tracing on each reflection standard horizon on the seismic profile, so as to determine all the positions of the reflection standard horizons on the migration profile.

[0062] Determining the position of each reflection standard layer on the migration profile is to ensure that during velocity picking, for each reflection standard layer position, velocity analysis is performed on the corresponding in-phase axis on the gather, and in subsequent steps, when statistically analyzing the velocity pick values of each CRP in-phase axis, it is ensured that they represent the same geological horizon (reflection standard layer).

[0063] The seismic horizon interpretation work assigns clear geological significance to the seismic reflection horizons based on the waveform characteristics and geological laws of the migration profile. The migration profile is the seismic response of the geological profile. In the seismic profile, a large amount of geological information is contained, and the geological phenomena involved in seismic reflection are all reflected in the seismic profile. However, the seismic profile also contains noise unrelated to geological phenomena, which has no geological significance. The selection of the reflection standard layer is to determine the reflection standard layer by combining the migration profile with known drilling, logging, vertical seismic profile and other data. In addition, the reflection standard layer must also have the characteristics of wide distribution range, stable distribution, and relatively clear geological horizon. Generally, a reflection layer with good continuity and stable waveform in the migration profile should be selected.

[0064] Step S03: According to the velocity picking file, statistically analyze the corresponding velocity pick values on each of the reflection standard layers.

[0065] In the present disclosure, the method for statistically analyzing the corresponding velocity pick values on each of the reflection standard layers according to the velocity picking file includes: performing inverse normal moveout on the CRP gather according to the velocity picking file to obtain an inverse normal moveout gather; obtaining a velocity spectrum by performing maximum amplitude coherent energy scanning on the inverse normal moveout gather; and statistically analyzing the velocity pick values on the velocity spectrum by using the reflection standard layer as a constraint condition.

[0066] In the embodiment of the present disclosure, velocity values are picked on the velocity spectrum, and the reflection standard layer is combined as a constraint condition during the picking process, that is, the velocity values corresponding to each reflection standard layer are respectively picked on the velocity spectrum.

[0067] Combining the reflection standard layer as a constraint condition can ensure that each reflection standard layer is accurately identified, and the picked velocity values can truly reflect the velocity change of the reflection standard layer.

[0068] During picking, the velocity pick values corresponding to the velocity spectrum are statistically analyzed for each reflection standard layer in a time window. The velocity pick value is the optimal imaging velocity value of each strong reflection axis, denoted as v ti , that is, the velocity value at time window t and the i-th point.

[0069] The principle of the time window is that one time window contains a strong reflection axis on a gather, that is, a velocity picking point, so that it can be ensured that the velocity picking points within the time window are for a specific formation.

[0070] Step S04: Sort the velocity pick values statistically obtained for each reflection standard layer respectively, and calculate the difference between every two adjacent velocity pick values after sorting.

[0071] In the embodiments of the present disclosure, the velocity pick values are sorted using the line numbers in the work area, and the difference calculation and statistics work are carried out in the front-to-back order. Calculate the difference Δv ti between two adjacent velocity pick values v ti as Δv t(i+1) = v ti - v

[0072] That is, within the t time window, it is the velocity difference between the (i + 1)-th point and the i-th point.

[0073] In the present disclosure, the method for respectively establishing the velocity difference plan view for each reflection standard layer according to the difference includes: using the line number of the velocity pick value as the ordinate and the line number as the abscissa to establish the velocity difference plan view for each reflection standard layer.

[0074] In the embodiments of the present disclosure, with the line number as the ordinate and the line number as the abscissa, a velocity difference plan view is respectively formed for each reflection standard layer, and difference abnormal points are searched on the velocity difference plan view. Under normal circumstances, the velocity change of each reflection standard layer is a stable quantity changing along the structural trend. If there is an abnormal change in the difference on the velocity difference plan view of a certain reflection standard layer, it is considered that there is a certain error in the velocity pick at this point.

[0075] Step S06: Correct the velocity pick values corresponding to the difference abnormal points. After all the difference abnormal points are corrected, a velocity pick file optimized by quality control is obtained.

[0076] In the present disclosure, the method for correcting the velocity pick values corresponding to the difference abnormal points includes: referring to the overall structural change trend and velocity change trend of the reflection standard layer corresponding to the difference abnormal point, adjusting the velocity pick value corresponding to the difference abnormal point to obtain an adjusted velocity pick value; judging whether the adjusted velocity pick value satisfies the flattening of the common reflection point (CRP) gather in-phase axis. If so, determine the difference corresponding to the adjusted velocity pick value; use the velocity error calculation formula to determine the error of the difference corresponding to the adjusted velocity pick value; judge whether the magnitude of the error is within a predetermined percentage of the velocity average value of its corresponding reflection standard layer. If so, the adjusted velocity pick value is the corrected velocity pick value.

[0077] In the present disclosure, the velocity error calculation formula includes:

[0078] Among them,

[0079] in the formula, Ea is the velocity error amount, and Δv ti is the velocity difference at the i-th point within the time window t, and v av is the average value of all picked-up velocities within the time window t.

[0080] In the embodiments of the present disclosure, the difference anomaly points found on the difference plan view correspond to two adjacent velocity pick-up values.

[0081] Analyze the differences between the two adjacent velocity pick-up values corresponding to the reflection standard layer of the difference anomaly point and their front and rear velocity pick-up values, and then refer to the overall structural change trend and velocity change trend of the corresponding reflection standard layer to manually adjust the two or one of the velocity pick-up values corresponding to the difference anomaly point.

[0082] The adjusted velocity pick-up value needs to satisfy the flattening of the CRP gather after reaction. Under this condition, recalculate the difference between the two adjusted velocity pick-up values, and use the velocity error calculation formula to calculate the error Ea of this difference.

[0083] Among them, the predetermined percentage is: 3%. Determine whether the error Ea is less than or equal to 3% of the average velocity of the corresponding reflection standard layer. If so, the adjusted velocity pick-up value is the corrected velocity pick-up value; if not, the velocity pick-up value needs to be adjusted again until the difference error between the velocity pick-up value and the adjacent point is within 3% of the average velocity of the reflection standard layer where it is located under the condition of satisfying the CRP gather flattening.

[0084] For example, if the average velocity of a certain formation (reflection standard layer) is 3000 m / s, then its value range should be within 3000 * 3% = 90, that is, the change amount of the adjacent point velocity pick-up value does not exceed 90 m / s.

[0085] Figure 3 is a comparison diagram of the abnormal velocity change of a certain velocity pick-up point and the velocity values of its adjacent points; in Figure 3 , the point with the velocity pick-up value v = 3156 is an abnormal velocity pick-up point, and the two velocity pick-up values before and after it are 2993 and 3000. Figure 4 is the velocity slice corresponding to the abnormal velocity pick-up point before correction; Figure 5 In the right figure of [], it is the velocity slice corresponding to the abnormal velocity pick-up point after correction, and the left figure is the structural diagram of the reflection standard layer position corresponding to the abnormal velocity pick-up point. It can be seen from Figure 4 , 5 that the velocity pick-up value coincides with the structural trend in the velocity plane slice after correction and optimization.

[0086] After sequentially completing the calibration of the velocity pick-up values corresponding to all the difference anomaly points of each reflection standard layer, a velocity pick-up file optimized by quality control is obtained.

[0087] It can be understood that, without violating the principle logic, the above-mentioned various method embodiments mentioned in the present disclosure can be combined with each other to form a combined embodiment. Due to space limitations, the present disclosure will not elaborate further.

[0088] The execution subject of a method for efficiently optimizing quality control of time-domain RMS velocity picking can be a device for efficiently optimizing quality control of time-domain RMS velocity picking. For example, the method for efficiently optimizing quality control of time-domain RMS velocity picking can be executed by a terminal device, a server, or other processing devices. Among them, the terminal device can be a user equipment (UE), a mobile device, a user terminal, a terminal, a cellular phone, a cordless phone, a personal digital assistant (PDA), a handheld device, a computing device, a vehicle-mounted device, a wearable device, etc. In some possible implementation manners, the method for efficiently optimizing quality control of time-domain RMS velocity picking can be implemented by a processor invoking computer-readable instructions stored in a memory.

[0089] Those skilled in the art can understand that in the above method of the specific implementation manner, the writing order of each step does not mean a strict execution order that constitutes any limitation to the implementation process. The specific execution order of each step should be determined according to its function and possible internal logic.

[0090] The present disclosure provides a device for efficiently optimizing quality control of time-domain RMS velocity picking, including: an acquisition unit for acquiring a velocity pick-up file of a work area and determining a corresponding migration profile according to the velocity pick-up file; an interpreted horizon determination unit for performing horizon interpretation on the migration profile to determine the reflection standard horizons on the migration profile; a velocity pick-up value determination unit for, according to the velocity pick-up file, statistically calculating the corresponding velocity pick-up values on each reflection standard horizon; a difference determination unit for respectively sorting the velocity pick-up values statistically calculated for each reflection standard horizon and calculating the difference between every two adjacent velocity pick-up values after sorting; a difference plane graph determination unit for respectively establishing a velocity difference plane graph for each reflection standard horizon according to the difference and finding the difference anomaly points on the velocity difference plane graph; an outlier correction unit for correcting the velocity pick-up values corresponding to the difference anomaly points. After all the difference anomaly points are corrected, a velocity pick-up file optimized by quality control is obtained.

[0091] In some embodiments, the functions or modules included in the apparatus provided by the embodiments of the present disclosure can be used to execute the methods described in the above method embodiments. The specific implementation can refer to the description of the above method embodiments. For the sake of brevity, it will not be repeated here.

[0092] The core of the present disclosure lies in determining each reflection standard layer by performing horizon interpretation on the migration profile, and analyzing and picking up velocity values through the time window constraint of the reflection standard layer; the velocity points of each reflection standard layer represent the same underground reflection layer, and the change amount of its velocity value is controlled within a certain range; the accuracy and rationality of each velocity picking point are quality-controlled by using the front and back change amounts of the velocity, and finally an optimized velocity picking file with higher accuracy is obtained. The operation process of the RMS velocity picking quality control optimization method of the present disclosure is faster, simpler, and has higher accuracy. It is not necessary to perform final migration imaging processing to quality-control the accuracy of each velocity picking file, thereby effectively improving the seismic imaging accuracy.

[0093] The embodiments of the present disclosure have been described above. The above description is exemplary and not exhaustive, and is not limited to the disclosed embodiments. Many modifications and variations are obvious to those of ordinary skill in the art in the technical field without departing from the scope and spirit of the described embodiments. The selection of the terms used herein is intended to best explain the principles, practical applications, or improvements to the technology in the market of the embodiments, or to enable other ordinary skill in the art in the technical field to understand the disclosed embodiments.

Claims

1. An efficient quality control optimization method for RMS velocity picking in the time domain, characterized in that, Including: Obtain the velocity picking file of the work area, and determine the corresponding migration profile according to the velocity picking file. The method includes: performing interpolation smoothing on the velocity picking file to form an initial RMS velocity; performing prestack time migration operation using the RMS velocity and the CMP gather to obtain a CRP gather; performing stacking processing on the CRP gather to obtain a migration profile; Perform horizon interpretation on the migration profile to determine the reflection standard horizon on the migration profile; According to the velocity picking file, on each of the reflection standard horizons, count the corresponding velocity picking values; Sort the velocity picking values counted for each reflection standard horizon respectively, and calculate the difference between every two adjacent velocity picking values after sorting; According to the differences, establish a velocity difference planar map for each reflection standard horizon respectively, and find the difference anomaly points on the velocity difference planar map; Correct the velocity picking value corresponding to the difference anomaly point. The method includes: referring to the overall structural change trend and velocity change trend of the reflection standard horizon corresponding to the difference anomaly point, adjusting the velocity picking value corresponding to the difference anomaly point to obtain an adjusted velocity picking value; determining whether the adjusted velocity picking value satisfies the flattening of the common reflection point (CRP) gather in-phase axis. If so, determine the difference corresponding to the adjusted velocity picking value; using the velocity error calculation formula, determine the error of the difference corresponding to the adjusted velocity picking value; determining whether the magnitude of the error is within a predetermined percentage of the velocity average value of its corresponding reflection standard horizon. If so, the adjusted velocity picking value is the corrected velocity picking value; The velocity error calculation formula includes: Among them, ; where Ea is the velocity error, and △ v ti is the velocity difference at the i-th point within the time window t, and v av is the average value of all picked-up velocities within the time window t; After all the difference anomaly points are corrected, obtain a quality-controlled and optimized velocity picking file.

2. The efficient quality control optimization method for time-domain RMS velocity picking according to claim 1, characterized in that Before the step of obtaining the velocity picking file of the work area, it further includes: creating the velocity picking file. The method includes: Obtain the CMP gather of the work area; Perform maximum amplitude coherence energy scanning on the CMP gather to obtain a velocity spectrum; Perform velocity picking on the center of the energy cluster of the velocity spectrum to form a velocity picking file.

3. The efficient quality control optimization method for time-domain RMS velocity picking according to claim 1, wherein The method of counting the corresponding velocity picking values on each of the reflection standard horizons according to the velocity picking file includes: According to the velocity picking file, perform de-migration on the CRP gather to obtain a de-migrated gather; Obtain a velocity spectrum by performing maximum amplitude coherence energy scanning on the de-migrated gather; On the velocity spectrum, use the reflection standard horizon as a constraint condition to count the velocity picking values.

4. The high-efficiency quality control optimization method for time-domain RMS velocity picking according to any one of claims 1-3, characterized in that The method of establishing a velocity difference planar map for each reflection standard horizon according to the differences includes: Establish a velocity difference planar map for each reflection standard horizon with the trace number of the velocity picking value as the vertical coordinate and the line number as the horizontal coordinate.

5. An efficient quality control optimization device for time-domain RMS velocity picking, characterized in that, Including: An acquisition unit for acquiring a velocity picking file of a work area and determining a corresponding migration profile according to the velocity picking file, the method comprising: performing interpolation and smoothing processing on the velocity picking file to form an initial RMS velocity; using the RMS velocity and the CMP gather to perform prestack time migration operation to obtain a CRP gather; performing stacking processing on the CRP gather to obtain a migration profile; An interpretation horizon determination unit for performing horizon interpretation on the migration profile to determine a reflection standard horizon on the migration profile; A velocity pick value determination unit for statistically calculating corresponding velocity pick values on each of the reflection standard horizons according to the velocity picking file; A difference determination unit for sorting the velocity pick values statistically calculated for each reflection standard horizon respectively and calculating the difference between every two adjacent velocity pick values after sorting; A difference plan view determination unit for respectively establishing a velocity difference plan view for each reflection standard horizon according to the difference and finding difference anomaly points on the velocity difference plan view; An outlier correction unit for correcting the velocity pick value corresponding to the difference anomaly point, the method comprising: referring to the overall structural change trend and velocity change trend of the reflection standard horizon corresponding to the difference anomaly point to adjust the velocity pick value corresponding to the difference anomaly point to obtain an adjusted velocity pick value; judging whether the adjusted velocity pick value satisfies the flattening of the common reflection point gather; if so, determining the difference corresponding to the adjusted velocity pick value; using a velocity error calculation formula to determine the error of the difference corresponding to the adjusted velocity pick value; judging whether the magnitude of the error is within a predetermined percentage of the velocity average value of its corresponding reflection standard horizon; if so, the adjusted velocity pick value is the corrected velocity pick value; the velocity error calculation formula includes: Among them, ; where Ea is the velocity error amount, △ v ti is the velocity difference at the i-th point within the time window t, and v av is the average value of all picked-up velocities within the time window t; After all the difference anomaly points are corrected, a quality-controlled and optimized velocity picking file is obtained.

Images