CN113466938B - Integral method prestack depth migration method and device - Google Patents
Integral method prestack depth migration method and device Download PDFInfo
- Publication number
- CN113466938B CN113466938B CN202010243070.5A CN202010243070A CN113466938B CN 113466938 B CN113466938 B CN 113466938B CN 202010243070 A CN202010243070 A CN 202010243070A CN 113466938 B CN113466938 B CN 113466938B
- Authority
- CN
- China
- Prior art keywords
- depth
- travel time
- time table
- stratum
- step length
- 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
- 238000000034 method Methods 0.000 title claims abstract description 55
- 230000005012 migration Effects 0.000 title claims abstract description 53
- 238000013508 migration Methods 0.000 title claims abstract description 53
- 238000003384 imaging method Methods 0.000 claims abstract description 42
- 238000005070 sampling Methods 0.000 claims abstract description 38
- 238000004364 calculation method Methods 0.000 claims abstract description 36
- 238000012545 processing Methods 0.000 claims abstract description 17
- 230000015572 biosynthetic process Effects 0.000 claims description 53
- 238000004590 computer program Methods 0.000 claims description 16
- 230000006870 function Effects 0.000 claims description 10
- 230000010354 integration Effects 0.000 claims description 10
- 230000003750 conditioning effect Effects 0.000 claims description 9
- 238000003860 storage Methods 0.000 claims description 6
- 230000008859 change Effects 0.000 abstract description 6
- 238000010586 diagram Methods 0.000 description 15
- 230000008569 process Effects 0.000 description 7
- 230000004048 modification Effects 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 238000012952 Resampling Methods 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000009776 industrial production Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/301—Analysis for determining seismic cross-sections or geostructures
- G01V1/302—Analysis for determining seismic cross-sections or geostructures in 3D data cubes
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/50—Corrections or adjustments related to wave propagation
- G01V2210/51—Migration
- G01V2210/512—Pre-stack
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention provides a method and a device for prestack depth migration by an integral method, wherein the method comprises the following steps: determining the travel time table step length corresponding to each stratum depth according to the seismic data and the calculation parameters; sampling in the depth direction of the stratum to obtain a travel time table of a target area; according to the travel time table step length and the travel time table corresponding to each stratum depth, carrying out linear interpolation on the sampling points, and carrying out integral processing to obtain the migration result of the seismic data of the target area; compared with the step length of a small and fixed travel time table in the prior art, the step length of the travel time table corresponding to each stratum depth is determined, namely the step length of the travel time table changes along with the change of the stratum depth, different step lengths are selected at different stratum depths, namely grid depth intervals of the travel time table are different at different stratum depths, so that the calculation amount of offset imaging and the I/O throughput of a hard disk for executing calculation are reduced on the premise of ensuring the imaging precision, and the offset efficiency is improved.
Description
Technical Field
The invention relates to the technical field of petroleum geophysical exploration, in particular to a method and a device for prestack depth migration by an integral method.
Background
Integral method prestack depth migration is a practical migration method in industrial production, and the realization process mainly comprises two parts, namely traveling calculation and superposition imaging processing along the traveling calculation. Therefore, the accuracy of offset imaging depends mainly on the calculation accuracy when traveling. However, in the integral prestack depth migration process, if the grid depth interval of the travel time table is selected too much, the deviation between interpolation travel time and real travel time is easy to be too large, and the imaging precision is affected; in order to meet imaging accuracy in the prior art, step lengths of a travel time table are set to be the same value which is small enough, namely grid depth intervals of the travel time table are small enough, so that calculation amount of offset imaging and I/O throughput of a hard disk for performing calculation are huge, and offset efficiency is low.
Disclosure of Invention
The embodiment of the invention provides an integral prestack depth migration method, which is used for reducing the calculation amount of migration imaging and the I/O throughput of a hard disk executing calculation on the premise of ensuring the imaging precision and improving the migration efficiency and comprises the following steps:
acquiring seismic data of a target area;
determining the travel time table step length corresponding to each stratum depth according to the seismic data and the calculation parameters;
sampling in the depth direction of the stratum according to the step length of the travel time table corresponding to each stratum depth to obtain the travel time table of the target area;
according to the travel time table step length corresponding to each stratum depth and the travel time table, carrying out linear interpolation on sampling points, and carrying out integral processing to obtain an offset result of seismic data of a target area;
the calculating parameters at least comprise: initial step length of a travel time table, stratum depth and effective frequency band;
determining the travel time table step length corresponding to each stratum depth according to the seismic data and the calculation parameters, wherein the steps comprise:
obtaining the maximum effective frequency of the seismic data of the target area according to the effective frequency band;
determining a minimum value of the velocity model at each formation depth according to the velocity model and the formation depth;
determining an adjustment parameter at each formation depth according to the formation depth;
determining a travel time table step length corresponding to each stratum depth according to the maximum effective frequency of the seismic data of the target area, the minimum value of the velocity model at each stratum depth, the depth interval of the imaging points, the adjusting parameters at each stratum depth and the initial step length of the travel time table;
determining the travel time table step length corresponding to each stratum depth according to the maximum effective frequency of the seismic data of the target area, the minimum value of the velocity model at each stratum depth, the depth interval of the imaging points, the adjusting parameters at each stratum depth and the initial step length of the travel time table according to the following formula:
wherein L represents the formation depth; TMINC [ L ]]Representing the step length of a travel time table corresponding to the stratum depth L; INT represents a rounding function; v. of min [L]Representing the minimum of the velocity model at the formation depth L; para [ L]Represents a conditioning parameter at the formation depth L; f. of max Representing a maximum effective frequency of seismic data for the target area; dz represents the depth interval of the imaging points and is a preset value; TMINC0 represents the initial step size of the travel time table, which is a preset value.
An embodiment of the present invention further provides an integral prestack depth migration apparatus, configured to reduce computation amount of migration imaging and I/O throughput of a hard disk that performs computation on the premise of ensuring imaging accuracy, and improve migration efficiency, where the apparatus includes:
the data acquisition module is used for acquiring seismic data of a target area;
the step length determining module is used for determining the step length of the travel time table corresponding to each stratum depth according to the seismic data and the calculation parameters;
the sampling module is used for sampling in the depth direction of the stratum according to the step length of the travel time table corresponding to each stratum depth to obtain the travel time table of the target area;
the data migration module is used for performing linear interpolation on sampling points according to the travel time table step length corresponding to each stratum depth and the travel time table, and performing integral processing to obtain a migration result of the seismic data of the target area;
the calculation parameters at least comprise: initial step length of a travel time table, stratum depth and effective frequency band;
the step length determining module is specifically configured to:
obtaining the maximum effective frequency of the seismic data of the target area according to the effective frequency band;
determining a minimum value of the velocity model at each formation depth according to a velocity model and the formation depths;
determining an adjustment parameter at each formation depth according to the formation depth;
determining a travel time table step length corresponding to each stratum depth according to the maximum effective frequency of the seismic data of the target area, the minimum value of the velocity model at each stratum depth, the depth interval of the imaging points, the adjusting parameters at each stratum depth and the initial step length of the travel time table;
determining the travel time table step length corresponding to each stratum depth according to the maximum effective frequency of the seismic data of the target area, the minimum value of the velocity model at each stratum depth, the depth interval of the imaging points, the adjusting parameters at each stratum depth and the initial step length of the travel time table according to the following formula:
wherein L represents the formation depth; TMINC [ L ]]Representing the step length of a travel time table corresponding to the stratum depth L; INT represents an integer function; v. of min [L]Represents the minimum of the velocity model at the formation depth L; para [ L]Represents a conditioning parameter at the formation depth L; f. of max Representing a maximum effective frequency of seismic data for the target area; dz represents the depth interval of the imaging points and is a preset value; TMINC0 represents the initial step size of the travel time table, and is a preset value.
The embodiment of the present invention further provides a computer device, which includes a memory, a processor, and a computer program stored in the memory and executable on the processor, and when the processor executes the computer program, the integration prestack depth migration method is implemented.
Embodiments of the present invention also provide a computer-readable storage medium storing a computer program for executing the integration prestack depth migration method.
In the embodiment of the invention, the seismic data of a target area are obtained; determining the travel time table step length corresponding to each stratum depth according to the seismic data and the calculation parameters; sampling in the depth direction of the stratum according to the step length of the travel time table corresponding to each stratum depth to obtain the travel time table of the target area; according to the travel time table step length and the travel time table corresponding to each stratum depth, carrying out linear interpolation on the sampling points, and carrying out integral processing to obtain the migration result of the seismic data of the target area; compared with the step length of a small and fixed travel time table in the prior art, the step length of the travel time table corresponding to each stratum depth is determined, namely the step length of the travel time table changes along with the change of the stratum depth, different step lengths are selected at different stratum depths, namely grid depth intervals of the travel time table are different at different stratum depths, so that the calculation amount of offset imaging and the I/O throughput of a hard disk for executing calculation are reduced on the premise of ensuring the imaging precision, and the offset efficiency is improved.
Drawings
In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the drawings required to be used in the description of the embodiments are briefly introduced below, and it is obvious that the drawings in the description below are only some embodiments of the present invention, and it is obvious for those skilled in the art that other drawings can be obtained according to the drawings without creative efforts.
Fig. 1 is a schematic diagram of an integration method prestack depth migration method according to an embodiment of the present invention.
FIG. 2 is a schematic diagram illustrating a method for determining a step size of a travel time table corresponding to each formation depth according to an embodiment of the present invention.
FIG. 3 (a) is a diagram of an offset common reflection point gather in a fixed step travel time table in an embodiment of the present invention.
FIG. 3 (b) is a schematic diagram of an offset common reflection point gather of a variable step travel time table in an embodiment of the present invention.
FIG. 3 (c) is a cross-sectional view of the residual error between FIG. 3 (a) and FIG. 3 (b) in an embodiment of the present invention.
FIG. 4 (a) is a schematic diagram of an offset common offset gather for a fixed step travel time table in an embodiment of the present invention.
FIG. 4 (b) is a schematic diagram of an offset common offset gather for a variable stride travel time table in an embodiment of the present invention.
FIG. 4 (c) is a cross-sectional view of the residual error between FIG. 4 (a) and FIG. 4 (b) in an embodiment of the present invention.
FIG. 5 is a block diagram of an integrated prestack depth migration apparatus according to an embodiment of the present invention.
Detailed Description
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 obtained by a person skilled in the art without making any creative effort based on the embodiments in the present invention, belong to the protection scope of the present invention.
An embodiment of the present invention provides an integral method prestack depth migration method, which is used to reduce the calculation amount of migration imaging and the I/O throughput of a hard disk executing calculation on the premise of ensuring imaging accuracy, and improve migration efficiency, and as shown in fig. 1, the method includes:
step 101: acquiring seismic data of a target area;
step 102: determining the travel time table step length corresponding to each stratum depth according to the seismic data and the calculation parameters;
step 103: sampling in the depth direction of the stratum according to the step length of the travel time table corresponding to each stratum depth to obtain the travel time table of the target area;
step 104: and performing linear interpolation on the sampling points according to the travel time table step length and the travel time table corresponding to each stratum depth, and performing integral processing to obtain the migration result of the seismic data of the target area.
As can be known from the process shown in FIG. 1, in the embodiment of the present invention, the seismic data of the target area is obtained; determining the travel time table step length corresponding to each stratum depth according to the seismic data and the calculation parameters; sampling in the depth direction of the stratum according to the step length of the travel time table corresponding to each stratum depth to obtain the travel time table of the target area; according to the travel time table step length and the travel time table corresponding to each stratum depth, carrying out linear interpolation on the sampling points, and carrying out integral processing to obtain the migration result of the seismic data of the target area; compared with the step length of a small and fixed travel-time table in the prior art, the step length of the travel-time table corresponding to each stratum depth is determined, namely the step length of the travel-time table changes along with the change of the stratum depth, different step lengths are selected at different stratum depths, namely the grid depth intervals of the travel-time table are different at different stratum depths, so that under the premise of ensuring the imaging precision, the calculation amount of offset imaging and the I/O throughput of a hard disk for executing calculation are reduced, and the offset efficiency is improved.
In specific implementation, firstly, the seismic data of a target area are acquired.
And after the seismic data are obtained, determining the travel time table step length corresponding to each stratum depth according to the seismic data and the calculation parameters. Wherein, calculating the parameters at least comprises: the travel time table initial step size, the formation depth and the effective frequency band. The specific process is shown in fig. 2, and comprises the following steps:
step 201: obtaining the maximum effective frequency of the seismic data of the target area according to the effective frequency band;
step 202: determining the minimum value of the velocity model at each stratum depth according to the velocity model and the stratum depth;
step 203: determining an adjustment parameter at each formation depth according to the formation depth;
step 204: and determining the travel time table step length corresponding to each stratum depth according to the maximum effective frequency of the seismic data of the target area, the minimum value of the velocity model at each stratum depth, the depth interval of the imaging points, the adjusting parameters at each stratum depth and the initial step length of the travel time table.
The above calculation parameters and the speed model may be calculated according to actual engineering, or may be initially given by a technician.
In this embodiment of the present invention, step 204 may determine the step size of the travel time table corresponding to each formation depth according to the following formula:
wherein L represents the formation depth; TMINC [ L ]]Representing the step length of a travel time table corresponding to the stratum depth L; INT represents an integer function; v. of min [L]Representing the minimum of the velocity model at the formation depth L; para [ L]Represents a conditioning parameter at the formation depth L; f. of max Representing a maximum effective frequency of seismic data for the target area; dz represents the imaging point depth interval, which is a preset value; TMINC0 represents the initial step length of the travel time table, which is a preset value, and the engineering generally takes 4 to 8 times of dz value.
In the embodiment, para [ L ] is generally greater than or equal to 2, the larger the adjusting parameter is, the more severe the condition for increasing the step length is, and according to the repeated experiments of the applicant on a large amount of actual seismic data, the following formula is obtained, and the adjusting parameter para [ L ] at each stratum depth is determined according to the stratum depth L:
different adjusting parameter values are obtained according to different stratum depths, when the depth is shallow, the larger adjusting parameter value is selected to ensure the imaging precision of the shallow layer, and when the depth is larger, the smaller adjusting parameter value is selected to ensure the imaging efficiency of the deep layer.
And after the step length of the travel time table corresponding to each stratum depth is determined, sampling in the stratum depth direction according to the step length of the travel time table corresponding to each stratum depth, and obtaining the travel time table of the target area. The specific implementation process may be, for example: according to the initial step length of the travel time table, the travel time table of the target area corresponding to the initial step length is obtained through ray tracing, resampling is carried out according to the travel time table step length corresponding to each stratum depth in the travel time table of the target area corresponding to the initial step length, at the moment, the number of sample points of the travel time table after resampling in the depth direction is greatly reduced, the depth distance between the sample points is not equal any more and is changed into dz multiplied by TMINC [ L ], the travel time table step length TMINC [ L ] is changed along with the change of the stratum depth L, and the depth distance between the sample points is changed along with the change of the stratum depth, so that the travel time table with the variable step length is obtained.
And then, according to the step length of the travel time table corresponding to each stratum depth and the variable step length travel time table, carrying out linear interpolation on sampling points, and carrying out integral processing to obtain a prestack depth migration result of the seismic data of the target area. In the specific embodiment, when linear interpolation is carried out on the sampling points, the number of interpolation sampling points between any two sampling points is TMINC [ L ] -2. In a specific embodiment, a Kirchhoff integral shift algorithm is used when the integration method is performed.
It should be understood that the above calculation formula is only an example, and the above formula may be modified or other formulas or methods may be adopted in the implementation, and those skilled in the art may understand that the above formula is only an example, and the modification of the formula or other solving methods all fall into the protection scope of the present invention, and are not described in detail in the embodiments.
A specific example is given below to illustrate how the integration method prestack depth migration is performed according to an embodiment of the present invention. The method is applied to three-dimensional actual seismic data of a certain area to prove that the method improves the calculation efficiency on the basis of ensuring the migration effect.
The applicant finds that the requirement on the sampling rate is gradually relaxed due to the gradual decrease of the main frequency of the seismic signals with the increase of the depth of the stratum, and the travel time difference of the same depth interval is gradually reduced due to the continuous increase of the speed of the underground medium. Therefore, the depth spacing of the travel time table at shallow layers needs to be small enough to ensure imaging accuracy; but high accuracy imaging can be obtained at deep levels even with large step sizes. Therefore, the embodiment of the invention adopts the travel time table with variable step length, avoids redundant calculation caused by adopting the same small grid, and can obtain the imaging result with similar precision.
Performing Kirchhoff integration method prestack depth migration twice on three-dimensional seismic data of a certain exploration area, and respectively setting the travel time table depth interval as a fixed step length (20 m) and a variable step length (the depth interval is gradually increased from 20m to 120m along with the increase of the depth) in the application to obtain a common reflection point gather as shown in fig. 3 (a) -3 (b) and a common offset gather as shown in fig. 4 (a) -4 (b). Fig. 3 (a) and 4 (a) show the fixed-step integration method prestack depth migration as a reference profile, fig. 3 (b) and 4 (b) show the variable-step migration results, and the two are subtracted to obtain a residual profile, as shown in fig. 3 (c) and 4 (c).
As can be seen from the offset results and the residual profiles shown in fig. 3 (c) and fig. 4 (c), the step-variable travel time table can obtain imaging results with accuracy equivalent to that of a fixed step, and the error is small and almost negligible.
Compare the computational efficiency of the two shifts:
in this test case, the throughput of the variable step travel time table is reduced by 73% compared to the fixed step travel time table. The offset time of the fixed step travel table is 8896s, the offset time of the variable step offset is 4302s, and the offset efficiency is improved by 107%.
When the integral method of the variable-step travel time meter is used for prestack depth migration, the number of stored travel time meter sample points is small, the I/O pressure and migration calculation amount of the travel time meter in the migration process are reduced, and the error of the sample point amplitude value obtained by interpolation is small, so that the method can prove that the migration efficiency can be greatly improved while the migration precision is ensured.
The implementation of the above specific application is only an example, and the rest of the embodiments are not described in detail.
Based on the same inventive concept, embodiments of the present invention further provide an integral method prestack depth migration apparatus, and because the principle of the problem solved by the integral method prestack depth migration apparatus is similar to that of the integral method prestack depth migration method, the apparatus may be implemented by referring to the implementation of the integral method prestack depth migration method, and the repeated parts are not described again, and the specific structure is shown in fig. 5:
a data acquisition module 501, configured to acquire seismic data of a target area;
a step length determining module 502, configured to determine a travel time table step length corresponding to each formation depth according to the seismic data and the calculation parameters;
the sampling module 503 is configured to sample in the depth direction of the stratum according to the step length of the travel time table corresponding to each depth of the stratum to obtain the travel time table of the target area;
and the data migration module 504 is configured to perform linear interpolation on the sampling points according to the travel time table step length and the travel time table corresponding to each stratum depth, and perform integral processing to obtain a migration result of the seismic data of the target area.
In a specific embodiment, calculating the parameters includes at least: travel time table initial step size, formation depth and effective frequency band. The step length determining module 502 is specifically configured to:
obtaining the maximum effective frequency of the seismic data of the target area according to the effective frequency band;
determining the minimum value of the velocity model at each stratum depth according to the velocity model and the stratum depth;
determining an adjustment parameter at each formation depth according to the formation depth;
and determining the travel time table step length corresponding to each stratum depth according to the maximum effective frequency of the seismic data of the target area, the minimum value of the velocity model at each stratum depth, the depth interval of the imaging points, the adjusting parameters at each stratum depth and the initial step length of the travel time table.
In specific implementation, according to the following formula, the travel time table step length corresponding to each stratum depth is determined according to the maximum effective frequency of the seismic data of the target area, the minimum value of the velocity model at each stratum depth, the imaging point depth interval, the adjusting parameter at each stratum depth and the travel time table initial step length:
wherein L represents the formation depth; TMINC [ L ]]Representing the step length of a travel time table corresponding to the stratum depth L; INT represents an integer function; v. of min [L]Represents the minimum of the velocity model at the formation depth L; para [ L]Represents a conditioning parameter at the formation depth L; f. of max Representing a maximum effective frequency of seismic data for the target area; dz represents the imaging point depth interval, which is a preset value; TMINC0 represents the initial step size of the travel time table, which is a preset value.
Determining the adjusting parameters at each stratum depth according to the stratum depths according to the following formula:
wherein para [ L ] represents the conditioning parameter at the formation depth L; l represents the formation depth.
In a specific embodiment, the data shifting module 504 is specifically configured to: when linear interpolation is carried out on the sampling points, the number of interpolation sampling points between any two sampling points is TMINC [ L ] -2.
The embodiment of the present invention further provides a computer device, which includes a memory, a processor, and a computer program stored in the memory and executable on the processor, and when the processor executes the computer program, the integration prestack depth migration method is implemented.
An embodiment of the present invention further provides a computer-readable storage medium, which stores a computer program for executing the integration method pre-stack depth migration method.
In summary, the method and the apparatus for integrating prestack depth migration according to the embodiments of the present invention have the following advantages:
acquiring seismic data of a target area; determining the travel time table step length corresponding to each stratum depth according to the seismic data and the calculation parameters; sampling in the depth direction of the stratum according to the step length of the travel time table corresponding to each stratum depth to obtain the travel time table of the target area; according to the travel time table step length and the travel time table corresponding to each stratum depth, carrying out linear interpolation on the sampling points, and carrying out integral processing to obtain the migration result of the seismic data of the target area; compared with the step length of a small and fixed travel time table in the prior art, the step length of the travel time table corresponding to each stratum depth is determined, namely the step length of the travel time table changes along with the change of the stratum depth, different step lengths are selected at different stratum depths, namely grid depth intervals of the travel time table are different at different stratum depths, so that the calculation amount of offset imaging and the I/O throughput of a hard disk for executing calculation are reduced on the premise of ensuring the imaging precision, and the offset efficiency is improved.
As will be appreciated by one skilled in the art, embodiments of the present invention may be provided as a method, apparatus, 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 is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus 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 above description is only a preferred embodiment of the present invention, and is not intended to limit the present invention, and various modifications and changes may be made to the embodiment of the present invention by those skilled in the art. Any modification, equivalent replacement, or improvement made within the spirit and principle of the present invention should be included in the protection scope of the present invention.
Claims (8)
1. An integral prestack depth migration method, comprising:
acquiring seismic data of a target area;
determining the travel time table step length corresponding to each stratum depth according to the seismic data and the calculation parameters;
sampling in the depth direction of the stratum according to the step length of the travel time table corresponding to each stratum depth to obtain the travel time table of the target area;
according to the travel time table step length corresponding to each stratum depth and the travel time table, performing linear interpolation on sampling points, and performing integral processing to obtain an offset result of seismic data of a target area;
the calculating parameters at least comprise: the initial step length, the stratum depth and the effective frequency band of the travel time table;
determining the travel time table step length corresponding to each stratum depth according to the seismic data and the calculation parameters, wherein the steps comprise:
obtaining the maximum effective frequency of the seismic data of the target area according to the effective frequency band;
determining a minimum value of the velocity model at each formation depth according to a velocity model and the formation depths;
determining an adjustment parameter at each formation depth based on the formation depths;
determining the travel time table step length corresponding to each stratum depth according to the maximum effective frequency of the seismic data of the target area, the minimum value of the velocity model at each stratum depth, the depth interval of the imaging points, the adjusting parameters at each stratum depth and the initial step length of the travel time table;
determining the travel time table step length corresponding to each stratum depth according to the maximum effective frequency of the seismic data of the target area, the minimum value of the velocity model at each stratum depth, the imaging point depth interval, the adjusting parameter at each stratum depth and the travel time table initial step length according to the following formula:
wherein L represents the formation depth; TMINC [ L ]]Representing the step length of a travel time table corresponding to the stratum depth L; INT represents an integer function; v. of min [L]Represents the minimum of the velocity model at the formation depth L; para [ L ]]Represents a conditioning parameter at the formation depth L; f. of max Representing a maximum effective frequency of seismic data for the target area; dz represents the depth interval of the imaging points and is a preset value; TMINC0 represents the initial step length of the travel time table, and is a preset value.
2. The method of claim 1, wherein the adjustment parameter at each formation depth is determined based on the formation depth according to the formula:
wherein para [ L ] represents the conditioning parameter at the formation depth L; l represents the formation depth.
3. The method as claimed in claim 1, wherein according to the travel time table step length corresponding to each stratum depth and the travel time table, performing linear interpolation on sampling points and performing integral processing to obtain the migration result of the seismic data of the target area, the method comprises:
when linear interpolation is carried out on the sampling points, the number of interpolation sampling points between any two sampling points is TMINC [ L ] -2.
4. An integration prestack depth migration apparatus, comprising:
the data acquisition module is used for acquiring seismic data of a target area;
the step length determining module is used for determining the step length of the travel time table corresponding to each stratum depth according to the seismic data and the calculation parameters;
the sampling module is used for sampling in the depth direction of the stratum according to the step length of the travel time table corresponding to each stratum depth to obtain the travel time table of the target area;
the data migration module is used for performing linear interpolation on sampling points according to the travel time table step length corresponding to each stratum depth and the travel time table, and performing integral processing to obtain a migration result of the seismic data of the target area;
the calculating parameters at least comprise: the initial step length, the stratum depth and the effective frequency band of the travel time table;
the step length determining module is specifically configured to:
obtaining the maximum effective frequency of the seismic data of the target area according to the effective frequency band;
determining a minimum value of the velocity model at each formation depth according to the velocity model and the formation depth;
determining an adjustment parameter at each formation depth according to the formation depth;
determining the travel time table step length corresponding to each stratum depth according to the maximum effective frequency of the seismic data of the target area, the minimum value of the velocity model at each stratum depth, the depth interval of the imaging points, the adjusting parameters at each stratum depth and the initial step length of the travel time table;
determining the travel time table step length corresponding to each stratum depth according to the maximum effective frequency of the seismic data of the target area, the minimum value of the velocity model at each stratum depth, the depth interval of the imaging points, the adjusting parameters at each stratum depth and the initial step length of the travel time table according to the following formula:
wherein L represents the formation depth; TMINC [ L ]]Representing the step length of a travel time table corresponding to the stratum depth L; INT represents a rounding function; v. of min [L]Represents the minimum of the velocity model at the formation depth L; para [ L ]]Represents a conditioning parameter at the formation depth L; f. of max Representing a maximum effective frequency of seismic data for the target area; dz represents the imaging point depth interval, which is a preset value; TMINC0 represents the initial step size of the travel time table, which is a preset value.
5. The apparatus of claim 4, wherein the adjustment parameter at each formation depth is determined based on the formation depth according to the following formula:
wherein para [ L ] represents the conditioning parameter at the formation depth L; l represents the formation depth.
6. The apparatus of claim 4, wherein the data migration module is specifically configured to:
when linear interpolation is carried out on the sampling points, the number of interpolation sampling points between any two sampling points is TMINC [ L ] -2.
7. A computer 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 of claims 1 to 3 when executing the computer program.
8. A computer-readable storage medium, characterized in that it stores a computer program for executing the method of any one of claims 1 to 3.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010243070.5A CN113466938B (en) | 2020-03-31 | 2020-03-31 | Integral method prestack depth migration method and device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010243070.5A CN113466938B (en) | 2020-03-31 | 2020-03-31 | Integral method prestack depth migration method and device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113466938A CN113466938A (en) | 2021-10-01 |
| CN113466938B true CN113466938B (en) | 2023-02-28 |
Family
ID=77865403
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010243070.5A Active CN113466938B (en) | 2020-03-31 | 2020-03-31 | Integral method prestack depth migration method and device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113466938B (en) |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104849751A (en) * | 2015-05-15 | 2015-08-19 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Prestack seismic data imaging method |
| CN106842304A (en) * | 2017-01-03 | 2017-06-13 | 中国石油天然气集团公司 | A kind of prestack depth migration method and device |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6519532B2 (en) * | 2001-01-31 | 2003-02-11 | Phillips Petroleum Company | Method and apparatus for 3D depth migration |
| US6687617B2 (en) * | 2001-06-28 | 2004-02-03 | Pgs America, Inc. | Method and system for migration of seismic data |
| US7388808B2 (en) * | 2004-09-23 | 2008-06-17 | Pgs Americas, Inc. | Method for depth migrating seismic data using pre-stack time migration, demigration, and post-stack depth migration |
| CN109085648B (en) * | 2018-07-16 | 2019-09-27 | 中国科学院地质与地球物理研究所 | Prestack depth migration method and device |
-
2020
- 2020-03-31 CN CN202010243070.5A patent/CN113466938B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104849751A (en) * | 2015-05-15 | 2015-08-19 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Prestack seismic data imaging method |
| CN106842304A (en) * | 2017-01-03 | 2017-06-13 | 中国石油天然气集团公司 | A kind of prestack depth migration method and device |
Non-Patent Citations (3)
| Title |
|---|
| Kirchhoff Pre-Stack Depth Migration: effective tool for depth imaging;Rahul Jain et al.;《9th Biennial lnternational Conference & Exposition on Petroleum Geophysics》;20121231;第1-8页 * |
| PS-wave Kirchhoff depth migration and its application to imaging gas clouds;Yue Yubo et al.;《SEG Houston 2013 Annual Meeting》;20131231;第1699-1702页 * |
| 波动方程深度偏移的频率相关变步长延拓方法;张剑锋等;《地球物理学报》;20080131;第51卷(第1期);第221-227页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113466938A (en) | 2021-10-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN102890290B (en) | Pre-stack depth migration method under condition of undulating surface | |
| CN105277978B (en) | A kind of method and device for determining near-surface velocity model | |
| CN103376466B (en) | A kind of multiple reflection drawing method | |
| CN104977607B (en) | The time-domain full waveform inversion method simulated using variable step grid acoustic wavefield | |
| CN109459789B (en) | Time-domain full waveform inversion method based on amplitude decaying and linear interpolation | |
| CN105388520A (en) | Seismic data pre-stack reverse time migration imaging method | |
| CN109946742A (en) | The pure rolling land qP shakes digital simulation method in a kind of TTI medium | |
| CN107765308A (en) | Reconstruct low-frequency data frequency domain full waveform inversion method based on convolution thought Yu accurate focus | |
| CN102540251B (en) | Two-dimensional transverse anisotropy HTI (transversely isotropicmedia with a horizontal symmetry axis) prestack depth migration modeling method and device | |
| CN111948708B (en) | Seismic wave field forward modeling method for dipping in undulating surface of boundary | |
| CN108957553B (en) | Method and device for dynamic correction without stretching distortion by recursion correction of dynamic correction value | |
| CN113466938B (en) | Integral method prestack depth migration method and device | |
| CN102162858A (en) | Method for analyzing dynamic correction velocity by using asymmetric travel time | |
| CN104123449A (en) | Subregion local variable-density non-equidistant dual mesh division method for complex mountainous region | |
| CN113866827B (en) | Interpretation velocity modeling seismic imaging method, system, medium and equipment | |
| CN110161565A (en) | A kind of Reconstruction of seismic data method | |
| CN110376642B (en) | Three-dimensional seismic velocity inversion method based on conical surface waves | |
| CN111929726B (en) | Seismic coherent data volume processing method and device | |
| CN112419493B (en) | Shale reservoir three-dimensional attribute model building method and device | |
| CN102749649B (en) | Determination method of optimal weathered layer speed | |
| CN106353801A (en) | Simulation method and device for 3D Laplace domain acoustic wave equation value | |
| CN111435172A (en) | Method and device for chromatographic static correction | |
| CN112363222A (en) | Post-stack adaptive broadband constraint wave impedance inversion method and device | |
| CN111999767B (en) | Offset imaging method and device for undulating surface | |
| CN111856557A (en) | Method and device for making depth domain synthetic seismic record |
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 |