CN-122017976-A - Mixed production seismic data separation method and device

CN122017976ACN 122017976 ACN122017976 ACN 122017976ACN-122017976-A

Abstract

The embodiment of the application discloses a method and a device for separating mixed mining seismic data. The application provides a three-dimensional FKK-domain grid-variable sparse inversion technology, which dynamically adjusts the scale of data grid blocks in an iterative process, matches the physical rule of signal extraction, is more in line with the essence of signal extraction from low frequency to high frequency from strong axis to weak axis during inversion, thereby obtaining better separation effect, provides a self-adaptive air channel compensation technology, fills air channels in a sparse transformation domain through regular interpolation, avoids abnormal attenuation of energy, effectively eliminates the phenomenon of whitening of processed data caused by changing the viewing defect channels, improves the effect, simultaneously selects a three-dimensional FKK transformation domain, combines signal coherence identification, transformation reversibility and calculation efficiency, finally realizes the improvement of separation precision and efficiency, obviously improves field production operation efficiency, shortens the operation period and reduces the exploration cost.

Inventors

ZHANG PENG
WU XUGUANG
WANG HAIKUN
Chen bang
WANG WEI

Assignees

中海油田服务股份有限公司

Dates

Publication Date: 20260512
Application Date: 20260311

Claims (12)

1. A method of commingled seismic data separation comprising: acquiring seismic data to be separated and mixed mining, and performing air channel detection on the seismic data to be separated and mixed mining, wherein the seismic data to be separated and mixed mining are synchronous source mixed mining data in a time domain; Transforming the to-be-separated mixed mining seismic data to FKK domains, performing channel region interpolation processing on the channels according to the FKK domain transformed seismic data corresponding to the effective channels in the preset range of each channel, and performing FKK inverse transformation processing on the to-be-separated mixed mining seismic data subjected to interpolation processing to obtain to-be-separated mixed mining seismic data subjected to channel compensation; Carrying out FKK-domain sparse transform on the spatial compensated to-be-separated mixed seismic data, and constructing an initial model without aliasing data according to the FKK-domain sparse transformed to-be-separated mixed seismic data; performing phased FKK-domain grid-change sparse inversion iterative processing on the initial model without the aliasing data, and updating the current model without the aliasing data and a sparse shrinkage threshold value after each iterative processing until a preset stop condition is met; And performing FKK inverse transformation processing on the model data corresponding to the non-aliased data model meeting the stop condition to obtain a seismic data separation result.
2. The method of claim 1, wherein the channel detection of the to-be-separated commingled seismic data further comprises: Carrying out amplitude statistics in a preset local time window on each seismic trace in the common detection point trace set to obtain root mean square amplitude of the corresponding seismic trace; and judging whether the root mean square amplitude is smaller than a preset root mean square amplitude threshold value, if so, determining the seismic channel as an empty channel and recording empty channel information.
3. The method of claim 2, wherein the performing the channel region interpolation on the channel according to the FKK-domain transformed seismic data corresponding to the effective channel within the preset range of each channel further comprises: Locating FKK empty channels in the seismic data to be separated and mixed according to the empty channel information; For any channel, calculating an interpolation weight coefficient according to the distance from each effective channel to the channel in the preset channel range, and carrying out channel region interpolation processing on seismic data after FKK-domain transformation corresponding to each effective channel in the preset channel range according to the interpolation weight coefficient.
4. The method of any of claims 1-3, wherein the constructing an aliased data-free initial model from FKK-domain sparsely transformed to-be-separated aliased seismic data further comprises: Calculating an initial sparse shrinkage threshold according to the FKK-domain sparse transformed to-be-separated mixed mining seismic data; And constructing an aliased data model according to the to-be-separated aliased seismic data subjected to FKK-domain sparse transformation, and performing pseudo-unmixed calculation on the aliased data model based on an initial sparse shrinkage threshold value to obtain an initial model without aliased data.
5. The method of any of claims 1-3, wherein performing a phased FKK-domain grid-change sparse inversion iterative process on the non-aliased data initial model, updating the current non-aliased data model and a sparse shrink threshold after each iterative process until a preset stop condition is met further comprises: The method comprises the steps of carrying out multi-stage FKK-domain sparse inversion iteration processing on an initial model of the non-aliasing data, dividing the iteration process into a plurality of iteration stages according to a change threshold of iteration residual errors, adopting a grid scale which is reduced step by step in each iteration stage to respectively adapt to seismic signal sparse extraction of corresponding frequency bands and amplitude characteristics, adopting a sparse shrinkage threshold which is reduced step by step in each iteration stage, and updating the current non-aliasing data model by each iteration based on the non-aliasing data model obtained by the previous iteration until a preset stop condition is met.
6. The method of claim 5, wherein the method further comprises calculating a current iteration objective function after the current non-aliased data model is updated; and if the current iteration objective function is smaller than the last iteration objective function, updating the aliasing noise estimation value based on the current aliasing-free data model.
7. The method of any of claims 1-3, wherein the updating the sparsity contraction threshold further comprises: After each iteration of each phase, the sparse shrink threshold is updated according to the following formula: Where i is the current iteration number, k is the attenuation sparsity, T i is the sparse shrinkage threshold corresponding to the ith iteration of the corresponding stage, and T i-1 is the sparse shrinkage threshold corresponding to the i-1 th iteration of the corresponding stage.
8. A method according to any one of claims 1-3, wherein prior to the channel detection of the to-be-separated commingled seismic data, the method further comprises: And denoising pretreatment is carried out on the to-be-separated mixed mining seismic data to suppress background low-frequency strong-amplitude surge noise, so as to obtain pretreated to-be-separated mixed mining seismic data.
9. A commingled seismic data separation device, comprising: The acquisition module is suitable for acquiring the to-be-separated mixed mining seismic data, wherein the to-be-separated mixed mining seismic data is synchronous source mixed mining data in a time domain; The detection module is suitable for carrying out space channel detection on the seismic data to be separated and mixed; The interpolation processing module is suitable for transforming the to-be-separated mixed mining seismic data to FKK domains, carrying out channel region interpolation processing on the channels according to FKK domain transformed seismic data corresponding to the effective channels in the preset range of each channel, and carrying out FKK inverse transformation processing on the to-be-separated mixed mining seismic data subjected to interpolation processing to obtain to-be-separated mixed mining seismic data subjected to channel compensation; The construction module is suitable for carrying out FKK-domain sparse conversion treatment on the channel-compensated to-be-separated mixed mining seismic data, and constructing an initial model without the mixed mining data according to the FKK-domain sparse converted to-be-separated mixed mining seismic data; The variable grid sparse inversion processing module is suitable for executing phased FKK-domain variable grid sparse inversion iterative processing on the non-aliased data initial model, and updating the current non-aliased data model and a sparse shrinkage threshold value after each iterative processing until a preset stop condition is met; And the inverse transformation processing module is suitable for executing FKK inverse transformation processing on the model data corresponding to the non-aliased data model meeting the stop condition to obtain a seismic data separation result.
10. A computing device comprising a processor, a memory, a communication interface, and a communication bus, the processor, the memory, and the communication interface completing communication with each other over the communication bus; The memory is configured to store at least one executable instruction that causes the processor to perform operations corresponding to the method for producing seismic data separation of any one of claims 1-8.
11. A computer storage medium having stored therein at least one executable instruction for causing a processor to perform operations corresponding to the method of commingled seismic data separation of any one of claims 1-8.
12. A computer program product comprising at least one executable instruction that causes a processor to perform operations corresponding to the method of commingled seismic data separation of any one of claims 1-8.

