CN112379413A - Irregular seismic source characterization method and device based on energy spectrum equivalence - Google Patents
Irregular seismic source characterization method and device based on energy spectrum equivalence Download PDFInfo
- Publication number
- CN112379413A CN112379413A CN202011173079.XA CN202011173079A CN112379413A CN 112379413 A CN112379413 A CN 112379413A CN 202011173079 A CN202011173079 A CN 202011173079A CN 112379413 A CN112379413 A CN 112379413A
- Authority
- CN
- China
- Prior art keywords
- seismic source
- grid
- regular
- nearest
- projection
- 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.)
- Granted
Links
- 238000001228 spectrum Methods 0.000 title claims abstract description 154
- 230000001788 irregular Effects 0.000 title claims abstract description 80
- 238000012512 characterization method Methods 0.000 title claims abstract description 70
- 238000000034 method Methods 0.000 claims abstract description 40
- 230000006870 function Effects 0.000 claims description 172
- 238000004590 computer program Methods 0.000 claims description 16
- 238000003860 storage Methods 0.000 claims description 10
- 238000010586 diagram Methods 0.000 description 18
- 238000012545 processing Methods 0.000 description 12
- 230000005284 excitation Effects 0.000 description 10
- 238000005286 illumination Methods 0.000 description 10
- 238000004458 analytical method Methods 0.000 description 9
- 238000003384 imaging method Methods 0.000 description 9
- 238000012986 modification Methods 0.000 description 8
- 230000004048 modification Effects 0.000 description 8
- 230000003595 spectral effect Effects 0.000 description 8
- 230000002411 adverse Effects 0.000 description 7
- 230000000694 effects Effects 0.000 description 7
- 230000014509 gene expression Effects 0.000 description 7
- 238000004088 simulation Methods 0.000 description 6
- 238000004364 calculation method Methods 0.000 description 5
- 230000006872 improvement Effects 0.000 description 3
- 238000013461 design Methods 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000012795 verification Methods 0.000 description 2
- 238000003491 array Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000013213 extrapolation Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
- 230000009466 transformation Effects 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/282—Application of seismic models, synthetic seismograms
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 an irregular seismic source characterization method and device based on energy spectrum equivalence, wherein the method comprises the following steps: establishing a regular gridding model, acquiring plane projection parameters of the seismic source, determining the most adjacent regular projection grid point index of the seismic source and the relative deviation between the actual seismic source plane position and the most adjacent grid projection direction, inputting the half-width parameter of the minimum approximate grid, determining the index range of the regular grid points of equivalent approximation of the seismic source and the weight function of the regular grid points of equivalent approximation of the seismic source, determining the two-dimensional space Fourier spectrum of the weight function of the regular grid points by combining the discretized wave number array, comparing the energy error of the space unit pulse function two-dimensional Fourier spectrum corresponding to the actual seismic source position, judging whether the energy error meets a preset energy error threshold value, and if the minimum approximate grid half-width parameter does not meet the preset energy error threshold, modifying the minimum approximate grid half-width parameter for iteration, and determining the grid node index and the corresponding node weight function value participating in the rule discretization representation. Spatial positioning errors can be reduced.
Description
Technical Field
The invention relates to the technical field of seismic data processing, in particular to an irregular seismic source characterization method and device based on energy spectrum equivalence.
Background
This section is intended to provide a background or context to the embodiments of the invention that are recited in the claims. The description herein is not admitted to be prior art by inclusion in this section.
Conventional seismic data acquisition designs tend to use relatively regular temporal and spatial sampling for signal excitation and reception. Driven by the inherent requirements of the natural geographic conditions, safety and the like of the earth surface and the improvement of underground illumination, excitation points and receiving points in the actual data acquisition process are often distributed in a spatially irregular mode, the irregular characteristic of the excitation points is particularly obvious, and in seismic data processing, a regular and equidistant discretization grid is required to be used for data operation and storage, so that the purpose of improving the calculation efficiency is achieved. The irregular characteristic of the observation system may be further aggravated by the seismic data acquisition design scheme based on compressed sensing and sparse transformation which is gradually developed in recent years.
The seismic wave field simulation is a necessary core link of seismic imaging, illumination analysis and full waveform inversion. The higher the simulation precision of the seismic wave field is, the more beneficial the improvement of the seismic imaging, illumination analysis and full waveform inversion precision is. Common seismic wavefield simulation methods include two categories, ray approximation solution and wave equation direct numerical solution. The wave equation direct numerical solution simulation precision is high, wave field extrapolation and wave field snapshot storage are required to be performed by using a regular and equidistant discretization grid, and the purpose of improving the calculation efficiency is achieved. When the actual seismic source position distribution is in a spatial irregular characteristic, all seismic source points cannot be accurately projected to a preset regular grid point, and a large spatial positioning error may be caused when the actual regular grid point is sparse, so that adverse effects are brought to subsequent imaging, illumination analysis and full waveform inversion.
The current solutions to the above problems are mainly two: (1) regular grids are further encrypted by reducing the grid intervals, and the space error caused by positioning of a single grid is reduced, so that the requirements of grid density, storage and calculation amount are rapidly expanded, the calculation efficiency is low, and the practicability is difficult; (2) and carrying out spatial regularization interpolation and reconstruction processing on the irregular excitation and receiving data to ensure that the processed seismic data excitation and receiving points are consistent with the regular grid. For simple geological structures and lithologic stratums, the change of the seismic wave homomorphic axis is regular, and the method can obtain reasonable results. For the conditions of sparse seismic source distribution or complex structure and lithologic stratum wave fields, the physical stability and the geological rationality of the reconstruction method cannot be verified.
Therefore, how to provide a new solution, which can solve the above technical problems, is a technical problem to be solved in the art.
Disclosure of Invention
The embodiment of the invention provides an irregular seismic source characterization method based on energy spectrum equivalence, which can better reduce larger space positioning errors and reduce adverse effects brought by subsequent imaging, illumination analysis and full waveform inversion, and comprises the following steps:
establishing a regular gridding model, and acquiring a seismic source plane projection parameter;
determining the nearest regular projection grid point index of the seismic source and the relative deviation of the actual seismic source plane position and the nearest grid projection direction according to the seismic source plane projection parameters;
inputting a minimum approximate grid half-width parameter, and determining a regular grid point index range of equivalent approximation of the seismic source and a regular grid point weight function of equivalent approximation of the seismic source according to the index of the nearest regular projection grid point of the seismic source and the relative deviation of the actual plane position of the seismic source and the nearest grid projection direction;
determining a two-dimensional space Fourier spectrum of the regular grid point weight function according to the discretization wave number array, the regular grid point index range of the seismic source equivalent approximation and the regular grid point weight function of the seismic source equivalent approximation;
comparing the energy error of the space unit pulse function two-dimensional Fourier spectrum corresponding to the actual seismic source position with the energy error of the regular grid point weight function two-dimensional space Fourier spectrum, and judging whether the energy error meets a preset energy error threshold value;
and when the energy error does not meet a preset energy error threshold, modifying the minimum approximate grid half-width parameter for iteration, and determining a grid node index and a corresponding node weight function value participating in the rule discretization representation until the energy error meets the preset energy error threshold.
The embodiment of the invention also provides an irregular seismic source characterization device based on energy spectrum equivalence, which comprises the following steps:
the seismic source plane projection parameter acquisition module is used for establishing a regular gridding model and acquiring seismic source plane projection parameters;
the seismic source nearest regular projection grid point index and the relative deviation between the actual seismic source plane position and the nearest grid projection direction are determined by the seismic source nearest regular projection grid point index and the relative deviation between the actual seismic source plane position and the nearest grid projection direction are determined according to the seismic source plane projection parameters;
the seismic source equivalent approximate regular grid point index range and the seismic source equivalent approximate regular grid point weight function determining module is used for inputting a minimum approximate grid half-width parameter, and determining the seismic source equivalent approximate regular grid point index range and the seismic source equivalent approximate regular grid point weight function according to the relative deviation of the seismic source nearest regular projection grid point index and the actual seismic source plane position and the nearest grid projection direction;
the regular grid point weight function two-dimensional space Fourier spectrum determining module is used for determining a regular grid point weight function two-dimensional space Fourier spectrum according to the discretization wave number array, the regular grid point index range of the seismic source equivalent approximation and the regular grid point weight function of the seismic source equivalent approximation;
the energy error judgment module is used for comparing the energy error of the space unit pulse function two-dimensional Fourier spectrum corresponding to the actual seismic source position with the energy error of the regular grid point weight function two-dimensional space Fourier spectrum and judging whether the energy error meets a preset energy error threshold value;
and the grid node index and corresponding node weight function value determination module is used for correcting the minimum approximate grid half-width parameter for iteration when the energy error does not meet the preset energy error threshold value, and determining the grid node index and corresponding node weight function value participating in the rule discretization representation until the energy error meets the preset energy error threshold value.
The embodiment of the invention also provides computer equipment which comprises a memory, a processor and a computer program stored on the memory and capable of running on the processor, wherein the processor executes the computer program to realize the irregular seismic source characterization method based on energy spectrum equivalence.
Embodiments of the present invention further provide a computer-readable storage medium storing a computer program for executing the irregular seismic source characterization method based on energy spectrum equivalence.
The embodiment of the invention provides a method and a device for representing an irregular seismic source based on energy spectrum equivalence.A regular gridding model is established to obtain a seismic source plane projection parameter; then, according to the plane projection parameters of the seismic source, determining the nearest regular projection grid point index of the seismic source and the relative deviation between the actual plane position of the seismic source and the nearest grid projection direction; then inputting a minimum approximate grid half-width parameter, and determining a regular grid point index range of equivalent approximation of the seismic source and a regular grid point weight function of equivalent approximation of the seismic source according to the index of the nearest regular projection grid point of the seismic source and the relative deviation of the actual plane position of the seismic source and the nearest grid projection direction; determining a two-dimensional space Fourier spectrum of the regular grid point weight function according to the discretization wave number array, the regular grid point index range of the seismic source equivalent approximation and the regular grid point weight function of the seismic source equivalent approximation; then comparing the energy error of the space unit pulse function two-dimensional Fourier spectrum corresponding to the actual seismic source position with the energy error of the regular grid point weight function two-dimensional space Fourier spectrum, and judging whether the energy error meets a preset energy error threshold value; and finally, when the energy error does not meet a preset energy error threshold, modifying the minimum approximate grid half-width parameter for iteration, and determining a grid node index and a corresponding node weight function value participating in the rule discretization representation until the energy error meets the preset energy error threshold. The embodiment of the invention establishes a regular discretization equivalent characterization method of an irregular excitation seismic source based on an energy spectrum equivalent principle, utilizes a method for characterizing the irregular seismic source by utilizing a regular discretization grid combined seismic source, finally generates a regular grid index and grid point weight which can be used for approximating the equivalent characterization, and can better reduce space positioning errors. For actual local sparse points, larger space positioning errors possibly caused are reduced, and adverse effects caused by subsequent imaging, illumination analysis and full waveform inversion are reduced. The method can calculate the weighting function of the seismic source combination in the corresponding rule discretization format under the condition of the given rule grid parameters and the plane coordinates of any seismic source, so that the error between the spectral energy of the corresponding equivalent seismic source combination and the actual seismic source spectral energy meets the given precision requirement.
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 or the prior art will be briefly described below, it is obvious that the drawings in the following description are only some embodiments of the present invention, and for those skilled in the art, other drawings can be obtained according to the drawings without creative efforts. In the drawings:
FIG. 1 is a schematic diagram of an irregular seismic source characterization method based on energy spectrum equivalence according to an embodiment of the present invention.
Fig. 2 is a graph showing a relationship between a regular grid model and an irregular seismic source planar projection of an irregular seismic source characterization method based on energy spectrum equivalence according to an embodiment of the present invention.
FIG. 3 is a graph of the index ranges of regular grid points of a seismic source equivalent approximation based on an irregular seismic source characterization method of energy spectrum equivalence according to an embodiment of the present invention.
FIG. 4 is a schematic diagram of values of a regular grid point weight function of a seismic source equivalent approximation based on an irregular seismic source characterization method of energy spectrum equivalence according to an embodiment of the present invention.
Fig. 5 is an energy error graph of a two-dimensional space fourier spectrum of a space unit impulse function and a two-dimensional space fourier spectrum of a weight function of a regular grid point corresponding to an actual seismic source position in an irregular seismic source characterization method based on energy spectrum equivalence according to an embodiment of the present invention.
FIG. 6 is a schematic diagram of a computer device running a method for irregular seismic source characterization based on energy spectrum equivalence implemented in the present invention.
FIG. 7 is a schematic diagram of an irregular seismic source characterization apparatus based on energy spectrum equivalence according to an embodiment of the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the embodiments of the present invention more apparent, the embodiments of the present invention are further described in detail below with reference to the accompanying drawings. The exemplary embodiments and descriptions of the present invention are provided to explain the present invention, but not to limit the present invention.
Fig. 1 is a schematic diagram of an irregular seismic source characterization method based on energy spectrum equivalence according to an embodiment of the present invention, and as shown in fig. 1, an irregular seismic source characterization method based on energy spectrum equivalence according to an embodiment of the present invention is provided, which can better reduce a large spatial positioning error and reduce adverse effects caused by subsequent imaging, illumination analysis, and full waveform inversion, and includes:
step 101: establishing a regular gridding model, and acquiring a seismic source plane projection parameter;
step 102: determining the nearest regular projection grid point index of the seismic source and the relative deviation of the actual seismic source plane position and the nearest grid projection direction according to the seismic source plane projection parameters;
step 103: inputting a minimum approximate grid half-width parameter, and determining a regular grid point index range of equivalent approximation of the seismic source and a regular grid point weight function of equivalent approximation of the seismic source according to the index of the nearest regular projection grid point of the seismic source and the relative deviation of the actual plane position of the seismic source and the nearest grid projection direction;
step 104: determining a two-dimensional space Fourier spectrum of the regular grid point weight function according to the discretization wave number array, the regular grid point index range of the seismic source equivalent approximation and the regular grid point weight function of the seismic source equivalent approximation;
step 105: comparing the energy error of the space unit pulse function two-dimensional Fourier spectrum corresponding to the actual seismic source position with the energy error of the regular grid point weight function two-dimensional space Fourier spectrum, and judging whether the energy error meets a preset energy error threshold value;
step 106: and when the energy error does not meet a preset energy error threshold, modifying the minimum approximate grid half-width parameter for iteration, and determining a grid node index and a corresponding node weight function value participating in the rule discretization representation until the energy error meets the preset energy error threshold.
The embodiment of the invention provides an irregular seismic source characterization method based on energy spectrum equivalence, which comprises the steps of firstly establishing a regular gridding model and obtaining a seismic source plane projection parameter; then, according to the plane projection parameters of the seismic source, determining the nearest regular projection grid point index of the seismic source and the relative deviation between the actual plane position of the seismic source and the nearest grid projection direction; then inputting a minimum approximate grid half-width parameter, and determining a regular grid point index range of equivalent approximation of the seismic source and a regular grid point weight function of equivalent approximation of the seismic source according to the index of the nearest regular projection grid point of the seismic source and the relative deviation of the actual plane position of the seismic source and the nearest grid projection direction; determining a two-dimensional space Fourier spectrum of the regular grid point weight function according to the discretization wave number array, the regular grid point index range of the seismic source equivalent approximation and the regular grid point weight function of the seismic source equivalent approximation; then comparing the energy error of the space unit pulse function two-dimensional Fourier spectrum corresponding to the actual seismic source position with the energy error of the regular grid point weight function two-dimensional space Fourier spectrum, and judging whether the energy error meets a preset energy error threshold value; and finally, when the energy error does not meet a preset energy error threshold, modifying the minimum approximate grid half-width parameter for iteration, and determining a grid node index and a corresponding node weight function value participating in the rule discretization representation until the energy error meets the preset energy error threshold. The embodiment of the invention establishes a regular discretization equivalent characterization method of an irregular excitation seismic source based on an energy spectrum equivalent principle, utilizes a method for characterizing the irregular seismic source by utilizing a regular discretization grid combined seismic source, finally generates a regular grid index and grid point weight which can be used for approximating the equivalent characterization, and can better reduce space positioning errors. For actual local sparse points, larger space positioning errors possibly caused are reduced, and adverse effects caused by subsequent imaging, illumination analysis and full waveform inversion are reduced. The method can calculate the weighting function of the seismic source combination in the corresponding rule discretization format under the condition of the given rule grid parameters and the plane coordinates of any seismic source, so that the error between the spectral energy of the corresponding equivalent seismic source combination and the actual seismic source spectral energy meets the given precision requirement.
The invention relates to an irregular seismic source characterization method in seismic wave field forward modeling, in particular to a discretization equivalent seismic source characterization method under the condition that the plane projection position of a seismic source of an actual observation system is not coincident with regular grid nodes. The embodiment of the invention is based on an energy spectrum equivalence principle, a method for representing an irregular seismic source by utilizing a regular discretization grid combined seismic source is used for finally generating regular grid indexes and grid point weights which can be used for approximate equivalent representation, and when the irregular seismic source representing method based on energy spectrum equivalence provided by the embodiment of the invention is specifically implemented, the method can comprise the following steps:
establishing a regular gridding model, and acquiring a seismic source plane projection parameter; determining the nearest regular projection grid point index of the seismic source and the relative deviation of the actual seismic source plane position and the nearest grid projection direction according to the seismic source plane projection parameters; inputting a minimum approximate grid half-width parameter, and determining a regular grid point index range of equivalent approximation of the seismic source and a regular grid point weight function of equivalent approximation of the seismic source according to the index of the nearest regular projection grid point of the seismic source and the relative deviation of the actual plane position of the seismic source and the nearest grid projection direction; determining a two-dimensional space Fourier spectrum of the regular grid point weight function according to the discretization wave number array, the regular grid point index range of the seismic source equivalent approximation and the regular grid point weight function of the seismic source equivalent approximation; comparing the energy error of the space unit pulse function two-dimensional Fourier spectrum corresponding to the actual seismic source position with the energy error of the regular grid point weight function two-dimensional space Fourier spectrum, and judging whether the energy error meets a preset energy error threshold value; and when the energy error does not meet a preset energy error threshold, modifying the minimum approximate grid half-width parameter for iteration, and determining a grid node index and a corresponding node weight function value participating in the rule discretization representation until the energy error meets the preset energy error threshold.
In a specific implementation of the irregular seismic source characterization method based on energy spectrum equivalence provided by the embodiment of the present invention, in an embodiment, the establishing a regular grid model to obtain the seismic source planar projection parameters includes:
establishing a regular gridding model according to the geological task, the seismic resolution requirement and the model space range; wherein the model space range comprises all actual seismic source positions to be calculated and all discretized equivalent seismic source positions;
and acquiring the plane projection parameters of the seismic source according to the regular gridding model.
In an embodiment, a regular gridding model is first built that satisfies the geological mission, seismic resolution requirements, and model spatial range, which should be large enough to encompass all actual source locations to be computed and all discretized equivalent source locations.
Fig. 2 is a graph showing a relationship between a regular grid model and a planar projection of an irregular seismic source based on an energy spectrum equivalent irregular seismic source characterization method according to an embodiment of the present invention, as shown in fig. 2, the relationship between the regular grid model and the planar projection of the irregular seismic source is shown, a hollow circle represents a regular grid point position, and an asterisk represents an actual irregular seismic source position.
In a specific implementation of the irregular seismic source characterization method based on energy spectrum equivalence provided by the embodiment of the present invention, in an embodiment, the obtaining of the seismic source planar projection parameters according to the regular grid model includes:
determining a horizontal projection plane of the regular gridding model according to the regular gridding model;
and acquiring the reference origin coordinate, the plane grid distance, the grid number in the x and y directions, the projection sequence number and the seismic source plane projection coordinate of the horizontal projection plane according to the horizontal projection plane of the regular gridding model.
In an embodiment, according to the established regular grid model, the horizontal projection plane xoy of the regular grid model can be determined, and the reference origin coordinate (x) thereof is obtained0,y0) The plane grid spacing (Deltax, Deltay), the number of grids in the x and y directions are NxAnd NyThe projection sequence numbers are (1,2, …, N)x) And (1,2, …, N)y) Corresponding to the actual single shot forward acting source plane projection coordinates(s)x,sy)。
In an embodiment of the invention, when the method for characterizing an irregular seismic source based on energy spectrum equivalence provided by the embodiment of the present invention is implemented specifically, the determining, according to the seismic source plane projection parameters, an index of a nearest regular projection grid point of the seismic source and a relative deviation between an actual seismic source plane position and a nearest grid projection direction includes:
determining the nearest regular projection grid point index of the seismic source according to the plane projection coordinate of the seismic source, the reference origin coordinate and the plane grid spacing;
and determining the relative deviation of the actual plane position of the seismic source and the projection direction of the nearest grid according to the index of the nearest regular projection grid point of the seismic source, the plane projection coordinate of the seismic source, the reference origin coordinate and the plane grid spacing.
In an embodiment, the source plane projection coordinates(s) are acquired according to the foregoingx,sy) Reference origin coordinate (x)0,y0) And plane grid spacing (Δ x, Δ y), calculating the index (i) of the nearest regular projection grid point of the seismic source in the two coordinate axis directions of the projection planex0,iy0)。
In a specific implementation of the irregular seismic source characterization method based on energy spectrum equivalence provided by the embodiment of the present invention, in an embodiment, the nearest neighbor regular projection grid point index of the seismic source is determined as follows:
ix0=NINT[(sx-x0)/Δx)]+1 ①
iy0=NINT[(sy-y0)/Δy)]+1 ②
wherein (i)x0,iy0) Projecting a grid point index for the nearest rule of the seismic source; (s)x,sy) Projecting coordinates for the seismic source plane; (x)0,y0) As reference origin coordinates; (Δ x, Δ y) is the planar grid spacing; NINT is rounded according to rounding.
The foregoing expressions for determining the index of the nearest regular projected grid points of the seismic source are given as examples, and it will be understood by those skilled in the art that the above equations may be modified in some form and other parameters or data may be added as needed, or other specific equations may be provided, and such modifications are intended to fall within the scope of the present invention.
Obtaining the nearest regular projection grid point index (i) of the seismic source in the calculationx0,iy0) Then, combining the plane projection coordinates of the seismic source, the coordinates of the reference origin and the plane grid spacing, calculating the relative deviation (e) between the actual plane position of the seismic source and the projection direction of the nearest gridx,ey)。
In an embodiment of the invention, when the irregular seismic source characterization method based on energy spectrum equivalence is implemented, the relative deviation between the actual seismic source plane position and the projection direction of the nearest grid is determined as follows:
wherein (e)x,ey) Relative deviation between the actual seismic source plane position and the projection direction of the nearest grid is obtained; (i)x0,iy0) Projecting a grid point index for the nearest rule of the seismic source; (s)x,sy) Projecting coordinates for the seismic source plane; (x)0,y0) As reference origin coordinates; (Δ x, Δ y) is the planar grid spacing.
While the foregoing expressions for determining the relative deviation of the actual source plane position from the projection direction of the nearest neighbor grid are provided as examples, those skilled in the art will appreciate that the above equations may be modified in some manner and other parameters or data may be added as needed, or other specific equations may be provided, and such modifications are intended to fall within the scope of the present invention.
In a specific implementation of the method for characterizing an irregular seismic source based on energy spectrum equivalence provided by the embodiment of the present invention, in an embodiment, the determining a range of indexes of regular grid points of equivalent approximation of the seismic source and a weight function of the regular grid points of equivalent approximation of the seismic source according to a maximum half-width parameter of a minimum approximate grid and a relative deviation between an actual plane position of the seismic source and a projection direction of a maximum approximate grid includes:
inputting a minimum approximate grid half-width parameter, and determining a regular grid point index range of equivalent approximation of the seismic source according to the nearest regular projection grid point index of the seismic source;
determining the unidirectional grid weight of the regular grid points of the seismic source equivalent approximation according to the relative deviation of the actual seismic source plane position and the projection direction of the nearest grid, the index range of the regular grid points of the seismic source equivalent approximation, the index of the nearest regular projection grid points of the seismic source and the half-width parameter of the minimum approximation grid;
and combining the unidirectional grid weights of the regular grid points of the equivalent approximation of the seismic source to determine a weight function of the regular grid points of the equivalent approximation of the seismic source.
In an embodiment, the minimum approximate grid half-width parameter d for the inputwRespectively calculating the index ranges (i) of all regular grid points participating in the equivalent approximation of the seismic source in the x direction and the y direction according to the index of the nearest regular projection grid point of the seismic sourcex,iy) (ii) a In an embodiment of the invention, when the irregular seismic source characterization method based on energy spectrum equivalence is implemented, the regular grid point index range of the equivalent approximation of the seismic source is determined as followsEnclosing:
ix0-dw≤ix≤ix0+d w ⑤
iy0-dw≤iy≤iy0+dw ⑥
wherein (i)x,iy) Index ranges in x-and y-directions for regular grid points, i, which are equivalent approximations of the seismic sourcexAnd iyIs a positive integer; dwIs a minimum approximate grid half-width parameter; (i)x0,iy0) The grid point index is projected for the nearest neighbor rule of the seismic source.
The foregoing expressions for determining the range of indices of regular grid points of equivalent approximation of a seismic source are merely exemplary, and it will be understood by those skilled in the art that the above equations may be modified in some manner and other parameters or data may be added as needed, or other specific equations may be provided, and such modifications are intended to fall within the scope of the present invention.
FIG. 3 is an index range diagram of regular grid points of equivalent approximation of a seismic source based on an irregular seismic source characterization method of energy spectrum equivalence according to an embodiment of the present invention, as shown in FIG. 3, nearest regular grid points calculated according to the half-width of the minimum approximation grid are represented by black dots; and then the relative deviation (e) from the nearest grid projection direction according to the actual source plane positionx,ey) Regular grid point index range (i) of equivalent approximation to seismic sourcex,iy) Source nearest regular projection grid point index (i)x0,iy0) And a minimum approximate grid half-width parameter dwRespectively calculating all the unidirectional grid weights C of the regular grid points participating in equivalent approximation of the seismic sourcex(ix) And Cy(iy) (ii) a Regular grid points participating in equivalent approximation in the model are represented by black dots; in a specific implementation of the irregular seismic source characterization method based on energy spectrum equivalence provided by the embodiment of the present invention, in an embodiment, the unidirectional grid weights of the regular grid points of the equivalent approximation of the seismic source are determined as follows:
wherein, Cx(ix) And Cy(iy) A one-way grid weight of a regular grid point which is equivalent and approximate to the seismic source; (i)x,iy) Index ranges in x-and y-directions for regular grid points, i, which are equivalent approximations of the seismic sourcexAnd iyIs a positive integer; dwIs a minimum approximate grid half-width parameter; (i)x0,iy0) Projecting a grid point index for the nearest rule of the seismic source; (e)x,ey) Relative deviation between the actual seismic source plane position and the projection direction of the nearest grid is obtained; sinc is a sine function notation; cos is the cosine function notation.
The foregoing expressions for determining the unidirectional grid weights of regular grid points of an equivalent approximation of a seismic source are merely exemplary, and it will be understood by those skilled in the art that the above equations may be modified in some manner and other parameters or data may be added as needed, or other specific equations may be provided, and such modifications are intended to fall within the scope of the present invention.
After the unidirectional grid weights of the regular grid points of the equivalent approximation of the seismic source are obtained, combining the unidirectional grid weights to obtain all the weight functions of the regular grid points participating in the equivalent approximation of the seismic source;
FIG. 4 is a schematic diagram showing values of regular grid point weight functions of seismic source equivalent approximation based on an irregular seismic source characterization method of energy spectrum equivalence, and a weight function C of all regular grid points participating in equivalent approximation is calculatedxy(ix,iy) In particular, when the irregular seismic source characterization method based on energy spectrum equivalence provided by the embodiment of the present invention is implemented, in an embodiment, a regular grid point weight function of seismic source equivalent approximation is determined as follows:
Cxy(ix,iy)=Cx(ix)Cy(iy) ⑨
wherein, Cxy(ix,iy) A regular grid point weight function that is an equivalent approximation of the seismic source; cx(ix) And Cy(iy) Unidirectional grid weights for regular grid points that are equivalent approximations of the seismic source.
The foregoing expressions of the weighting function for regular grid points for determining an equivalent approximation of a seismic source are merely exemplary, and it will be understood by those skilled in the art that the above equations may be modified in some manner and other parameters or data may be added or other specific equations may be provided as desired, and such modifications are intended to fall within the scope of the present invention.
According to the input discretization wave number arrays (kx, ky) in the x-axis direction and the y-axis direction, the dimensions are p and q respectively, and the subscripts of corresponding array elements are represented by l and m respectively. Computing a regular grid point weight function Cxy(ix,iy) At each pair (k)xl,kym) Two-dimensional Fourier spectrum F (kx) of (A)l,kym);
In a specific implementation of the irregular seismic source characterization method based on energy spectrum equivalence provided by the embodiment of the present invention, in an embodiment, a two-dimensional spatial fourier spectrum of a weight function of a regular grid point is determined as follows:
wherein, F (kx)l,kym) Two-dimensional spatial Fourier spectrum of a weight function of a regular grid point; cxy(ix,iy) A regular grid point weight function that is an equivalent approximation of the seismic source; (i)x,iy) Index ranges in x-and y-directions for regular grid points, i, which are equivalent approximations of the seismic sourcexAnd iyIs a positive integer; (i)x0,iy0) Projecting a grid point index for the nearest rule of the seismic source; (e)x,ey) For the actual source plane position and the nearest grid projection directionRelative deviation; (kx, ky) is a discretization wave number array in the x-axis direction and the y-axis direction, the dimensions are p and q respectively, and subscripts of corresponding array elements are represented by l and m respectively; j is a unit imaginary number; exp is an exponential function.
The aforementioned expression for determining the two-dimensional spatial fourier spectrum of the regular grid point weight function is an example, and those skilled in the art will understand that, in implementation, the above formula may be modified in some form and other parameters or data may be added as needed, or other specific formulas may be provided, and these modifications are all within the scope of the present invention.
Fig. 5 is an energy error diagram of a two-dimensional fourier spectrum of a space unit pulse function and a two-dimensional fourier spectrum of a weight function of a regular grid point corresponding to an actual seismic source position of an irregular seismic source characterization method based on energy spectrum equivalence, as shown in fig. 5, next, an energy error between a two-dimensional fourier spectrum energy of the weight function and an impulse spectrum is calculated, and the two-dimensional spatial fourier spectrum error of the equivalent seismic source spectrum is compared with the energy error between the impulse seismic source spectrum, which specifically includes: comparing the energy error of the space unit pulse function two-dimensional Fourier spectrum corresponding to the actual seismic source position with the energy error of the regular grid point weight function two-dimensional space Fourier spectrum; in a specific implementation of the irregular seismic source characterization method based on energy spectrum equivalence provided by the embodiment of the present invention, in an embodiment, the energy error of the two-dimensional fourier spectrum of the spatial unit impulse function corresponding to the actual seismic source position and the two-dimensional spatial fourier spectrum of the regular grid point weight function is compared as follows:
wherein, F (kx)l,kym) Two-dimensional spatial Fourier spectrum of a weight function of a regular grid point; e (kx)l,kym) A two-dimensional fourier spectrum of the unit impulse function in space corresponding to the actual source location.
The aforementioned expression comparing the energy error of the two-dimensional space fourier spectrum of the pulse function in space unit corresponding to the actual source position with the two-dimensional space fourier spectrum of the weighting function of the regular grid points is an example, and it can be understood by those skilled in the art that, in implementation, the above formula may be modified in some form and other parameters or data may be added as needed, or other specific formulas may be provided, and these modifications all fall within the protection scope of the present invention.
After the energy error of the space unit impulse function two-dimensional Fourier spectrum corresponding to the actual seismic source position and the regular grid point weight function two-dimensional space Fourier spectrum is obtained,
judging whether all the discretization Fourier spectrum energy errors meet a preset energy error threshold value or not;
and when the energy error meets a preset energy error threshold, finishing the estimation of the equivalent approximate weight function, and directly outputting the grid node index participating in the rule discretization representation and the corresponding node weight function value.
When the energy error does not meet the preset energy error threshold, modifying the minimum approximate grid half-width parameter for iterative repetition until the energy error meets the preset energy error threshold, and determining a grid node index (i) participating in the rule discretization characterizationx,iy) And corresponding node weight function value Cxy(ix,iy)。
The embodiment of the invention also provides an implementation step of the irregular seismic source characterization method based on energy spectrum equivalence, which mainly comprises the following steps:
(1) establishing a regular grid model, and acquiring seismic source plane projection parameters;
(2) determining the nearest regular grid point index of the seismic source and the relative deviation of the actual coordinates, mainly comprising the following steps: determining the nearest regular projection grid point index of the seismic source and the relative deviation of the actual seismic source plane position and the nearest grid projection direction according to the seismic source plane projection parameters;
(3) determining a regular grid point index range and a grid weight function for equivalent approximation, mainly comprising: inputting a minimum approximate grid half-width parameter, and determining a regular grid point index range of equivalent approximation of the seismic source and a regular grid point weight function of equivalent approximation of the seismic source according to the index of the nearest regular projection grid point of the seismic source and the relative deviation of the actual plane position of the seismic source and the nearest grid projection direction;
(4) calculating a two-dimensional spatial Fourier spectrum of a weight function, which mainly comprises the following steps: determining a two-dimensional space Fourier spectrum of the regular grid point weight function according to the discretization wave number array, the regular grid point index range of the seismic source equivalent approximation and the regular grid point weight function of the seismic source equivalent approximation;
(5) comparing the equivalent seismic source spectrum two-dimensional space Fourier spectrum with the pulse seismic source spectrum error, mainly comprising: comparing the energy error of the space unit pulse function two-dimensional Fourier spectrum corresponding to the actual seismic source position with the energy error of the regular grid point weight function two-dimensional space Fourier spectrum;
(6) judging whether all the discretization Fourier spectrum energies meet the precision requirement or not according to a preset energy error threshold, and finishing the estimation of the equivalent approximate weight function if the discretization Fourier spectrum energies meet the precision requirement; and when the energy error does not meet the preset energy error threshold, modifying the minimum approximate grid half-width parameter for iteration, and performing the processing from (3) to (5) again until the discretization Fourier spectrum energy error meets the preset energy error threshold requirement.
(7) Outputting the grid node index (i) meeting the precision requirement and participating in the rule discretization characterizationx,iy) And corresponding node weight function value Cxy(ix,iy)。
Table 1 shows the grid node index and the corresponding node weight function of the rule discretization representation that the output meets the requirements.
TABLE 1
The embodiment of the invention establishes a regular discretization equivalent characterization method of an irregular excitation seismic source based on an energy spectrum equivalent principle, utilizes a method for characterizing the irregular seismic source by utilizing a regular discretization grid combined seismic source, finally generates a regular grid index and grid point weight which can be used for approximating the equivalent characterization, and can better reduce space positioning errors. For actual local sparse points, larger space positioning errors possibly caused are reduced, and adverse effects caused by subsequent imaging, illumination analysis and full waveform inversion are reduced. The method can calculate the weighting function of the seismic source combination in the corresponding rule discretization format under the condition of the given rule grid parameters and the plane coordinates of any seismic source, so that the error between the spectral energy of the corresponding equivalent seismic source combination and the actual seismic source spectral energy meets the given precision requirement. The method has wide application prospect because the regularization method is an important basic step in the current data processing in the seismic data processing. The feasibility of the method is proved through theoretical data simulation verification.
Fig. 6 is a schematic diagram of a computer device for executing an irregular seismic source characterization method based on energy spectrum equivalence implemented by the present invention, and as shown in fig. 6, an embodiment of the present invention further provides a computer device including a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor implements the irregular seismic source characterization method based on energy spectrum equivalence when executing the computer program.
Embodiments of the present invention further provide a computer-readable storage medium storing a computer program for implementing the irregular seismic source characterization method based on energy spectrum equivalence.
The embodiment of the invention also provides an irregular seismic source characterization device based on energy spectrum equivalence, and the irregular seismic source characterization device is described in the following embodiments. Because the principle of solving the problems of the device is similar to the irregular seismic source characterization method based on energy spectrum equivalence, the implementation of the device can refer to the implementation of the irregular seismic source characterization method based on energy spectrum equivalence, and repeated details are omitted.
Fig. 7 is a schematic diagram of an irregular seismic source characterization device based on energy spectrum equivalence according to an embodiment of the present invention, and as shown in fig. 7, an irregular seismic source characterization device based on energy spectrum equivalence is further provided in an embodiment of the present invention, and may include:
the seismic source plane projection parameter acquisition module 701 is used for establishing a regular gridding model and acquiring seismic source plane projection parameters;
a seismic source nearest regular projection grid point index and actual seismic source plane position relative deviation from nearest grid projection direction determination module 702, configured to determine, according to the seismic source plane projection parameters, a seismic source nearest regular projection grid point index and actual seismic source plane position relative deviation from nearest grid projection direction;
a seismic source equivalent approximate regular grid point index range and seismic source equivalent approximate regular grid point weight function determination module 703, configured to input a minimum approximate grid half-width parameter, and determine a seismic source equivalent approximate regular grid point index range and a seismic source equivalent approximate regular grid point weight function according to a relative deviation between a seismic source nearest regular projection grid point index and an actual seismic source plane position and a nearest grid projection direction;
a regular grid point weight function two-dimensional space fourier spectrum determination module 704, configured to determine a regular grid point weight function two-dimensional space fourier spectrum according to the discretized wave number array, the regular grid point index range of the seismic source equivalent approximation, and the regular grid point weight function of the seismic source equivalent approximation;
the energy error judgment module 705 is configured to compare an energy error of a two-dimensional space fourier spectrum of a spatial unit pulse function corresponding to an actual seismic source position with a two-dimensional space fourier spectrum of a regular grid point weight function, and judge whether the energy error meets a preset energy error threshold;
and a grid node index and corresponding node weight function value determining module 706, configured to modify the minimum approximate grid half-width parameter for iteration when the energy error does not meet the preset energy error threshold, and determine the grid node index and corresponding node weight function value participating in the rule discretization characterization until the energy error meets the preset energy error threshold.
In an embodiment of the invention, when the irregular seismic source characterization device based on energy spectrum equivalence is implemented specifically, the seismic source planar projection parameter obtaining module is specifically configured to:
establishing a regular gridding model according to the geological task, the seismic resolution requirement and the model space range; wherein the model space range comprises all actual seismic source positions to be calculated and all discretized equivalent seismic source positions;
and acquiring the plane projection parameters of the seismic source according to the regular gridding model.
In an embodiment of the invention, when the irregular seismic source characterization device based on energy spectrum equivalence is implemented specifically, the seismic source planar projection parameter obtaining module is further configured to:
determining a horizontal projection plane of the regular gridding model according to the regular gridding model;
and acquiring the reference origin coordinate, the plane grid distance, the grid number in the x and y directions, the projection sequence number and the seismic source plane projection coordinate of the horizontal projection plane according to the horizontal projection plane of the regular gridding model.
In an embodiment of the invention, when the irregular seismic source characterization device based on energy spectrum equivalence provided by the embodiment of the present invention is implemented specifically, the module for determining the nearest regular projection grid point index of the seismic source and the relative deviation between the actual seismic source plane position and the nearest grid projection direction is specifically configured to:
determining the nearest regular projection grid point index of the seismic source according to the plane projection coordinate of the seismic source, the reference origin coordinate and the plane grid spacing;
and determining the relative deviation of the actual plane position of the seismic source and the projection direction of the nearest grid according to the index of the nearest regular projection grid point of the seismic source, the plane projection coordinate of the seismic source, the reference origin coordinate and the plane grid spacing.
In an embodiment, when the irregular seismic source characterization device based on energy spectrum equivalence provided by the embodiment of the present invention is implemented specifically, the module for determining the nearest-to-regular projection grid point index of the seismic source and the relative deviation between the actual seismic source plane position and the nearest-to-regular projection grid point direction is further configured to determine the nearest-to-regular projection grid point index of the seismic source as follows:
ix0=NINT[(sx-x0)/Δx)]+1
iy0=NINT[(sy-y0)/Δy)]+1
wherein (i)x0,iy0) Projecting a grid point index for the nearest rule of the seismic source; (s)x,sy) Projecting coordinates for the seismic source plane; (x)0,y0) As reference origin coordinates; (Δ x, Δ y) is the planar grid spacing; NINT is rounded according to rounding.
In an embodiment, when the irregular seismic source characterization device based on energy spectrum equivalence provided by the embodiment of the present invention is implemented specifically, the module for determining the index of the nearest regular projection grid point of the seismic source and the relative deviation between the actual seismic source plane position and the nearest projection direction of the grid is further configured to determine the relative deviation between the actual seismic source plane position and the nearest projection direction of the grid as follows:
wherein (e)x,ey) Relative deviation between the actual seismic source plane position and the projection direction of the nearest grid is obtained; (i)x0,iy0) Projecting a grid point index for the nearest rule of the seismic source; (s)x,sy) Projecting coordinates for the seismic source plane; (x)0,y0) As reference origin coordinates; (Δ x, Δ y) is the planar grid spacing.
In an embodiment of the invention, when the irregular seismic source characterization device based on energy spectrum equivalence provided in the embodiment of the present invention is implemented specifically, the module for determining the index range of the regular grid points of the equivalent approximation of the seismic source and the weight function of the regular grid points of the equivalent approximation of the seismic source is specifically configured to:
inputting a minimum approximate grid half-width parameter, and determining a regular grid point index range of equivalent approximation of the seismic source according to the nearest regular projection grid point index of the seismic source;
determining the unidirectional grid weight of the regular grid points of the seismic source equivalent approximation according to the relative deviation of the actual seismic source plane position and the projection direction of the nearest grid, the index range of the regular grid points of the seismic source equivalent approximation, the index of the nearest regular projection grid points of the seismic source and the half-width parameter of the minimum approximation grid;
and combining the unidirectional grid weights of the regular grid points of the equivalent approximation of the seismic source to determine a weight function of the regular grid points of the equivalent approximation of the seismic source.
In an embodiment of the invention, when the irregular seismic source characterization device based on energy spectrum equivalence provided by the embodiment of the present invention is implemented specifically, the module for determining the index range of the regular grid points of the equivalent approximation of the seismic source and the weight function of the regular grid points of the equivalent approximation of the seismic source is further configured to determine the index range of the regular grid points of the equivalent approximation of the seismic source according to the following manner:
ix0-dw≤ix≤ix0+dw
iy0-dw≤iy≤iy0+dw
wherein (i)x,iy) Index ranges in x-and y-directions for regular grid points, i, which are equivalent approximations of the seismic sourcexAnd iyIs a positive integer; dwIs a minimum approximate grid half-width parameter; (i)x0,iy0) The grid point index is projected for the nearest neighbor rule of the seismic source.
In an embodiment of the invention, when the irregular seismic source characterization device based on energy spectrum equivalence provided by the embodiment of the present invention is implemented specifically, the module for determining the index range of the regular grid points of the seismic source equivalent approximation and the weight function of the regular grid points of the seismic source equivalent approximation is further configured to determine the unidirectional grid weight of the regular grid points of the seismic source equivalent approximation according to the following manner:
wherein, Cx(ix) And Cy(iy) A one-way grid weight of a regular grid point which is equivalent and approximate to the seismic source; (i)x,iy) Index ranges in x-and y-directions for regular grid points, i, which are equivalent approximations of the seismic sourcexAnd iyIs a positive integer; dwIs a minimum approximate grid half-width parameter; (i)x0,iy0) Projecting a grid point index for the nearest rule of the seismic source; (e)x,ey) Relative deviation between the actual seismic source plane position and the projection direction of the nearest grid is obtained; sinc is a sine function notation; cos is the cosine function notation.
In an embodiment of the invention, when the irregular seismic source characterization device based on energy spectrum equivalence provided by the embodiment of the present invention is implemented specifically, the aforementioned regular grid point index range of seismic source equivalent approximation and the regular grid point weight function determination module of seismic source equivalent approximation are further configured to determine the regular grid point weight function of seismic source equivalent approximation according to the following manner:
Cxy(ix,iy)=Cx(ix)Cy(iy)
wherein, Cxy(ix,iy) A regular grid point weight function that is an equivalent approximation of the seismic source; cx(ix) And Cy(iy) Unidirectional grid weights for regular grid points that are equivalent approximations of the seismic source.
In an embodiment of the invention, when the irregular seismic source characterization device based on energy spectrum equivalence is implemented specifically, the module for determining a regular grid point weight function two-dimensional space fourier spectrum is specifically configured to determine a regular grid point weight function two-dimensional space fourier spectrum according to the following manner:
wherein, F (kx)l,kym) Two-dimensional spatial Fourier spectrum of a weight function of a regular grid point; cxy(ix,iy) A regular grid point weight function that is an equivalent approximation of the seismic source; (i)x,iy) Index ranges in x-and y-directions for regular grid points, i, which are equivalent approximations of the seismic sourcexAnd iyIs a positive integer; (i)x0,iy0) Projecting a grid point index for the nearest rule of the seismic source; (e)x,ey) Relative deviation between the actual seismic source plane position and the projection direction of the nearest grid is obtained; (kx, ky) is a discretization wave number array in the x-axis direction and the y-axis direction, the dimensions are p and q respectively, and subscripts of corresponding array elements are represented by l and m respectively; j is a unit imaginary number; exp is an exponential function.
In an embodiment of the invention, when the irregular seismic source characterization device based on energy spectrum equivalence provided in the embodiment of the present invention is implemented specifically, the energy error determination module is specifically configured to compare an energy error of a two-dimensional fourier spectrum of a unit impulse function in space corresponding to an actual seismic source position with an energy error of a two-dimensional fourier spectrum of a weight function of a regular grid point according to the following manner:
E(kxl,kym)=|F(kxl,kym)|2-1
wherein, F (kx)l,kym) Two-dimensional spatial Fourier spectrum of a weight function of a regular grid point; e (kx)l,kym) A two-dimensional fourier spectrum of the unit impulse function in space corresponding to the actual source location.
To sum up, the irregular seismic source characterization method and device based on energy spectrum equivalence provided by the embodiment of the invention firstly establish a regular gridding model and obtain the plane projection parameters of the seismic source; then, according to the plane projection parameters of the seismic source, determining the nearest regular projection grid point index of the seismic source and the relative deviation between the actual plane position of the seismic source and the nearest grid projection direction; then inputting a minimum approximate grid half-width parameter, and determining a regular grid point index range of equivalent approximation of the seismic source and a regular grid point weight function of equivalent approximation of the seismic source according to the index of the nearest regular projection grid point of the seismic source and the relative deviation of the actual plane position of the seismic source and the nearest grid projection direction; determining a two-dimensional space Fourier spectrum of the regular grid point weight function according to the discretization wave number array, the regular grid point index range of the seismic source equivalent approximation and the regular grid point weight function of the seismic source equivalent approximation; then comparing the energy error of the space unit pulse function two-dimensional Fourier spectrum corresponding to the actual seismic source position with the energy error of the regular grid point weight function two-dimensional space Fourier spectrum, and judging whether the energy error meets a preset energy error threshold value; and finally, when the energy error does not meet a preset energy error threshold, modifying the minimum approximate grid half-width parameter for iteration, and determining a grid node index and a corresponding node weight function value participating in the rule discretization representation until the energy error meets the preset energy error threshold. The embodiment of the invention establishes a regular discretization equivalent characterization method of an irregular excitation seismic source based on an energy spectrum equivalent principle, utilizes a method for characterizing the irregular seismic source by utilizing a regular discretization grid combined seismic source, finally generates a regular grid index and grid point weight which can be used for approximating the equivalent characterization, and can better reduce space positioning errors. For actual local sparse points, larger space positioning errors possibly caused are reduced, and adverse effects caused by subsequent imaging, illumination analysis and full waveform inversion are reduced. The method can calculate the weighting function of the seismic source combination in the corresponding rule discretization format under the condition of the given rule grid parameters and the plane coordinates of any seismic source, so that the error between the spectral energy of the corresponding equivalent seismic source combination and the actual seismic source spectral energy meets the given precision requirement.
The irregular seismic source characterization method and device based on energy spectrum equivalence have wide application prospect because the irregular seismic source characterization method and device based on energy spectrum equivalence is an important basic step in the current data processing in the seismic data processing, which is caused by the deep application of more and more nodes receiving, pile number-free excitation and compressed sensing. The feasibility of the method and the device is proved through theoretical data simulation verification.
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 is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the invention. It will be understood that each flow and/or block of the flow diagrams and/or block diagrams, and combinations of flows and/or blocks in the flow diagrams and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
The above-mentioned embodiments are intended to illustrate the objects, technical solutions and advantages of the present invention in further detail, and it should be understood that the above-mentioned embodiments are only exemplary embodiments of the present invention, and are not intended to limit the scope of the present invention, and any modifications, equivalent substitutions, improvements and the like made within the spirit and principle of the present invention should be included in the scope of the present invention.
Claims (26)
1. An irregular seismic source characterization method based on energy spectrum equivalence is characterized by comprising the following steps:
establishing a regular gridding model, and acquiring a seismic source plane projection parameter;
determining the nearest regular projection grid point index of the seismic source and the relative deviation of the actual seismic source plane position and the nearest grid projection direction according to the seismic source plane projection parameters;
inputting a minimum approximate grid half-width parameter, and determining a regular grid point index range of equivalent approximation of the seismic source and a regular grid point weight function of equivalent approximation of the seismic source according to the index of the nearest regular projection grid point of the seismic source and the relative deviation of the actual plane position of the seismic source and the nearest grid projection direction;
determining a two-dimensional space Fourier spectrum of the regular grid point weight function according to the discretization wave number array, the regular grid point index range of the seismic source equivalent approximation and the regular grid point weight function of the seismic source equivalent approximation;
comparing the energy error of the space unit pulse function two-dimensional Fourier spectrum corresponding to the actual seismic source position with the energy error of the regular grid point weight function two-dimensional space Fourier spectrum, and judging whether the energy error meets a preset energy error threshold value;
and when the energy error does not meet a preset energy error threshold, modifying the minimum approximate grid half-width parameter for iteration, and determining a grid node index and a corresponding node weight function value participating in the rule discretization representation until the energy error meets the preset energy error threshold.
2. The method of claim 1, wherein establishing a regular gridding model to obtain source plane projection parameters comprises:
establishing a regular gridding model according to the geological task, the seismic resolution requirement and the model space range; wherein the model space range comprises all actual seismic source positions to be calculated and all discretized equivalent seismic source positions;
and acquiring the plane projection parameters of the seismic source according to the regular gridding model.
3. The method of claim 2, wherein obtaining the source plane projection parameters according to a regular gridding model comprises:
determining a horizontal projection plane of the regular gridding model according to the regular gridding model;
and acquiring the reference origin coordinate, the plane grid distance, the grid number in the x and y directions, the projection sequence number and the seismic source plane projection coordinate of the horizontal projection plane according to the horizontal projection plane of the regular gridding model.
4. The method of claim 3, wherein determining, from the source plane projection parameters, a source nearest neighbor regular projection grid point index and a relative deviation of an actual source plane location from a nearest neighbor grid projection direction comprises:
determining the nearest regular projection grid point index of the seismic source according to the plane projection coordinate of the seismic source, the reference origin coordinate and the plane grid spacing;
and determining the relative deviation of the actual plane position of the seismic source and the projection direction of the nearest grid according to the index of the nearest regular projection grid point of the seismic source, the plane projection coordinate of the seismic source, the reference origin coordinate and the plane grid spacing.
5. The method of claim 4, wherein the source nearest-neighbor regular projection grid point index is determined as follows:
ix0=NINT[(sx-x0)/Δx)]+1
iy0=NINT[(sy-y0)/Δy)]+1
wherein (i)x0,iy0) Is the nearest neighbor of the seismic sourceRegularly projecting grid point indexes; (s)x,sy) Projecting coordinates for the seismic source plane; (x)0,y0) As reference origin coordinates; (Δ x, Δ y) is the planar grid spacing; NINT is rounded according to rounding.
6. The method of claim 4, wherein the relative deviation of the actual source plane location from the nearest grid projection direction is determined as follows:
wherein (e)x,ey) Relative deviation between the actual seismic source plane position and the projection direction of the nearest grid is obtained; (i)x0,iy0) Projecting a grid point index for the nearest rule of the seismic source; (s)x,sy) Projecting coordinates for the seismic source plane; (x)0,y0) As reference origin coordinates; (Δ x, Δ y) is the planar grid spacing.
7. The method of claim 1, wherein inputting a minimum approximation grid half-width parameter, determining a range of regular grid point indices for the source equivalent approximation and a regular grid point weight function for the source equivalent approximation based on the source nearest regular projection grid point index and a relative deviation of an actual source plane position from a nearest grid projection direction comprises:
inputting a minimum approximate grid half-width parameter, and determining a regular grid point index range of equivalent approximation of the seismic source according to the nearest regular projection grid point index of the seismic source;
determining the unidirectional grid weight of the regular grid points of the seismic source equivalent approximation according to the relative deviation of the actual seismic source plane position and the projection direction of the nearest grid, the index range of the regular grid points of the seismic source equivalent approximation, the index of the nearest regular projection grid points of the seismic source and the half-width parameter of the minimum approximation grid;
and combining the unidirectional grid weights of the regular grid points of the equivalent approximation of the seismic source to determine a weight function of the regular grid points of the equivalent approximation of the seismic source.
8. The method of claim 7, wherein the range of regular grid point indices for the equivalent approximation of the seismic source is determined as follows:
ix0-dw≤ix≤ix0+dw
iy0-dw≤iy≤iy0+dw
wherein (i)x,iy) Index ranges in x-and y-directions for regular grid points, i, which are equivalent approximations of the seismic sourcexAnd iyIs a positive integer; dwIs a minimum approximate grid half-width parameter; (i)x0,iy0) The grid point index is projected for the nearest neighbor rule of the seismic source.
9. The method of claim 7, wherein the unidirectional grid weights for the regular grid points of the equivalent approximation of the seismic source are determined as follows:
wherein, Cx(ix) And Cy(iy) A one-way grid weight of a regular grid point which is equivalent and approximate to the seismic source; (i)x,iy) Index ranges in x-and y-directions for regular grid points, i, which are equivalent approximations of the seismic sourcexAnd iyIs a positive integer; dwIs a minimum approximate grid half-width parameter; (i)x0,iy0) Is a vibrationProjecting grid point indexes by using a source nearest rule; (e)x,ey) Relative deviation between the actual seismic source plane position and the projection direction of the nearest grid is obtained; sinc is a sine function notation; cos is the cosine function notation.
10. The method of claim 7, wherein the regular grid point weight function for the equivalent approximation of the seismic source is determined as follows:
Cxy(ix,iy)=Cx(ix)Cy(iy)
wherein, Cxy(ix,iy) A regular grid point weight function that is an equivalent approximation of the seismic source; cx(ix) And Cy(iy) Unidirectional grid weights for regular grid points that are equivalent approximations of the seismic source.
11. The method of claim 1, characterized in that the regular grid point weight function two-dimensional spatial fourier spectrum is determined in the following manner:
wherein, F (kx)l,kym) Two-dimensional spatial Fourier spectrum of a weight function of a regular grid point; cxy(ix,iy) A regular grid point weight function that is an equivalent approximation of the seismic source; (i)x,iy) Index ranges in x-and y-directions for regular grid points, i, which are equivalent approximations of the seismic sourcexAnd iyIs a positive integer; (i)x0,iy0) Projecting a grid point index for the nearest rule of the seismic source; (e)x,ey) Relative deviation between the actual seismic source plane position and the projection direction of the nearest grid is obtained; (kx, ky) is a discretization wave number array in the x-axis direction and the y-axis direction, the dimensions are p and q respectively, and subscripts of corresponding array elements are represented by l and m respectively; j is a unit imaginary number; exp is an exponential function.
12. The method of claim 1, wherein the energy error of the two-dimensional space fourier spectrum of the unit impulse function in space corresponding to the actual source location is compared to the two-dimensional space fourier spectrum of the weight function of the regular grid points by:
E(kxl,kym)=|F(kxl,kym)|2-1
wherein, F (kx)l,kym) Two-dimensional spatial Fourier spectrum of a weight function of a regular grid point; e (kx)l,kym) A two-dimensional fourier spectrum of the unit impulse function in space corresponding to the actual source location.
13. An irregular seismic source characterization device based on energy spectrum equivalence, which is characterized by comprising:
the seismic source plane projection parameter acquisition module is used for establishing a regular gridding model and acquiring seismic source plane projection parameters;
the seismic source nearest regular projection grid point index and the relative deviation between the actual seismic source plane position and the nearest grid projection direction are determined by the seismic source nearest regular projection grid point index and the relative deviation between the actual seismic source plane position and the nearest grid projection direction are determined according to the seismic source plane projection parameters;
the seismic source equivalent approximate regular grid point index range and the seismic source equivalent approximate regular grid point weight function determining module is used for inputting a minimum approximate grid half-width parameter, and determining the seismic source equivalent approximate regular grid point index range and the seismic source equivalent approximate regular grid point weight function according to the relative deviation of the seismic source nearest regular projection grid point index and the actual seismic source plane position and the nearest grid projection direction;
the regular grid point weight function two-dimensional space Fourier spectrum determining module is used for determining a regular grid point weight function two-dimensional space Fourier spectrum according to the discretization wave number array, the regular grid point index range of the seismic source equivalent approximation and the regular grid point weight function of the seismic source equivalent approximation;
the energy error judgment module is used for comparing the energy error of the space unit pulse function two-dimensional Fourier spectrum corresponding to the actual seismic source position with the energy error of the regular grid point weight function two-dimensional space Fourier spectrum and judging whether the energy error meets a preset energy error threshold value;
and the grid node index and corresponding node weight function value determination module is used for correcting the minimum approximate grid half-width parameter for iteration when the energy error does not meet the preset energy error threshold value, and determining the grid node index and corresponding node weight function value participating in the rule discretization representation until the energy error meets the preset energy error threshold value.
14. The apparatus of claim 13, wherein the source plane projection parameters acquisition module is specifically configured to:
establishing a regular gridding model according to the geological task, the seismic resolution requirement and the model space range; wherein the model space range comprises all actual seismic source positions to be calculated and all discretized equivalent seismic source positions;
and acquiring the plane projection parameters of the seismic source according to the regular gridding model.
15. The apparatus of claim 14, wherein the source plane projection parameters acquisition module is further to:
determining a horizontal projection plane of the regular gridding model according to the regular gridding model;
and acquiring the reference origin coordinate, the plane grid distance, the grid number in the x and y directions, the projection sequence number and the seismic source plane projection coordinate of the horizontal projection plane according to the horizontal projection plane of the regular gridding model.
16. The apparatus of claim 15, wherein the source nearest regular projected grid point index and the relative deviation of the actual source plane position from the nearest grid projection direction determination module are specifically configured to:
determining the nearest regular projection grid point index of the seismic source according to the plane projection coordinate of the seismic source, the reference origin coordinate and the plane grid spacing;
and determining the relative deviation of the actual plane position of the seismic source and the projection direction of the nearest grid according to the index of the nearest regular projection grid point of the seismic source, the plane projection coordinate of the seismic source, the reference origin coordinate and the plane grid spacing.
17. The apparatus of claim 16, wherein the source nearest regular projected grid point index and the relative deviation of the actual source plane location from the nearest grid projection direction determination module are further configured to determine the source nearest regular projected grid point index as follows:
ix0=NINT[(sx-x0)/Δx)]+1
iy0=NINT[(sy-y0)/Δy)]+1
wherein (i)x0,iy0) Projecting a grid point index for the nearest rule of the seismic source; (s)x,sy) Projecting coordinates for the seismic source plane; (x)0,y0) As reference origin coordinates; (Δ x, Δ y) is the planar grid spacing; NINT is rounded according to rounding.
18. The apparatus of claim 16, wherein the source nearest regular projected grid point index and actual source plane position relative deviation from nearest grid projection direction determination module is further configured to determine the actual source plane position relative deviation from nearest grid projection direction as follows:
wherein (e)x,ey) Relative deviation between the actual seismic source plane position and the projection direction of the nearest grid is obtained; (i)x0,iy0) Projecting a grid point index for the nearest rule of the seismic source; (s)x,sy) Projecting coordinates for the seismic source plane; (x)0,y0) As reference origin coordinates; (Δ x, Δ y) is the planar grid spacing.
19. The apparatus of claim 13, wherein the index range of regular grid points for the source equivalent approximation and the weight function determination module for the regular grid points for the source equivalent approximation are specifically configured to:
inputting a minimum approximate grid half-width parameter, and determining a regular grid point index range of equivalent approximation of the seismic source according to the nearest regular projection grid point index of the seismic source;
determining the unidirectional grid weight of the regular grid points of the seismic source equivalent approximation according to the relative deviation of the actual seismic source plane position and the projection direction of the nearest grid, the index range of the regular grid points of the seismic source equivalent approximation, the index of the nearest regular projection grid points of the seismic source and the half-width parameter of the minimum approximation grid;
and combining the unidirectional grid weights of the regular grid points of the equivalent approximation of the seismic source to determine a weight function of the regular grid points of the equivalent approximation of the seismic source.
20. The apparatus of claim 19, wherein the regular grid point index range for the source equivalent approximation and the regular grid point weight function determination module for the source equivalent approximation are further configured to determine the regular grid point index range for the source equivalent approximation as follows:
ix0-dw≤ix≤ix0+dw
iy0-dw≤iy≤iy0+dw
wherein (i)x,iy) Index ranges in x-and y-directions for regular grid points, i, which are equivalent approximations of the seismic sourcexAnd iyIs a positive integer; dwIs a minimum approximate grid half-width parameter; (i)x0,iy0) The grid point index is projected for the nearest neighbor rule of the seismic source.
21. The apparatus of claim 19, wherein the index range of regular grid points for the source equivalent approximation and the function determination module for weights of regular grid points for the source equivalent approximation are further configured to determine unidirectional grid weights for regular grid points for the source equivalent approximation as follows:
wherein, Cx(ix) And Cy(iy) A one-way grid weight of a regular grid point which is equivalent and approximate to the seismic source; (i)x,iy) Index ranges in x-and y-directions for regular grid points, i, which are equivalent approximations of the seismic sourcexAnd iyIs a positive integer; dwIs a minimum approximate grid half-width parameter; (i)x0,iy0) Projecting a grid point index for the nearest rule of the seismic source; (e)x,ey) Relative deviation between the actual seismic source plane position and the projection direction of the nearest grid is obtained; sinc is a sine function notation; cos is the cosine function notation.
22. The apparatus of claim 19, wherein the index range of regular grid points for the source equivalent approximation and the weight function determination module for regular grid points for the source equivalent approximation are further configured to determine the weight function for regular grid points for the source equivalent approximation as follows:
Cxy(ix,iy)=Cx(ix)Cy(iy)
wherein, Cxy(ix,iy) A regular grid point weight function that is an equivalent approximation of the seismic source; cx(ix) And Cy(iy) Regular grid point list for equivalent approximation of seismic sourceAnd weighting to the grid.
23. The apparatus according to claim 13, wherein the regular grid point weight function two-dimensional spatial fourier spectrum determining module is specifically configured to determine the regular grid point weight function two-dimensional spatial fourier spectrum in the manner:
wherein, F (kx)l,kym) Two-dimensional spatial Fourier spectrum of a weight function of a regular grid point; cxy(ix,iy) A regular grid point weight function that is an equivalent approximation of the seismic source; (i)x,iy) Index ranges in x-and y-directions for regular grid points, i, which are equivalent approximations of the seismic sourcexAnd iyIs a positive integer; (i)x0,iy0) Projecting a grid point index for the nearest rule of the seismic source; (e)x,ey) Relative deviation between the actual seismic source plane position and the projection direction of the nearest grid is obtained; (kx, ky) is a discretization wave number array in the x-axis direction and the y-axis direction, the dimensions are p and q respectively, and subscripts of corresponding array elements are represented by l and m respectively; j is a unit imaginary number; exp is an exponential function.
24. The apparatus according to claim 13, wherein the energy error determination module is specifically configured to compare the energy error of the two-dimensional space fourier spectrum of the unit impulse function in space corresponding to the actual source location with the two-dimensional space fourier spectrum of the weight function of the regular grid points, as follows:
E(kxl,kym)=|F(kxl,kym)|2-1
wherein, F (kx)l,kym) Two-dimensional spatial Fourier spectrum of a weight function of a regular grid point; e (kx)l,kym) A two-dimensional fourier spectrum of the unit impulse function in space corresponding to the actual source location.
25. A computer device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, characterized in that the processor implements the method of any of claims 1 to 12 when executing the computer program.
26. A computer-readable storage medium, characterized in that the computer-readable storage medium stores a computer program for executing a method according to any one of claims 1 to 12.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011173079.XA CN112379413B (en) | 2020-10-28 | Irregular seismic source characterization method and device based on energy spectrum equivalence |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011173079.XA CN112379413B (en) | 2020-10-28 | Irregular seismic source characterization method and device based on energy spectrum equivalence |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112379413A true CN112379413A (en) | 2021-02-19 |
CN112379413B CN112379413B (en) | 2024-07-26 |
Family
ID=
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114114395A (en) * | 2021-11-09 | 2022-03-01 | 同济大学 | Irregular grid aliasing seismic data separation method and system based on interpolation operator |
CN114545491A (en) * | 2021-12-30 | 2022-05-27 | 中国石油天然气股份有限公司 | Optimization method and device for arrangement parameters in seismic data acquisition |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20060098529A1 (en) * | 2004-11-08 | 2006-05-11 | Exxonmobil Upstream Research Company | Method for data regulariization for shot domain processing |
CN1797038A (en) * | 2004-12-29 | 2006-07-05 | 中国石油天然气集团公司 | Method for shifting depth before superposition in seismic data process of undulating the earth's surface |
US20090290449A1 (en) * | 2008-05-25 | 2009-11-26 | Smith Patrick J | Processing Seismic Data Using Combined Regularization and 4D Binning |
CN104977607A (en) * | 2014-04-09 | 2015-10-14 | 中国石油集团东方地球物理勘探有限责任公司 | Time-domain full waveform inversion method by utilizing variable-step-length grid acoustic wave field simulation |
CN104977608A (en) * | 2014-04-09 | 2015-10-14 | 中国石油集团东方地球物理勘探有限责任公司 | Time-domain full waveform inversion method by utilizing fixed grid acoustic wave field simulation |
CN105093265A (en) * | 2014-05-09 | 2015-11-25 | 中国石油化工股份有限公司 | Method of simulating the transmission rules of seismic waves in a TI medium |
CN105549068A (en) * | 2015-12-09 | 2016-05-04 | 中国科学院地质与地球物理研究所 | 3D anisotropic micro seismic interference inverse positioning method and 3D anisotropic micro seismic interference inverse positioning system |
KR101635791B1 (en) * | 2015-02-27 | 2016-07-04 | 전남대학교산학협력단 | Determination method for location and origin time of earthquake using arrival time of primary wave |
CN108957538A (en) * | 2018-06-21 | 2018-12-07 | 成都启泰智联信息科技有限公司 | A kind of virtual focus two dimension wavefront construction seimic travel time calculation method |
CN111694051A (en) * | 2020-06-17 | 2020-09-22 | 电子科技大学 | Gaussian beam-based viscoacoustic medium seismic wave forward modeling method |
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20060098529A1 (en) * | 2004-11-08 | 2006-05-11 | Exxonmobil Upstream Research Company | Method for data regulariization for shot domain processing |
CN1797038A (en) * | 2004-12-29 | 2006-07-05 | 中国石油天然气集团公司 | Method for shifting depth before superposition in seismic data process of undulating the earth's surface |
US20090290449A1 (en) * | 2008-05-25 | 2009-11-26 | Smith Patrick J | Processing Seismic Data Using Combined Regularization and 4D Binning |
CN104977607A (en) * | 2014-04-09 | 2015-10-14 | 中国石油集团东方地球物理勘探有限责任公司 | Time-domain full waveform inversion method by utilizing variable-step-length grid acoustic wave field simulation |
CN104977608A (en) * | 2014-04-09 | 2015-10-14 | 中国石油集团东方地球物理勘探有限责任公司 | Time-domain full waveform inversion method by utilizing fixed grid acoustic wave field simulation |
CN105093265A (en) * | 2014-05-09 | 2015-11-25 | 中国石油化工股份有限公司 | Method of simulating the transmission rules of seismic waves in a TI medium |
KR101635791B1 (en) * | 2015-02-27 | 2016-07-04 | 전남대학교산학협력단 | Determination method for location and origin time of earthquake using arrival time of primary wave |
CN105549068A (en) * | 2015-12-09 | 2016-05-04 | 中国科学院地质与地球物理研究所 | 3D anisotropic micro seismic interference inverse positioning method and 3D anisotropic micro seismic interference inverse positioning system |
CN108957538A (en) * | 2018-06-21 | 2018-12-07 | 成都启泰智联信息科技有限公司 | A kind of virtual focus two dimension wavefront construction seimic travel time calculation method |
CN111694051A (en) * | 2020-06-17 | 2020-09-22 | 电子科技大学 | Gaussian beam-based viscoacoustic medium seismic wave forward modeling method |
Non-Patent Citations (6)
Title |
---|
周松 等: "《基于压缩感知的非规则地震勘探观测***设计与数据重建》", 《石油物探》, vol. 56, no. 5, pages 617 - 625 * |
孙婧 等: "《基于双程波动方程的多次波照明分析方法研究》", 《地球物理学进展》, vol. 33, no. 1, pages 243 - 249 * |
张钋 等: "《一种基于Kirchhoff积分的快速三维非均匀弹性介质地震波场模拟方法》", 《石油物探》, vol. 46, no. 6, pages 621 - 627 * |
景月岭 等: "《多次各向异性散射模式对S波能量密度包络曲线的影响》", 《地球物理学报》, vol. 55, no. 6, pages 1993 - 2003 * |
王庆 等: "《时间域地震全波形反演方法进展》", 《地球物理学进展》, vol. 30, no. 6, 31 December 2015 (2015-12-31), pages 2797 - 2806 * |
陈浩 等: "<标量波Laplace-Fourier域有限差分正演模拟方法>", 《地球物理学进展》, vol. 32, no. 2, pages 584 - 560 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114114395A (en) * | 2021-11-09 | 2022-03-01 | 同济大学 | Irregular grid aliasing seismic data separation method and system based on interpolation operator |
CN114114395B (en) * | 2021-11-09 | 2022-10-25 | 同济大学 | Irregular grid aliasing seismic data separation method and system based on interpolation operator |
CN114545491A (en) * | 2021-12-30 | 2022-05-27 | 中国石油天然气股份有限公司 | Optimization method and device for arrangement parameters in seismic data acquisition |
CN114545491B (en) * | 2021-12-30 | 2024-06-07 | 中国石油天然气股份有限公司 | Optimization method and device for arrangement parameters in seismic data acquisition |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110031896B (en) | Seismic random inversion method and device based on multi-point geostatistics prior information | |
Sraj et al. | Uncertainty quantification and inference of Manning’s friction coefficients using DART buoy data during the Tōhoku tsunami | |
EP1877834A1 (en) | 3d fast fault restoration | |
Gao et al. | An Efficient Multiscale Finite‐Element Method for Frequency‐Domain Seismic Wave Propagation | |
Branchini et al. | Reconstructing Positions and Peculiar Velocities of Galaxy Clusters within 20000 km/sec. I: The Cluster 3D Dipole | |
De Curtis et al. | Collision integrals for cosmological phase transitions | |
Qu et al. | 3-D Least-Squares Reverse Time Migration in Curvilinear-τ Domain | |
Siahkoohi et al. | Velocity continuation with Fourier neural operators for accelerated uncertainty quantification | |
CN104081437A (en) | Bayesian approach for gas concentration reconstruction based on tunable diode laser absorption spectroscopy | |
CN112379413B (en) | Irregular seismic source characterization method and device based on energy spectrum equivalence | |
CN112379413A (en) | Irregular seismic source characterization method and device based on energy spectrum equivalence | |
Gu et al. | Analysis and quantification of data assimilation based on sequential Monte Carlo methods for wildfire spread simulation | |
He et al. | Accelerating full waveform inversion using HSS solver and limited memory conjugate gradient method | |
Mao | Towards front-tracking based on conservation in two space dimensions II, tracking discontinuities in capturing fashion | |
Yin et al. | A fast 3D gravity forward algorithm based on circular convolution | |
Picasso | Numerical study of an anisotropic error estimator in the L^2(H^1) norm for the finite element discretization of the wave equation | |
Baker | Anisotropic Adaptivity and Subgrid Scale Modelling for the Solution of the Neutron Transport Equation with an Emphasis on Shielding Applications | |
Coulibaly et al. | Monte Carlo and quasi-Monte Carlo algorithms for a linear integro-differential equation | |
Hiltebrand et al. | Efficient preconditioners for a shock capturing space-time discontinuous Galerkin method for systems of conservation laws | |
Zhou et al. | New high-fidelity terrain modeling method constrained by terrain semanteme | |
Bashkardin et al. | Phase-space computation of multi-arrival traveltimes, Part I: Theory and concepts | |
Li et al. | A low order non linear transport acceleration scheme for the method of the characteristic | |
Zhou et al. | A topography-dependent eikonal solver for accurate and efficient computation of traveltimes and their derivatives in 3D heterogeneous media | |
Yan et al. | Modeling seismic wave propagation in heterogeneous medium using overlap domain pseudospectral method | |
Gao et al. | LATTE: Open-source, high-performance acoustic and elastic traveltime computation, tomography, and source location |
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 |