CN109031410B - The bilateral beam synthetic method in big gun inspection domain and system of common offset field result constraint - Google Patents
The bilateral beam synthetic method in big gun inspection domain and system of common offset field result constraint Download PDFInfo
- Publication number
- CN109031410B CN109031410B CN201810776052.6A CN201810776052A CN109031410B CN 109031410 B CN109031410 B CN 109031410B CN 201810776052 A CN201810776052 A CN 201810776052A CN 109031410 B CN109031410 B CN 109031410B
- Authority
- CN
- China
- Prior art keywords
- data
- local plane
- shot
- plane wave
- seismic
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 230000002146 bilateral effect Effects 0.000 title claims abstract description 59
- 238000007689 inspection Methods 0.000 title abstract description 6
- 238000010189 synthetic method Methods 0.000 title abstract description 5
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 71
- 238000003786 synthesis reaction Methods 0.000 claims abstract description 71
- 238000000034 method Methods 0.000 claims abstract description 33
- 238000001514 detection method Methods 0.000 claims abstract description 24
- 230000002194 synthesizing effect Effects 0.000 claims description 22
- 238000013507 mapping Methods 0.000 claims description 20
- 238000001308 synthesis method Methods 0.000 claims description 19
- 238000005070 sampling Methods 0.000 claims description 15
- 238000012937 correction Methods 0.000 claims description 14
- 239000000523 sample Substances 0.000 claims description 12
- 239000011159 matrix material Substances 0.000 claims description 8
- 230000009466 transformation Effects 0.000 claims description 7
- 238000012163 sequencing technique Methods 0.000 claims description 5
- 238000013075 data extraction Methods 0.000 claims description 4
- 238000003384 imaging method Methods 0.000 abstract description 11
- 238000007906 compression Methods 0.000 abstract description 4
- 230000006835 compression Effects 0.000 abstract description 4
- 239000000284 extract Substances 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 11
- 238000013508 migration Methods 0.000 description 9
- 230000005012 migration Effects 0.000 description 9
- 150000003839 salts Chemical class 0.000 description 4
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000013144 data compression Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 230000000452 restraining effect Effects 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- HOWHQWFXSLOJEF-MGZLOUMQSA-N systemin Chemical compound NCCCC[C@H](N)C(=O)N[C@@H](CCSC)C(=O)N[C@@H](CCC(N)=O)C(=O)N[C@@H]([C@@H](C)O)C(=O)N[C@@H](CC(O)=O)C(=O)OC(=O)[C@@H]1CCCN1C(=O)[C@H]1N(C(=O)[C@H](CC(O)=O)NC(=O)[C@H](CCCN=C(N)N)NC(=O)[C@H](CCCCN)NC(=O)[C@H](CO)NC(=O)[C@H]2N(CCC2)C(=O)[C@H]2N(CCC2)C(=O)[C@H](CCCCN)NC(=O)[C@H](CO)NC(=O)[C@H](CCC(N)=O)NC(=O)[C@@H](NC(=O)[C@H](C)N)C(C)C)CCC1 HOWHQWFXSLOJEF-MGZLOUMQSA-N 0.000 description 1
- 108010050014 systemin Proteins 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
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 discloses the bilateral beam synthetic methods in big gun inspection domain and system of a kind of constraint of common offset field result, and this method includes that the big gun detection point range that will be counted in advance is divided into multiple rectangle face elements and label；The face element data directory structure comprising face element number is established according to the earthquake record data of acquisition；The unilateral Local plane wave of common offset of multiple single side members is synthesized according to the Local plane wave synthetic parameters of face element data directory structure and input and extracts effective lineups data of this plane wave, this effective lineups data is mapped to big gun inspection domain from common offset, the bilateral Local plane wave in domain is then examined according to the big gun that the bilateral Shu Hecheng inversion problem formula in big gun inspection domain under the constraint of common offset field result and the unilateral Local plane wave of common offset after map operation synthesize multiple single side members.Earthquake beam energy out-focus problem in bilateral synthesis beam is able to solve using method provided by the invention or system, improves the compression ratio of earthquake beam data, imaging efficiency and precision.
Description
Technical Field
The invention relates to the technical field of seismic data processing of petroleum geophysical exploration, in particular to a shot-check domain bilateral beam synthesis method and system based on common offset domain result constraint.
Background
The Gaussian beam migration method is an important method in the field of seismic exploration data processing migration imaging, is between a wave equation method and a high-frequency approximate ray method, and has the precision close to that of the wave equation method and the efficiency close to that of the ray migration method. Beam shifting is a class of methods that encompasses and is simplified over gaussian beams, which in the narrow sense refers to simplified gaussian beam shifting imaging methods, and has been extensively studied and developed in academia and industry. Beam migration in the narrow sense should be more efficient than the kirchhoff method and the traditional Gaussian beam migration, wherein one of the methods for improving the efficiency is to determine the ray parameters of the shot point incidence and the ray parameters of the wave detection point emergence through the plane spectrum, so as to avoid summing ray bundles in all directions, and therefore bilateral synthesis of the seismic beam is required to be performed in the shot detection domain.
For bilateral beam synthesis of actual field data, firstly, the data volume of five-dimensional local plane waves output by a single surface element is huge, and secondly, in many cases, shot points are sparse, so that local plane wave synthesis (seismic beam synthesis) of a common-detection-point gather cannot be focused, the ray direction of the shot points cannot be determined, and the compression and rapid imaging of the local plane wave data are not facilitated.
Disclosure of Invention
The invention aims to solve the problem that the energy of a seismic beam in a bilateral synthetic beam is not focused, and provides a shot-geophone domain bilateral beam synthetic method with result constraint in a common offset domain.
In order to achieve the purpose, the invention provides the following scheme:
a shot domain bilateral beam synthesis method based on common offset domain result constraint comprises the following steps:
acquiring seismic record data and local plane wave synthesis parameters; the seismic record data is seismic exploration field acquisition data or laboratory synthesis data; the local plane wave synthesis parameters comprise the size of a surface element, the angle interval of the direction of the local plane wave, the number of the direction angles of the local plane wave, the number of effective homophase axis data pick-up during output and a seismic velocity model; the standard format of the seismic record data is SEGY format, and the SEGY format generally comprises three parts, namely an EBCDIC file header, a binary file header and a seismic channel; each seismic channel comprises channel head information and seismic channel data; the line number, the track number, the sampling point number, the geodetic coordinate, the shot point coordinate and the demodulator probe coordinate corresponding to the seismic channel are stored in the channel head information;
uniformly dividing a shot point range and a receiver point range which are counted in advance by using the bin size parameters to obtain a plurality of rectangular bins, sequencing the rectangular bins in sequence, recording bin numbers of the rectangular bins, and determining a central point shot point coordinate and a central point receiver point coordinate which correspond to each rectangular bin;
establishing a surface element data index structure containing surface element numbers according to all seismic channel data in the seismic record data; the header file of the bin data index structure comprises shot offset distances of rectangular bin center coordinates, and the bin data index structure further comprises offset distances of each seismic channel data; the offset of the central coordinate of the rectangular bin is the difference between the coordinate of the central point geophone point and the coordinate of the central point shot point, and the offset of the seismic channel data is the difference between the coordinate of the geophone point and the coordinate of the shot point stored in the channel head information corresponding to the seismic channel data;
synthesizing the seismic channel data of each surface element according to the surface element data index structure, the local plane wave direction angle interval and the local plane wave direction angle number to obtain a common offset range unilateral local plane wave of a plurality of single surface elements;
extracting effective in-phase axis data of the common offset range unilateral local plane wave of each surface element according to the number of the effective in-phase axis data during output; the effective in-phase axis data is in-phase axis data of effective local plane waves;
and mapping the effective homophase axis data of the common offset range single-sided local plane waves from the common offset range to a shot-geophone range, and synthesizing the shot-geophone range double-sided local plane waves of a plurality of single-surface elements according to a shot-geophone range double-sided beam synthesis inverse problem formula under the common offset range result constraint and the common offset range single-sided local plane waves after mapping operation.
Optionally, the shot domain bilateral beam synthesis method further includes: and extracting effective homophase axis data of bilateral local plane waves of the shot domain of each surface element and compressing the effective homophase axis data.
Optionally, the synthesizing the seismic channel data of each surface element according to the surface element data index structure, the local plane wave direction angle interval, and the local plane wave direction angle number to obtain a common offset range single-sided local plane wave of a plurality of single surface elements specifically includes:
searching all seismic channel data corresponding to each bin number from the bin data index structure;
correcting time difference caused by the fact that the offset distance of all seismic channel data corresponding to each bin number is different from the offset distance of shot-geophone of the rectangular bin center coordinate corresponding to the bin number to obtain corrected seismic channel data corresponding to each bin number;
and synthesizing the common offset domain unilateral local plane waves of the single surface elements according to the corrected seismic channel data corresponding to each surface element number, the local plane wave direction angle interval and the local plane wave direction angle number.
Optionally, the formula for correcting the seismic trace data isWherein Δ t is the travel time difference between the seismic traces within the bin and the seismic trace at the center point of the bin; Δ h is the difference in offset between the seismic traces in the bin and the seismic trace at the center point of the bin; h ismIs the offset, t, of the seismic trace at the center point of the bin0The travel time of the seismic channel with zero offset distance of the surface element center point is shown; v. ofnIs the motion correction speed.
Optionally, a synthesis formula of the common offset domain single-sided local plane wave is as follows:
wherein,is the coordinate of the center point and is, the data are single-sided local plane wave data of a common offset range of a single surface element;is the direction of the plane wave at the center point,the method is determined by the angle interval of the direction of the local plane wave, the number of the angle of the direction of the local plane wave, the number of the picked effective homophase axis data during output and a seismic velocity model; ω is the angular frequency; a transformation matrix, AHIs the conjugate transpose of A; d is the seismic trace data contained in the bin,is the seismic trace coordinates within the bin of the geophone point,and the coordinates of the seismic channels in the shot surface element.
Optionally, the extracting, according to the number of the effective event axis data pickups during the outputting, effective event axis data of the common offset range single-sided local plane wave of each bin specifically includes:
and extracting and compressing the effective in-phase axis data of the common offset range unilateral local plane wave corresponding to each surface element by adopting a threshold value method according to the number of the effective in-phase axis data during output.
Optionally, the mapping the effective in-phase axis data of the common offset domain single-sided local plane wave from the common offset domain to the shot domain specifically includes: adopting the following formula to make the common offset range single-side local planeMapping the effective in-phase axis data of the wave from the common offset domain to a shot domain; the formula isWherein,respectively the plane wave directions of the central point, the shot point and the wave detection point,the method is determined by the local plane wave direction angle interval, the local plane wave direction angle number, the effective homophase axis data pick-up number during output and a seismic velocity model.
Optionally, the shot-geophone domain bilateral beam synthesis inverse problem formula under the common offset distance domain result constraint is as follows:
the shot-geophone domain bilateral beam synthesis inverse problem formula under the constraint of the common offset domain result is as follows:
wherein E is an error function or an objective function; ω is the angular frequency;andrespectively representing the seismic channel coordinates in a shot point bin and the seismic channel coordinates in a detector point bin; a, converting a matrix, wherein lambda is a hyper-parameter coefficient and is a constant; d is all seismic trace data contained in the singleton;the local plane wave data of the shot domain and the two sides of the shot domain of the single surface element are obtained;is a common offset field of a single binSingle-sided local plane wave data;respectively the plane wave directions of the central point, the shot point and the wave detection point,the method is determined by the local plane wave direction angle interval, the local plane wave direction angle number, the effective homophase axis data pick-up number during output and a seismic velocity model.
Optionally, a synthesis formula of the bilateral local plane wave in the shot and geophone domain is as follows:
the invention also provides a shot-geophone domain bilateral beam synthesis system with result constraint of common offset range, which comprises:
the seismic record data and local plane wave synthesis parameter acquisition module is used for acquiring seismic record data and local plane wave synthesis parameters; the seismic record data is seismic exploration field acquisition data or laboratory synthesis data; the local plane wave synthesis parameters comprise the size of a surface element, the angle interval of the direction of the local plane wave, the number of the direction angles of the local plane wave, the number of effective homophase axis data pick-up during output and a seismic velocity model; the standard format of the seismic record data is SEGY format, and the SEGY format generally comprises three parts, namely an EBCDIC file header, a binary file header and a seismic channel; each seismic channel comprises channel head information and seismic channel data; the line number, the track number, the sampling point number, the geodetic coordinate, the shot point coordinate and the demodulator probe coordinate corresponding to the seismic channel are stored in the channel head information;
the rectangular surface element and surface element number obtaining module is used for uniformly dividing a shot point range and a receiver point range which are counted in advance by adopting surface element size parameters to obtain a plurality of rectangular surface elements, sequencing the rectangular surface elements in sequence, recording surface element numbers of the rectangular surface elements, and determining a shot point coordinate of a central point and a receiver point coordinate of the central point corresponding to each rectangular surface element;
the bin data index structure establishing module is used for establishing a bin data index structure containing bin numbers according to all seismic channel data in the seismic record data; the header file of the bin data index structure comprises shot offset distances of rectangular bin center coordinates, and the bin data index structure further comprises offset distances of each seismic channel data; the offset of the central coordinate of the rectangular bin is the difference between the coordinate of the central point geophone point and the coordinate of the central point shot point, and the offset of the seismic channel data is the difference between the coordinate of the geophone point and the coordinate of the shot point stored in the channel head information corresponding to the seismic channel data;
the common offset range single-sided local plane wave synthesis module is used for synthesizing the seismic channel data of each surface element according to the surface element data index structure, the local plane wave direction angle interval and the local plane wave direction angle number to obtain common offset range single-sided local plane waves of a plurality of single surface elements;
the first effective in-phase axis data extraction module is used for extracting effective in-phase axis data of the common offset range unilateral local plane wave of each surface element according to the number of the effective in-phase axis data during output; the effective in-phase axis data is in-phase axis data of effective local plane waves;
and the shot domain bilateral local plane wave synthesis module is used for mapping the effective homophase axis data of the common offset domain unilateral local plane wave to a shot domain from the common offset domain, synthesizing an inversion problem formula according to a shot domain bilateral beam under the common offset domain result constraint and the common offset domain bilateral local plane wave after mapping operation, and synthesizing a plurality of single-surface element shot domain bilateral local plane waves.
According to the specific embodiment provided by the invention, the invention discloses the following technical effects:
the invention provides a shot and examine field bilateral beam synthetic method and system of common offset range result constraint, the method includes obtaining earthquake record data, surface element size, local plane wave direction angle interval, local plane wave direction angle number, effective homophase axis data pick-up number and earthquake velocity model when outputting, and adopt surface element size parameter to divide the shot point range and examine wave point range counted in advance into a plurality of rectangular surface elements, and sort the rectangular surface elements in turn, record the surface element number of the rectangular surface element; establishing a surface element data index structure containing surface element numbers according to all seismic channel data in the seismic record data; synthesizing seismic channel data of each surface element according to a surface element data index structure, local plane wave direction angle intervals and the number of local plane wave direction angles to obtain common offset range unilateral local plane waves of a plurality of single surface elements; extracting effective in-phase axis data of the common offset range unilateral local plane wave of each surface element according to the number of the effective in-phase axis data during output; and mapping effective homophase axis data of the common offset range unilateral local plane waves from the common offset range to a shot-geophone field, and synthesizing an inverse problem formula and the common offset range unilateral local plane waves after mapping operation according to the shot-geophone field bilateral beam under the common offset range result constraint and synthesizing the shot-geophone field bilateral local plane waves of a plurality of single-surface elements. By applying the method or the system provided by the invention, the shot-geophone domain bilateral seismic beam synthesis is restricted through the common offset distance domain unilateral seismic beam synthesis result so as to solve the problem that the seismic beam energy in the bilateral synthetic beam is not focused, and the compression rate, the imaging efficiency and the precision of seismic beam data are improved.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings needed in the embodiments will be briefly described below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and it is obvious for those skilled in the art to obtain other drawings without creative efforts.
FIG. 1 is a schematic flowchart of a shot-geophone domain bilateral beam synthesis method with common offset domain result constraint according to an embodiment of the present invention;
FIG. 2 is a schematic structural diagram of a shot-geophone domain bilateral beam synthesis system with common offset domain result constraint according to an embodiment of the present invention;
FIG. 3 is a velocity model diagram of a three-dimensional salt dome model crossline 246 in accordance with the present invention;
FIG. 4 is a velocity model diagram of a three-dimensional salt dome model inline 376 of the present invention;
FIG. 5 is a diagram of seismic synthetic recording data according to the present invention;
FIG. 6 is a schematic diagram of a local plane wave without partial motion correction moveout correction for the 300 th surface element of the common offset domain according to the present invention;
FIG. 7 is a schematic diagram of a local plane wave obtained by partially correcting the moveout difference of the 300 th surface element in the common offset range according to the present invention;
FIG. 8 is a schematic diagram of an effective plane wave obtained by a threshold method through 300 th bin common offset domain single-sided local plane wave transformation according to the present invention;
FIG. 9 is a result graph showing three-dimensional local plane waves of the edge of a wave detection point when the azimuth angle and the inclination angle of the shot edge obtained by mapping the common offset local plane waves to the shot detection domain are both zero according to the present invention;
FIG. 10 is a result graph showing three-dimensional local plane waves of shot boundaries when azimuth angles and inclination angles of the common offset local plane waves are mapped to shot detection domains and the receiver point edges are zero according to the present invention;
FIG. 11 is a three-dimensional local plane wave result diagram of a displayed wave detection point edge when the azimuth angle and the inclination angle of the shot-detection domain local plane wave are both zero under the constraint of the common offset result of the invention;
FIG. 12 is a three-dimensional local plane wave result diagram of a displayed wave detection point edge when the azimuth angle and the inclination angle of the shot-detection domain local plane wave are both zero under the constraint of the common offset result of the invention;
FIG. 13 is a schematic representation of the seismic imaging results of the present invention obtained by beam migration for Inline number 246.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
The current industry selects a single-side beam synthesis method to perform beam migration imaging, so that the problem that double-side synthesized seismic beam data is huge and low in efficiency is avoided, and the problem that shot point sparsity is not beneficial to beam synthesis at the shot edge is avoided. If the problem that the synthetic energy of the seismic beam is not focused due to sparse cannons or insufficient data is solved, the data compression rate, the imaging efficiency and the precision of beam migration can be improved, and the application of a beam migration method of bilateral beam synthesis in the industry can be promoted.
Therefore, the invention aims to solve the problem that the energy of a seismic beam in a bilateral synthetic beam is not focused, and provides a shot-geophone domain bilateral beam synthetic method and system with common offset domain result constraint. According to the method, the common offset range domain single-side seismic beam synthesis result is used for restraining shot-geophone domain double-side seismic beam synthesis, and the problem of energy unfocusing is finally solved through iteration of the common offset range domain beam synthesis and the shot-geophone domain double-side synthesis, so that synthetic data in a shot-geophone domain high-dimensional space is sparse, and the compression rate, imaging efficiency and accuracy of seismic beam data are improved.
In order to make the aforementioned objects, features and advantages of the present invention more comprehensible, the present invention is described in detail with reference to the accompanying drawings and the detailed description thereof.
The beam synthesis algorithm is an optimization problem, and an effective way for solving the problem of non-focusing of local plane wave transformation caused by shot point sparsity (sparse gather of the wave detection points) from the constraint angle is provided.
The shot domain and the central point offset domain are coordinate spaces which can be mutually converted, although the coordinate spaces are in a linear transformation relation, the missing performance of actual data in the two spaces is different, and therefore the precision of the double-side beam synthesis method of the shot domain can be improved by using the advantages of each space through a cross-coordinate space constraint method.
Fig. 1 is a schematic flowchart of a shot-geophone domain bilateral beam synthesis method based on result constraint of a common offset range domain according to an embodiment of the present invention, and as shown in fig. 1, the shot-geophone domain bilateral beam synthesis method according to the embodiment of the present invention includes the following steps.
Step 101: acquiring seismic record data and local plane wave synthesis parameters; the seismic record data is seismic exploration field acquisition data or laboratory synthesis data; the local plane wave synthesis parameters mainly comprise the size of a surface element, the angle interval of the direction of the local plane wave, the number of the direction angles of the local plane wave, the number of the effective homophase axis data during output and a seismic velocity model.
The standard format of the seismic record data is SEGY format, the SEGY format generally comprises three parts, the first part is an EBCDIC file header (3200 bytes) and consists of 40 cards (for example, 80 characters per line by 40 lines) for storing some information describing a seismic data body; the second part is a binary file header (400 bytes) used for storing some key information describing the SEGY file, including some information such as the data format, the sampling point number, the sampling interval, the measurement unit and the like of the SEGY file, and the information is generally stored at a fixed position of the binary file header; the third part is the actual seismic traces, each containing 240 bytes of trace header information and seismic trace data. The trace header information generally stores information such as line number, trace number, sampling point number, geodetic coordinates, shot point coordinates, and demodulator probe coordinates corresponding to the seismic trace, but some key parameter positions (such as the line number and the position of the trace number in the trace header) are not fixed. The seismic trace data is obtained by sampling the waveform of the seismic signal at a certain time interval delta t and recording the series of discrete amplitude values in a certain mode.
Step 102: the shot point range and the receiver point range which are counted in advance are uniformly divided by the bin size parameters to obtain a plurality of rectangular bins, the rectangular bins are sequentially sorted and recorded, and bin numbers of the rectangular bins are determined to determine each of the shot point coordinates and the receiver point coordinates of the central point, which correspond to the rectangular bins.
Step 103: establishing a surface element data index structure containing surface element numbers according to all seismic channel data in the seismic record data; the header file of the bin data index structure comprises shot offset distances of rectangular bin center coordinates, and the bin data index structure further comprises offset distances of each seismic channel data; the offset of the central coordinate of the rectangular surface element is the difference between the central point demodulator probe coordinate and the central point shot coordinate, and the offset of the seismic channel data is the difference between the demodulator probe coordinate and the shot coordinate stored in the channel head information corresponding to the seismic channel data.
Step 104: and synthesizing the seismic channel data of each surface element according to the surface element data index structure, the local plane wave direction angle interval and the local plane wave direction angle number to obtain the common offset range unilateral local plane waves of a plurality of single surface elements.
Step 105: extracting effective in-phase axis data of the common offset range unilateral local plane wave of each surface element according to the number of the effective in-phase axis data during output; the effective in-phase axis data is in-phase axis data of the effective local plane wave.
Step 106: and mapping the effective homophase axis data of the common offset range single-sided local plane waves to a shot-geophone domain from the common offset range, synthesizing an inverse problem formula according to shot-geophone domain double-sided beams under the result constraint of the common offset range, and synthesizing the common offset range single-sided local plane waves after mapping operation into a plurality of single-surface element shot-geophone domain double-sided local plane waves.
Step 107: and extracting effective homophase axis data of bilateral local plane waves of the shot domain of each surface element and compressing the effective homophase axis data.
And stopping until all the surface elements are subjected to the operations, and then checking all the data and the indexes, and checking the consistency of the two files. All data here are generated local plane wave data and associated index file data,
the main contents of the steps 101 to 103 are the bilateral rectangular binning of the shot domain and the binning index establishment of input data.
The binning is according to the conventional method, namely, the binning size parameters are read in from the input local plane wave synthesis parameters, then according to the shot point range and the demodulator probe range of statistics, the binning size parameters are used for uniformly dividing the shot point range and the demodulator probe range of statistics to obtain a plurality of rectangular bins, the rectangular bins are sequentially marked in sequence, the bin numbers are recorded, and then a bin data index structure comprising the bin numbers is established for all seismic channel data of the input seismic record data; the method is mainly different from a non-offset domain result constraint method in that shot-geophone offset of rectangular bin center coordinates is added into a header file of a bin data index structure, and offset information of each seismic channel data is added into a built face metadata index structure.
Step 104 mainly comprises single-surface element seismic channel search, partial dynamic correction of common offset range seismic channels and local plane wave synthesis of common offset range.
Seismic trace search for a single bin: searching all seismic channel data corresponding to each bin number from the bin data index structure, specifically, firstly sequencing the head files of the bin data index structure, and then searching.
Portions of common offset domain seismic tracesDynamic correction: correcting time difference caused by the fact that the offset distance of all seismic channel data corresponding to each bin number is different from the offset distance of shot-geophone of the rectangular bin center coordinate corresponding to the bin number to obtain seismic channel data after partial dynamic correction; the correction formula is as follows:
where Δ t is the travel time difference between the seismic traces in the bin and the bin center point seismic trace, Δ h is the offset difference between the seismic traces in the bin and the bin center point seismic trace, and h is the offset difference between the seismic traces in the bin and the bin center point seismic tracemIs the offset, t, of the seismic trace at the center point of the bin0Is the travel time of the seismic trace with zero offset of the surface element center point, vnIs the motion correction speed.
Local plane wave synthesis in the common offset domain: and synthesizing the common offset domain unilateral local plane waves of the single surface elements according to the seismic channel data corresponding to each surface element after partial dynamic correction, the local plane wave direction angle interval in the local plane wave synthesis parameters and the local plane wave direction angle number.
In order to ensure the efficiency of practical application, a simple plane wave synthesis method is used, and the synthesis formula is as follows:wherein,is the coordinate of the center point and is, the data are single-sided local plane wave data of a common offset range of a single surface element;is the direction of the plane wave at the center point,the method is determined by the angle interval of the direction of the local plane wave, the number of the angle of the direction of the local plane wave, the number of the picked effective homophase axis data during output and a seismic velocity model; ω is the angular frequency; a transformation matrix, AHIs the conjugate transpose of A; d is the seismic trace data contained in the bin,is the seismic trace coordinates within the bin of the geophone point,and the coordinates of the seismic channels in the shot surface element.
And extracting and compressing the effective homophase axis data of the unilateral local plane wave in the common offset range corresponding to each surface element by adopting a threshold value method according to the number of the effective homophase axis data during output in the local plane wave synthesis parameters.
The unilateral local plane synthesis has no problem of huge data volume, and the effective homophase axis data extracted by a threshold value method can ensure higher precision. And (3) after single-side local plane wave synthesis of the single-side element common offset range is carried out by the formula (2), extracting effective in-phase axis data by using a defined threshold value.
The shot and examine regional local plane wave synthesis of result constraint of common offset distance region includes three main steps: projecting effective homophase axis data to a shot domain bilateral plane wave space; synthesizing local plane waves on two sides of a single-surface shot domain; and (4) picking and storing effective in-phase axis data of bilateral local plane waves of the single-surface element shot domain.
The shot-geophone domain bilateral beam synthesis inverse problem formula under the constraint of the common offset domain result is as follows:
wherein E is an error function orThis is called the objective function; ω is the angular frequency;andrespectively representing the seismic channel coordinates in a shot point bin and the seismic channel coordinates in a detector point bin; a is a conversion matrix, and lambda is a super parameter coefficient and is a small constant value; d is all seismic trace data contained in the singleton;the local plane wave data of the shot domain and the two sides of the shot domain of the single surface element are obtained;the data are single-sided local plane wave data of a common offset range of a single surface element;respectively the plane wave directions of the central point, the shot point and the wave detection point,the parameters are input from the outside: the direction angle interval of the local plane wave, the direction angle number of the local plane wave, the picking number of effective homophase axis data during output and the determination of a seismic velocity model.
Common offset range single-sided local plane wave data of single surface element in formula (3)Is in a coordinate systemIn the coordinate system, the shot domain bilateral local plane wave data of the single surface elementIn case of useAs the synthetic constraint condition of bilateral local plane waves of the shot domain of the single surface element, the synthetic constraint condition needs to be satisfiedIs mapped into
The effective homophase axis data of the common offset range unilateral local plane wave of each surface element is stored in a vector list of a user-defined data structure (containing position and amplitude information), and the effective homophase axis data is obtained according to the ray direction of the central point of the plane waveAnd shot inspection point ray directionRelationships, e.g. equation (4), take valid in-phase axis data from the common offset domainMapping to shot domain
The mapping formula is
Mapping the single-sided local plane wave of the common offset range of the single surface element after the operationAnd (3) carrying out local plane wave synthesis in the shot domain, and then picking up effective homophase axis data of the bilateral local plane waves of each shot domain by using a homophase axis picking method and compressing and storing the effective homophase axis data.
Minimizing the objective function E in equation (3) can result,
the above formula is further simplified, and the Hessian matrix term (A) is removedH+A+λI)-1To obtain
Wherein A is a transformation matrix, AHIs the conjugate transpose of a.
In order to achieve the purpose, the invention also discloses a shot-geophone domain bilateral beam synthesis system with result constraint in a common offset range domain.
Fig. 2 is a schematic structural diagram of a shot-geophone domain bilateral beam synthesis system based on result constraint of a common offset range domain in an embodiment of the present invention, and as shown in fig. 2, the shot-geophone domain bilateral beam synthesis system provided in an embodiment of the present invention includes:
a module 100 for acquiring seismic recording data and local plane wave synthesis parameters, configured to acquire the seismic recording data and the local plane wave synthesis parameters; the seismic record data is seismic exploration field acquisition data or laboratory synthesis data; the local plane wave synthesis parameters comprise the size of a surface element, the angle interval of the direction of the local plane wave, the number of the direction angles of the local plane wave, the number of effective homophase axis data pick-up during output and a seismic velocity model; the standard format of the seismic record data is SEGY format, and the SEGY format generally comprises three parts, namely an EBCDIC file header, a binary file header and a seismic channel; each seismic channel comprises channel head information and seismic channel data; and the line number, the track number, the sampling point number, the geodetic coordinate, the shot point coordinate and the demodulator probe coordinate corresponding to the seismic channel are stored in the channel head information.
The rectangular surface element and the surface element number obtain a module 200 for uniformly dividing the shot point range and the receiver point range counted in advance by adopting the surface element size parameters to obtain a plurality of rectangular surface elements, and right the rectangular surface elements are sequentially sorted and recorded the surface element number of the rectangular surface element determines each of the shot point coordinates and the receiver point coordinates of the central point corresponding to the rectangular surface element.
A surface element data index structure establishing module 300, configured to establish a surface element data index structure including surface elements according to all seismic channel data in the seismic record data; the header file of the bin data index structure comprises shot offset distances of rectangular bin center coordinates, and the bin data index structure further comprises offset distances of each seismic channel data; the offset of the central coordinate of the rectangular surface element is the difference between the central point demodulator probe coordinate and the central point shot coordinate, and the offset of the seismic channel data is the difference between the demodulator probe coordinate and the shot coordinate stored in the channel head information corresponding to the seismic channel data.
And the common offset range single-sided local plane wave synthesis module 400 is configured to synthesize seismic channel data of each surface element according to the surface element data index structure, the local plane wave direction angle interval, and the local plane wave direction angle number, so as to obtain a plurality of single-sided common offset range single-sided local plane waves of the surface elements.
A first effective in-phase axis data extraction module 500, configured to extract effective in-phase axis data of the common offset range single-sided local plane wave of each bin according to the number of effective in-phase axis data picks up during the output; the effective in-phase axis data is in-phase axis data of the effective local plane wave.
And the shot-check domain bilateral local plane wave synthesis module 600 is configured to map effective homophase axis data of the common-offset-range unilateral local plane wave from the common offset range to the shot-check range, and synthesize an inversion problem formula and a common-offset-range unilateral local plane wave after mapping operation according to a shot-check-range bilateral beam combination under the common offset range result constraint, so as to synthesize a plurality of single-surface-element shot-check-range bilateral local plane waves.
And the second effective event data extraction module 700 is configured to extract and compress effective event data of the local plane wave at both sides of the shot and examine region of each bin.
In order to illustrate the beneficial effects of the invention, a three-dimensional salt dome model test is used for illustration.
As shown in fig. 3 and 4, the seismic velocity model of the three-dimensional salt dome model has a 676x676x210 sampling grid, which represents the number of sampling points in the y direction, the x direction, and the z direction, respectively, and the sampling interval is 20mx20mx20 m. The recorded data is synthesized, as shown in fig. 5, the shot number 47, the longitudinal time sampling interval is 15ms, one shot line (line number is 246), 96 shots in total, the shot interval is 100m, the first shot position (2000m, 4900m), 4 lines (line numbers 241, 246, 251, 256, and the line interval of the detection points is 100m) are arranged on the detection points, the distance between the detection points is 25m, the first transverse position is 800m, and 97 detection points are in total. The imaging point range Inline number ranges from 241 to 256, the Xline number ranges from 71 to 606, and therefore the final imaging result is detected by using only the Inline number 246 line.
The local plane wave synthesis parameters are bin sizes of 300mx300m, the total number of the local plane wave original bilateral angle sampling total number of one bin is 961, the common offset range single-sided local plane wave synthesis has 570 effective bins, the maximum inclination angle is 60 degrees, 4-degree intervals are provided, and the three-dimensional angle sampling shot detection edges are the same. FIG. 6 is a schematic diagram of a local plane wave of the 300 th surface element of the common offset range without partial motion correction moveout correction, as shown in FIG. 6, with the abscissa beingDevelopment of two-dimensional coordinates of inclination and azimuthFrom-180 deg. to 180 deg. with 12 deg. intervals, one azimuth angle comprising 31 dip samples, from-60 deg. to 60 deg. with 4 deg. intervals, and the abscissa of the following figures is so set. FIG. 7 is partial plane wave data, more focused than FIG. 6, obtained with partial motion correction moveout correctionEspecially the shallow layer is more focused than in fig. 6; then, the result of fig. 7 is subjected to a threshold method to obtain the effective local plane wave data of fig. 8; further, the data in fig. 8 are mapped to the two sides of the shot domain to obtain results in fig. 9 and fig. 10, fig. 9 shows that the azimuth angle and the inclination angle of the shot edge are both zero, and the three-dimensional local plane wave result of the edge of the wave detection point is displayed, fig. 10 shows that the azimuth angle and the inclination angle of the edge of the wave detection point are zero, and the three-dimensional local plane wave result of the shot edge is displayed, and fig. 9 and fig. 10 are symmetrical to each other.
The local plane wave synthesis in the shot domain is constrained by the formula (3) using the results of fig. 9 and 10, so that the local plane wave result with a relatively focused shot domain can be obtained, and becomes very sparse and easy to compress as shown in fig. 11 and 12. The results of the three-dimensional study Inline246 were obtained by beam shifting the data shown in fig. 11 and 12, and as shown in fig. 13, it was demonstrated that the sparse compressed data contained sufficient data information to obtain a clear subsurface structure.
The principles and embodiments of the present invention have been described herein using specific examples, which are provided only to help understand the method and the core concept of the present invention; meanwhile, for a person skilled in the art, according to the idea of the present invention, the specific embodiments and the application range may be changed. In view of the above, the present disclosure should not be construed as limiting the invention.
Claims (10)
1. A shot domain bilateral beam synthesis method based on common offset domain result constraint is characterized by comprising the following steps:
acquiring seismic record data and local plane wave synthesis parameters; the seismic record data is seismic exploration field acquisition data or laboratory synthesis data; the local plane wave synthesis parameters comprise the size of a surface element, the angle interval of the direction of the local plane wave, the number of the direction angles of the local plane wave, the number of effective homophase axis data pick-up during output and a seismic velocity model; the standard format of the seismic record data is SEGY format, and the SEGY format generally comprises three parts, namely an EBCDIC file header, a binary file header and a seismic channel; each seismic channel comprises channel head information and seismic channel data; the line number, the track number, the sampling point number, the geodetic coordinate, the shot point coordinate and the demodulator probe coordinate corresponding to the seismic channel are stored in the channel head information;
uniformly dividing a shot point range and a receiver point range which are counted in advance by using the bin size parameters to obtain a plurality of rectangular bins, sequencing the rectangular bins in sequence, recording bin numbers of the rectangular bins, and determining a central point shot point coordinate and a central point receiver point coordinate which correspond to each rectangular bin;
establishing a surface element data index structure containing surface element numbers according to all seismic channel data in the seismic record data; the header file of the bin data index structure comprises shot offset distances of rectangular bin center coordinates, and the bin data index structure further comprises offset distances of each seismic channel data; the offset of the central coordinate of the rectangular bin is the difference between the coordinate of the central point geophone point and the coordinate of the central point shot point, and the offset of the seismic channel data is the difference between the coordinate of the geophone point and the coordinate of the shot point stored in the channel head information corresponding to the seismic channel data;
synthesizing the seismic channel data of each surface element according to the surface element data index structure, the local plane wave direction angle interval and the local plane wave direction angle number to obtain a common offset range unilateral local plane wave of a plurality of single surface elements;
extracting effective in-phase axis data of the common offset range unilateral local plane wave of each surface element according to the number of the effective in-phase axis data during output; the effective in-phase axis data is in-phase axis data of effective local plane waves;
and mapping the effective homophase axis data of the common offset range single-sided local plane waves from the common offset range to a shot-geophone range, and synthesizing the shot-geophone range double-sided local plane waves of a plurality of single-surface elements according to a shot-geophone range double-sided beam synthesis inverse problem formula under the common offset range result constraint and the common offset range single-sided local plane waves after mapping operation.
2. The shot domain double-side beam synthesis method according to claim 1, further comprising: and extracting effective homophase axis data of bilateral local plane waves of the shot domain of each surface element and compressing the effective homophase axis data.
3. The shot-geophone domain double-edge beam synthesis method according to claim 1, wherein said synthesizing seismic trace data of each bin according to said bin data index structure, said local plane wave direction angle interval, and said local plane wave direction angle number to obtain single-edge local plane waves of common offset range of a plurality of single bins specifically comprises:
searching all seismic channel data corresponding to each bin number from the bin data index structure;
correcting time difference caused by the fact that the offset distance of all seismic channel data corresponding to each bin number is different from the offset distance of shot-geophone of the rectangular bin center coordinate corresponding to the bin number to obtain corrected seismic channel data corresponding to each bin number;
and synthesizing the common offset domain unilateral local plane waves of the single surface elements according to the corrected seismic channel data corresponding to each surface element number, the local plane wave direction angle interval and the local plane wave direction angle number.
4. The shot domain double-side beam synthesis method of claim 3, wherein the correction formula of the seismic trace data isWherein Δ t is the travel time difference between the seismic traces within the bin and the seismic trace at the center point of the bin; Δ h is the difference in offset between the seismic traces in the bin and the seismic trace at the center point of the bin; h ismIs the offset, t, of the seismic trace at the center point of the bin0Travel time of seismic channel with zero offset of surface element center point；vnIs the motion correction speed.
5. The shot-geophone domain double-side beam synthesis method according to claim 3, wherein said common offset domain single-side local plane wave synthesis formula is:
wherein,is the coordinate of the center point and is, the data are single-sided local plane wave data of a common offset range of a single surface element;is the direction of the plane wave at the center point,the method is determined by the angle interval of the direction of the local plane wave, the number of the angle of the direction of the local plane wave, the number of the picked effective homophase axis data during output and a seismic velocity model; ω is the angular frequency; a transformation matrix, AHIs the conjugate transpose of A; d is the seismic trace data contained in the bin,is the seismic trace coordinates within the bin of the geophone point,and the coordinates of the seismic channels in the shot surface element.
6. The shot-examine domain double-side beam synthesis method according to claim 1, wherein the extracting effective in-phase axis data of the common offset domain single-side local plane wave of each bin according to the number of the effective in-phase axis data picks up at the time of the outputting specifically comprises:
and extracting and compressing the effective in-phase axis data of the common offset range unilateral local plane wave corresponding to each surface element by adopting a threshold value method according to the number of the effective in-phase axis data during output.
7. The shot domain double-side beam synthesis method according to claim 5, wherein the mapping the effective in-phase axis data of the common offset domain single-side local plane wave from the common offset domain to the shot domain specifically comprises: mapping effective homophase axis data of the common offset range single-sided local plane wave to a shot detection range from the common offset range by adopting the following formula; the formula isWherein,respectively the plane wave directions of the central point, the shot point and the wave detection point,the method is determined by the local plane wave direction angle interval, the local plane wave direction angle number, the effective homophase axis data pick-up number during output and a seismic velocity model.
8. The shot-geophone domain double-sided beam synthesis method according to claim 7, wherein said shot-geophone domain double-sided beam synthesis inverse problem formula under common offset range result constraint is:
the shot-geophone domain bilateral beam synthesis inverse problem formula under the constraint of the common offset domain result is as follows:
wherein E is an error function or an objective function; ω is the angular frequency;andrespectively representing the seismic channel coordinates in a shot point bin and the seismic channel coordinates in a detector point bin; a, converting a matrix, wherein lambda is a hyper-parameter coefficient and is a constant; d is all seismic trace data contained in the singleton;the local plane wave data of the shot domain and the two sides of the shot domain of the single surface element are obtained;the data are single-sided local plane wave data of a common offset range of a single surface element;respectively the plane wave directions of the central point, the shot point and the wave detection point,the method is determined by the local plane wave direction angle interval, the local plane wave direction angle number, the effective homophase axis data pick-up number during output and a seismic velocity model.
9. The shot domain bilateral beam synthesis method of claim 8, wherein the synthesis formula of the shot domain bilateral local plane wave is as follows:
10. a shot domain bilateral beam synthesis system with common offset domain result constraint is characterized by comprising:
the seismic record data and local plane wave synthesis parameter acquisition module is used for acquiring seismic record data and local plane wave synthesis parameters; the seismic record data is seismic exploration field acquisition data or laboratory synthesis data; the local plane wave synthesis parameters comprise the size of a surface element, the angle interval of the direction of the local plane wave, the number of the direction angles of the local plane wave, the number of effective homophase axis data pick-up during output and a seismic velocity model; the standard format of the seismic record data is SEGY format, and the SEGY format generally comprises three parts, namely an EBCDIC file header, a binary file header and a seismic channel; each seismic channel comprises channel head information and seismic channel data; the line number, the track number, the sampling point number, the geodetic coordinate, the shot point coordinate and the demodulator probe coordinate corresponding to the seismic channel are stored in the channel head information;
the rectangular surface element and surface element number obtaining module is used for uniformly dividing a shot point range and a receiver point range which are counted in advance by adopting surface element size parameters to obtain a plurality of rectangular surface elements, sequencing the rectangular surface elements in sequence, recording surface element numbers of the rectangular surface elements, and determining a shot point coordinate of a central point and a receiver point coordinate of the central point corresponding to each rectangular surface element;
the bin data index structure establishing module is used for establishing a bin data index structure containing bin numbers according to all seismic channel data in the seismic record data; the header file of the bin data index structure comprises shot offset distances of rectangular bin center coordinates, and the bin data index structure further comprises offset distances of each seismic channel data; the offset of the central coordinate of the rectangular bin is the difference between the coordinate of the central point geophone point and the coordinate of the central point shot point, and the offset of the seismic channel data is the difference between the coordinate of the geophone point and the coordinate of the shot point stored in the channel head information corresponding to the seismic channel data;
the common offset range single-sided local plane wave synthesis module is used for synthesizing the seismic channel data of each surface element according to the surface element data index structure, the local plane wave direction angle interval and the local plane wave direction angle number to obtain common offset range single-sided local plane waves of a plurality of single surface elements;
the first effective in-phase axis data extraction module is used for extracting effective in-phase axis data of the common offset range unilateral local plane wave of each surface element according to the number of the effective in-phase axis data during output; the effective in-phase axis data is in-phase axis data of effective local plane waves;
and the shot domain bilateral local plane wave synthesis module is used for mapping the effective homophase axis data of the common offset domain unilateral local plane wave to a shot domain from the common offset domain, synthesizing an inversion problem formula according to a shot domain bilateral beam under the common offset domain result constraint and the common offset domain bilateral local plane wave after mapping operation, and synthesizing a plurality of single-surface element shot domain bilateral local plane waves.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810776052.6A CN109031410B (en) | 2018-07-16 | 2018-07-16 | The bilateral beam synthetic method in big gun inspection domain and system of common offset field result constraint |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810776052.6A CN109031410B (en) | 2018-07-16 | 2018-07-16 | The bilateral beam synthetic method in big gun inspection domain and system of common offset field result constraint |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109031410A CN109031410A (en) | 2018-12-18 |
| CN109031410B true CN109031410B (en) | 2019-07-02 |
Family
ID=64642646
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810776052.6A Active CN109031410B (en) | 2018-07-16 | 2018-07-16 | The bilateral beam synthetic method in big gun inspection domain and system of common offset field result constraint |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109031410B (en) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102854532A (en) * | 2011-06-30 | 2013-01-02 | 中国石油天然气集团公司 | Three-dimensional pre-stack offset stochastic noise suppression method |
| CN102866421A (en) * | 2012-09-04 | 2013-01-09 | 中国科学院地质与地球物理研究所 | Scattered wave pre-stack imaging method for identifying small-fault throw breakpoints |
| CN102914796A (en) * | 2012-08-21 | 2013-02-06 | 北京多分量地震技术研究院 | Control method for acquiring offset speeds of longitudinal and transverse waves based on Gaussian beam |
| CN103744115A (en) * | 2013-12-31 | 2014-04-23 | 中国石油天然气股份有限公司 | Method and system for carrying out data regularization processing in common offset distance and common azimuth angle domain |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US9435903B2 (en) * | 2013-10-11 | 2016-09-06 | Chevron U.S.A. Inc. | System and method for regularizing seismic data |
-
2018
- 2018-07-16 CN CN201810776052.6A patent/CN109031410B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102854532A (en) * | 2011-06-30 | 2013-01-02 | 中国石油天然气集团公司 | Three-dimensional pre-stack offset stochastic noise suppression method |
| CN102914796A (en) * | 2012-08-21 | 2013-02-06 | 北京多分量地震技术研究院 | Control method for acquiring offset speeds of longitudinal and transverse waves based on Gaussian beam |
| CN102866421A (en) * | 2012-09-04 | 2013-01-09 | 中国科学院地质与地球物理研究所 | Scattered wave pre-stack imaging method for identifying small-fault throw breakpoints |
| CN103744115A (en) * | 2013-12-31 | 2014-04-23 | 中国石油天然气股份有限公司 | Method and system for carrying out data regularization processing in common offset distance and common azimuth angle domain |
Non-Patent Citations (1)
| Title |
|---|
| 五维规则化技术研究与应用;王伟 等;《石油地球物理勘探》;20171130;第52卷;第28-33页 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109031410A (en) | 2018-12-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Schneider | Developments in seismic data processing and analysis (1968–1970) | |
| CN101907725B (en) | Crack prediction method and device | |
| EP2920615A1 (en) | Process for creating image gathers | |
| CN101609160A (en) | A kind of automatic seismic data first-arrival wave pickup method of strong background noise resistance | |
| CN108828668B (en) | Pre-stack time migration data processing method and device | |
| US11397273B2 (en) | Full waveform inversion in the midpoint-offset domain | |
| CN110187383A (en) | A kind of quick method for separating of sea wide-azimuth seismic data COV trace gather | |
| CN112180433B (en) | Method and device for picking up first arrival wave of earthquake | |
| CN112394414B (en) | Two-step seismic diffraction wave field prestack separation method | |
| Sopher | Converting scanned images of seismic reflection data into SEG-Y format | |
| CN108226997A (en) | Seismic facies division method based on pre-stack seismic data | |
| US7369944B2 (en) | Method of processing seismic data for AVO or AVOA characterization | |
| CN113885077A (en) | Multi-seismic-source seismic data separation method based on deep learning | |
| CN105510975B (en) | Improve the method and device of geological data signal to noise ratio | |
| CN109031410B (en) | The bilateral beam synthetic method in big gun inspection domain and system of common offset field result constraint | |
| CA2606203A1 (en) | Seismic data processing method for rmo picking | |
| CA2533011C (en) | Method for bispectral picking of anelliptical nmo correction parameters | |
| CN114545495A (en) | OVT trace gather processing method and device for PS wave seismic trace data and electronic equipment | |
| CN112630822B (en) | Method and system for processing OVT domain pre-stack seismic data | |
| CN109031417B (en) | The bilateral synthesis beam offset imaging method and system picked up based on lineups | |
| AU2002333947B2 (en) | Improvements to surface picking techniques | |
| CN106019373B (en) | A kind of common conversion point gathering method and device | |
| GB2149503A (en) | Improvements relating to seismic surveying | |
| CN112698394B (en) | Method and device for determining size of design surface element of two-dimensional observation system | |
| CN111722274B (en) | Method and device for picking up seismic wave time-distance curve |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |