CN108107485B - The appraisal procedure and device of seismic data multiple wave pollution level - Google Patents
The appraisal procedure and device of seismic data multiple wave pollution level Download PDFInfo
- Publication number
- CN108107485B CN108107485B CN201711144555.3A CN201711144555A CN108107485B CN 108107485 B CN108107485 B CN 108107485B CN 201711144555 A CN201711144555 A CN 201711144555A CN 108107485 B CN108107485 B CN 108107485B
- Authority
- CN
- China
- Prior art keywords
- vsp
- data
- deconvolution
- depth
- wave
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V13/00—Manufacturing, calibrating, cleaning, or repairing instruments or devices covered by groups G01V1/00 – G01V11/00
Landscapes
- Engineering & Computer Science (AREA)
- Manufacturing & Machinery (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The embodiment of the present application provides the appraisal procedure and device of a kind of seismic data multiple wave pollution level, this method comprises: obtaining in work area before the deconvolution of VSP data VSP upgoing wave data after VSP down going wave data and deconvolution；It determines and once involves multiple wave before the deconvolution in VSP down going wave data at any time or the multiple wave of change in depth pollution discretization curve, and determine after the deconvolution VSP upgoing wave data at any time or the bed boundary response sequence of change in depth；Pollute discretization curve and the bed boundary response sequence according to the multiple wave, obtain the work area seismic data at any time or the multiple wave pollution level of change in depth.The embodiment of the present application can realize the qualitative assessment to seismic data multiple wave pollution level.
Description
Technical Field
The application relates to the technical field, in particular to a method and a device for evaluating the multiple pollution degree of seismic data.
Background
At present, more analysis and suppression methods are used for seismic data multiples, and the representative methods are as follows: in 1973, Cassano et al proposed an optimal filtering superposition method, which solved the filtering factors of each superposition channel by using a least square method to make the superposition reach the optimal suppression multiple so as to optimally approximate the primary wave. The deconvolution suppression method proposed by Lokshtanov et al is to calculate a deconvolution operator in frequency and slowness domains on the basis of a one-dimensional and two-dimensional reflection model, and further perform deconvolution suppression on multiples. Doicin et al propose a specific pegleg attenuation method, which uses a spatial matrix filtering method in the f-x domain to obtain a model of multiples.
The degree of multiple pollution in the seismic data directly affects the quality of the seismic data. However, although various methods of multiple analysis and suppression are currently available, as described above, there are rare solutions for quantitative analysis and evaluation of the degree of multiple contamination in seismic data. Therefore, a technical scheme for quantitatively evaluating the multiple pollution degree of seismic data is needed.
Disclosure of Invention
The embodiment of the application aims to provide a method and a device for evaluating the multiple pollution degree of seismic data, so as to realize quantitative evaluation of the multiple pollution degree of the seismic data.
In order to achieve the above object, in one aspect, an embodiment of the present application provides a method for evaluating a degree of seismic data multiple contamination, including:
acquiring VSP (vertical seismic profile) downlink wave data before deconvolution and VSP uplink wave data after deconvolution of VSP data in a work area;
determining a multiple pollution discretization curve of the first wave and the multiple changing along with time or depth in the VSP downlink wave data before deconvolution, and determining a stratum interface response sequence of the VSP uplink wave data after deconvolution along with time or depth;
and acquiring the multiple pollution degree of the seismic data of the work area along with the time or depth change according to the multiple pollution discretization curve and the stratum interface response sequence.
Preferably, the acquiring VSP downlink wave data before deconvolution and VSP uplink wave data after deconvolution of VSP data in the work area includes:
decoding VSP data in the work area;
carrying out polarization rotation on the decoded VSP data to obtain a longitudinal wave field component and a tangential wave field component;
performing wave field separation processing on the longitudinal wave field component to obtain VSP downlink wave data before deconvolution of the VSP data; and carrying out deconvolution and wave field separation processing on the tangential wave field component in sequence to obtain the deconvoluted VSP upgoing wave data of the VSP data.
Preferably, the determining a discretization curve of contamination of multiples of the first and the multiple with time or depth in the pre-deconvolution VSP downgoing data includes:
carrying out absolute value processing on the VSP downlink wave data before deconvolution;
determining the amplitude and the peak position of VSP downlink wave data before deconvolution after absolute value processing;
and normalizing the peak amplitude, and fitting the peak position and the normalized peak amplitude to obtain a multiple pollution discretization curve of multiple changes along with time or depth in the VSP downlink wave data before deconvolution.
Preferably, the determining a formation interface response sequence of the post-deconvolution VSP upgoing data with time or depth includes:
carrying out absolute value processing on the VSP uplink wave data after deconvolution;
and extracting the amplitude and the peak position of the deconvoluted VSP uplink wave data after the absolute value processing to obtain a stratum interface response sequence of the deconvoluted VSP uplink wave data changing with time or depth.
Preferably, the obtaining of the multiple contamination degree of the seismic data of the work area along with the time or depth change according to the multiple contamination discretization curve and the stratum interface response sequence includes:
and multiplying each stratum interface response value in the stratum interface response sequence by the multiple pollution discretization curve respectively, and correspondingly placing the obtained result at the position of the corresponding stratum interface response value and below the corresponding stratum interface response value, thereby obtaining the multiple pollution degree of the seismic data of the work area along with the change of time or depth.
On the other hand, this application embodiment still provides a seismic data multiple pollution degree's evaluation device, includes:
the uplink and downlink data acquisition module is used for acquiring VSP downlink wave data before deconvolution and VSP uplink wave data after deconvolution of the VSP data in the work area;
a curve and sequence acquisition module for determining a multiple pollution discretization curve of primary waves and multiple changes along with time or depth in the VSP downlink wave data before deconvolution and determining a stratum interface response sequence of the VSP uplink wave data after deconvolution along with time or depth changes;
and the pollution degree determining module is used for acquiring the multiple pollution degree of the seismic data of the work area along with the change of time or depth according to the multiple pollution discretization curve and the stratum interface response sequence.
Preferably, the acquiring VSP downlink wave data before deconvolution and VSP uplink wave data after deconvolution of VSP data in the work area includes:
decoding VSP data in the work area;
carrying out polarization rotation on the decoded VSP data to obtain a longitudinal wave field component and a tangential wave field component;
performing wave field separation processing on the longitudinal wave field component to obtain VSP downlink wave data before deconvolution of the VSP data; and carrying out deconvolution and wave field separation processing on the tangential wave field component in sequence to obtain the deconvoluted VSP upgoing wave data of the VSP data.
Preferably, the determining a discretization curve of contamination of multiples of the first and the multiple with time or depth in the pre-deconvolution VSP downgoing data includes:
carrying out absolute value processing on the VSP downlink wave data before deconvolution;
determining the amplitude and the peak position of VSP downlink wave data before deconvolution after absolute value processing;
and normalizing the peak amplitude, and fitting the peak position and the normalized peak amplitude to obtain a multiple pollution discretization curve of multiple changes along with time or depth in the VSP downlink wave data before deconvolution.
Preferably, the determining a formation interface response sequence of the post-deconvolution VSP upgoing data with time or depth includes:
carrying out absolute value processing on the VSP uplink wave data after deconvolution;
and extracting the amplitude and the peak position of the deconvoluted VSP uplink wave data after the absolute value processing to obtain a stratum interface response sequence of the deconvoluted VSP uplink wave data changing with time or depth.
Preferably, the obtaining of the multiple contamination degree of the seismic data of the work area along with the time or depth change according to the multiple contamination discretization curve and the stratum interface response sequence includes:
and multiplying each stratum interface response value in the stratum interface response sequence by the multiple pollution discretization curve respectively, and correspondingly placing the obtained result at the position of the corresponding stratum interface response value and below the corresponding stratum interface response value, thereby obtaining the multiple pollution degree of the seismic data of the work area along with the change of time or depth.
In another aspect, an embodiment of the present application further provides another apparatus for evaluating a degree of seismic data multiple contamination, including a memory, a processor, and a computer program stored on the memory, where the computer program, when executed by the processor, performs the following steps:
acquiring VSP downlink wave data before deconvolution and VSP uplink wave data after deconvolution of the VSP data in a work area;
determining a multiple pollution discretization curve of the first wave and the multiple changing along with time or depth in the VSP downlink wave data before deconvolution, and determining a stratum interface response sequence of the VSP uplink wave data after deconvolution along with time or depth;
and acquiring the multiple pollution degree of the seismic data of the work area along with the time or depth change according to the multiple pollution discretization curve and the stratum interface response sequence.
According to the technical scheme provided by the embodiment of the application, firstly, VSP downlink wave data before deconvolution and VSP uplink wave data after deconvolution of VSP data in a work area are obtained; then determining a multiple pollution discretization curve of the primary wave and the multiple changing with time or depth in the VSP downlink wave data before deconvolution, and determining a stratum interface response sequence of the VSP uplink wave data after deconvolution changing with time or depth; and finally, acquiring the multiple pollution degree of the seismic data of the work area along with the change of time or depth according to the multiple pollution discretization curve and the stratum interface response sequence. Therefore, the multi-wave pollution degree of seismic data in a work area at different depths (or different time) is calculated quantitatively by utilizing the characteristics that the multi-wave in the VSP data can be tracked in the whole process and a deconvolution operator is known, and the method has great practical value in the aspects of quality control and the like of the seismic data.
Drawings
In order to more clearly illustrate the embodiments of the present application or the technical solutions in the prior art, the drawings needed to be used in the description of the embodiments or the prior art will be briefly introduced below, it is obvious that the drawings in the following description are only some embodiments described in the present application, and for those skilled in the art, other drawings can be obtained according to the drawings without any creative effort. In the drawings:
FIG. 1 is a flow chart of a method for evaluating a degree of seismic data multiple contamination according to an embodiment of the present disclosure;
FIG. 2 is a schematic diagram of pre-deconvolution VSP downlink data obtained in an embodiment of the present application;
FIG. 3 is a schematic diagram of a multiple contamination discretization curve obtained in an embodiment of the subject application;
FIG. 4 is a schematic diagram of post-deconvolution VSP upgoing data acquired in an embodiment of the present application;
FIG. 5 is a schematic illustration of a sequence of formation boundary responses obtained in an embodiment of the present application;
fig. 6 is a schematic diagram of a multiple contamination discretization curve obtained in an embodiment of the present application and a stratum boundary response sequence operation (an upper diagram is a conventional convolution calculation method, and a lower diagram is a calculation method of the present embodiment);
FIG. 7 is a schematic diagram of the evaluation results of the contamination level of VSP multiples obtained in one embodiment of the present application (dark color is the effective wave component and light color is the multiple component);
FIG. 8 is a block diagram of an apparatus for evaluating a multiple contamination level of seismic data according to an embodiment of the present disclosure;
fig. 9 is a block diagram showing a structure of an apparatus for evaluating a degree of contamination of multiples of seismic data according to another embodiment of the present invention.
Detailed Description
In order to make those skilled in the art better understand the technical solutions in the present application, the technical solutions in the embodiments of the present application will be clearly and completely described below with reference to the drawings in the embodiments of the present application, and it is obvious that the described embodiments are only a part of the embodiments of the present application, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present application.
VSP is a geophysical technique that is intermediate between seismic and well logging methods, and is known to have a higher signal-to-noise ratio and resolution relative to conventional seismic observation methods. In the process of realizing the method, the inventor finds that for the multiples, a down-going wave and an up-going wave can be observed simultaneously by adopting a VSP seismic observation method, and the multiples in a VSP data wave field have the characteristics of being tracked in the whole process and known by a deconvolution operator (down-going wave), so that the method has unique advantages for tracking the occurrence and propagation conditions of the multiples. In view of this, VSP data in the work area can be obtained by performing VSP seismic observation in the work area, and then the multiple wave pollution degree of the seismic data in the work area at different depths (or different times) can be quantitatively calculated by using the characteristics that multiple waves in the VSP data can be completely tracked and known by a deconvolution operator. The method has great practical value in the aspects of quality control and the like of the seismic data.
On the basis of the above theory, referring to fig. 1, an evaluation method for a multiple contamination degree of seismic data according to an embodiment of the present application may include the following steps:
s101, VSP downlink wave data before deconvolution and VSP uplink wave data after deconvolution of the VSP data in the work area are obtained.
In some embodiments of the present application, the acquiring the pre-deconvolution VSP downlink data and the post-deconvolution VSP uplink data of the VSP data in the work area may include the following steps:
1) decoding the VSP data in the work area, namely performing data decoding on the VSP data in the work area;
2) carrying out polarization rotation on the decoded VSP data to obtain a longitudinal wave field component and a tangential wave field component;
3) performing wave field separation processing on the longitudinal wave field component to obtain pre-deconvolution VSP downlink wave data of the VSP data, for example, as shown in FIG. 2; and performing deconvolution and wave field separation processing on the tangential wave field component in sequence to obtain deconvoluted VSP upgoing wave data of the VSP data, for example, as shown in fig. 4. The deconvolution processing is to extract a deconvolution operator by using VSP (vertical seismic profiling) downlink wave data before deconvolution, perform deconvolution processing on the tangential wave field component according to the deconvolution operator, and the data after deconvolution processing does not contain a multiple wave component any more.
S102, determining a multiple pollution discretization curve of the primary wave and the multiple changing along with time or depth in the VSP downlink wave data before deconvolution, and determining a stratum interface response sequence of the VSP uplink wave data after deconvolution along with time or depth.
In some embodiments of the present application, the determining a discretization curve of a pollution of multiples of a primary wave and a multiple varying with time or depth in the pre-deconvolution VSP downgoing data may include the following steps:
1) carrying out absolute value processing on the VSP downlink wave data before deconvolution;
2) determining the amplitude and the peak position of VSP downlink wave data before deconvolution after absolute value processing; these amplitude and peak position information directly reflect the relative relationship between primary and multiple waves in seismic data.
3) And normalizing the peak amplitude, and fitting the peak position and the normalized peak amplitude to obtain a multiple pollution discretization curve of multiple changes with time or depth in the VSP downlink wave data before deconvolution, for example, as shown in FIG. 3. The multiple pollution discretization curve may be recorded in a discrete data sequence, the data interval may be VSP data time sampling interval, and the data value may be the fitted amplitude value.
In some embodiments of the present application, the determining a formation interface response sequence of the post-deconvolution VSP upgoing data over time or depth may include the following steps:
1) carrying out absolute value processing on the VSP uplink wave data after deconvolution;
2) and extracting the amplitude and the peak position of the deconvoluted VSP upgoing wave data after the absolute value processing to obtain a formation interface response sequence of the deconvoluted VSP upgoing wave data changing with time or depth, for example, as shown in FIG. 5. Wherein, the stratum interface response sequence is different from the seismic reflection coefficient sequence, is a comprehensive response effect, has no positive and negative difference, and the data interval of the sequence can be the VSP data time sampling interval.
S103, acquiring the multiple pollution degree of the seismic data of the work area along with the time or depth change according to the multiple pollution discretization curve and the stratum interface response sequence.
In some embodiments of the present application, the obtaining a multiple contamination degree of the seismic data of the work area over time or depth according to the multiple contamination discretization curve and the formation interface response sequence includes:
and multiplying each stratum interface response value in the stratum interface response sequence by the multiple pollution discretization curve respectively, and correspondingly placing the obtained result at the position of the corresponding stratum interface response value and below the corresponding stratum interface response value, so as to obtain the multiple pollution degree of the seismic data of the work area along with the change of time or depth, for example, as shown in fig. 7.
In an exemplary embodiment of the present application, as shown in fig. 6, the multiple contamination level calculation process is similar to the convolution calculation, in which: the upper graph is a conventional convolution calculation; the following diagram is the operation applied in this example, and in the diagram, the following diagram sequence 1 represents a multiple pollution discretization curve, and the following diagram sequence 2 represents a formation interface response sequence. It can be seen that the main difference from the convolution operation is that the first sample of sequence 1 is aligned with the sample to be calculated for sequence 2, i.e. it corresponds to only half the sequence of sequence 1 in the convolution operation.
While the process flows described above include operations that occur in a particular order, it should be appreciated that the processes may include more or less operations that are performed sequentially or in parallel (e.g., using parallel processors or a multi-threaded environment).
Referring to fig. 8, an apparatus for evaluating a degree of seismic data multiple contamination according to an embodiment of the present disclosure may include:
the uplink and downlink data acquisition module 81 may be configured to acquire VSP downlink wave data before deconvolution and VSP uplink wave data after deconvolution of the VSP data in the work area;
a curve and sequence obtaining module 82, configured to determine a multiple pollution discretization curve of a primary wave and a multiple that change with time or depth in the VSP downlink data before deconvolution, and determine a formation interface response sequence of the VSP uplink data that changes with time or depth after deconvolution;
and the pollution degree determining module 83 may be configured to obtain the multiple pollution degree of the seismic data of the work area, which changes with time or depth, according to the multiple pollution discretization curve and the stratum interface response sequence.
Referring to fig. 9, another apparatus for evaluating a degree of seismic data multiple contamination according to an embodiment of the present application may include a memory, a processor, and a computer program stored in the memory, where the computer program is executed by the processor to perform the following steps:
acquiring VSP downlink wave data before deconvolution and VSP uplink wave data after deconvolution of the VSP data in a work area;
determining a multiple pollution discretization curve of the first wave and the multiple changing along with time or depth in the VSP downlink wave data before deconvolution, and determining a stratum interface response sequence of the VSP uplink wave data after deconvolution along with time or depth;
and acquiring the multiple pollution degree of the seismic data of the work area along with the time or depth change according to the multiple pollution discretization curve and the stratum interface response sequence.
The apparatus of the above embodiment of the present application corresponds to the method of the above embodiment of the present application, and therefore, for details about the apparatus of the above embodiment of the present application, please refer to the method of the above embodiment of the present application, which is not described herein again.
For convenience of description, the above devices are described as being divided into various units by function, and are described separately. Of course, the functionality of the units may be implemented in one or more software and/or hardware when implementing the present application.
The present invention is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the invention. It will be understood that each flow and/or block of the flow diagrams and/or block diagrams, and combinations of flows and/or blocks in the flow diagrams and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
In a typical configuration, a computing device includes one or more processors (CPUs), input/output interfaces, network interfaces, and memory.
The memory may include forms of volatile memory in a computer readable medium, Random Access Memory (RAM) and/or non-volatile memory, such as Read Only Memory (ROM) or flash memory (flash RAM). Memory is an example of a computer-readable medium.
Computer-readable media, including both non-transitory and non-transitory, removable and non-removable media, may implement information storage by any method or technology. The information may be computer readable instructions, data structures, modules of a program, or other data. Examples of computer storage media include, but are not limited to, phase change memory (PRAM), Static Random Access Memory (SRAM), Dynamic Random Access Memory (DRAM), other types of Random Access Memory (RAM), Read Only Memory (ROM), Electrically Erasable Programmable Read Only Memory (EEPROM), flash memory or other memory technology, compact disc read only memory (CD-ROM), Digital Versatile Discs (DVD) or other optical storage, magnetic cassettes, magnetic tape magnetic disk storage or other magnetic storage devices, or any other non-transmission medium that can be used to store information that can be accessed by a computing device. As defined herein, a computer readable medium does not include a transitory computer readable medium such as a modulated data signal and a carrier wave.
It should also be noted that the terms "comprises," "comprising," or any other variation thereof, are intended to cover a non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements does not include only those elements but may include other elements not expressly listed or inherent to such process, method, article, or apparatus. Without further limitation, an element defined by the phrase "comprising an … …" does not exclude the presence of other like elements in a process, method, article, or apparatus that comprises the element.
As will be appreciated by one skilled in the art, embodiments of the present application may be provided as a method, system, or computer program product. Accordingly, the present application may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present application may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.
The application may be described in the general context of computer-executable instructions, such as program modules, being executed by a computer. Generally, program modules include routines, programs, objects, components, data structures, etc. that perform particular tasks or implement particular abstract data types. The application may also be practiced in distributed computing environments where tasks are performed by remote processing devices that are linked through a communications network. In a distributed computing environment, program modules may be located in both local and remote computer storage media including memory storage devices.
The embodiments in the present specification are described in a progressive manner, and the same and similar parts among the embodiments are referred to each other, and each embodiment focuses on the differences from the other embodiments. In particular, for the system embodiment, since it is substantially similar to the method embodiment, the description is simple, and for the relevant points, reference may be made to the partial description of the method embodiment.
The above description is only an example of the present application and is not intended to limit the present application. Various modifications and changes may occur to those skilled in the art. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present application should be included in the scope of the claims of the present application.
Claims (5)
1. A method for evaluating the degree of multiple contamination of seismic data, comprising:
acquiring VSP downlink wave data before deconvolution and VSP uplink wave data after deconvolution of the VSP data in a work area;
determining a multiple pollution discretization curve of the first wave and the multiple changing along with time or depth in the VSP downlink wave data before deconvolution, and determining a stratum interface response sequence of the VSP uplink wave data after deconvolution along with time or depth;
acquiring the multiple pollution degree of the seismic data of the work area along with the time or depth change according to the multiple pollution discretization curve and the stratum interface response sequence; wherein,
the determining a multiple pollution discretization curve of the first wave and the multiple changing with time or depth in the VSP downlink wave data before deconvolution comprises:
carrying out absolute value processing on the VSP downlink wave data before deconvolution;
determining the amplitude and the peak position of VSP downlink wave data before deconvolution after absolute value processing;
normalizing the peak amplitude, and fitting the peak position and the normalized peak amplitude to obtain a multiple pollution discretization curve of multiple changes along with time or depth in the VSP downlink wave data before deconvolution;
the determining a formation interface response sequence of the post-deconvolution VSP upgoing data changing with time or depth comprises:
carrying out absolute value processing on the VSP uplink wave data after deconvolution;
extracting the amplitude and the peak position of the deconvoluted VSP uplink wave data after absolute value processing to obtain a stratum interface response sequence of the deconvoluted VSP uplink wave data changing with time or depth;
the acquiring the multiple pollution degree of the seismic data of the work area along with the change of time or depth according to the multiple pollution discretization curve and the stratum interface response sequence comprises the following steps:
and multiplying each stratum interface response value in the stratum interface response sequence by the multiple pollution discretization curve respectively, and correspondingly placing the obtained result at the position of the corresponding stratum interface response value and below the corresponding stratum interface response value, thereby obtaining the multiple pollution degree of the seismic data of the work area along with the change of time or depth.
2. The method of evaluating the degree of seismic data multiples contamination of claim 1, wherein said obtaining pre-deconvolution VSP downgoing data and post-deconvolution VSP upgoing data for VSP data in a work area comprises:
decoding VSP data in the work area;
carrying out polarization rotation on the decoded VSP data to obtain a longitudinal wave field component and a tangential wave field component;
performing wave field separation processing on the longitudinal wave field component to obtain VSP downlink wave data before deconvolution of the VSP data; and carrying out deconvolution and wave field separation processing on the tangential wave field component in sequence to obtain the deconvoluted VSP upgoing wave data of the VSP data.
3. An apparatus for evaluating a degree of contamination of seismic data with multiples, comprising:
the uplink and downlink data acquisition module is used for acquiring VSP downlink wave data before deconvolution and VSP uplink wave data after deconvolution of the VSP data in the work area;
a curve and sequence acquisition module for determining a multiple pollution discretization curve of primary waves and multiple changes along with time or depth in the VSP downlink wave data before deconvolution and determining a stratum interface response sequence of the VSP uplink wave data after deconvolution along with time or depth changes;
the pollution degree determining module is used for acquiring the multiple pollution degree of the seismic data of the work area along with the time or depth change according to the multiple pollution discretization curve and the stratum interface response sequence; wherein,
the determining a multiple pollution discretization curve of the first wave and the multiple changing with time or depth in the VSP downlink wave data before deconvolution comprises:
carrying out absolute value processing on the VSP downlink wave data before deconvolution;
determining the amplitude and the peak position of VSP downlink wave data before deconvolution after absolute value processing;
normalizing the peak amplitude, and fitting the peak position and the normalized peak amplitude to obtain a multiple pollution discretization curve of multiple changes along with time or depth in the VSP downlink wave data before deconvolution;
the determining a formation interface response sequence of the post-deconvolution VSP upgoing data changing with time or depth comprises:
carrying out absolute value processing on the VSP uplink wave data after deconvolution;
extracting the amplitude and the peak position of the deconvoluted VSP uplink wave data after absolute value processing to obtain a stratum interface response sequence of the deconvoluted VSP uplink wave data changing with time or depth;
the acquiring the multiple pollution degree of the seismic data of the work area along with the change of time or depth according to the multiple pollution discretization curve and the stratum interface response sequence comprises the following steps:
and multiplying each stratum interface response value in the stratum interface response sequence by the multiple pollution discretization curve respectively, and correspondingly placing the obtained result at the position of the corresponding stratum interface response value and below the corresponding stratum interface response value, thereby obtaining the multiple pollution degree of the seismic data of the work area along with the change of time or depth.
4. The apparatus for evaluating the degree of seismic data multiples contamination of claim 3, wherein said obtaining pre-deconvolution VSP down-going data and post-deconvolution VSP up-going data for VSP data in a work area comprises:
decoding VSP data in the work area;
carrying out polarization rotation on the decoded VSP data to obtain a longitudinal wave field component and a tangential wave field component;
performing wave field separation processing on the longitudinal wave field component to obtain VSP downlink wave data before deconvolution of the VSP data; and carrying out deconvolution and wave field separation processing on the tangential wave field component in sequence to obtain the deconvoluted VSP upgoing wave data of the VSP data.
5. An apparatus for assessing the degree of multiple contamination of seismic data, comprising a memory, a processor, and a computer program stored on said memory, wherein said computer program when executed by said processor performs the steps of:
acquiring VSP downlink wave data before deconvolution and VSP uplink wave data after deconvolution of the VSP data in a work area;
determining a multiple pollution discretization curve of the first wave and the multiple changing along with time or depth in the VSP downlink wave data before deconvolution, and determining a stratum interface response sequence of the VSP uplink wave data after deconvolution along with time or depth;
acquiring the multiple pollution degree of the seismic data of the work area along with the time or depth change according to the multiple pollution discretization curve and the stratum interface response sequence; wherein,
the determining a multiple pollution discretization curve of the first wave and the multiple changing with time or depth in the VSP downlink wave data before deconvolution comprises:
carrying out absolute value processing on the VSP downlink wave data before deconvolution;
determining the amplitude and the peak position of VSP downlink wave data before deconvolution after absolute value processing;
normalizing the peak amplitude, and fitting the peak position and the normalized peak amplitude to obtain a multiple pollution discretization curve of multiple changes along with time or depth in the VSP downlink wave data before deconvolution;
the determining a formation interface response sequence of the post-deconvolution VSP upgoing data changing with time or depth comprises:
carrying out absolute value processing on the VSP uplink wave data after deconvolution;
extracting the amplitude and the peak position of the deconvoluted VSP uplink wave data after absolute value processing to obtain a stratum interface response sequence of the deconvoluted VSP uplink wave data changing with time or depth;
the acquiring the multiple pollution degree of the seismic data of the work area along with the change of time or depth according to the multiple pollution discretization curve and the stratum interface response sequence comprises the following steps:
and multiplying each stratum interface response value in the stratum interface response sequence by the multiple pollution discretization curve respectively, and correspondingly placing the obtained result at the position of the corresponding stratum interface response value and below the corresponding stratum interface response value, thereby obtaining the multiple pollution degree of the seismic data of the work area along with the change of time or depth.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711144555.3A CN108107485B (en) | 2017-11-17 | 2017-11-17 | The appraisal procedure and device of seismic data multiple wave pollution level |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711144555.3A CN108107485B (en) | 2017-11-17 | 2017-11-17 | The appraisal procedure and device of seismic data multiple wave pollution level |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108107485A CN108107485A (en) | 2018-06-01 |
| CN108107485B true CN108107485B (en) | 2019-11-08 |
Family
ID=62207672
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201711144555.3A Active CN108107485B (en) | 2017-11-17 | 2017-11-17 | The appraisal procedure and device of seismic data multiple wave pollution level |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108107485B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109471162A (en) * | 2018-10-08 | 2019-03-15 | 中国石油天然气集团有限公司 | Interbed multiple processing method, system, electronic equipment and readable medium |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| GB0425022D0 (en) * | 2003-11-14 | 2004-12-15 | Schlumberger Holdings | High-frequency processing of seismic vibrator data |
| CN101071175A (en) * | 2006-05-11 | 2007-11-14 | 中国石油集团东方地球物理勘探有限责任公司 | Zero hypocentral distance vertical seismic section compressional-shear wave data depth field corridor stacked section processing method |
| CN104199107A (en) * | 2014-08-06 | 2014-12-10 | 中国石油天然气集团公司 | Vertical seismic multi-wave data based predication method and system for depth before drilling |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US9823369B2 (en) * | 2013-02-28 | 2017-11-21 | Cgg Services Sas | System and method for correcting near surface statics by using internal multiples prediction |
-
2017
- 2017-11-17 CN CN201711144555.3A patent/CN108107485B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| GB0425022D0 (en) * | 2003-11-14 | 2004-12-15 | Schlumberger Holdings | High-frequency processing of seismic vibrator data |
| CN101071175A (en) * | 2006-05-11 | 2007-11-14 | 中国石油集团东方地球物理勘探有限责任公司 | Zero hypocentral distance vertical seismic section compressional-shear wave data depth field corridor stacked section processing method |
| CN104199107A (en) * | 2014-08-06 | 2014-12-10 | 中国石油天然气集团公司 | Vertical seismic multi-wave data based predication method and system for depth before drilling |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108107485A (en) | 2018-06-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Shephard et al. | On the consistency of seismically imaged lower mantle slabs | |
| Edwards et al. | Determination of site amplification from regional seismicity: application to the Swiss National Seismic Networks | |
| Górszczyk et al. | Application of curvelet denoising to 2D and 3D seismic data—Practical considerations | |
| Bindi et al. | Site amplifications observed in the Gubbio Basin, Central Italy: hints for lateral propagation effects | |
| Goutorbe et al. | Rayleigh wave group velocities at periods of 6–23 s across Brazil from ambient noise tomography | |
| Oth et al. | Spectral analysis of K-NET and KiK-net data in Japan, Part I: Database compilation and peculiarities | |
| Zhao et al. | Separating prestack diffractions with SVMF in the flattened shot domain | |
| CN108107485B (en) | The appraisal procedure and device of seismic data multiple wave pollution level | |
| CN108445532B (en) | A kind of Depth Domain inverse migration method and device | |
| Drouet et al. | Attenuation, seismic moments, and site effects for weak-motion events: application to the Pyrenees | |
| Zhang et al. | Improving receiver function imaging with high-resolution Radon transform | |
| CN107422377B (en) | Utilize the method and device of seismic data prediction hydrocarbon source rock total content of organic carbon | |
| Muller et al. | Complete identification of complex salt geometries from inaccurate migrated subsurface offset gathers using deep learning | |
| RU2573125C1 (en) | Phase modulation and noise minimisation for simultaneous vibroseis acquisition | |
| CN111751870A (en) | Method and device for suppressing multiple waves between stacked layers | |
| US10330809B2 (en) | Device and method for optimization of 4D and 3D seismic data | |
| Picozzi et al. | Deblurring of frequency–wavenumber images from small-scale seismic arrays | |
| CN111596358A (en) | Multiple suppression method, device and system | |
| CN106291697A (en) | A kind of method and system of the value of the Q of interval quality factors definitely | |
| CN112049623B (en) | Method and device for judging water outlet of medium-base volcanic bound water | |
| Zare | Recent development of the earthquake strong motion-intensity catalog and intensity prediction equations for Iran | |
| Oren et al. | An overview of reproducible 3D seismic data processing and imaging using Madagascar | |
| US20180059275A1 (en) | System and method for mapping horizons in seismic images | |
| Wang et al. | High dimensional multistep deblending using supervised training and transfer learning | |
| Falconer et al. | Attribute-driven footprint suppression |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |