Processing method and device for trace gather record and storage medium

[0020] The embodiments of this specification provide a processing method, device and storage medium for gather recording.

[0021] In order to enable those skilled in the art to better understand the technical solutions in this specification, the following will clearly and completely describe the technical solutions in the embodiments of this specification with reference to the drawings in the embodiments of this specification. Obviously, the described The embodiments are only a part of the embodiments in this specification, rather than all the embodiments. Based on the implementation of this specification, all other implementations obtained by those of ordinary skill in the art without creative work shall fall within the protection scope of this specification.

[0022] The processing method of gather recording provided by the embodiment of this specification will be described in detail below with reference to the accompanying drawings. Although this specification provides the method operation steps as described in the following embodiments or flowcharts, more or less operation steps can be included in the method based on routine or without creative labor, and the order of execution between the operation steps is not limit.

[0023] See figure 1 The processing method for gather records provided by the embodiment of this specification may include the following steps.

[0024] Step S10: Perform dynamic correction on the gather record according to a seismic wave velocity value.

[0025] In this embodiment, the seismic wave velocity may include the propagation velocity of the seismic wave in the ground. In the case where the seismic wave velocity is the same as the real velocity of the formation, the image generated after the gather record motion correction is a horizontal straight line. The seismic wave velocity value can be used as a processing parameter for motion correction, migration, time-depth conversion and other processing.

[0026] In this embodiment, the seismic wave velocity value may be randomly selected or selected according to actual engineering needs. For example, a certain velocity value interval is determined, and a seismic wave velocity value is selected from the velocity value interval.

[0027] In this embodiment, the dynamic correction may include a process of correcting the seismic wave reflection time of different offsets to the seismic wave reflection time of zero offsets. The dynamic correction amount can be calculated according to the seismic wave time, and the dynamic correction processing is performed on the gather record according to the dynamic correction amount. The selected seismic wave velocity is different, the calculated motion correction amount is also different, and the motion correction effect brought by different motion correction amount is also different.

[0028] In this embodiment, the step of dynamically correcting the gather record according to a seismic wave velocity value may include, according to the reflection time of the seismic wave in the case of zero offset and the known offset of the seismic trace record in the set record And the seismic wave velocity value, according to the formula, The seismic wave reflection time recorded by the seismic trace can be calculated; where t i The seismic wave reflection time recorded for the seismic trace, t 0 Is the reflection time of the zero offset seismic wave, x i Is the offset recorded by the seismic trace, and v is the seismic wave velocity value. The seismic wave reflection time recorded by the seismic trace is subtracted from the zero offset seismic wave reflection time to calculate the motion correction amount of the seismic trace record. According to the formula, Calculated. Perform dynamic correction processing on the gather record according to the dynamic correction amount.

[0029] Step S12: Calculate the average amplitude, similarity coefficient and amplitude variance of the seismic trace records in the gather record after the dynamic correction; wherein the similarity coefficient is used to indicate the similarity of the waveforms of the seismic trace records in the gather record Degree; wherein the gather record includes at least one seismic trace record.

[0030] In this embodiment, after performing dynamic correction according to the seismic wave velocity, the seismic wave recorded by the seismic trace can be corrected to a position close to the same phase. If the seismic wave recorded by the seismic trace is corrected to the in-phase position through the dynamic correction processing, the superimposed energy of the gather record can reach the maximum.

[0031] In this embodiment, the amplitude may include the magnitude of the deviation of the proton from the equilibrium position in the seismic wave. For different seismic wave reflection times, the amplitude can be different. The seismic wave reflection time may correspond to a sample point in the seismic trace record, and the sample point value of the sample point may be the amplitude of the seismic trace record under the seismic wave reflection time.

[0032] In this embodiment, after the dynamic correction processing is performed on the gather record according to the seismic wave velocity. A time window can be selected or randomly selected to obtain the amplitude of the seismic trace recorded in the window. Specifically, for example, in the gather record, a time interval of seismic wave reflection time is randomly selected, and the amplitude of the seismic trace record within 250ms to 500ms is selected; in the gather record, one of the seismic traces is selected The seismic wave reflection time corresponding to the connection line of the recorded seismic wave with the same phase, the reflection time is 300ms, and the amplitude of the seismic trace recorded within 200ms to 400ms is selected with 300ms as the center.