Description

Mixed production seismic data separation method and device Technical Field The embodiment of the application relates to the technical field of seismic exploration, in particular to a method and a device for separating mixed mining seismic data. Background Along with the increasing development degree of offshore oil exploration, the technology of 'two wide and one high' seismic acquisition of the submarine nodes is a necessary trend, the submarine node acquisition generally needs high gun channel density to realize high coverage and wide azimuth acquisition, and the traditional sequential gun excitation acquisition mode is low in construction efficiency and high in operation cost. The synchronous source efficient aliasing acquisition is a great assistance for accelerating the promotion of the large-scale application of the submarine nodes, and the existing application examples show that better quality can be obtained under the condition of the same cost, and the cost is less under the condition of the same quality. The mixed data separation technology is one of core technologies for realizing mixed acquisition, and greatly determines the selection of a field operation mode and the quality of seismic data. At present, a mixed mining data separation scheme based on sparse inversion and denoising is mainly adopted for data separation, however, after separation based on the method, the problem of poor separation effect exists, and therefore, a mixed mining seismic data separation scheme capable of effectively improving the separation effect is needed. Disclosure of Invention The present application has been made in view of the above problems, and it is therefore an object of the present application to provide a method, apparatus, computing device, computer storage medium and computer program product for commingled seismic data separation that overcomes or at least partially solves the above problems. According to an aspect of an embodiment of the present application, there is provided a method for separating seismic data from a mixed production, including: Acquiring seismic data to be separated and mixed mining, and performing air channel detection on the seismic data to be separated and mixed mining, wherein the seismic data to be separated and mixed mining are synchronous source mixed mining data in a time domain; Transforming the to-be-separated mixed mining seismic data to FKK domains, carrying out channel region interpolation processing on the channels according to the FKK domain transformed seismic data corresponding to the effective channels in the preset range of each channel, and carrying out FKK inverse transformation processing on the to-be-separated mixed mining seismic data subjected to interpolation processing to obtain the to-be-separated mixed mining seismic data subjected to channel compensation; Carrying out FKK-domain sparse transform on the spatial compensated to-be-separated mixed seismic data, and constructing an initial model without aliasing data according to the FKK-domain sparse transformed to-be-separated mixed seismic data; Performing phased FKK domain grid-change sparse inversion iterative processing on the initial model without the aliasing data, and updating the current model without the aliasing data and a sparse shrinkage threshold value after each iterative processing until a preset stop condition is met; And performing FKK inverse transformation processing on the model data corresponding to the non-aliased data model meeting the stop condition to obtain a seismic data separation result. Further, performing the space-channel detection on the seismic data to be separated and mixed with the mining further comprises: Carrying out amplitude statistics in a preset local time window on each seismic trace in the common detection point trace set to obtain root mean square amplitude of the corresponding seismic trace; and judging whether the root mean square amplitude is smaller than a preset root mean square amplitude threshold value, if so, determining the seismic channel as an empty channel and recording empty channel information. Further, performing the channel region interpolation processing on the channels according to the FKK-domain transformed seismic data corresponding to the effective channels in the preset range of each channel further comprises: Locating FKK empty channels in the seismic data to be separated and mixed according to the empty channel information; For any channel, calculating an interpolation weight coefficient according to the distance from each effective channel to the channel in a preset channel range, and carrying out channel region interpolation processing on seismic data after FKK-domain transformation corresponding to each effective channel in the preset channel range according to the interpolation weight coefficient. Further, constructing an initial model of the aliasing-free data according to the FKK-domain sparsely transformed to-be-separated aliased s