CN115857007A - Full-signal imaging method and device - Google Patents
Full-signal imaging method and device Download PDFInfo
- Publication number
- CN115857007A CN115857007A CN202211424416.7A CN202211424416A CN115857007A CN 115857007 A CN115857007 A CN 115857007A CN 202211424416 A CN202211424416 A CN 202211424416A CN 115857007 A CN115857007 A CN 115857007A
- Authority
- CN
- China
- Prior art keywords
- imaging
- seismic
- point
- data
- processing
- 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.)
- Pending
Links
Images
Landscapes
- Image Processing (AREA)
Abstract
The invention provides a full-signal imaging method and a device, wherein the method comprises the following steps: carrying out noise suppression processing on the acquired pre-stack seismic data, and carrying out seismic data development processing on the pre-stack seismic data after noise suppression processing to obtain a seismic section; determining imaging point data corresponding to each central point according to a plurality of central points on the seismic profile, and carrying out field angle arrangement processing on the imaging point data to obtain a central point field angle gather; and determining an imaging area by using the central point field angle gather, and performing imaging processing on each imaging point in the imaging area to obtain a seismic imaging section. According to the invention, the seismic data is suppressed and denoised, and the processing of each imaging point in the imaging area is combined to obtain the seismic imaging section, so that the frequency bands participating in imaging are expanded, and the signals participating in imaging are richer and more complete, thereby solving the problems that the existing conventional imaging method cannot effectively utilize the data which is in the frequency suppression band and has a lower signal-to-noise ratio, and the problem of severe noise amplification caused by the processing of a plurality of broadband signals.
Description
Technical Field
The invention relates to the technical field of seismic data processing, in particular to a full-signal imaging method and a full-signal imaging device.
Background
In the conventional imaging processing process of the reflection seismic data, only the seismic data in a frequency band range with a higher signal-to-noise ratio is generally extracted and imaged, and low-frequency and high-frequency data which are in a frequency suppression band and have a relatively lower signal-to-noise ratio are not processed, so that effective signals of seismic original data are not fully mined and utilized, the effective bandwidth of the seismic imaging data is reduced, and the precision and the accuracy of subsequent seismic exploration work are influenced. Generally, the wider the imaging frequency band is, the more the number of signals participating in imaging is, the more easily the problem of severe noise amplification is generated during imaging, and the industrial industry generally adopts the methods of limiting frequency band width, partial signal stacking, post-stack filtering and the like to improve the signal-to-noise ratio of a seismic section, but such methods lose and damage seismic effective signals, are not seismic full-signal processing, and are not favorable for seismic data amplitude-preserving fidelity.
In summary, the conventional imaging method has the disadvantages that data in a frequency suppression band and with a low signal-to-noise ratio cannot be effectively utilized, and for the problem of noise amplification caused by broadband multi-signal imaging, the method commonly used in the industry has adverse effects on seismic effective signals. In view of the above problems, no effective solution has been proposed.
Disclosure of Invention
In view of the problems in the prior art, a primary object of embodiments of the present invention is to provide a full-signal imaging method and apparatus, which implement expanding an imaging frequency band, so that the imaging signals are richer and more complete.
In order to achieve the above object, an embodiment of the present invention provides a full-signal imaging method, including:
carrying out noise suppression processing on the acquired pre-stack seismic data, and carrying out seismic data development processing on the pre-stack seismic data subjected to noise suppression processing to obtain a seismic section;
determining imaging point data corresponding to each central point according to a plurality of central points on the seismic profile, and carrying out field angle arrangement processing on the imaging point data to obtain a central point field angle gather;
and determining an imaging area by using the central opening angle gather, and performing imaging processing on each imaging point in the imaging area to obtain a seismic imaging section.
Optionally, in an embodiment of the present invention, determining, according to a plurality of central points on the seismic section, imaging point data corresponding to each central point includes:
and calculating the seismic wave travel time, the migration amplitude and the stratum opening angle corresponding to the imaging point at each central point position according to the central point position corresponding to each central point on the seismic section.
Optionally, in an embodiment of the present invention, the performing opening angle arrangement processing on the imaging point data to obtain a central opening angle gather includes:
and carrying out field angle arrangement and data superposition by using the size of the stratum field angle in the imaging point data to obtain a central point field angle gather corresponding to the central point position.
Optionally, in an embodiment of the present invention, the imaging processing on each imaging point in the imaging area to obtain a seismic imaging section includes:
and carrying out opening angle calculation on each imaging point in the imaging area, and if the opening angle obtained by calculation is positioned in the imaging area, calculating the walking time and amplitude of the imaging point to finish seismic imaging processing to obtain a seismic imaging section.
An embodiment of the present invention further provides a full-signal imaging device, including:
the seismic profile module is used for suppressing noise of the acquired pre-stack seismic data and performing seismic data development processing on the pre-stack seismic data after noise suppression processing to obtain a seismic profile;
the field angle gather module is used for determining imaging point data corresponding to each central point according to a plurality of central points on the seismic section, and carrying out field angle arrangement processing on the imaging point data to obtain a central point field angle gather;
and the seismic imaging module is used for determining an imaging area by using the central point field angle gather, and performing imaging processing on each imaging point in the imaging area to obtain a seismic imaging section.
Optionally, in an embodiment of the present invention, the field angle gather module is further configured to calculate, according to a central point position corresponding to each central point on the seismic profile, a seismic travel time, a migration amplitude value, and a formation field angle corresponding to the imaging point at each central point position.
Optionally, in an embodiment of the present invention, the aperture angle gather module is further configured to perform aperture angle arrangement and data superposition by using the size of the formation aperture angle in the imaging point data, so as to obtain a central aperture angle gather corresponding to the central point position.
Optionally, in an embodiment of the present invention, the seismic imaging module is further configured to perform opening angle calculation on each imaging point in the imaging area, and if the opening angle obtained through calculation is located in the imaging area, the imaging point is calculated with the amplitude when walking, so as to complete seismic imaging processing, and obtain the seismic imaging section.
The invention also provides an electronic device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, wherein the processor implements the method when executing the program.
The present invention also provides a computer-readable storage medium having stored thereon a computer program for executing the above method.
The invention also provides a computer program product comprising computer programs/instructions which, when executed by a processor, implement the steps of the above method.
The invention reduces the noise of the seismic data by pressing, combines the processing of each imaging point in the imaging area to obtain the seismic imaging section, expands the frequency band participating in imaging, and enables the signals participating in imaging to be richer and more complete, thereby solving the problems that the prior conventional imaging method can not effectively utilize the data which are in the frequency pressing band and have lower signal-to-noise ratio, and the problem of severe noise amplification caused by the processing of the wide-frequency multi-signal.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments will be briefly introduced below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and it is obvious for those skilled in the art that other drawings can be obtained based on these drawings without creative efforts.
FIG. 1 is a flow chart of a full signal imaging method according to an embodiment of the present invention;
FIG. 2 is a schematic cross-sectional view of a conventional deconvolution stack in an embodiment of the present invention;
FIG. 3 is a schematic cross-sectional view of pulse deconvolution stack in an embodiment of the present disclosure;
FIG. 4 is a schematic diagram of a comparison of conventional deconvolution and pulse deconvolution spectra in an embodiment of the present invention;
FIG. 5 is a schematic cross-sectional view of a full-frequency imaging method using a conventional post-stack denoising method according to an embodiment of the present invention;
FIG. 6 is a schematic cross-sectional view of a full-frequency imaging system using the method of the present invention in an embodiment of the present invention;
FIG. 7 is a schematic diagram illustrating a comparison between the frequency spectrums of FIGS. 5 and 6 according to an embodiment of the present invention;
FIG. 8 is a schematic structural diagram of a full-signal imaging apparatus according to an embodiment of the present invention;
fig. 9 is a schematic structural diagram of an electronic device according to an embodiment of the present invention.
Detailed Description
The embodiment of the invention provides a full-signal imaging method and device.
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, 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 invention.
Fig. 1 is a flowchart illustrating a full-signal imaging method according to an embodiment of the present invention, and an implementation subject of the full-signal imaging method according to the embodiment of the present invention includes, but is not limited to, a computer. The invention reduces the noise of the seismic data by pressing, combines the processing of each imaging point in the imaging area to obtain the seismic imaging section, expands the frequency band participating in imaging, and enables the signals participating in imaging to be richer and more complete, thereby solving the problems that the prior conventional imaging method can not effectively utilize the data which are in the frequency pressing band and have lower signal-to-noise ratio, and the problem of severe noise amplification caused by the processing of the wide-frequency multi-signal. The method shown in the figure comprises the following steps:
s1, performing noise suppression processing on the acquired pre-stack seismic data, and performing seismic data development processing on the pre-stack seismic data subjected to noise suppression processing to obtain a seismic section;
s2, determining imaging point data corresponding to each central point according to the plurality of central points on the seismic section, and carrying out field angle arrangement processing on the imaging point data to obtain a central point field angle gather;
and S3, determining an imaging area by using the central point field angle gather, and performing imaging processing on each imaging point in the imaging area to obtain a seismic imaging section.
The pre-stack seismic data acquired in advance is subjected to noise suppression processing, such as surface wave removal, industrial electrical interference, multiple waves and the like, so that the purpose of improving the signal-to-noise ratio of the seismic data is achieved. And (3) carrying out processing on the seismic data by using an impulse deconvolution technology to obtain a seismic section. And aiming at a plurality of common center point positions (marked as SAM-CDP) on the seismic section, solving the seismic wave travel time, the offset amplitude and the stratum opening angle of each imaging point corresponding to the common center point positions, arranging according to the opening angle, and partially overlapping to obtain an opening angle gather at the SAM-CDP position.
Further, an imaging area is determined according to an opening angle gather corresponding to the SAM-CDP position, and the imaging areas corresponding to the rest CDP positions are determined in an interpolation continuation mode by utilizing the determined imaging areas. And calculating the opening angle of each imaging point on all CDP positions of the data before folding, if the opening angle is positioned in the imaging area, continuously calculating the travel time and the amplitude of the imaging point to finish imaging, and if the opening angle is positioned outside the imaging area, repeating the operation of the step for the next imaging point. And finishing the calculation of all CDP position imaging points to obtain a final seismic imaging section.
As an embodiment of the present invention, determining imaging point data corresponding to each central point according to a plurality of central points on a seismic section includes: and calculating the seismic wave travel time, the migration amplitude and the stratum opening angle corresponding to the imaging point at each central point position according to the central point position corresponding to each central point on the seismic section.
In this embodiment, the processing of field angle arrangement on the imaging point data to obtain a central field angle gather includes: and carrying out field angle arrangement and data superposition by using the size of the stratum field angle in the imaging point data to obtain a central point field angle gather corresponding to the central point position.
Wherein, the positions of a plurality of concentric points on the seismic imaging section are determined as follows: the common center point position is selected according to different intervals on the conventional section according to the structure change condition, the interval is reduced when the structure is complex, and the interval can be increased when the structure is simple.
Further, the calculation of the seismic travel time, the offset amplitude and the opening angle is realized as follows: defining the horizontal coordinate of a shot point as alpha, the horizontal coordinate of a demodulator probe as beta, the lambda as an azimuth angle, the coordinate of an imaging point as (x, tau), and the offset speed of the imaging point as vel rms And seismic wave travel time from a shot point to an imaging point and then to a wave detection point:
further, integers are defined
Wherein, delta t is sampling interval of seismic data, amplitude values of four points of M-1, M +1 and M +2 on the seismic data are utilized to carry out quadratic curve interpolation, seismic amplitude value xi at the travel time t is solved, and an imaging weight coefficient is calculated as follows:
the magnitude of the shift at the imaging point (x, τ) is then:
as an embodiment of the present invention, the imaging processing of each imaging point in the imaging area to obtain a seismic imaging section includes: and carrying out opening angle calculation on each imaging point in the imaging area, and if the opening angle obtained by calculation is positioned in the imaging area, calculating the walking time and amplitude of the imaging point to finish seismic imaging processing to obtain a seismic imaging section.
Wherein, the field angle is obtained: the following variables, v 1 ＝α-v,v 2 α + β -2v; calculating the distance variable z 1 ＝v 1 cosφ；z 2 ＝v 1 v 2 +2(vel rms τ) 2 ；z 3 ＝v 2 cosφ；
Is the x-axis abscissa, then z 5 ＝(z 1 z 4 2 -z 2 z 3 )/(z 3 2 -z 4 2 ) Then the opening angle theta is arctan [ (v) 1 +z 5 cosλ)/(vel rms τ)]。
On the basis of the conventional frequency band, the invention excavates and recycles the data with low signal-to-noise ratio in the frequency suppression band of the seismic data, avoids the defects of insufficient excavation and effective signal loss of the conventional method to the data signal and achieves the purpose of full signal processing. The invention utilizes the field angle gather to determine the accurate imaging area, suppresses signal noise in the imaging process and solves the problem of severe noise amplification caused by wide-frequency multi-signal processing.
In a specific embodiment of the present invention, the imaging method in the present invention is used to solve the problem that the conventional imaging method cannot effectively utilize data within a frequency suppression band and with a low signal-to-noise ratio, and the problem of severe noise amplification caused by processing of multiple broadband signals. Compared with the conventional method, the method expands the frequency band participating in imaging, and the signals participating in imaging are richer and more complete.
In this embodiment, taking a seismic data of a certain oil field as an example, the method specifically includes the following steps:
1. pre-stack seismic data is subjected to noise suppression processing, wherein the data sampling interval is 1ms, the recording time of seismic signals is 5000ms, the CDP interval is 20m, the minimum offset is 300m, the maximum offset is 4800m, and the offset interval is 100m.
2. And (3) carrying out processing on the seismic data obtained in the step (1) by using an impulse deconvolution technology.
3. And aiming at a plurality of common center point positions (marked as SAM-CDP) on the seismic imaging section, solving the seismic wave travel time, the offset amplitude and the stratum opening angle of each imaging point corresponding to the common center point positions, arranging according to the opening angle, and partially overlapping to obtain an opening angle gather at the SAM-CDP position.
4. Determining imaging areas aiming at flare angle trace sets corresponding to SAM-CDP positions, and determining the imaging areas corresponding to the rest CDP positions by utilizing the determined imaging areas in an interpolation continuation mode;
5. and (3) solving the opening angle of each imaging point on all CDP positions of the pre-stack data, if the opening angle is positioned in the imaging area, continuously calculating the travel time and the amplitude of the imaging point to finish imaging, and if the opening angle is positioned outside the imaging area, repeating the operation of the step for the next imaging point.
6. And finishing the calculation of all CDP position imaging points to obtain a final seismic imaging section.
Wherein, fig. 2 is a conventional deconvolution cross section, fig. 3 is a deconvolution cross section obtained by the method of the present invention, it can be seen that the signal-to-noise ratio of the deep layer in fig. 3 is significantly low, and the high-frequency noise is severely developed, fig. 4 is a comparison of the conventional deconvolution with the deconvolution frequency spectrum obtained by the method of the present invention, the dotted line on the graph represents the conventional deconvolution frequency spectrum, the solid line represents the method of the present invention, and it can be seen that the high-frequency part of the frequency spectrum is significantly broadened after the pulse deconvolution processing. Fig. 5 is a full-frequency imaging section using a conventional post-stack denoising method, fig. 6 is a full-frequency imaging section using the method of the present invention, it can be seen that interlayer information on fig. 6 is richer, and details are depicted more clearly, fig. 7 is a comparison of the section spectra of fig. 5 and fig. 6, a dotted line on the graph represents the frequency spectrum of the conventional post-stack denoising method, and a solid line represents the frequency spectrum processed by the method of the present invention, it can be seen that the frequency spectrum obtained by the method of the present invention is significantly better than the conventional post-stack denoising method in terms of bandwidth, although the conventional post-stack denoising suppresses high-frequency noise, which increases the signal-to-noise ratio, but damages the high-frequency effective signal, and the method of the present invention does not damage the high-frequency effective signal in the process of increasing the signal-to-noise ratio.
The invention reduces the noise of the seismic data by pressing, combines the processing of each imaging point in the imaging area to obtain the seismic imaging section, expands the frequency band participating in imaging, and enables the signals participating in imaging to be richer and more complete, thereby solving the problems that the prior conventional imaging method can not effectively utilize the data which are in the frequency pressing band and have lower signal-to-noise ratio, and the problem of severe noise amplification caused by the processing of the wide-frequency multi-signal.
Fig. 8 is a schematic structural diagram of a full-signal imaging apparatus according to an embodiment of the present invention, where the apparatus includes:
the seismic section module 10 is configured to perform noise suppression processing on the acquired pre-stack seismic data, and perform seismic data development processing on the pre-stack seismic data after the noise suppression processing to obtain a seismic section;
the field angle gather module 20 is configured to determine imaging point data corresponding to each central point according to a plurality of central points on the seismic section, and perform field angle arrangement processing on the imaging point data to obtain a central point field angle gather;
and the seismic imaging module 30 is configured to determine an imaging area by using the central point field angle gather, and perform imaging processing on each imaging point in the imaging area to obtain a seismic imaging section.
As an embodiment of the present invention, the flare angle gather module 20 is further configured to calculate seismic travel time, offset amplitude and formation flare angle corresponding to the imaging point at each central point position according to the central point position corresponding to each central point on the seismic section.
In this embodiment, the aperture angle gather module 20 is further configured to perform aperture angle arrangement and data superposition by using the size of the formation aperture angle in the imaging point data, so as to obtain a central aperture angle gather corresponding to the central point position.
As an embodiment of the present invention, the seismic imaging module 30 is further configured to perform an opening angle calculation on each imaging point in the imaging area, and if the opening angle obtained through calculation is located in the imaging area, calculate the amplitude and the walking time of the imaging point to complete the seismic imaging processing, so as to obtain the seismic imaging profile.
Based on the same application concept as the full-signal imaging method, the invention also provides the full-signal imaging device. Because the principle of solving the problems of the full-signal imaging device is similar to that of the full-signal imaging method, the implementation of the full-signal imaging device can refer to the implementation of the full-signal imaging method, and repeated parts are not described again.
The invention reduces the noise of the seismic data by pressing, combines the processing of each imaging point in the imaging area to obtain the seismic imaging section, expands the frequency band participating in imaging, and enables the signals participating in imaging to be richer and more complete, thereby solving the problems that the prior conventional imaging method can not effectively utilize the data which are in the frequency pressing band and have lower signal-to-noise ratio, and the problem of severe noise amplification caused by the processing of the wide-frequency multi-signal.
The invention also provides an electronic device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, wherein the processor implements the method when executing the program.
The present invention also provides a computer program product comprising computer programs/instructions which, when executed by a processor, implement the steps of the above method.
The present invention also provides a computer-readable storage medium having stored thereon a computer program for executing the above method.
As shown in fig. 9, the electronic device 600 may further include: communication module 110, input unit 120, audio processor 130, display 160, power supply 170. It is noted that the electronic device 600 does not necessarily include all of the components shown in FIG. 9; furthermore, the electronic device 600 may also comprise components not shown in fig. 9, which may be referred to in the prior art.
As shown in fig. 9, the central processor 100, sometimes referred to as a controller or operational control, may include a microprocessor or other processor device and/or logic device, the central processor 100 receiving input and controlling the operation of the various components of the electronic device 600.
The memory 140 may be, for example, one or more of a buffer, a flash memory, a hard drive, a removable media, a volatile memory, a non-volatile memory, or other suitable device. The information relating to the failure may be stored, and a program for executing the information may be stored. And the central processing unit 100 may execute the program stored in the memory 140 to realize information storage or processing, etc.
The input unit 120 provides input to the cpu 100. The input unit 120 is, for example, a key or a touch input device. The power supply 170 is used to provide power to the electronic device 600. The display 160 is used to display an object to be displayed, such as an image or a character. The display may be, for example, an LCD display, but is not limited thereto.
The memory 140 may be a solid state memory such as Read Only Memory (ROM), random Access Memory (RAM), a SIM card, or the like. There may also be a memory that holds information even when power is off, can be selectively erased, and is provided with more data, an example of which is sometimes called an EPROM or the like. The memory 140 may also be some other type of device. Memory 140 includes buffer memory 141 (sometimes referred to as a buffer). The memory 140 may include an application/function storage section 142, and the application/function storage section 142 is used to store application programs and function programs or a flow for executing the operation of the electronic device 600 by the central processing unit 100.
The memory 140 may also include a data store 143, the data store 143 for storing data, such as contacts, digital data, pictures, sounds, and/or any other data used by the electronic device. The driver storage portion 144 of the memory 140 may include various drivers of the electronic device for communication functions and/or for performing other functions of the electronic device (e.g., messaging application, address book application, etc.).
The communication module 110 is a transmitter/receiver 110 that transmits and receives signals via an antenna 111. The communication module (transmitter/receiver) 110 is coupled to the central processor 100 to provide an input signal and receive an output signal, which may be the same as in the case of a conventional mobile communication terminal.
Based on different communication technologies, a plurality of communication modules 110, such as a cellular network module, a bluetooth module, and/or a wireless local area network module, may be provided in the same electronic device. The communication module (transmitter/receiver) 110 is also coupled to a speaker 131 and a microphone 132 via an audio processor 130 to provide audio output via the speaker 131 and receive audio input from the microphone 132 to implement general telecommunications functions. Audio processor 130 may include any suitable buffers, decoders, amplifiers and so forth. In addition, an audio processor 130 is also coupled to the central processor 100, so that recording on the local can be enabled through a microphone 132, and so that sound stored on the local can be played through a speaker 131.
As will be appreciated by one skilled in the art, embodiments of the present invention may be provided as a method, system, or computer program product. Accordingly, the present invention may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present invention 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 present invention has been 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.
The principle and the implementation mode of the invention are explained by applying specific embodiments in the invention, and the description of the embodiments is only used for helping to understand the method and the core idea of the invention; meanwhile, for a person skilled in the art, according to the idea of the present invention, there may be variations in the specific embodiments and the application scope, and in summary, the content of the present specification should not be construed as a limitation to the present invention.
Claims (11)
1. A method of full signal imaging, the method comprising:
carrying out noise suppression processing on the acquired pre-stack seismic data, and carrying out seismic data development processing on the pre-stack seismic data subjected to noise suppression processing to obtain a seismic section;
determining imaging point data corresponding to each central point according to the plurality of central points on the seismic section, and carrying out field angle arrangement processing on the imaging point data to obtain a central point field angle gather;
and determining an imaging area by using the central point field angle gather, and performing imaging processing on each imaging point in the imaging area to obtain a seismic imaging section.
2. The method of claim 1, wherein determining imaging point data corresponding to each of the plurality of centerpoints on the seismic section from the plurality of centerpoints comprises:
and calculating the seismic wave travel time, the migration amplitude and the size of the stratum opening angle corresponding to the imaging point at each central point position according to the central point position corresponding to each central point on the seismic section.
3. The method of claim 2, wherein the opening angle ranking the imaging point data to obtain a central opening angle gather comprises:
and carrying out field angle arrangement and data superposition by using the size of the stratum field angle in the imaging point data to obtain a central point field angle gather corresponding to the central point position.
4. The method of claim 1, wherein imaging each imaging point in the imaging zone to obtain a seismic imaging profile comprises:
and calculating the field angle of each imaging point in the imaging area, and if the field angle obtained by calculation is positioned in the imaging area, calculating the walking time and amplitude of the imaging point to finish seismic imaging processing to obtain the seismic imaging section.
5. A full signal imaging apparatus, characterized in that the apparatus comprises:
the seismic profile module is used for suppressing noise of the acquired pre-stack seismic data and performing seismic data development processing on the pre-stack seismic data after noise suppression processing to obtain a seismic profile;
the field angle gather module is used for determining imaging point data corresponding to each central point according to the plurality of central points on the seismic section, and carrying out field angle arrangement processing on the imaging point data to obtain a central point field angle gather;
and the seismic imaging module is used for determining an imaging area by using the central point field angle gather, and performing imaging processing on each imaging point in the imaging area to obtain a seismic imaging section.
6. The apparatus of claim 5, wherein the field angle gather module is further configured to calculate seismic travel time, migration amplitude and formation field angle corresponding to the imaging point at each central point position according to the central point position corresponding to each central point on the seismic section.
7. The apparatus of claim 6, wherein the aperture angle gather module is further configured to perform aperture angle arrangement and data superposition according to the size of the formation aperture angle in the imaging point data to obtain a central aperture angle gather corresponding to the central point position.
8. The device of claim 5, wherein the seismic imaging module is further configured to perform a flare angle calculation on each imaging point in the imaging area, and if the calculated flare angle is located in the imaging area, perform a walking and amplitude calculation on the imaging point to complete a seismic imaging process to obtain the seismic imaging profile.
9. An electronic device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, wherein the processor implements the method of any one of claims 1 to 4 when executing the computer program.
10. A computer-readable storage medium, characterized in that the computer-readable storage medium stores a computer program for executing the method of any one of claims 1 to 4.
11. A computer program product comprising computer programs/instructions, characterized in that the computer programs/instructions, when executed by a processor, implement the steps of the method of any of claims 1 to 4.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211424416.7A CN115857007A (en) | 2022-11-15 | 2022-11-15 | Full-signal imaging method and device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211424416.7A CN115857007A (en) | 2022-11-15 | 2022-11-15 | Full-signal imaging method and device |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN115857007A true CN115857007A (en) | 2023-03-28 |
Family
ID=85663411
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202211424416.7A Pending CN115857007A (en) | 2022-11-15 | 2022-11-15 | Full-signal imaging method and device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115857007A (en) |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20130182538A1 (en) * | 2012-01-12 | 2013-07-18 | Westerngeco L.L.C. | Reverse time migration model dip-guided imaging |
| US20140254320A1 (en) * | 2013-03-11 | 2014-09-11 | Saudi Arabian Oil Company | Vertical seismic profiling-common midpoint gathers (vsp-cmp) by wavefield extrapolation |
| CN110646840A (en) * | 2018-06-27 | 2020-01-03 | 中国石油化工股份有限公司 | Angle gather extraction method and system |
| CN110703327A (en) * | 2019-10-13 | 2020-01-17 | 东北石油大学 | Full-band imaging method |
| CN111538077A (en) * | 2020-05-07 | 2020-08-14 | 中国石油天然气集团有限公司 | Pre-stack depth migration method and device based on dip angle constraint |
-
2022
- 2022-11-15 CN CN202211424416.7A patent/CN115857007A/en active Pending
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20130182538A1 (en) * | 2012-01-12 | 2013-07-18 | Westerngeco L.L.C. | Reverse time migration model dip-guided imaging |
| US20140254320A1 (en) * | 2013-03-11 | 2014-09-11 | Saudi Arabian Oil Company | Vertical seismic profiling-common midpoint gathers (vsp-cmp) by wavefield extrapolation |
| CN110646840A (en) * | 2018-06-27 | 2020-01-03 | 中国石油化工股份有限公司 | Angle gather extraction method and system |
| CN110703327A (en) * | 2019-10-13 | 2020-01-17 | 东北石油大学 | Full-band imaging method |
| CN111538077A (en) * | 2020-05-07 | 2020-08-14 | 中国石油天然气集团有限公司 | Pre-stack depth migration method and device based on dip angle constraint |
Non-Patent Citations (1)
| Title |
|---|
| 陈生昌;马在田;WU RU-SHAN;: "波动方程角度域偏移成像", 计算物理, no. 02, pages 211 - 215 * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105355210B (en) | Preprocessing method and device for far-field speech recognition | |
| US10270412B2 (en) | Sound-field correction device, sound-field correction method, and sound-field correction program | |
| CN108717495A (en) | The method, apparatus and electronic equipment of multi-beam beam forming | |
| CN107369460A (en) | Speech sound enhancement device and method based on acoustics vector sensor space sharpening technique | |
| CN115857007A (en) | Full-signal imaging method and device | |
| CN113687305A (en) | Method, device and equipment for positioning sound source azimuth and computer readable storage medium | |
| CN110703327A (en) | Full-band imaging method | |
| CN113743193B (en) | Pre-stack seismic data linear interference suppression method and system | |
| CN114296134B (en) | Deep convolutional network seismic data unmixing method and system | |
| CN112565124B (en) | Distortion-free signal processing method based on overlapping windowing, storage medium and system | |
| CN111781644B (en) | Linear interference attenuation method and device for shallow and medium stratum seismic data | |
| CN112558004B (en) | Method and device for determining wave arrival direction of beam information and storage medium | |
| CN112147693B (en) | Method and device for suppressing residual multiples with near offset distance | |
| CN112014884B (en) | Method and device for suppressing near shot point strong energy noise | |
| CN112067927A (en) | Medium-high frequency oscillation detection method and device | |
| WO2020186434A1 (en) | Flexible differential microphone arrays with fractional order | |
| CN111366973B (en) | Frequency domain noise generation and addition method and device of forward model | |
| CN113687420B (en) | Method and device for removing surface wave noise of seismic data | |
| CN111724808A (en) | Audio signal processing method, device, terminal and storage medium | |
| CN112859160A (en) | Method and device for predicting thickness of thin interbed sand body | |
| CN110737019B (en) | Processing method and system for improving velocity spectrum precision | |
| CN114578418A (en) | Partition denoising method and device based on near-surface type | |
| CN111781646B (en) | Amplitude-preserving and fidelity-preserving linear coherent noise suppression method and processing terminal | |
| CN113514799A (en) | Sound source positioning method, device and equipment based on microphone array and storage medium | |
| CN113031071B (en) | Static correction method and device for long wavelength of earthquake wave |
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 |