[0033] In this embodiment, the average amplitude of the gather record may include, in the gather record, an average value of the amplitude of at least one seismic trace record. Calculating the average amplitude of the gather records may include adding the amplitudes of the at least one seismic trace record in-phase and dividing by the number of seismic trace records of the at least one seismic trace record to obtain the average amplitude. Specifically, the average amplitude may be based on the formula: It is calculated that the seismic wave reflection time recorded by the seismic trace t 0 Represents the seismic wave reflection time at zero offset, x i Represents the offset, v represents the seismic wave velocity, x i Is the offset of the i-th channel, N represents the number of channels, λ represents the length of the time window, u(t i +j,x i ) Represents the amplitude of the seismic trace record, t i +j is the jth sample point recorded by the seismic trace.

[0034] In this embodiment, the similarity coefficient may be the degree of similarity of the seismic waves recorded by the seismic traces in the gather record. When the seismic waves of each trace are exactly the same, the similarity coefficient may be 1; When it is random, the similarity coefficient approaches 0. When the selected seismic wave velocity is equal to the dynamic correction speed, the waveforms of the seismic traces are most similar. In the time window, the similarity coefficient is close to 1. Specifically, the similarity coefficient can be based on the formula It is calculated that the seismic wave reflection time recorded by the seismic trace t0 represents the seismic wave reflection time at zero offset, x i Represents the offset, v represents the seismic wave velocity, x i Is the offset of the i-th channel, N represents the number of channels, λ represents the length of the time window, u(t i +j,x i ) Represents the amplitude of the seismic trace record, t i +j is the jth sample point recorded by the seismic trace.

[0035] In this embodiment, the amplitude variance may be the degree of deviation between each sample value in each seismic trace record in the gather record and the average value. Specifically, the in-phase addition of the amplitudes of the seismic trace records in the gather record can be performed according to the formula The amplitude variance is calculated, where ti is the seismic wave reflection time recorded by the seismic trace, t 0 Represents zero offset time, v represents seismic wave velocity, x i Is the offset of the i-th channel, N represents the number of channels, λ represents the width of the time window, u(t i +j,x i ) Represents the amplitude of the seismic trace record, Indicates the expectation of the data, σ 2 It is a positive number with a small value.

[0036] In this embodiment, when different seismic wave velocities are selected, the calculated motion correction amount may be different; the motion correction amount is different, and the result of in-phase superposition of the seismic trace records after motion correction may also be different.

[0037] Step S14: calculating a velocity evaluation value according to the average amplitude, similarity coefficient and amplitude variance; wherein the velocity evaluation value is used to indicate the closeness of the seismic wave velocity value to the formation velocity.

[0038] In this embodiment, the target dynamic correction velocity may include the seismic wave velocity value that is closest to the formation velocity. After the target dynamic correction velocity is dynamically corrected with the gather record, the profile of the gather record can be leveled into one. Horizontal straight line.

[0039] In this embodiment, the formation velocity may include the velocity of seismic wave propagation in the layered formation. When the seismic wave velocity value is equal to the formation velocity, dynamic correction is performed on the gather record, and the profile of the gather record can be leveled into a horizontal straight line.

[0040] In this embodiment, when the average amplitude and the similarity coefficient are the largest, the correspondingly selected seismic wave velocity can be closest to the target motion correction velocity recorded by the gather, and the corresponding seismic wave velocity when the amplitude variance is the smallest The target motion correction speed recorded by the gather can be the closest.

[0041] In this embodiment, the speed evaluation value can be based on the formula Calculated, where A is the average amplitude, S is the similarity coefficient, V is the amplitude variance, σ 2 It can be a positive number with a small value to prevent the denominator from being zero.

[0042] In this embodiment, it can be seen from the formula that when the average amplitude and similarity coefficient are the largest, the seismic wave velocity value is closest to the formation velocity, and when the amplitude variance is the smallest, the seismic wave velocity is closest to the ground velocity. The formation velocity. Therefore, the greater the velocity evaluation value obtained by multiplying the average amplitude and the similarity coefficient and then by the reciprocal of the amplitude variance, the greater the seismic wave velocity is, the closer to the formation velocity.

[0043] Step S16: Changing the seismic wave velocity value to obtain the maximum value of the velocity evaluation value; wherein the seismic wave velocity value corresponding to the maximum value of the velocity evaluation value is the target motion correction velocity recorded by the gather.

[0044] In this embodiment, the target motion correction velocity may be the seismic wave velocity value closest to the formation velocity. Performing dynamic correction processing on the gather record according to the target dynamic correction speed can level the profile of the gather record into a horizontal straight line.

[0045] In this embodiment, to change the seismic wave velocity value, a velocity value can be randomly selected or the seismic wave velocity value can be changed according to a specified rule. A certain speed range can be selected according to actual needs, and a certain speed interval can be set. The value can be changed from large to small according to the size of the speed value, or the value can be changed from small to large, and the seismic wave velocity can be changed in any way .

[0046] In this embodiment, the calculation of the maximum value of the velocity evaluation value may be that every time the seismic wave velocity value is changed, the corresponding average amplitude, amplitude variance, and similarity coefficient are calculated to obtain a velocity evaluation value. The seismic wave velocity value is changed by the method, and the calculation of the velocity evaluation value is ended when no more changes are made, and all the velocity evaluation values ​​obtained are compared to obtain the maximum value of the velocity evaluation value. The corresponding seismic wave velocity can be close to the formation velocity. Specifically, the speed range is selected from 1000m/s to 5000m/s, with 100m/s as the speed interval, the speed range is divided into 50 speed sections, 100m/s is used to calculate the speed evaluation value, and 200m/s is used to calculate the speed Evaluation value, 300m/s calculate the speed evaluation value until the speed evaluation value corresponding to 5000m/s is calculated, compare the calculated 50 speed evaluation values ​​to obtain the maximum value of the speed evaluation value, and obtain The seismic wave velocity corresponding to the maximum value of the velocity evaluation value.

[0047] In this embodiment, the results of the in-phase addition of the seismic wave amplitudes of the seismic trace records obtained by changing different seismic wave speeds may be different; the results of the in-phase addition of the seismic wave amplitudes of the seismic trace records are different, and in the gather record The average amplitude, similarity coefficient, and amplitude variance of the seismic trace records may be different; the velocity evaluation value calculated for different average amplitude, similarity coefficient, and amplitude variance may also be different. When the maximum value of the velocity evaluation value is calculated, the seismic wave velocity corresponding to the velocity evaluation value may be used as the target motion correction velocity recorded in the gather.

[0048] In this embodiment, the average amplitude, correlation coefficient, and amplitude variance of the gather record are calculated using the seismic wave velocity as a variable to obtain the velocity evaluation value. The average amplitude, correlation coefficient, and amplitude variance are all related to the dynamic correction speed of the gather record. Relatedly, calculating the maximum value of the velocity evaluation value and using the seismic wave velocity corresponding to the maximum value of the objective function as the dynamic correction velocity can obtain a more accurate dynamic correction velocity and improve the imaging quality of the formation profile.

[0049] Step S18: Process the gather record according to the target motion correction speed.

[0050] In this embodiment, the target dynamic correction speed can be used to provide the speed parameters used for dynamic correction to determine the optimal superimposed speed; the formation velocity can be obtained according to the target dynamic correction speed; the formation speed can be obtained according to the target dynamic correction speed The superimposed profile of the gather record, etc.

[0051] In this embodiment, by changing the seismic wave velocity, the average amplitude, similarity coefficient, and amplitude variance of the gather records corresponding to different seismic wave velocities are calculated to obtain the velocity evaluation values ​​corresponding to the different seismic wave velocities. , Taking the seismic wave velocity corresponding to the maximum value of the velocity evaluation value as the correct motion correction velocity. Taking into account the individual difference of each seismic trace record in the gather record, and since the average amplitude, correlation coefficient and amplitude variance of the gather record are all related to the dynamic correction speed of the gather record, the seismic wave velocity evaluation is used. Processing the gather records with the seismic wave velocity corresponding to the maximum value can improve the imaging effect of the stratum, especially the stratum sensitive to velocity.

