CN108387926A - A kind of method and device of determining air-gun array far-field wavelet - Google Patents
A kind of method and device of determining air-gun array far-field wavelet Download PDFInfo
- Publication number
- CN108387926A CN108387926A CN201810104599.1A CN201810104599A CN108387926A CN 108387926 A CN108387926 A CN 108387926A CN 201810104599 A CN201810104599 A CN 201810104599A CN 108387926 A CN108387926 A CN 108387926A
- Authority
- CN
- China
- Prior art keywords
- field
- wavelet
- field wavelet
- air gun
- initial
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 50
- 238000001914 filtration Methods 0.000 claims abstract description 43
- 239000011159 matrix material Substances 0.000 claims description 36
- 238000001228 spectrum Methods 0.000 claims description 20
- 238000001514 detection method Methods 0.000 claims description 12
- 238000012545 processing Methods 0.000 claims description 9
- 230000003595 spectral effect Effects 0.000 claims description 7
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 6
- 238000010586 diagram Methods 0.000 description 11
- 238000003860 storage Methods 0.000 description 6
- 230000006870 function Effects 0.000 description 5
- 238000009826 distribution Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 4
- 238000005070 sampling Methods 0.000 description 4
- 238000003491 array Methods 0.000 description 2
- 238000004590 computer program Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000013016 damping Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000005055 memory storage Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000010979 ruby Substances 0.000 description 1
- 229910001750 ruby Inorganic materials 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 230000001052 transient effect Effects 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/02—Generating seismic energy
- G01V1/133—Generating seismic energy using fluidic driving means, e.g. highly pressurised fluids; using implosion
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/38—Seismology; Seismic or acoustic prospecting or detecting specially adapted for water-covered areas
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/10—Aspects of acoustic signal generation or detection
- G01V2210/12—Signal generation
- G01V2210/129—Source location
- G01V2210/1293—Sea
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/30—Noise handling
- G01V2210/32—Noise reduction
- G01V2210/324—Filtering
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Geology (AREA)
- Environmental & Geological Engineering (AREA)
- Acoustics & Sound (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Oceanography (AREA)
- Geophysics And Detection Of Objects (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
The embodiment of the present application discloses a kind of method and device of determining air-gun array far-field wavelet.The method is provided with the initial near-field wavelet for the receiving point for being excited by the air gun source in air-gun array and being recorded in the near field wave detector in the geophone group of near field；The method includes：According to the air gun source in the near field geophone group at a distance from the wave detector of near field and the virtual image of the air gun source at a distance from the wave detector of near field, determine filtering factor in the near field geophone group；The initial near-field wavelet is filtered using the filtering factor, obtains Target near field wavelet；According to the Target near field wavelet, the corresponding far-field wavelet of the air-gun array is determined.Technical solution provided by the embodiments of the present application can improve the precision of determined air-gun array far-field wavelet.
Description
Technical Field
The application relates to the technical fields of seismic exploration and acquisition of air gun sources, geological survey and the like in marine environments and inland lakes, in particular to a method and a device for determining far-field wavelets of an air gun array.
Background
The air gun seismic source has the characteristics of stable and reliable performance, safety, environmental protection and the like, and is the most widely applied seismic source in marine seismic exploration at present. When the air gun seismic source is excited, high-pressure gas in the air gun is instantaneously released and generates a large-volume bubble, and along with the damping type periodic motion of the bubble, far-field wavelets of the air gun array are formed and spread to the periphery. The designed air gun array is limited by field construction conditions and is difficult to achieve circular symmetry in a strict sense but has a certain length and width, so that the air gun array does not meet the assumed conditions of a point source, namely far-field wavelets of the air gun array do not spread to the periphery in a spherical mode, the energy of the far-field wavelets of the air gun array changes along with the changes of a horizontal azimuth angle and a vertical azimuth angle, and the far-field wavelets are generally called as the directional characteristics of the far-field wavelets of the air gun array. The directivity of the far-field wavelet of the air gun array is analyzed, the distribution condition of the far-field wavelet energy in different frequency ranges in a three-dimensional space can be determined, the air gun array for marine seismic exploration is designed and optimized according to the distribution condition, and obviously, the research is also significant for recovering the true amplitude in the processing of marine seismic exploration data.
At present, researchers usually directly determine the far-field wavelet of the air gun array by using the simulated wavelet and analyze the directivity of the far-field wavelet of the air gun array, however, the simulated wavelet is often different from the actual far-field wavelet of the air gun array, and therefore a method for determining the far-field wavelet of the air gun array with higher precision is urgently needed.
Disclosure of Invention
An object of the embodiments of the present application is to provide a method and an apparatus for determining far-field wavelets of an airgun array, so as to improve the accuracy of the determined far-field wavelets of the airgun array.
To solve the above technical problem, an embodiment of the present application provides a method and an apparatus for determining far-field wavelets of an airgun array, which are implemented as follows:
a method of determining a far-field wavelet for an airgun array, the method providing an initial near-field wavelet excited by an airgun source in the airgun array and recorded at a receiving point of a near-field detector in a near-field detector set; the method comprises the following steps:
determining a filter factor according to the distance between the air gun seismic source and a near-field detector in the near-field detector group and the distance between a virtual image of the air gun seismic source and the near-field detector in the near-field detector group;
filtering the initial near-field wavelet by using the filtering factor to obtain a target near-field wavelet;
and determining the far-field wavelet corresponding to the air gun array according to the target near-field wavelet.
In a preferred embodiment, the filtering the initial near-field wavelet by using the filtering factor to obtain a target near-field wavelet includes:
converting the initial near-field wavelet from a time domain to a frequency domain to obtain an initial near-field wavelet of the frequency domain, and establishing an initial near-field wavelet spectrum matrix according to the initial near-field wavelet of the frequency domain; the initial near-field wavelet spectrum matrix is used for representing an initial near-field wavelet set of a frequency domain recorded by the near-field detector group correspondingly;
and filtering the initial near-field wavelet set represented by the initial near-field wavelet spectrum matrix by using the filtering factor to obtain a filtered initial near-field wavelet set of the frequency domain, and taking the filtered initial near-field wavelet set of the frequency domain as the target near-field wavelet.
In a preferred embodiment, the filter factor is determined by the following formula:
wherein A represents the filter factor, rmnRepresents the distance between the receiving point of the mth near-field detector in the near-field detector group and the nth air gun seismic source in the air gun array，(rg)mnRepresenting the distance between the receiving point of the mth near-field detector in the near-field detector group and the virtual image of the nth air gun seismic source in the air gun array, C representing the speed of sound waves in water, omega representing angular frequency, i representing an imaginary number2＝-1；
Obtaining a target near-field wavelet spectrum matrix for characterizing the initial near-field wavelet set of the filtered frequency domain using the following formula:
X(ω)＝(ATA+λI)-1ATB
wherein X (ω) represents the target near-field wavelet spectral matrix, ATDenotes the transposed matrix of A, λ denotes the stability factor, Nm(ω) represents an initial near-field wavelet of a frequency domain recorded at a receiving point of an mth near-field detector in the near-field detector set, m represents the number of near-field detectors in the near-field detector set, and n represents the number of air gun vibration sources in an air gun array.
In a preferred embodiment, the value of the stability factor is 10-8。
In a preferred scheme, the far-field wavelet corresponding to the air gun array is determined by adopting the following formula:
wherein,the propagation distance of the far field sub-wave corresponding to the air gun array is RA vertical azimuth angle of theta and a horizontal azimuth angle of thetaAnd far-field wavelets with angular frequency omega, wherein the vertical azimuth angle theta represents an included angle between the propagation direction of the far-field wavelets in a Cartesian coordinate system and the z-axis of the Cartesian coordinate system, and the horizontal azimuth angle isRepresenting the included angle between the projection of the far-field wavelet propagation direction in the plane xoy of the Cartesian coordinate system and the x axis in the Cartesian coordinate system; r isn(r) represents the distance between the detection point of the far-field wavelet and the nth air gun source in the air gun arrayg)nRepresenting a distance between a detection point of the far-field wavelet and a virtual image of an nth airgun source in the airgun array.
In a preferred scheme, the distance r between the detection point of the far-field wavelet and the kth air gun seismic source in the air gun array is calculated by adopting the following formulak：
Wherein x isk、ykAnd zkRespectively representing the coordinates of the k-th air gun source in the directions of an x-axis, a y-axis and a z-axis in a Cartesian coordinate system.
In a preferred scheme, the distance (r) between the detection point of the far-field wavelet and the virtual image of the kth air gun seismic source in the air gun array is calculated by adopting the following formulag)k：
Wherein x isk、ykAnd zkRespectively representing the x-axis, y-axis and k-th air gun seismic source in a Cartesian coordinate systemCoordinates in the z-axis direction.
An apparatus for determining a far-field wavelet for an airgun array, said apparatus providing an initial near-field wavelet excited by an airgun source in the airgun array and recorded at a receiving point of a near-field detector in a near-field detector set; the device comprises: the device comprises a filtering factor determining module, a target near-field wavelet determining module and a far-field wavelet determining module; wherein,
the filter factor determination module is used for determining a filter factor according to the distance between the air gun seismic source and a near-field detector in the near-field detector group and the distance between a virtual image of the air gun seismic source and the near-field detector in the near-field detector group;
the target near-field wavelet determining module is used for filtering the initial near-field wavelet by using the filtering factor to obtain a target near-field wavelet;
and the far-field wavelet determining module is used for determining the far-field wavelet corresponding to the air gun array according to the target near-field wavelet.
In a preferred embodiment, the target near-field wavelet determining module is configured to convert the initial near-field wavelet from a time domain to a frequency domain to obtain an initial near-field wavelet of the frequency domain, and establish an initial near-field wavelet matrix according to the initial near-field wavelet of the frequency domain, where the initial near-field wavelet matrix is used to represent an initial near-field wavelet set of the frequency domain recorded by the near-field detector set, and the initial near-field wavelet set represented by the initial near-field wavelet matrix is filtered by using the filtering factor to obtain a filtered initial near-field wavelet set of the frequency domain, and the filtered initial near-field wavelet set of the frequency domain is used as the target near-field wavelet.
In a preferred embodiment, the filtering factor determining module is configured to determine the filtering factor by using the following formula:
wherein A represents the filter factor, rmnRepresenting the distance between the receiving point of the mth near-field detector in the near-field detector group and the nth air gun seismic source in the air gun array, (r)g)mnRepresenting the distance between the receiving point of the mth near-field detector in the near-field detector group and the virtual image of the nth air gun seismic source in the air gun array, C representing the speed of sound waves in water, omega representing angular frequency, i representing an imaginary number2＝-1；
The target near-field wavelet determining module is used for obtaining a target near-field wavelet spectrum matrix of the initial near-field wavelet set for representing the filtered frequency domain by adopting the following formula:
X(ω)＝(ATA+λI)-1ATB
wherein X (ω) represents the target near-field wavelet spectral matrix, ATDenotes the transposed matrix of A, λ denotes the stability factor, Nm(ω) represents an initial near-field wavelet of a frequency domain recorded at a receiving point of an mth near-field detector in the near-field detector set, m represents the number of near-field detectors in the near-field detector set, and n represents the number of air gun vibration sources in an air gun array.
According to the technical scheme provided by the embodiment of the application, the method and the device for determining the far-field wavelet of the air gun array in the embodiment of the application can determine the filter factor according to the distance between the air gun seismic source and the near-field detector in the near-field detector group and the distance between the virtual image of the air gun seismic source and the near-field detector in the near-field detector group; filtering the initial near-field wavelet by using the filtering factor to obtain a target near-field wavelet; the far-field wavelet corresponding to the airgun array may be determined based on the target near-field wavelet. Compared with the prior art, the method provided by the embodiment of the application utilizes the actually measured near-field wavelet to simulate the far-field wavelet of the air gun array, the simulation result has scientific basis, and meanwhile, the near-field wavelet is subjected to filtering processing through the determined filtering factor so as to eliminate the virtual reflection of the air gun seismic source and the influence of other air gun seismic sources in the air gun array except the air gun seismic source, and therefore the precision of the determined air gun array far-field wavelet is improved.
Drawings
In order to more clearly illustrate the embodiments of the present application or the technical solutions in the prior art, the drawings needed to be used in the description of the embodiments or the prior art will be briefly introduced below, it is obvious that the drawings in the following description are only some embodiments described in the present application, and for those skilled in the art, other drawings can be obtained according to the drawings without any creative effort.
FIG. 1 is a schematic diagram of the distribution of an air gun array and detector array in an embodiment of the present application;
FIG. 2 is a flow diagram of an embodiment of a method of determining far-field wavelets for an airgun array according to the present application;
FIG. 3 is a schematic diagram of an initial near-field wavelet recorded at a receiving point of each near-field detector in a near-field detector group in an embodiment of the present application;
FIG. 4 is a schematic diagram of far-field wavelets with different vertical azimuths in an embodiment of the present application;
FIG. 5 is a block diagram illustrating the components of one embodiment of the apparatus for determining far-field wavelets for an airgun array of the present application;
FIG. 6 is a block diagram illustrating the components of another embodiment of the apparatus for determining far-field wavelets for an airgun array according to the present application.
Detailed Description
The embodiment of the application provides a method and a device for determining far-field wavelets of an air gun array.
In order to make those skilled in the art better understand the technical solutions in the present application, the technical solutions in the embodiments of the present application will be clearly and completely described below with reference to the drawings in the embodiments of the present application, and it is obvious that the described embodiments are only a part of the embodiments of the present application, 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 application.
The embodiment of the application provides a method for determining far-field wavelets of an air gun array. The method of determining the far-field wavelet of an airgun array provides an initial near-field wavelet that is excited by an airgun source in the airgun array and recorded at a receiving point of a near-field detector in a near-field detector set.
In this embodiment, the airgun array may include a plurality of airgun sources, and the near-field detector may include a plurality of near-field detectors. Wherein the near-field detector may be disposed at a position 1 meter directly above the air gun seismic source. For example, FIG. 2 is a schematic diagram of the distribution of air gun arrays and detector arrays in an embodiment of the present application. N1, N2, N3, N4, N5, N6, N7, N8, N9, N10, N11, N12, N13, N14, N15, N16, and N17 in fig. 1 respectively denote near-field detectors in the near-field detector group. As shown in fig. 1, 29 air gun sources are included in the air gun array, and 17 near field detectors are included in the near field detector group.
In this embodiment, the initial near-field wavelet at the receiving point of each near-field detector in the near-field detector group may be a result of superposition of near-field wavelets excited by each air-gun source in the air-gun array and recorded at the receiving point of the near-field detector.
FIG. 2 is a flow chart of an embodiment of a method of determining far-field wavelets for an airgun array according to the present application. As shown in FIG. 2, the method for determining far-field wavelets of an airgun array comprises the following steps.
Step S101: and determining a filter factor according to the distance between the air gun seismic source and the near-field detector in the near-field detector group and the distance between a virtual image of the air gun seismic source and the near-field detector in the near-field detector group.
In this embodiment, the filter factor may be determined using the following equation:
wherein A represents the filter factor, rmnRepresenting the distance between the receiving point of the mth near-field detector in the near-field detector group and the nth air gun seismic source in the air gun array, (r)g)mnRepresenting the distance between the receiving point of the mth near-field detector in the near-field detector group and the virtual image of the nth air gun seismic source in the air gun array, C representing the speed of sound waves in water, omega representing angular frequency, i representing an imaginary number2＝-1。
Step S102: and carrying out filtering processing on the initial near-field wavelet by using the filtering factor to obtain a target near-field wavelet.
In this embodiment, the filtering processing is performed on the initial near-field wavelet by using the filtering factor to obtain a target near-field wavelet, and specifically, the method may include converting the initial near-field wavelet from a time domain to a frequency domain by using a fourier transform method to obtain an initial near-field wavelet of the frequency domain, and establishing an initial near-field wavelet spectrum matrix according to the initial near-field wavelet of the frequency domain; and the initial near-field wavelet spectrum matrix is used for representing an initial near-field wavelet set of a frequency domain recorded correspondingly by the near-field detector group. The initial near-field wavelet set represented by the initial near-field wavelet spectrum matrix can be filtered by using the filtering factor to obtain a filtered initial near-field wavelet set of the frequency domain, and the filtered initial near-field wavelet set of the frequency domain is used as the target near-field wavelet.
In this embodiment, the following formula may be used to obtain a target near-field wavelet spectrum matrix for characterizing the initial near-field wavelet set of the filtered frequency domain:
X(ω)＝(ATA+λI)-1ATB
wherein X (ω) represents the target near-field wavelet spectral matrix, A represents the filter factor, ATDenotes the transposed matrix of A, λ denotes the stability factor, Nm(ω) represents an initial near-field wavelet of a frequency domain recorded at a receiving point of an mth near-field detector in the near-field detector set, m represents the number of near-field detectors in the near-field detector set, and n represents the number of air gun vibration sources in an air gun array.
In this embodiment, the value of the stability factor may be 10-8。
Step S103: and determining the far-field wavelet corresponding to the air gun array according to the target near-field wavelet.
In this embodiment, the far-field wavelet corresponding to the airgun array may be determined using the following formula:
wherein,the propagation distance of the far field sub-wave corresponding to the air gun array is R, the vertical azimuth angle is theta, and the horizontal azimuth angle is thetaAnd far-field wavelets with angular frequency ω; the resulting far-field wavelet is now in the frequency domain. The vertical azimuth angle theta represents an included angle between the far-field wavelet propagation direction in the Cartesian coordinate system and the z axis of the Cartesian coordinate system, and the horizontal azimuth angle theta isRepresenting the included angle between the projection of the far-field wavelet propagation direction in the plane xoy of the Cartesian coordinate system and the x axis in the Cartesian coordinate system; r isn(r) represents the distance between the detection point of the far-field wavelet and the nth air gun source in the air gun arrayg)nRepresenting a distance between a detection point of the far-field wavelet and a virtual image of an nth airgun source in the airgun array.
In this embodiment, the following formula may be used to calculate the distance r between the detection point of the far-field wavelet and the kth air gun seismic source in the air gun arrayk：
Wherein x isk、ykAnd zkRespectively representing the coordinates of the k-th air gun source in the directions of an x-axis, a y-axis and a z-axis in a Cartesian coordinate system.
In this embodiment, the distance (r) between the detection point of the far-field wavelet and the virtual image of the kth air-gun seismic source in the air-gun array can be calculated by using the following formulag)k：
Wherein x isk、ykAnd zkRespectively representing the coordinates of the k-th air gun source in the directions of an x-axis, a y-axis and a z-axis in a Cartesian coordinate system.
In this embodiment, the propagation distance of the far-field sub-wave corresponding to the air gun array may be R, the vertical azimuth may be θ, and the horizontal azimuth may be θ by an inverse fourier transform methodFar-field wavelet with sum angular frequency omegaAnd converting from the frequency domain to the time domain to obtain the corresponding far-field wavelet of the time domain.
For example, FIG. 3 is a schematic illustration of an initial near-field wavelet recorded at the receiving point of each near-field detector in the near-field detector set of FIG. 1. The abscissa and ordinate in fig. 3 are the number and sampling time of the near-field detector, respectively, and the unit of the sampling time is milliseconds (ms). FIG. 4 is a schematic diagram of far-field wavelets with different vertical azimuth angles in an embodiment of the present application. The abscissa and the ordinate in fig. 3 are the vertical azimuth and the sampling time of the far-field wavelet, respectively, the value range of the vertical azimuth is 0 to 90 °, and the unit of the sampling time is millisecond (ms). FIG. 4 is a schematic diagram of far-field wavelets with different vertical azimuth angles obtained using the method of the present application based on the initial near-field wavelets in FIG. 3.
In the embodiment of the method for determining the far-field wavelet of the air gun array, the filter factor can be determined according to the distance between the air gun seismic source and the near-field detector in the near-field detector group and the distance between the virtual image of the air gun seismic source and the near-field detector in the near-field detector group; filtering the initial near-field wavelet by using the filtering factor to obtain a target near-field wavelet; the far-field wavelet corresponding to the airgun array may be determined based on the target near-field wavelet. Compared with the prior art, the method provided by the embodiment of the application utilizes the actually measured near-field wavelet to simulate the far-field wavelet of the air gun array, the simulation result has scientific basis, and meanwhile, the near-field wavelet is subjected to filtering processing through the determined filtering factor so as to eliminate the virtual reflection of the air gun seismic source and the influence of other air gun seismic sources in the air gun array except the air gun seismic source, and therefore the precision of the determined air gun array far-field wavelet is improved.
FIG. 5 is a block diagram illustrating the components of one embodiment of the apparatus for determining far-field wavelets for an airgun array according to the present application. The means for determining the far field wavelet of the airgun array provides a near field wavelet that is excited by an airgun source in the airgun array and recorded at a receiving point in a near field detector bank. As shown in fig. 5, the apparatus for determining far-field wavelets for an airgun array may comprise: a filter factor determination module 100, a target near-field wavelet determination module 200, and a far-field wavelet determination module 300.
The filter factor determination module 100 may be configured to determine a filter factor according to a distance between the air gun seismic source and a near-field detector in the near-field detector set, and a distance between a virtual image of the air gun seismic source and a near-field detector in the near-field detector set.
The target near-field wavelet determining module 200 may be configured to perform filtering processing on the initial near-field wavelet by using the filtering factor to obtain a target near-field wavelet.
The far-field wavelet determining module 300 may be configured to determine a far-field wavelet corresponding to the air gun array according to the target near-field wavelet.
In this embodiment, the target near-field wavelet determining module 200 may be configured to convert the initial near-field wavelet from a time domain to a frequency domain to obtain an initial near-field wavelet of the frequency domain, and establish an initial near-field wavelet spectrum matrix according to the initial near-field wavelet of the frequency domain, where the initial near-field wavelet spectrum matrix is used to represent an initial near-field wavelet set of the frequency domain recorded by the near-field detector set correspondingly, and perform filtering processing on the initial near-field wavelet set represented by the initial near-field wavelet spectrum matrix by using the filtering factor to obtain a filtered initial near-field wavelet set of the frequency domain, and use the filtered initial near-field wavelet set of the frequency domain as the target near-field wavelet.
In this embodiment, the filter factor determining module 100 may be configured to determine the filter factor by using the following formula:
wherein A represents the filter factor, rmnRepresenting the distance between the receiving point of the mth near-field detector in the near-field detector group and the nth air gun seismic source in the air gun array, (r)g)mnRepresenting the distance between the receiving point of the mth near-field detector in the near-field detector group and the virtual image of the nth air gun seismic source in the air gun array, C representing the speed of sound waves in water, omega representing angular frequency, i representing an imaginary number2＝-1；
The target near-field wavelet determination module 200 may be configured to derive a target near-field wavelet spectral matrix for characterizing the initial near-field wavelet set of the filtered frequency domain using the following equation:
X(ω)＝(ATA+λI)-1ATB
wherein X (ω) represents the target near-field wavelet spectral matrix, ATDenotes the transposed matrix of A, λ denotes the stability factor, Nm(ω) represents an initial near-field wavelet of a frequency domain recorded at a receiving point of an mth near-field detector in the near-field detector set, m represents the number of near-field detectors in the near-field detector set, and n represents the number of air gun vibration sources in an air gun array.
FIG. 6 is a block diagram illustrating the components of another embodiment of the apparatus for determining far-field wavelets for an airgun array according to the present application. As shown in fig. 6, the apparatus for determining far-field wavelets for an airgun array may include a memory, a processor, and a computer program stored on the memory. The memory also stores initial near-field wavelets excited by an air gun seismic source in the air gun array and recorded at the receiving points of near-field detectors in the near-field detector group; the computer program when executed by the processor performs the steps of:
determining a filter factor according to the distance between the air gun seismic source and a near-field detector in the near-field detector group and the distance between a virtual image of the air gun seismic source and the near-field detector in the near-field detector group;
filtering the initial near-field wavelet by using the filtering factor to obtain a target near-field wavelet;
and determining the far-field wavelet corresponding to the air gun array according to the target near-field wavelet.
The device embodiment for determining the far-field wavelet of the air gun array corresponds to the method embodiment for determining the far-field wavelet of the air gun array, the technical scheme of the method embodiment for determining the far-field wavelet of the air gun array can be realized, and the technical effect of the method embodiment can be obtained.
In the 90 s of the 20 th century, improvements in a technology could clearly distinguish between improvements in hardware (e.g., improvements in circuit structures such as diodes, transistors, switches, etc.) and improvements in software (improvements in process flow). However, as technology advances, many of today's process flow improvements have been seen as direct improvements in hardware circuit architecture. Designers almost always obtain the corresponding hardware circuit structure by programming an improved method flow into the hardware circuit. Thus, it cannot be said that an improvement in the process flow cannot be realized by hardware physical modules. For example, programmable logic devices
A Programmable Logic Device (PLD), such as a Field Programmable Gate Array (FPGA), is an integrated circuit whose Logic functions are determined by a user programming the Device. A digital system is "integrated" on a PLD by the designer's own programming without requiring the chip manufacturer to design and fabricate application-specific integrated circuit chips. Furthermore, instead of manually manufacturing the integrated circuit chip, the programming is mostly implemented by "logic compiler" software, which is similar to the software compiler used in the program development and writing, and the original code before compiling is written by a specific programming language, which is called hardware description language
(HDL), but HDL is not only one, but many, such as ABEL (advanced Boolean Expression Language), AHDL (Altera Hardware Description Language), Confluence, CUPL (Central University Programming Language), HDCal, JHDL (Java Hardware Description Language), Lava, Lola, MyHDL, LASPAM, RHDL (Ruby Hardware Description Language), etc., and VHDL (Very-High-speed-Language) and Verilog2 are currently most commonly used. It will also be apparent to those skilled in the art that hardware circuitry that implements the logical method flows can be readily obtained by merely slightly programming the method flows into an integrated circuit using the hardware description languages described above.
Those skilled in the art will also appreciate that, in addition to implementing the controller as pure computer readable program code, the same functionality can be implemented by logically programming method steps such that the controller is in the form of logic gates, switches, application specific integrated circuits, programmable logic controllers, embedded microcontrollers and the like. Such a controller may thus be considered a hardware component, and the means included therein for performing the various functions may also be considered as a structure within the hardware component. Or even means for performing the functions may be regarded as being both a software module for performing the method and a structure within a hardware component.
The apparatuses and modules illustrated in the above embodiments may be implemented by a computer chip or an entity, or by a product with certain functions.
For convenience of description, the above devices are described as being divided into various modules by functions, and are described separately. Of course, the functionality of the various modules may be implemented in the same one or more software and/or hardware implementations as the present application.
From the above description of the embodiments, it is clear to those skilled in the art that the present application can be implemented by software plus necessary general hardware platform. With this understanding in mind, the present solution, or portions thereof that contribute to the prior art, may be embodied in the form of a software product, which in a typical configuration includes one or more processors (CPUs), input/output interfaces, network interfaces, and memory. The computer software product may include instructions for causing a computing device (which may be a personal computer, a server, or a network device, etc.) to perform the methods described in the various embodiments or portions of embodiments of the present application. The computer software product may be stored in a memory, which may include forms of volatile memory in a computer readable medium, Random Access Memory (RAM) and/or non-volatile memory, such as Read Only Memory (ROM) or flash memory (flash RAM). Memory is an example of a computer-readable medium. Computer-readable 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 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, computer readable media does not include transitory computer readable media (transient media), such as modulated data signals and carrier waves.
The embodiments in the present specification are described in a progressive manner, and the same and similar parts among the embodiments are referred to each other, and each embodiment focuses on the differences from the other embodiments. In particular, as for the apparatus embodiment, since it is substantially similar to the method embodiment, the description is relatively simple, and for the relevant points, reference may be made to the partial description of the method embodiment.
The application is operational with numerous general purpose or special purpose computing system environments or configurations. For example: personal computers, server computers, hand-held or portable devices, tablet-type devices, multiprocessor systems, microprocessor-based systems, set top boxes, programmable consumer electronics, network PCs, minicomputers, mainframe computers, distributed computing environments that include any of the above systems or devices, and the like.
The application may be described in the general context of computer-executable instructions, such as program modules, being executed by a computer. Generally, program modules include routines, programs, objects, components, data structures, etc. that perform particular tasks or implement particular abstract data types. The application may also be practiced in distributed computing environments where tasks are performed by remote processing devices that are linked through a communications network. In a distributed computing environment, program modules may be located in both local and remote computer storage media including memory storage devices.
While the present application has been described with examples, those of ordinary skill in the art will appreciate that there are numerous variations and permutations of the present application without departing from the spirit of the application, and it is intended that the appended claims encompass such variations and permutations without departing from the spirit of the application.
Claims (10)
1. A method of determining far field wavelets for an airgun array, the method providing an initial near field wavelet excited by an airgun source in the airgun array and recorded at a receiving point of a near field detector in a near field detector bank; the method comprises the following steps:
determining a filter factor according to the distance between the air gun seismic source and a near-field detector in the near-field detector group and the distance between a virtual image of the air gun seismic source and the near-field detector in the near-field detector group;
filtering the initial near-field wavelet by using the filtering factor to obtain a target near-field wavelet;
and determining the far-field wavelet corresponding to the air gun array according to the target near-field wavelet.
2. The method of claim 1, wherein said filtering said initial near-field wavelet with said filter factor to obtain a target near-field wavelet comprises:
converting the initial near-field wavelet from a time domain to a frequency domain to obtain an initial near-field wavelet of the frequency domain, and establishing an initial near-field wavelet spectrum matrix according to the initial near-field wavelet of the frequency domain; the initial near-field wavelet spectrum matrix is used for representing an initial near-field wavelet set of a frequency domain recorded by the near-field detector group correspondingly;
and filtering the initial near-field wavelet set represented by the initial near-field wavelet spectrum matrix by using the filtering factor to obtain a filtered initial near-field wavelet set of the frequency domain, and taking the filtered initial near-field wavelet set of the frequency domain as the target near-field wavelet.
3. The method of claim 2, wherein the filter factor is determined using the following equation:
wherein A represents the filter factor, rmnRepresenting the distance between the receiving point of the mth near-field detector in the near-field detector group and the nth air gun seismic source in the air gun array, (r)g)mnRepresenting the distance between the receiving point of the mth near-field detector in the near-field detector group and the virtual image of the nth air gun seismic source in the air gun array, C representing the speed of sound waves in water, omega representing angular frequency, i representing an imaginary number2＝-1；
Obtaining a target near-field wavelet spectrum matrix for characterizing the initial near-field wavelet set of the filtered frequency domain using the following formula:
X(ω)＝(ATA+λI)-1ATB
wherein X (ω) represents the target near-field wavelet spectral matrix, ATDenotes the transposed matrix of A, λ denotes the stability factor, Nm(ω) represents an initial near-field wavelet of a frequency domain recorded at a receiving point of an mth near-field detector in the near-field detector set, m represents the number of near-field detectors in the near-field detector set, and n represents the number of air gun vibration sources in an air gun array.
4. The method of claim 3, wherein the stability factor has a value of 10-8。
5. The method of claim 3, wherein the far-field wavelet corresponding to the airgun array is determined using the following equation:
wherein,the propagation distance of the far field sub-wave corresponding to the air gun array is R, the vertical azimuth angle is theta, and the horizontal azimuth angle is thetaAnd a far-field wavelet with an angular frequency omega, the vertical azimuth angle theta representing the propagation direction of the far-field wavelet in a Cartesian coordinate system and the z-axis of the Cartesian coordinate systemAngle of horizontal azimuth ofRepresenting the included angle between the projection of the far-field wavelet propagation direction in the plane xoy of the Cartesian coordinate system and the x axis in the Cartesian coordinate system; r isn(r) represents the distance between the detection point of the far-field wavelet and the nth air gun source in the air gun arrayg)nRepresenting a distance between a detection point of the far-field wavelet and a virtual image of an nth airgun source in the airgun array.
6. The method of claim 5, wherein the distance r between the detection point of the far-field wavelet and the kth air-gun source in the air-gun array is calculated using the following equationk：
Wherein x isk、ykAnd zkRespectively representing the coordinates of the k-th air gun source in the directions of an x-axis, a y-axis and a z-axis in a Cartesian coordinate system.
7. The method of claim 5, wherein the distance (r) between the detection point of the far-field wavelet and the virtual image of the kth air-gun source in the air-gun array is calculated using the following equationg)k：
Wherein x isk、ykAnd zkRespectively representing the coordinates of the k-th air gun source in the directions of an x-axis, a y-axis and a z-axis in a Cartesian coordinate system.
8. An apparatus for determining far field wavelets for an airgun array, said apparatus providing an initial near field wavelet excited by an airgun source in the airgun array and recorded at a receiving point of a near field detector in a near field detector bank; the device comprises: the device comprises a filtering factor determining module, a target near-field wavelet determining module and a far-field wavelet determining module; wherein,
the filter factor determination module is used for determining a filter factor according to the distance between the air gun seismic source and a near-field detector in the near-field detector group and the distance between a virtual image of the air gun seismic source and the near-field detector in the near-field detector group;
the target near-field wavelet determining module is used for filtering the initial near-field wavelet by using the filtering factor to obtain a target near-field wavelet;
and the far-field wavelet determining module is used for determining the far-field wavelet corresponding to the air gun array according to the target near-field wavelet.
9. The apparatus of claim 8, wherein the target near-field wavelet determining module is configured to convert the initial near-field wavelet from a time domain to a frequency domain to obtain an initial near-field wavelet of the frequency domain, and establish an initial near-field wavelet spectrum matrix according to the initial near-field wavelet of the frequency domain, where the initial near-field wavelet spectrum matrix is used to represent an initial near-field wavelet set of the frequency domain recorded by the near-field detector set, and perform filtering processing on the initial near-field wavelet set represented by the initial near-field wavelet spectrum matrix by using the filtering factor to obtain a filtered initial near-field wavelet set of the frequency domain, and use the filtered initial near-field wavelet set of the frequency domain as the target near-field wavelet.
10. The apparatus of claim 9, wherein the filter factor determining module is configured to determine the filter factor using the following equation:
wherein A represents the filter factor, rmnRepresenting the distance between the receiving point of the mth near-field detector in the near-field detector group and the nth air gun seismic source in the air gun array, (r)g)mnRepresenting the distance between the receiving point of the mth near-field detector in the near-field detector group and the virtual image of the nth air gun seismic source in the air gun array, C representing the speed of sound waves in water, omega representing angular frequency, i representing an imaginary number2＝-1；
The target near-field wavelet determining module is used for obtaining a target near-field wavelet spectrum matrix of the initial near-field wavelet set for representing the filtered frequency domain by adopting the following formula:
X(ω)＝(ATA+λI)-1ATB
wherein X (ω) represents the target near-field wavelet spectral matrix, ATDenotes the transposed matrix of A, λ denotes the stability factor, Nm(ω) represents an initial near-field wavelet of a frequency domain recorded at a receiving point of an mth near-field detector in the near-field detector set, m represents the number of near-field detectors in the near-field detector set, and n represents the number of air gun vibration sources in an air gun array.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810104599.1A CN108387926B (en) | 2018-02-02 | 2018-02-02 | Method and device for determining far-field wavelets of air gun array |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810104599.1A CN108387926B (en) | 2018-02-02 | 2018-02-02 | Method and device for determining far-field wavelets of air gun array |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108387926A true CN108387926A (en) | 2018-08-10 |
| CN108387926B CN108387926B (en) | 2020-01-07 |
Family
ID=63074944
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810104599.1A Active CN108387926B (en) | 2018-02-02 | 2018-02-02 | Method and device for determining far-field wavelets of air gun array |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108387926B (en) |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110568479A (en) * | 2019-08-26 | 2019-12-13 | 广州海洋地质调查局 | Method for determining far-field wavelet of marine air gun seismic source |
| CN111538089A (en) * | 2020-05-08 | 2020-08-14 | 中国石油天然气集团有限公司 | Multi-ship multi-source seismic source array synchronous control method and device for marine node exploration |
| CN112114355A (en) * | 2019-06-21 | 2020-12-22 | 中国石油天然气集团有限公司 | Air gun array energy center determination method and device |
| CN112162333A (en) * | 2020-09-16 | 2021-01-01 | 中海油田服务股份有限公司 | Method and device for detecting synchronism of air gun array |
| CN112558181A (en) * | 2019-09-26 | 2021-03-26 | 中国石油天然气集团有限公司 | Sensitivity calibration method and device for near-field detector of marine air gun |
| CN112649875A (en) * | 2019-10-11 | 2021-04-13 | 中国石油天然气集团有限公司 | Method and device for positioning submarine node geophone |
| CN114636971A (en) * | 2022-04-26 | 2022-06-17 | 海南浙江大学研究院 | Hydrophone array data far-field signal separation method and device |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101241192A (en) * | 2007-02-06 | 2008-08-13 | 中国石油集团东方地球物理勘探有限责任公司 | Method for eliminating pneumatic gun near-field wavelet imaginary reaction |
| WO2015044207A1 (en) * | 2013-09-26 | 2015-04-02 | Cgg Services Sa | Systems and methods for far field signature reconstruction using data from near field, mid field, and surface field sensors |
| CN104849747A (en) * | 2015-05-19 | 2015-08-19 | 中国海洋石油总公司 | Air gun array optimization method and device |
| CN107561586A (en) * | 2017-07-31 | 2018-01-09 | 中国石油天然气集团公司 | A kind of method and apparatus of bubble compacting |
-
2018
- 2018-02-02 CN CN201810104599.1A patent/CN108387926B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101241192A (en) * | 2007-02-06 | 2008-08-13 | 中国石油集团东方地球物理勘探有限责任公司 | Method for eliminating pneumatic gun near-field wavelet imaginary reaction |
| WO2015044207A1 (en) * | 2013-09-26 | 2015-04-02 | Cgg Services Sa | Systems and methods for far field signature reconstruction using data from near field, mid field, and surface field sensors |
| CN104849747A (en) * | 2015-05-19 | 2015-08-19 | 中国海洋石油总公司 | Air gun array optimization method and device |
| CN107561586A (en) * | 2017-07-31 | 2018-01-09 | 中国石油天然气集团公司 | A kind of method and apparatus of bubble compacting |
Non-Patent Citations (2)
| Title |
|---|
| SONG JIANGUO ET AL.: "The simulation of far-field wavelets using frequency-domain air-gun array near-field wavelets", 《APPLIED GEOPHYSICS》 * |
| 杨博: "气枪震源机理与子波处理", 《中国优秀硕士学位论文全文数据库 基础科学辑》 * |
Cited By (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112114355A (en) * | 2019-06-21 | 2020-12-22 | 中国石油天然气集团有限公司 | Air gun array energy center determination method and device |
| CN112114355B (en) * | 2019-06-21 | 2024-03-01 | 中国石油天然气集团有限公司 | Air gun array energy center determining method and device |
| CN110568479A (en) * | 2019-08-26 | 2019-12-13 | 广州海洋地质调查局 | Method for determining far-field wavelet of marine air gun seismic source |
| CN110568479B (en) * | 2019-08-26 | 2021-04-06 | 广州海洋地质调查局 | Method for determining far-field wavelet of marine air gun seismic source |
| CN112558181A (en) * | 2019-09-26 | 2021-03-26 | 中国石油天然气集团有限公司 | Sensitivity calibration method and device for near-field detector of marine air gun |
| CN112649875A (en) * | 2019-10-11 | 2021-04-13 | 中国石油天然气集团有限公司 | Method and device for positioning submarine node geophone |
| CN111538089A (en) * | 2020-05-08 | 2020-08-14 | 中国石油天然气集团有限公司 | Multi-ship multi-source seismic source array synchronous control method and device for marine node exploration |
| CN112162333A (en) * | 2020-09-16 | 2021-01-01 | 中海油田服务股份有限公司 | Method and device for detecting synchronism of air gun array |
| CN112162333B (en) * | 2020-09-16 | 2022-10-18 | 中海油田服务股份有限公司 | Method and device for detecting synchronism of air gun array |
| CN114636971A (en) * | 2022-04-26 | 2022-06-17 | 海南浙江大学研究院 | Hydrophone array data far-field signal separation method and device |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108387926B (en) | 2020-01-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108387926B (en) | Method and device for determining far-field wavelets of air gun array | |
| Martí et al. | Numerical relativistic hydrodynamics: Local characteristic approach | |
| Yang et al. | A nearly analytic discrete method for acoustic and elastic wave equations in anisotropic media | |
| Florinski et al. | Galactic cosmic ray transport in the global heliosphere | |
| CN106842304B (en) | A kind of prestack depth migration method and device | |
| CN108983285B (en) | moment tensor-based multi-seismic source wave field simulation method and device | |
| Fabien-Ouellet et al. | Time-domain seismic modeling in viscoelastic media for full waveform inversion on heterogeneous computing platforms with OpenCL | |
| CN107193045B (en) | A kind of seismic data processing technique and device | |
| Zhang et al. | A discontinuous Galerkin method for seismic wave propagation in coupled elastic and poroelastic media | |
| Tu et al. | Dependence of the proton beam drift velocity on the proton core plasma beta in the solar wind | |
| CN107765298A (en) | A kind of method and device for determining near-surface velocity model | |
| CN106054251A (en) | First break picking method and device | |
| CN108897040B (en) | Method and device for determining coupling between land detection and seabed | |
| Oudompheng et al. | A theoretical and experimental comparison of the iterative equivalent source method and the generalized inverse beamforming | |
| Abdalla et al. | The dark sector cosmology | |
| CN107179546A (en) | Method and device for determining velocity field | |
| Qiao et al. | Error analysis and filtering methods for absolute ocean gravity data | |
| CN106291691A (en) | A kind of seismic migration imaging method and device | |
| CN104122584A (en) | Method and device for determining directionality according to seismic data | |
| CN107561588A (en) | Seismic data noise suppression method and device | |
| CN110618463A (en) | Nuclear magnetic resonance data inversion method, device, storage medium and equipment | |
| CN111239832A (en) | 3D3C-VSP imaging processing method, device and equipment | |
| CN111142157B (en) | Method, device and equipment for processing three-dimensional inhomogeneous dielectric elastic wave | |
| Li et al. | A 3D reflection ray-tracing method based on linear traveltime perturbation interpolation | |
| CN108983289B (en) | Method and device for determining seismic wave travel time |
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 |