CN112698394B - Method and device for determining size of design surface element of two-dimensional observation system - Google Patents
Method and device for determining size of design surface element of two-dimensional observation system Download PDFInfo
- Publication number
- CN112698394B CN112698394B CN201911004066.7A CN201911004066A CN112698394B CN 112698394 B CN112698394 B CN 112698394B CN 201911004066 A CN201911004066 A CN 201911004066A CN 112698394 B CN112698394 B CN 112698394B
- Authority
- CN
- China
- Prior art keywords
- seismic waves
- interpolation
- seismic
- field
- gradient field
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 31
- 238000005070 sampling Methods 0.000 claims abstract description 58
- 230000006870 function Effects 0.000 claims description 14
- 238000004590 computer program Methods 0.000 claims description 6
- 238000011084 recovery Methods 0.000 abstract description 28
- 230000000694 effects Effects 0.000 abstract description 6
- 238000003491 array Methods 0.000 abstract 1
- 230000002349 favourable effect Effects 0.000 abstract 1
- 238000005457 optimization Methods 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 10
- 239000011295 pitch Substances 0.000 description 9
- 238000005516 engineering process Methods 0.000 description 7
- 238000004458 analytical method Methods 0.000 description 5
- 230000008569 process Effects 0.000 description 5
- 238000004422 calculation algorithm Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000003068 static effect 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/30—Analysis
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/16—Receiving elements for seismic signals; Arrangements or adaptations of receiving elements
- G01V1/20—Arrangements of receiving elements, e.g. geophone pattern
-
- 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
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Remote Sensing (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 embodiment of the invention provides a method and a device for determining the size of a design surface element of a two-dimensional observation system, wherein the method comprises the following steps: acquiring a direction field of seismic waves according to the actual propagation speed of the seismic waves sampled in the area to be sampled; acquiring a gradient field of the seismic wave according to the direction field of the seismic wave; carrying out transverse sampling interpolation on the shot gather data of the seismic waves according to the gradient field of the seismic waves; and determining an interpolation channel with the frequency and the signal-to-noise ratio according with the preset identification degree, and taking the determined interpolation distance of the interpolation channel as the size of the design surface element of the two-dimensional observation system. The scheme can realize data recovery capability and acquisition cost, has good guiding effect on selection of field acquisition parameters, sampling encryption and optimization of acquisition design schemes, has non-negligible effect on development of high-density seismic data acquisition design with ultra-long and ultra-wide arrays, and is favorable for having wide application prospect.
Description
Technical Field
The invention relates to the technical field of geophysical exploration, in particular to a method and a device for determining the size of a design surface element of a two-dimensional observation system.
Background
Aiming at a two-dimensional observation system in the geophysical exploration process, the two-dimensional observation system is used for acquiring seismic data, the sparsity degree of the acquired data depends on the size of a design surface element of the two-dimensional observation system, the larger the size of the design surface element is, the more sparse the acquired data is, the smaller the size of the design surface element is, and the denser the acquired data is, however, the two-dimensional observation system has the following defects at present: the too large size of the design surface element causes the sparsity of sparsely sampled data, and the data recovery capability is reduced; or the size of the designed surface element is too small, so that the acquired data is too dense, and the cost of acquiring the data is increased.
Disclosure of Invention
The embodiment of the invention provides a method for determining the size of a design surface element of a two-dimensional observation system, which aims to solve the technical problem that the two-dimensional observation system in the prior art cannot give consideration to both data recovery capability and cost in the data acquisition process. The method comprises the following steps:
acquiring a direction field of seismic waves according to the actual propagation speed of the seismic waves sampled in the area to be sampled;
acquiring a gradient field of the seismic wave according to the direction field of the seismic wave;
carrying out transverse sampling interpolation on the shot gather data of the seismic waves according to the gradient field of the seismic waves;
and taking the interpolation distance corresponding to the interpolation channel with the frequency and the signal-to-noise ratio conforming to the preset identification degree as the size of the design surface element of the two-dimensional observation system.
The embodiment of the invention also provides a device for determining the size of the design surface element of the two-dimensional observation system, so as to solve the technical problem that the two-dimensional observation system in the prior art cannot give consideration to both the data recovery capability and the cost in the data acquisition process. The device includes:
the directional field acquisition module is used for acquiring a directional field of seismic waves according to the actual propagation speed of the seismic waves sampled in the area to be sampled;
the gradient field acquisition module is used for acquiring a gradient field of the seismic wave according to the direction field of the seismic wave;
the interpolation module is used for carrying out transverse sampling interpolation on the shot gather data according to the gradient field of the seismic waves;
and the determining module is used for taking the interpolation distance corresponding to the interpolation channel with the frequency and the signal-to-noise ratio conforming to the preset identification degree as the size of the two-dimensional observation system design surface element.
In the embodiment of the invention, the direction field of the seismic waves is obtained through the actual propagation speed of the seismic waves sampled in the area to be sampled, the gradient field of the seismic waves is further obtained, then the transverse sampling interpolation is carried out on the shot gather data of the seismic waves according to the gradient field of the seismic waves, finally, the interpolation channel which accords with the preset identification degree is determined according to the frequency and the signal to noise ratio of the interpolation channel after interpolation, the interpolation distance of the interpolation channel is taken as the size of the design bin of the two-dimensional observation system, namely, the direction field of the seismic waves is determined, the data recovery capability of the seismic waves along different directions for sparse sampling can be researched by utilizing vector wave field analysis technology based on the direction field, the transverse sampling interpolation is carried out after the gradient field is determined based on the direction field with higher data recovery capability, finally, whether the interpolated data meet the data recovery requirement or not is judged according to the preset identification degree through whether the frequency and the signal to noise ratio or not, the interpolation distance which corresponds to the interpolation channel which accords with the identification degree is taken as the design bin size, the design bin size which is determined based on the direction field with higher data recovery capability, the design, the super-wide application of the super-wide data acquisition and the super-wide sampling recovery requirement can be met.
Drawings
The accompanying drawings, which are included to provide a further understanding of the invention and are incorporated in and constitute a part of this application, illustrate embodiment(s) of the invention and together with the description serve to explain the principles of the invention. In the drawings:
fig. 1 is a flowchart of a method for determining a design bin size of a two-dimensional observation system according to an embodiment of the present invention;
fig. 2 is a schematic diagram of aliasing and signal distortion of sampled data according to an embodiment of the present invention;
fig. 3 is a schematic diagram of data after interpolation recovery based on directionality according to an embodiment of the present invention;
FIG. 4 is a schematic diagram of a programming interface of a proper wavelength interpolation algorithm according to an embodiment of the present invention;
FIG. 5 is a schematic diagram of a forward simulation theoretical model according to an embodiment of the present invention;
FIG. 6 is a first schematic diagram illustrating comparison between sparsely sampled data and recovered data according to an embodiment of the present invention;
FIG. 7 is a second schematic diagram illustrating comparison between sparsely sampled data and recovered data according to an embodiment of the present invention;
FIG. 8 is a schematic diagram of comparing sparsely sampled data and recovered data according to a third embodiment of the present invention;
FIG. 9 is a block diagram of a computer device according to an embodiment of the present invention;
fig. 10 is a block diagram of a device for determining a design bin size of a two-dimensional observation system according to an embodiment of the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the present invention will be described in further detail with reference to the following embodiments and 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.
In an embodiment of the present invention, a method for determining a design bin size of a two-dimensional observation system is provided, as shown in fig. 1, the method includes:
step 102: acquiring a direction field of seismic waves according to the actual propagation speed of the seismic waves sampled in the area to be sampled;
step 104: acquiring a gradient field of the seismic wave according to the direction field of the seismic wave;
step 106: carrying out transverse sampling interpolation on the shot gather data of the seismic waves according to the gradient field of the seismic waves;
step 108: and taking the interpolation distance corresponding to the interpolation channel with the frequency and the signal-to-noise ratio conforming to the preset identification degree as the size of the design surface element of the two-dimensional observation system.
It can be known from the process shown in fig. 1 that, in the embodiment of the present invention, a direction field of seismic waves is obtained through an actual propagation speed of seismic waves sampled in a region to be sampled, and then a gradient field of the seismic waves is obtained, then transverse sampling interpolation is performed on shot gather data of the seismic waves according to the gradient field of the seismic waves, finally, an interpolation channel conforming to a preset identification degree is determined according to a frequency and a signal-to-noise ratio of the interpolation channel after interpolation, and an interpolation distance of the interpolation channel is taken as a size of a two-dimensional observation system design bin, that is, determination of the direction field of the seismic waves is achieved, a data recovery capability of sparse sampling of the seismic waves in different directions can be studied by using a wave field analysis technology based on a vector of the direction field, and then transverse sampling interpolation is performed after the gradient field is determined based on a direction field with a high data recovery capability, and finally, whether the interpolated data meeting a data recovery requirement of the data recovery requirement is measured by judging whether the frequency and the signal-to-noise ratio conforms to the preset identification degree or not, and further, and the interpolation distance corresponding to the data acquisition requirement of the design bin size can be taken into consideration when the super-wide sampling application of super-wide data acquisition and super-sampling data acquisition is achieved.
In specific implementation, the inventors of the present application found that sampling data obtained by sparse sampling may cause severe aliasing and signal distortion, as shown in fig. 2 (curve 1 represents an original signal, and curve 2 represents an acquired signal). If the motion direction information of each sampling point is obtained, a wave field signal without spurious frequency can be reconstructed through a vector interpolation technology, and the fitting effect of the seismic signal with directionality after interpolation reconstruction is shown in figure 3 (the inclined direction of a thick line at each black point in figure 3 represents the motion direction special direction of the sampling point).
In specific implementation, the directional field of the sampling point can be calculated by the following formula:
ATAN (Z component sample value/X component sample value)/pi X180
Aiming at seismic waves collected by using a two-component velocity detector, the seismic wave data comprise data parallel to a receiving direction and data vertical to the receiving direction, an X component is a component parallel to the receiving line direction, and the value of the X component is positive and negative characteristics; the Z component is a component perpendicular to the receive line direction and has positive and negative characteristics. And (3) calculating positive and negative values and changes in the result, namely calculating the directional characteristics of the target sampling point, wherein the numerical value set of the positive and negative values and changes is the numerical characteristics of the target directional field.
In specific implementation, after a direction field of seismic waves is acquired, data recovery capability of sparse sampling of the seismic waves in different directions can be analyzed by adopting a vector wave field analysis technology, and in order to improve the data recovery capability and reduce acquisition cost, the inventor of the present application finds that data parallel to a receiving direction has data recovery capability better than data perpendicular to the receiving direction, and therefore, proposes to perform interpolation data recovery based on the direction field parallel to the receiving direction, for example, obtain the direction field parallel to the receiving direction in the seismic waves according to the actual propagation speed of the seismic waves sampled in a region to be sampled, further obtain a gradient field based on the direction field parallel to the receiving direction, perform transverse sampling interpolation on shot gather data of the seismic waves based on the gradient field, and finally determine an interpolation channel according with a preset identification degree according to frequency and a signal-to-noise ratio, and use the interpolation distance of the interpolation channel as the size of a two-dimensional observation system design bin as a corresponding acquisition parameter for a geological structure of the region to be sampled for optimal design of an observation system. Because the data parallel to the receiving direction has stronger data recovery capability relatively, the method calculates the direction field and the gradient field based on the data parallel to the receiving direction, determines the interpolation channel according with the preset identification degree based on the frequency and the signal-to-noise ratio, and the like, so that the data recovery capability can be further ensured by adopting the determined interpolation distance as the size of the design surface element of the two-dimensional observation system, the acquired data can be relatively prevented from being too dense, and the data acquisition cost is favorably reduced.
In specific implementation, the gradient field of the seismic waves can be obtained according to a direction field parallel to the receiving direction in the seismic waves by the following formula:
wherein f (x) is seismic data parallel to the receiving direction in the seismic waves;g(x;θ j ) Is a parameter set theta j The basis function of (decomposing a single channel of data into a plurality of single frequency signals as the basis function); g x (x;θ j ) Is a basis function g (x; theta.theta. j ) The spatial gradient of (a); f. of x (x) A gradient field corresponding to seismic data parallel to the receive direction;to representDefinitional formulas (e.g., meaning "noted", "defined", "equivalent" and the like); x is the sampling frequency and j is the sampling dominant frequency.
In specific implementation, in order to implement interpolation based on seismic data parallel to the receiving direction to recover data and further determine the size of a two-dimensional observation system design surface element, in this embodiment, performing lateral sampling interpolation on shot gather data of seismic waves according to the gradient field of the seismic waves includes:
determining an interpolation distance according to the number of interpolation channels and a gradient field corresponding to the sampling main frequency in the seismic data parallel to the receiving direction; wherein, the gradient field corresponding to the sampling main frequency in the seismic data parallel to the receiving direction is a passing formulaA determined gradient field; and then carrying out transverse sampling interpolation on the shot gather data of the seismic waves according to the interpolation distance.
During specific implementation, in the process of performing transverse sampling interpolation on seismic wave shot gather data, transverse sampling interpolation can be performed by adopting an interpolation distance, whether the frequency and the signal-to-noise ratio of the interpolated recovered data meet the preset identification degree or not is judged, namely whether the frequency and the signal-to-noise ratio are greater than the preset identification degree or not is judged, and if yes, the interpolation distance can be used as the size of a design surface element of the two-dimensional observation system; and transverse sampling interpolation can be carried out by adopting a plurality of different interpolation distances, whether the frequency and the signal-to-noise ratio of the recovered data after each interpolation meet the preset identification degree or not is judged, namely whether the frequency and the signal-to-noise ratio are greater than the preset identification degree or not is judged, an interpolation channel with the frequency and the signal-to-noise ratio greater than the preset identification degree is determined, and the interpolation distance corresponding to the interpolation channel can be used as the size of a two-dimensional observation system design surface element.
In specific implementation, the preset identification degree can be a preset signal-to-noise ratio curve or amplitude change, and can also be an identifiable preset seismic wave hyperbolic characteristic.
In specific implementation, the method for determining the size of the design surface element of the two-dimensional observation system may be implemented by a program, an interface implemented by the program is shown in fig. 4, and the method for determining the size of the design surface element of the two-dimensional observation system may include the following steps:
step 4, sampling main frequency (namely basis function g (x; theta) in the seismic data parallel to the receiving direction according to the number of interpolation channels j ) Corresponding gradient field (i.e. g) x (x;θ j ) Determining an interpolation distance, and then performing interpolation calculation and wave field data recovery based on the interpolation distance;
step 5, analyzing interpolation effects of all basis functions of different track pitches, and performing frequency and signal-to-noise ratio related analysis;
and 6, taking the interpolation track pitch corresponding to the interpolation track with the frequency and the signal-to-noise ratio larger than the preset identification degree (specifically, the preset identification degree can be set according to the signal-to-noise ratio condition of the basis function autocorrelation data) as the size of the design surface element.
The method for determining the size of the design bin of the two-dimensional observation system is described below with reference to a specific example, where the region to be sampled is an example of a forward modeling theoretical model shown in fig. 5:
1) Extracting a three-component single shot data sample from test data which are acquired by an observation system and have seismic wave data of 100m track pitch, 20m shot pitch and 3S3L1000T wide line design, and identifying and extracting data of three single components X, Y, Z by referring to data track header word information;
2) Selecting horizontal X component data (as shown in a left graph in fig. 8) parallel to the receiving line direction, and respectively thinning the data with the track pitch of 200m (as shown in the left graph in fig. 7) and the track pitch of 400m (as shown in the left graph in fig. 6) for the horizontal X component data to simulate two kinds of sparse sampling data;
3) Recovering the 200 m-track sparse sampling data generated in the step 2 into 100 m-track (as shown in the right graph in fig. 8) sparse sampling data by using a vector wave field interpolation algorithm in the determination method for designing the bin size of the two-dimensional observation system, and recovering the 400 m-track sparse sampling data generated in the step 2 into 200 m-track (as shown in the right graph in fig. 6 and fig. 7) sparse sampling data by using a 100m interpolation distance;
4) The effect comparison is carried out on the original data of the track pitch of 100m and the three recovery data of the track pitch of 100m after the simulated thinning-vector interpolation in sequence, so that the obvious result that the recovery data of the track pitch of 100m after the thinning is 200m still has good identification degree in the aspect of reflected wave energy.
5) Accordingly, the design scheme that the channel distance of 200m is the size of the design surface element of the two-dimensional observation system can be optimized by considering the relevant requirements of controlling the acquisition cost and guaranteeing the acquisition quality.
In this embodiment, a computer device is provided, as shown in fig. 9, and includes a memory 902, a processor 904, and a computer program stored on the memory and executable on the processor, and the processor executes the computer program to implement any of the above-mentioned methods for determining a design bin size of a two-dimensional observation system.
In particular, the computer device may be a computer terminal, a server or a similar computing device.
In the present embodiment, there is provided a computer-readable storage medium storing a computer program for executing the above-described method for determining a design bin size of a two-dimensional observation system.
In particular, computer-readable storage media, including both non-transitory and non-transitory, removable and non-removable media, may implement information storage by any method or technology. The information may be computer readable instructions, data structures, modules of a program, or other data. Examples of computer-readable storage media include, but are not limited to, phase change memory (PRAM), static Random Access Memory (SRAM), dynamic Random Access Memory (DRAM), other types of Random Access Memory (RAM), read Only Memory (ROM), electrically Erasable Programmable Read Only Memory (EEPROM), flash memory or other memory technology, compact disc read only memory (CD-ROM), digital Versatile Discs (DVD) or other optical storage, magnetic cassettes, magnetic tape magnetic disk storage or other magnetic storage devices, or any other non-transmission medium that can be used to store information that can be accessed by a computing device. As defined herein, a computer readable storage medium does not include a transitory computer readable medium such as a modulated data signal and a carrier wave.
Based on the same inventive concept, the embodiment of the present invention further provides a device for determining a design bin size of a two-dimensional observation system, as described in the following embodiments. Because the principle of solving the problem of the device for determining the size of the design surface element of the two-dimensional observation system is similar to the method for determining the size of the design surface element of the two-dimensional observation system, the implementation of the device for determining the size of the design surface element of the two-dimensional observation system can refer to the implementation of the method for determining the size of the design surface element of the two-dimensional observation system, and repeated parts are not described again. As used hereinafter, the term "unit" or "module" may be a combination of software and/or hardware that implements a predetermined function. Although the means described in the embodiments below are preferably implemented in software, an implementation in hardware or a combination of software and hardware is also possible and contemplated.
Fig. 10 is a block diagram of an arrangement for determining a design bin size of a two-dimensional observation system according to an embodiment of the present invention, and as shown in fig. 10, the arrangement includes:
the directional field obtaining module 1002 is configured to obtain a directional field of a seismic wave according to an actual propagation velocity of the seismic wave sampled in a region to be sampled;
a gradient field obtaining module 1004, configured to obtain a gradient field of the seismic wave according to a directional field of the seismic wave;
an interpolation module 1006, configured to perform lateral sampling interpolation on shot gather data according to the gradient field of the seismic wave;
the determining module 1008 is configured to use an interpolation distance corresponding to an interpolation channel with a frequency and a signal-to-noise ratio that meet a preset identification degree as a size of a two-dimensional observation system design bin.
In one embodiment, the direction field obtaining module is specifically configured to obtain a direction field parallel to a receiving direction in the seismic waves according to an actual propagation speed of the seismic waves sampled in the region to be sampled;
the gradient field acquisition module is specifically used for acquiring a gradient field of the seismic waves according to a direction field parallel to the receiving direction in the seismic waves.
In one embodiment, the gradient field acquisition module acquires the gradient field of the seismic waves according to the direction field parallel to the receiving direction in the seismic waves by:
wherein f (x) is seismic data parallel to the receiving direction in the seismic waves;g(x;θ j ) Is a parameter set theta j A basis function of (a);g x (x;θ j ) Is the basis function g (x; theta j ) The spatial gradient of (a); f. of x (x) A gradient field corresponding to seismic data parallel to the receive direction;represents a defined formula; x is the sampling frequency and j is the sampling dominant frequency.
In an embodiment, the interpolation module is specifically configured to determine an interpolation distance according to the number of interpolation channels and a gradient field corresponding to a sampling main frequency in the seismic data parallel to the receiving direction; and carrying out transverse sampling interpolation on the shot gather data of the seismic waves according to the interpolation distance.
The embodiment of the invention realizes the following technical effects: in the embodiment of the invention, the direction field of the seismic waves is obtained through the actual propagation speed of the seismic waves sampled in the area to be sampled, the gradient field of the seismic waves is further obtained, then the transverse sampling interpolation is carried out on the shot gather data of the seismic waves according to the gradient field of the seismic waves, finally, the interpolation channel which accords with the preset identification degree is determined according to the frequency and the signal to noise ratio of the interpolation channel after interpolation, the interpolation distance of the interpolation channel is taken as the size of the design bin of the two-dimensional observation system, namely, the direction field of the seismic waves is determined, the data recovery capability of the seismic waves along different directions for sparse sampling can be researched by utilizing vector wave field analysis technology based on the direction field, the transverse sampling interpolation is carried out after the gradient field is determined based on the direction field with higher data recovery capability, finally, whether the interpolated data meet the data recovery requirement or not is judged according to the preset identification degree through whether the frequency and the signal to noise ratio or not, the interpolation distance which corresponds to the interpolation channel which accords with the identification degree is taken as the design bin size, the design bin size which is determined based on the direction field with higher data recovery capability, the design, the super-wide application of the super-wide data acquisition and the super-wide sampling recovery requirement can be met.
It will be apparent to those skilled in the art that the modules or steps of the embodiments of the invention described above may be implemented by a general purpose computing device, they may be centralized in a single computing device or distributed across a network of multiple computing devices, and alternatively, they may be implemented by program code executable by a computing device, such that it may be stored in a memory device and executed by a computing device, and in some cases, the steps shown or described may be executed out of order, or separately as individual integrated circuit modules, or multiple modules or steps may be implemented as a single integrated circuit module. Thus, embodiments of the invention are not limited to any specific combination of hardware and software.
The above description is only a preferred embodiment of the present invention, and is not intended to limit the present invention, and various modifications and changes may be made to the embodiment of the present invention by those skilled in the art. Any modification, equivalent replacement, or improvement made within the spirit and principle of the present invention should be included in the protection scope of the present invention.
Claims (10)
1. A method for determining the size of a design bin of a two-dimensional observation system is characterized by comprising the following steps:
acquiring a direction field of seismic waves according to the actual propagation speed of the seismic waves sampled in the area to be sampled;
acquiring a gradient field of the seismic wave according to the direction field of the seismic wave;
carrying out transverse sampling interpolation on the shot gather data of the seismic waves according to the gradient field of the seismic waves;
and taking the interpolation distance corresponding to the interpolation channel with the frequency and the signal-to-noise ratio conforming to the preset identification degree as the size of the design surface element of the two-dimensional observation system.
2. A method for determining a design bin size for a two dimensional observation system as claimed in claim 1 wherein obtaining the directional field of seismic waves based on the actual propagation velocity of seismic waves sampled in the area to be sampled comprises:
acquiring a direction field parallel to a receiving direction in the seismic waves according to the actual propagation speed of the seismic waves sampled in the area to be sampled;
acquiring a gradient field of the seismic wave according to the direction field of the seismic wave, wherein the gradient field of the seismic wave comprises the following steps:
and acquiring a gradient field of the seismic waves according to a direction field parallel to the receiving direction in the seismic waves.
3. The method for determining the design bin size of the two-dimensional observation system according to claim 2, wherein the step of obtaining the gradient field of the seismic waves according to the directional field parallel to the receiving direction in the seismic waves is realized by the following formula:
wherein f (x) is seismic data parallel to the receiving direction in the seismic waves;g(x;θ j ) Is a parameter set theta j A basis function of (a);g x (x;θ j ) Is the basis function g (x; theta j ) The spatial gradient of (a); f. of x (x) A gradient field corresponding to seismic data parallel to the receive direction;represents a definitional formula; x is the sampling frequency and j is the sampling dominant frequency.
4. A method for determining a design bin size for a two-dimensional observation system according to any one of claims 1 to 3, wherein the transversely sampling interpolation of the shot gather data of the seismic waves according to the gradient field of the seismic waves comprises:
determining an interpolation distance according to the number of interpolation channels and a gradient field corresponding to the sampling main frequency in the seismic data parallel to the receiving direction;
and carrying out transverse sampling interpolation on the shot gather data of the seismic waves according to the interpolation distance.
5. An apparatus for determining a design bin size for a two-dimensional vision system, comprising:
the directional field acquisition module is used for acquiring a directional field of seismic waves according to the actual propagation speed of the seismic waves sampled in the area to be sampled;
the gradient field acquisition module is used for acquiring a gradient field of the seismic wave according to the direction field of the seismic wave;
the interpolation module is used for carrying out transverse sampling interpolation on the shot gather data according to the gradient field of the seismic waves;
and the determining module is used for taking the interpolation distance corresponding to the interpolation channel with the frequency and the signal-to-noise ratio conforming to the preset identification degree as the size of the two-dimensional observation system design surface element.
6. The two-dimensional observation system design bin size determining apparatus of claim 5,
the directional field acquisition module is specifically used for acquiring a directional field parallel to the receiving direction in the seismic waves according to the actual propagation speed of the seismic waves sampled in the area to be sampled;
the gradient field acquisition module is specifically used for acquiring a gradient field of the seismic waves according to a direction field parallel to the receiving direction in the seismic waves.
7. The apparatus for determining the design bin size of a two-dimensional observation system according to claim 6, wherein the gradient field obtaining module obtains the gradient field of the seismic waves according to the direction field parallel to the receiving direction in the seismic waves by:
wherein f (x) is seismic data parallel to the receiving direction in the seismic waves;g(x;θ j ) Is a parameter set theta j A basis function of (a);g x (x;θ j ) Is the basis function g (x; theta.theta. j ) The spatial gradient of (a); f. of x (x) A gradient field corresponding to seismic data parallel to the receive direction;represents a definitional formula; x is the sampling frequency and j is the sampling main frequency.
8. The apparatus for determining a design bin size of a two-dimensional observation system according to any one of claims 5 to 7, wherein the interpolation module is specifically configured to determine the interpolation distance according to the number of interpolation channels and a gradient field corresponding to the sampling main frequency in the seismic data parallel to the receiving direction; and carrying out transverse sampling interpolation on the shot gather data of the seismic waves according to the interpolation distance.
9. A computer device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, wherein the processor when executing the computer program implements the method for determining a design bin size for a two-dimensional vision system as set forth in any one of claims 1 to 4.
10. A computer-readable storage medium, characterized in that the computer-readable storage medium stores a computer program for executing the method for determining a design bin size of a two-dimensional vision system as claimed in any one of claims 1 to 4.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911004066.7A CN112698394B (en) | 2019-10-22 | 2019-10-22 | Method and device for determining size of design surface element of two-dimensional observation system |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911004066.7A CN112698394B (en) | 2019-10-22 | 2019-10-22 | Method and device for determining size of design surface element of two-dimensional observation system |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112698394A CN112698394A (en) | 2021-04-23 |
CN112698394B true CN112698394B (en) | 2022-11-04 |
Family
ID=75504793
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911004066.7A Active CN112698394B (en) | 2019-10-22 | 2019-10-22 | Method and device for determining size of design surface element of two-dimensional observation system |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112698394B (en) |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5787051A (en) * | 1996-05-08 | 1998-07-28 | Pancanadian Petroleum Limited | Three-dimensional seismic acquisition |
CN103698807A (en) * | 2012-09-28 | 2014-04-02 | 中国石油天然气集团公司 | Scalariform two-dimensional wide-band observation system design method |
CN104422957A (en) * | 2013-08-30 | 2015-03-18 | 中国石油天然气集团公司 | Optimizing design method of observing system |
CN104536041A (en) * | 2014-12-17 | 2015-04-22 | 中国石油天然气集团公司 | Optimization method of seismological observation system parameters |
CN106501841A (en) * | 2016-09-19 | 2017-03-15 | 中国石油天然气集团公司 | A kind of Optimization Design of two-dimentional wode line seismic observation system and device |
CN106855635A (en) * | 2016-12-09 | 2017-06-16 | 中国石油天然气集团公司 | A kind of method and apparatus of the bin attributes information for determining orthogonal observation system |
CN109655937A (en) * | 2017-10-10 | 2019-04-19 | 中国石油化工股份有限公司 | A kind of evaluation method of earthquake data before superposition regularization observation system |
-
2019
- 2019-10-22 CN CN201911004066.7A patent/CN112698394B/en active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5787051A (en) * | 1996-05-08 | 1998-07-28 | Pancanadian Petroleum Limited | Three-dimensional seismic acquisition |
CN103698807A (en) * | 2012-09-28 | 2014-04-02 | 中国石油天然气集团公司 | Scalariform two-dimensional wide-band observation system design method |
CN104422957A (en) * | 2013-08-30 | 2015-03-18 | 中国石油天然气集团公司 | Optimizing design method of observing system |
CN104536041A (en) * | 2014-12-17 | 2015-04-22 | 中国石油天然气集团公司 | Optimization method of seismological observation system parameters |
CN106501841A (en) * | 2016-09-19 | 2017-03-15 | 中国石油天然气集团公司 | A kind of Optimization Design of two-dimentional wode line seismic observation system and device |
CN106855635A (en) * | 2016-12-09 | 2017-06-16 | 中国石油天然气集团公司 | A kind of method and apparatus of the bin attributes information for determining orthogonal observation system |
CN109655937A (en) * | 2017-10-10 | 2019-04-19 | 中国石油化工股份有限公司 | A kind of evaluation method of earthquake data before superposition regularization observation system |
Non-Patent Citations (3)
Title |
---|
关于面元计算和观测***设计的思考;高银波 等;《石油地球物理勘探》;20080815;第43卷(第04期);383-386,396 * |
观测***对高密度地震采集资料的影响;王海 等;《石油地球物理勘探》;20090415;第44卷(第02期);131-135 * |
高密度地震数据采集中参数选取方法探讨;王梅生 等;《勘探地球物理进展》;20091225;第32卷(第06期);404-408 * |
Also Published As
Publication number | Publication date |
---|---|
CN112698394A (en) | 2021-04-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US9915742B2 (en) | Method and system for geophysical modeling of subsurface volumes based on label propagation | |
CN108828668B (en) | Pre-stack time migration data processing method and device | |
US9529115B2 (en) | Geophysical modeling of subsurface volumes based on horizon extraction | |
CN115240149A (en) | Three-dimensional point cloud detection and identification method and device, electronic equipment and storage medium | |
CN114839673B (en) | Separation method, separation system and computer equipment for multi-seismic-source efficient acquisition wave field | |
CN108845350B (en) | Method and device for inverting two-dimensional velocity model | |
WO2014102616A2 (en) | Efficient wavefield extrapolation in anisotropic media | |
WO2021127382A1 (en) | Full waveform inversion in the midpoint-offset domain | |
CN112529841A (en) | Method and system for processing seabed gas plume in multi-beam water column data and application | |
CN112698394B (en) | Method and device for determining size of design surface element of two-dimensional observation system | |
CN113557528A (en) | Generating point cloud completion network and processing point cloud data | |
CN115586575B (en) | 3D passive source simulated cannon set calculation method, system, equipment and storage medium | |
CN114114421B (en) | Deep learning-based guided self-learning seismic data denoising method and device | |
CN108828655B (en) | Method and device for processing trace gather recording and stacking profiles and computer storage medium | |
CN114002742B (en) | Euler Gaussian beam offset imaging method and device | |
CN112630840B (en) | Random inversion method based on statistical characteristic parameters and processor | |
CN113424220B (en) | Processing for generating point cloud completion network and point cloud data | |
CN111965705B (en) | Method, device, equipment and medium for calibrating geological horizon recorded by earthquake single shot | |
Li et al. | Seismic full-waveform inversion based on superresolution for high-precision prediction of reservoir parameters | |
CN111914609A (en) | Well-seismic combined prestack geostatistics elastic parameter inversion method and device | |
CN113589364A (en) | Seismic data regularization processing method based on Zornia equation constraint | |
CN112083490A (en) | Seismic data noise attenuation method and device | |
CN113126163B (en) | Five-dimensional seismic data noise attenuation method and device | |
CN104704393A (en) | System and method for producing local images of subsurface targets | |
CN115293980B (en) | Small-size dynamic noise filtering method and device based on historical information |
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 |