[0052] This manual provides a scenario example, please refer to figure 2 , A gather record with a common center point in the work area, the number of shallow traces is small and it is more sensitive to velocity. If the selected seismic wave velocity and the formation velocity have a small error, some offsets cannot be leveled, and non-in-phase superposition occurs. Due to the small number of channels, the non-in-phase stacking of this part of the offset will have a greater impact on the stacking result.

[0053] In this scenario example, see image 3 Set the seismic wave velocity range from 0m/s to 6500m/s, the velocity interval is 500m/s, calculate the seismic wave superimposed velocity recorded by the gather, and generate velocity spectrum, according to the velocity spectrum, the seismic wave velocity is 2600m/s, corresponding When the in-phase axis time is 2.5s, the corresponding superimposed amplitude is the largest, and the initial dynamic correction speed is 2600m/s.

[0054] In this scenario example, see Figure 4 , Is the result of dynamic correction of the common center point gather record according to the initial dynamic correction speed.

[0055] In this scenario example, see Figure 5 , In order to obtain the initial dynamic correction speed of each gather record of the entire survey line by the FK interpolation method according to the initial dynamic correction speed and the in-phase axis time corresponding to the initial dynamic correction speed. The collection records are superimposed horizontally to obtain the superimposed profile. The event axis in the ellipse is messy, and the event axis in the rectangular frame is messy and jitters.

[0056] In this scenario example, according to the speed spectrum and the initial dynamic correction speed, the speed adjustment range is within ±5% of the initial dynamic correction speed, that is, the speed adjustment range is 2470m/s to 2730m/s, and the speed adjustment interval is 5m/s. s.

[0057] In this scenario example, the seismic wave velocity is selected according to the speed adjustment range and the speed adjustment interval to perform dynamic correction processing on the gather record. After the dynamic correction processing, the seismic trace in the gather record is processed. The recorded amplitudes are added in phase to calculate the average amplitude, similarity coefficient and amplitude variance of the seismic trace records in the gather record, and the average amplitude is multiplied by the similarity coefficient and then multiplied by the derivative of the amplitude variance to obtain the corresponding seismic wave velocity The velocity evaluation value is compared with the obtained velocity evaluation value, and the seismic wave velocity corresponding to the maximum value of the velocity evaluation value is taken as the correct motion correction velocity.

[0058] In this scenario example, it also includes calculating the correct motion correction speed in other time windows of the gather record, please refer to Image 6 , The seismic wave velocities around 1s, 1.7s, 2.3s, 3.5s, and 4.5s have all been corrected to a certain extent. The maximum value of the speed evaluation value in each time window length of the gather record is compared to obtain the maximum value of the maximum value of the speed evaluation value.

[0059] In this scenario example, see Figure 7 , The result of performing dynamic correction processing on the gather record according to the seismic wave velocity corresponding to the maximum value of the maximum value of the velocity evaluation value and the in-phase axis time corresponding to the seismic wave velocity.

[0060] In this scenario example, see Figure 8 According to the dynamic correction speed and the in-phase axis time corresponding to the dynamic correction speed, the dynamic correction speed of each gather record of the entire survey line is obtained through the spline interpolation method, and each gather record is horizontally superimposed according to the dynamic correction speed of each gather record Obtain the superimposed profile, and the event axis within the ellipse and rectangle changes naturally.

[0061] In this scenario example, after the speed is adjusted, the maximum value of the objective function is calculated, and the seismic wave velocity value corresponding to the maximum value of the objective function is used as the dynamic correction speed. The result of the dynamic correction processing on the gather record is compared with the initial dynamic correction speed The result of dynamic correction processing on gather records has achieved better results. After performing horizontal stacking separately, the result of stacking according to the initial dynamic correction speed finds that the event axis is messy and the event axis jitters; the result of horizontal stacking according to the seismic wave velocity corresponding to the maximum value of the objective function shows that the event axis changes naturally. A better stacking effect.

[0062] In one embodiment, the method further includes calculating the maximum value of the speed evaluation value in each time window in the gather record. Wherein, the time window is used to indicate the time length of a certain time interval in the gather record.

[0063] In this embodiment, the gather record may include many reflected waves from different interfaces, and the gather record may include multiple in-phase axes, and the in-phase axes may correspond to the reflected waves at zero offset. Reflection time. Therefore, the entire gather record can be divided into many time windows at a certain time interval, and these time windows can be selected according to the time corresponding to the in-phase axis at a certain time interval.

[0064] In this embodiment, calculating the maximum value of the objective function in each time window length may include dividing the gather record into many time windows at a certain time interval, and each time window may correspond to an in-phase axis, and By selecting and changing the seismic wave velocity value, the maximum value of the velocity evaluation value in each time window can be obtained.

[0065] In this embodiment, by calculating the maximum value of the objective function in each time window length, the correct dynamic correction speed of each reflected wave from the shallow layer to the deep layer in the entire gather record can be obtained.

[0066] In one embodiment, the method further includes comparing the maximum value of the speed evaluation value in each time window; wherein, the largest value among the maximum value of the speed evaluation value in each time window corresponds to The seismic wave velocity is the target motion correction velocity recorded in the gather.

[0067] In this embodiment, the comparing the maximum value of the objective function in the length of each time window may include, for different reflection layers, the time window selected in the gather record and the selected seismic wave velocity are both Variable, the size of the objective function can be changed with the selection of the window and the change of the seismic wave velocity, the maximum value of the objective function in each time window length is compared, and the maximum value of the maximum value of the objective function in each time window length is The corresponding seismic wave velocity may be the target motion correction velocity recorded by the gather.

[0068] In this embodiment, by comparing the maximum value of the velocity evaluation value in each time window, the seismic wave velocity corresponding to the largest value of the maximum value is used as the target motion correction velocity of the gather record. Since the greater the velocity evaluation value, the closer the corresponding velocity is to the true formation velocity. The seismic wave velocity value corresponding to the largest value among the maximum values ​​in each time window can be processed to the gather record to obtain higher Formation imaging effect.

[0069] In one embodiment, the method further includes: generating a velocity spectrum of the gather record according to the set seismic wave velocity range; wherein the velocity value corresponding to the maximum superimposed energy of the gather record in the velocity spectrum The initial motion correction speed recorded for the gather; the speed adjustment range is obtained according to the initial motion correction speed; the seismic wave velocity value is changed according to the speed adjustment range.

[0070] In this embodiment, the set seismic wave velocity range and velocity interval can be set according to the gather records and actual engineering needs. The larger the seismic wave velocity range is and the smaller the velocity interval is. The seismic wave superimposed velocity is more accurate.

[0071] In this embodiment, the step of generating the velocity spectrum of the gather record according to the set seismic wave velocity range may include selecting a seismic wave velocity, and dynamically correcting the gather record according to the seismic wave velocity, After dynamic correction, calculate the superimposed energy or average amplitude of the gather record; select the seismic wave velocity again according to the set seismic wave speed range and speed interval, and calculate the superimposed energy or average amplitude of the gather record accordingly, so that it can be generated The relationship curve between seismic wave velocity and the superimposed energy or average amplitude. Such a relationship curve can be the velocity spectrum recorded by the gather; in addition, it can also be based on the seismic wave velocity and the superimposed energy or average amplitude and the corresponding in-phase The axis time generates a three-dimensional data image, and the three-dimensional data image may be the velocity spectrum.

[0072] In this embodiment, the seismic wave velocity corresponding to the maximum superimposed energy or average amplitude in the velocity spectrum may be the initial motion correction velocity recorded by the gather.

[0073] In this embodiment, the initial motion correction speed may include that since the set seismic wave velocity range and velocity interval may have a certain error, the obtained motion correction speed may also have a certain error. The speed can be further adjusted to obtain a more accurate dynamic correction speed, so the speed value obtained from the velocity spectrum can be called the initial dynamic correction speed.

[0074] In this embodiment, the speed adjustment range may include: the initial dynamic correction speed and the energy cluster of the speed spectrum can be obtained according to the speed spectrum, and the energy cluster according to the initial dynamic correction speed and the speed spectrum The divergence of, the error range between the initial motion correction speed and the correct motion correction speed can be obtained, and the error range may be the speed adjustment range.

[0075] In this embodiment, changing the seismic wave velocity according to the velocity adjustment range may include, within the velocity adjustment range, selecting the seismic wave velocity according to engineering needs or randomly; within the velocity adjustment range, specifying The speed interval is to change the seismic wave speed according to the speed interval; within the speed adjustment range, a designated change rule is established, and the seismic wave speed is changed according to the rule. Specifically, for example, the speed interval is set to 10m/s, the speed adjustment range is 30m/s to 50m/s, 30m/s is selected, the first time is changed to 40m/s, and the second time is changed to 50m/s ; Set the previous change to 10m/s less than the next one. The speed adjustment range is from 0m/s to 60m/s. Select 0m/s, the first change is 10m/s, the second change is 30m/s, The third change is 60m/s.

[0076] In this embodiment, on the basis of obtaining the initial dynamic correction speed, further speed adjustment is performed, and the maximum value of the speed evaluation value is calculated to obtain the correct dynamic correction speed, which can not only improve the imaging quality of speed-sensitive formations, but also Increased computational efficiency.

[0077] In one embodiment, the step of changing the seismic wave speed value according to the speed adjustment range includes: setting a speed adjustment interval; and changing the seismic wave speed value according to the speed adjustment interval.

[0078] In this embodiment, the speed adjustment interval may include, within the speed adjustment range, taking a certain seismic wave velocity value to calculate the velocity evaluation value once, and then adding a fixed value for each seismic wave velocity value taken. Speed ​​value, this fixed speed value can adjust the interval for the speed. The speed adjustment interval may increase or decrease a certain speed according to the initial motion correction speed.

[0079] In this embodiment, changing the seismic wave velocity value according to the velocity adjustment range and the velocity adjustment interval can be within the velocity adjustment range, and can be in accordance with the magnitude of the velocity value from small to large or from large to small or in any manner. The seismic wave velocity value is changed according to the velocity adjustment interval, and the maximum value of the objective function is calculated. The seismic wave velocity value corresponding to the maximum value of the objective function may be the motion correction velocity recorded by the gather.

[0080] In this embodiment, on the basis of obtaining the initial dynamic correction velocity, further velocity adjustment is performed, the maximum value of the objective function is calculated, and the velocity adjustment interval is set to change the seismic wave velocity within the multiple adjustment range to obtain the correct dynamic correction. Speed ​​can not only improve the imaging quality of speed-sensitive formations, but also improve calculation efficiency.

[0081] In one embodiment, the gather record is dynamically corrected according to the target motion correction speed and the in-phase axis time corresponding to the target motion correction speed; wherein, the in-phase axis time is used to represent the seismic trace record The time corresponding to the connection of the same phase of the seismic wave.

[0082] In this embodiment, the in-phase axis time may include the time corresponding to the connection of the seismic wave in the same phase in the seismic trace record. The in-phase axis time may correspond to the time elapsed by a shot point in the case of zero offset for a reflective layer, and the seismic wave is received by the detector after the seismic wave is reflected by the reflective layer.

[0083] In this embodiment, the step of performing dynamic correction on the gather record according to the target motion correction speed and the in-phase axis time corresponding to the target motion correction speed may include, according to the target motion correction speed and the corresponding The in-phase axis time and the known offset of the seismic trace in the gather record are calculated to calculate the dynamic correction amount, and the dynamic correction processing can be performed on the gather record according to the dynamic correction amount.

[0084] In this embodiment, an accurate dynamic correction amount can be calculated according to the seismic wave velocity corresponding to the maximum value of the objective function and the event axis time, and the dynamic correction processing can be performed with the dynamic correction amount to improve the imaging quality of the formation.

[0085] In one embodiment, the method further includes, according to the target motion correction speed and the in-phase axis time corresponding to the target motion correction speed, obtaining the correct motion correction of other gather records on the entire survey line through interpolation The velocity and the in-phase axis time corresponding to the target motion correction speed; the seismic wave superimposed section of the survey line is obtained according to the target motion correction speed recorded in the gather and the in-phase axis time corresponding to the target motion correction speed.

[0086] In this embodiment, the interpolation method may include spline interpolation, FK interpolation, or weighted interpolation according to distance to perform interpolation operations. You can choose to obtain the dynamic correction speed of several gather records on the entire survey line. According to the interpolation operation, you can obtain the dynamic correction speed of other gather records for which the dynamic correction speed is not calculated.

[0087] In this embodiment, according to the motion correction speed corresponding to each gather record and the corresponding in-phase axis time, the motion correction value of each gather record in the work area can be calculated; the corresponding motion correction value can be calculated according to the motion correction value. Collect records for dynamic correction processing. After the dynamic correction processing, the gather records in the work area can be superimposed to obtain the superimposed track records of each gather record, and then the seismic wave superposition section of the entire work area can be obtained.

[0088] In this embodiment, the correct motion correction speed recorded in all gathers on the entire survey line can be obtained by interpolation, and it is not necessary to calculate the maximum value of the speed evaluation value recorded in all gathers. Obtain the correct dynamic correction speed of all gather records in the entire survey line, avoid troublesome calculations, and achieve the effect of quickly obtaining the correct dynamic correction speed of gather records in the work area.

[0089] See Picture 9 The embodiment of this specification also provides a processing device for gather records, a dynamic correction module; for dynamic correction of gather records according to a seismic wave velocity value; a first calculation module; for calculating the dynamic correction of the trace The average amplitude, similarity coefficient and amplitude variance of the seismic trace records in the gather record; wherein the similarity coefficient is used to indicate the similarity of the waveforms of the seismic trace records in the gather record; wherein, the gather record includes At least one seismic trace record; a second calculation module; used to calculate a velocity evaluation value based on the average amplitude, similarity coefficient, and amplitude variance; the velocity evaluation value is used to indicate how close the seismic wave velocity value is to the formation velocity; The third calculation module; for changing the seismic wave velocity value to obtain the maximum value of the velocity evaluation value; wherein the seismic wave velocity value corresponding to the maximum value of the velocity evaluation value is the target motion recorded by the gather Correction speed; processing module; used to process the gather record according to the target motion correction speed.

[0090] The devices or units explained in the above embodiments may be specifically implemented by computer chips or entities, or products with certain functions. For the convenience of description, when describing the above device, the functions are divided into various modules and described separately. Of course, when implementing this specification, the functions of each module can be implemented in the same one or more software and/or hardware.

[0091] It can be known from the description of the above embodiments that those skilled in the art can also understand that the various illustrative logic blocks, modules, and steps listed in the embodiments of this specification can be implemented by hardware, software or a combination of both. Whether it is implemented by hardware or software depends on the specific application and the design requirements of the entire system. Those skilled in the art can use various methods to implement the described functions for each specific application, but such implementation should not be construed as going beyond the protection scope of the present invention.

[0092] The various illustrative modules described in the embodiments of this specification can be implemented by general-purpose processors, digital signal processors, application specific integrated circuits, field programmable gate arrays or other programmable logic devices, discrete hardware components, or any combination of the foregoing The design to achieve or operate the described function. The general-purpose processor may be a microprocessor. Optionally, the general-purpose processor may also be any conventional processor, controller, microcontroller, or state machine. The processor may also be implemented by a combination of computing devices, such as a digital signal processor, multiple microprocessors, one or more microprocessors combined with a digital signal processor core, or other similar configurations.

[0093] The embodiment of the present specification also provides a storage medium, the storage medium stores computer program instructions, when the computer program instructions are executed: perform dynamic correction of gather records according to a seismic wave velocity value; calculate the dynamic correction The average amplitude, similarity coefficient and amplitude variance of the seismic trace records in the gather record; wherein the similarity coefficient is used to indicate the degree of similarity of the waveforms of the seismic trace records in the gather record; wherein, the gather The record includes at least one seismic trace record; the velocity evaluation value is calculated according to the average amplitude, similarity coefficient and amplitude variance; wherein, the velocity evaluation value is used to indicate the proximity of the seismic wave velocity value to the formation velocity; The seismic wave velocity value is used to obtain the maximum value of the velocity evaluation value; wherein the seismic wave velocity value corresponding to the maximum value of the velocity evaluation value is the target motion correction speed recorded by the gather; Speed ​​processing of the gather records.

[0094] In this embodiment, the storage medium includes, but is not limited to, random access memory (Random Access Memory, RAM), read-only memory (Read-Only Memory, ROM), cache (Cache), hard disk (Hard Disk Drive, HDD), or Memory Card (Memory Card).

[0095] For the storage medium provided in this embodiment, the functions and effects achieved when the program instructions are executed can be explained by reference to other embodiments.

[0096] The implementation of this specification provides a processing method, device and storage medium for gather recording. A seismic wave velocity is used to dynamically correct the gather records. After the dynamic correction, the average amplitude, amplitude variance and similarity coefficient of the seismic trace records in the gather record are calculated, and the said average amplitude, amplitude variance and similarity coefficient are calculated according to the average amplitude, amplitude variance and similarity coefficient. The velocity evaluation value is calculated by changing the seismic wave velocity to calculate the maximum value of the velocity evaluation value, the seismic wave velocity corresponding to the maximum value of the velocity evaluation value is regarded as the correct dynamic correction velocity, and the correct dynamic correction velocity Gather records are processed. The seismic wave velocity value corresponding to the maximum value of the value is the target motion correction velocity of the gather record to correct the gather record. Taking into account the individual differences of each seismic track record in the gather record, and since the average amplitude, correlation coefficient and amplitude variance of the gather record are all related to the dynamic correction speed of the gather record, the calculations calculated by the above method The target dynamic correction speed can be obtained closer to the true speed of the formation, and processing the gather record by using the seismic wave speed corresponding to the seismic wave speed evaluation value can improve the imaging effect of the formation, especially the formation that is more sensitive to speed.

[0097] The steps of the methods or algorithms described in the embodiments of this specification can be directly embedded in hardware, software and hardware modules executed by a processor, or a combination of the two. The software module can be stored in RAM memory, flash memory, ROM memory, EPROM memory, registers, hard disk, removable disk, CD-ROM or any other storage medium in the field.

[0098] The aforementioned functions described in the embodiments of this specification can be implemented in hardware, software, firmware, or any combination of the three. If implemented in software, these functions can be stored on a computer-readable medium or transmitted on the computer-readable medium in the form of one or more instructions or codes. Computer-readable media include computer storage media and communication media that facilitate the transfer of computer programs from one place to another. The storage medium can be any available medium that can be accessed by a general-purpose or special computer. For example, such computer-readable media may include but are not limited to RAM, ROM, EPROM, CD-ROM or other optical disk storage, magnetic disk storage or other magnetic storage devices, or any other device that can be used to carry or store instructions or data structures and Other program code media that can be read by general or special computers, or general or special processors. In addition, any connection can be appropriately defined as a computer-readable medium, for example, if the software is from a website site, server, or other remote source through a coaxial cable, fiber optic cable, twisted pair, digital subscriber line (DSL) Or wireless transmission such as infrared, wireless and microwave is also included in the defined computer readable medium. The discs and magnetic discs include compressed discs, laser discs, optical discs, DVDs, floppy discs, and Blu-ray discs. Discs usually copy data magnetically, while discs usually copy data optically by laser. The above combination can also be contained in a computer readable medium.

[0099] The various embodiments in this specification are described in a progressive manner, and the same or similar parts between the various embodiments can be referred to each other, and each embodiment focuses on the differences from other embodiments. In particular, for the device embodiment, since it is basically similar to the method embodiment, the description is relatively simple, and the relevant part can refer to the description of the method embodiment.

[0100] Although this specification has been described through implementations, those of ordinary skill in the art know that there are many variations and changes in this specification without departing from the spirit of this specification, and it is hoped that the appended claims include these variations and changes without departing from the spirit of this specification.