CN113608261A - Diffracted wave imaging method and device and electronic equipment - Google Patents
Diffracted wave imaging method and device and electronic equipment Download PDFInfo
- Publication number
- CN113608261A CN113608261A CN202110878253.9A CN202110878253A CN113608261A CN 113608261 A CN113608261 A CN 113608261A CN 202110878253 A CN202110878253 A CN 202110878253A CN 113608261 A CN113608261 A CN 113608261A
- Authority
- CN
- China
- Prior art keywords
- preset
- weight function
- diffracted wave
- amplitude
- order weight
- 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
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. analysis, for interpretation, for correction
- G01V1/282—Application of seismic models, synthetic seismograms
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V8/00—Prospecting or detecting by optical means
- G01V8/10—Detecting, e.g. by using light barriers
Abstract
The invention provides a diffracted wave imaging method, a diffracted wave imaging device and electronic equipment, which relate to the technical field of geological exploration and comprise the following steps: acquiring a dip angle gather of a region to be imaged; determining an amplitude dataset of a target imaging point based on the dip gather; determining a diffracted wave imaging result of the target imaging point by combining a preset double-order weight function and the amplitude data set of the target imaging point; and determining a comprehensive diffracted wave imaging result of the area to be imaged based on diffracted wave imaging results of all imaging points in the area to be imaged. In the method, the preset double-order weight function value corresponding to the diffracted wave amplitude is larger than the first preset threshold, the preset double-order weight function value corresponding to the reflected wave amplitude is smaller than the second preset threshold, and the first preset threshold is larger than the second preset threshold, so that the method can effectively suppress the reflected wave energy, and the diffracted wave imaging is more focused, thereby realizing the high-precision diffracted wave imaging and effectively solving the technical problem of low imaging result accuracy of the diffracted wave imaging method in the prior art.
Description
Technical Field
The invention relates to the technical field of geological exploration, in particular to a diffracted wave imaging method, a diffracted wave imaging device and electronic equipment.
Background
Discontinuous geologic bodies such as fault breakpoints, collapse columns, stratum sharp vanishing points and the like are closely related to oil and gas migration and coal mining safety, and accurate identification of the discontinuous geologic bodies is beneficial to improving oil and gas mining efficiency, reducing mining cost and reducing geological hazards in the coal mining process. The geological discontinuous information is represented as diffracted wave characteristics in the seismic wave field, so that the acquisition of the geological discontinuous information can be realized by utilizing a diffracted wave imaging technology, and the seismic exploration precision is improved.
The traditional diffracted wave imaging method usually depends on formation angle information or needs to provide filter window information, and when the method for calculating the formation angle information or the filter window information has limitations and/or the data quality signal-to-noise ratio is low, the information with higher precision is difficult to obtain, so that the accuracy of the diffracted wave imaging result is low.
Disclosure of Invention
The invention aims to provide a diffracted wave imaging method, a diffracted wave imaging device and electronic equipment, so as to solve the technical problem that the imaging result accuracy is low in the diffracted wave imaging method in the prior art.
In a first aspect, the present invention provides a diffracted wave imaging method, comprising: acquiring a dip angle gather of a region to be imaged; determining an amplitude dataset of a target imaging point based on the dip gather; the target imaging point is any one of all imaging points in the region to be imaged; determining a diffracted wave imaging result of the target imaging point by combining a preset double-order weight function and the amplitude data set of the target imaging point; wherein each amplitude data in the set of amplitude data comprises: diffraction wave amplitude and reflected wave amplitude; a preset double-order weight function value corresponding to the amplitude of the diffracted wave is larger than a first preset threshold value, a preset double-order weight function value corresponding to the amplitude of the reflected wave is smaller than a second preset threshold value, and the first preset threshold value is larger than the second preset threshold value; and determining a comprehensive diffracted wave imaging result of the to-be-imaged area based on the diffracted wave imaging results of all imaging points in the to-be-imaged area.
In an optional embodiment, determining a diffracted wave imaging result of the target imaging point by combining a preset biquadratic weight function and the amplitude data set of the target imaging point includes: processing the amplitude data set by using the preset double-order weight function to obtain a preset double-order weight function value array of the target imaging point; and determining a diffracted wave imaging result of the target imaging point based on the amplitude data set of the target imaging point and the preset double-order weight function value array.
In an optional embodiment, determining a diffracted wave imaging result of the target imaging point based on the amplitude data set of the target imaging point and the preset bivariate weight function value array includes: equation of utilization
Determining a target imaging point m0The imaging result of the diffracted wave; wherein, I (m)0) Representing the target imaging point m0The diffracted wave imaging result of (1), M represents a target imaging point M0Corresponding total number of observed inclinations, xiThe ith amplitude data, w, in the amplitude data set X representing the target imaging pointd(xi) Represents the ith amplitude data xiAnd the corresponding preset double-order weight function value.
In an alternative embodiment, the method further comprises: acquiring a preset zeroth-order weight function and a preset first-order weight function; wherein the predetermined zeroth order weight function is expressed as
The preset first-order weight function is expressed as
Wherein x isiI-th amplitude data in the amplitude data set representing the target imaging point, u represents an amplitude data mean value of the amplitude data set, and epsilon represents a preset minimum value, x'iDenotes xiU' represents the mean value of the absolute values of the first derivatives of the amplitude data in the amplitude data set, and ω represents a preset attenuation threshold; and constructing the preset double-order weight function based on the preset zero-order weight function and the preset first-order weight function.
In an alternative embodiment, the preset biquadratic weight function is expressed as
Wherein, wd(xi) Represents the ith amplitude data xiCorresponding predetermined two-order weight function value, s1Representing a preset reflection suppression threshold, s2Indicating a predetermined diffraction damage threshold.
In a second aspect, the present invention provides a diffracted wave imaging apparatus, comprising: the first acquisition module is used for acquiring a dip angle gather of a region to be imaged; a first determination module to determine an amplitude dataset of a target imaging point based on the dip gather; the target imaging point is any one of all imaging points in the region to be imaged; the second determining module is used for determining a diffracted wave imaging result of the target imaging point by combining a preset double-order weight function and the amplitude data set of the target imaging point; wherein each amplitude data in the set of amplitude data comprises: diffraction wave amplitude and reflected wave amplitude; a preset double-order weight function value corresponding to the amplitude of the diffracted wave is larger than a first preset threshold value, a preset double-order weight function value corresponding to the amplitude of the reflected wave is smaller than a second preset threshold value, and the first preset threshold value is larger than the second preset threshold value; and the third determining module is used for determining a comprehensive diffracted wave imaging result of the to-be-imaged area based on the diffracted wave imaging results of all imaging points in the to-be-imaged area.
In an alternative embodiment, the second determining module comprises: the processing unit is used for processing the amplitude data set by using the preset double-order weight function to obtain a preset double-order weight function value array of the target imaging point; and the determining unit is used for determining the diffracted wave imaging result of the target imaging point based on the amplitude data set of the target imaging point and the preset double-order weight function value array.
In an optional embodiment, the determining unit is specifically configured to: equation of utilization
Determining a diffracted wave imaging result of the target imaging point m 0; wherein, I (m)0) Representing the target imaging point m0Is wound aroundThe result of the radio wave imaging, M represents the target imaging point M0Corresponding total number of observed inclinations, xiThe ith amplitude data, w, in the amplitude data set X representing the target imaging pointd(xi) Represents the ith amplitude data xiAnd the corresponding preset double-order weight function value.
In a third aspect, the present invention provides an electronic device, comprising a memory and a processor, wherein the memory stores a computer program operable on the processor, and the processor executes the computer program to implement the steps of the method according to any of the foregoing embodiments.
In a fourth aspect, the invention provides a computer readable medium having non-volatile program code executable by a processor, the program code causing the processor to perform the method of any of the preceding embodiments.
According to the diffracted wave imaging method, after an inclination angle gather of a region to be imaged is obtained, an amplitude data set of a target imaging point is determined based on the inclination angle gather; the target imaging point is any one of all imaging points in the region to be imaged; then determining a diffracted wave imaging result of the target imaging point by combining a preset double-order weight function and the amplitude data set of the target imaging point; and finally, determining a comprehensive diffracted wave imaging result of the to-be-imaged area based on diffracted wave imaging results of all imaging points in the to-be-imaged area. In the method of the present invention, each amplitude data in the amplitude data set comprises: diffraction wave amplitude and reflected wave amplitude; the preset double-order weight function value corresponding to the diffracted wave amplitude is larger than a first preset threshold, the preset double-order weight function value corresponding to the reflected wave amplitude is smaller than a second preset threshold, and the first preset threshold is larger than the second preset threshold, so that the method can effectively suppress the energy of reflected waves, the diffracted wave imaging is focused, the diffracted wave high-precision imaging is realized, and the technical problem of low imaging result accuracy in the diffracted wave imaging method in the prior art is effectively solved.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described below, and it is obvious that the drawings in the following description are some embodiments of the present invention, and other drawings can be obtained by those skilled in the art without creative efforts.
FIG. 1 is a flowchart of a diffracted wave imaging method according to an embodiment of the present invention;
fig. 2 is a flowchart of determining a diffracted wave imaging result of a target imaging point by combining a preset biquadratic weight function and an amplitude data set of the target imaging point according to an embodiment of the present invention;
FIG. 3 is a functional block diagram of a diffracted wave imaging apparatus according to an embodiment of the present invention;
fig. 4 is a schematic diagram of an electronic device according to an embodiment of the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are some, but not all, embodiments of the present invention. The components of embodiments of the present invention generally described and illustrated in the figures herein may be arranged and designed in a wide variety of different configurations.
Thus, the following detailed description of the embodiments of the present invention, presented in the figures, is not intended to limit the scope of the invention, as claimed, but is merely representative of selected embodiments of the invention. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
Some embodiments of the invention are described in detail below with reference to the accompanying drawings. The embodiments described below and the features of the embodiments can be combined with each other without conflict.
According to the Huygens principle, when an underground geologic body is taken as a seismic source to excite a seismic wave to the earth surface, the continuous geologic body, namely a reflector can form a wave front surface in a plane wave form; the discontinuous geologic body, i.e. the diffraction body, can form spherical wave diffusion. Discontinuous geologic bodies such as fault breakpoints, collapse columns, stratum sharp vanishing points and the like are closely related to oil and gas migration and coal mining safety, and accurate identification of the discontinuous geologic bodies is beneficial to improving oil and gas mining efficiency, reducing mining cost and reducing geological hazards in the coal mining process.
The geological discontinuous information is represented as diffracted wave characteristics in the seismic wave field, so that the acquisition of the geological discontinuous information can be realized by utilizing a diffracted wave imaging technology, and the seismic exploration precision is improved. Diffracted wave energy is weak (small in amplitude) and is often disturbed by strong reflected waves (large in amplitude), so that a key problem in diffracted wave imaging is weak diffracted wave signal extraction in the background of strong reflected waves.
The traditional diffracted wave imaging method usually depends on formation angle information or needs to provide filter window information, and when the method for calculating the formation angle information or the filter window information has limitations and/or the data quality signal-to-noise ratio is low, the information with higher precision is difficult to obtain, so that the accuracy of the diffracted wave imaging result is low. Embodiments of the present invention provide a diffracted wave imaging method to alleviate the above-mentioned technical problems.
Example one
Fig. 1 is a flowchart of a diffracted wave imaging method according to an embodiment of the present invention, and as shown in fig. 1, the method specifically includes the following steps:
and S102, acquiring a dip angle gather of the to-be-imaged area.
Specifically, diffracted wave imaging is performed on an area to be imaged, firstly, seismic shot gather data U (S, r, t) of the area are acquired, wherein S represents shot point information, r represents demodulator probe information, and t represents recording time, then the seismic shot gather data are processed by a kirchhoff pre-stack migration method, and an inclination gather S (theta, m) of the area to be imaged can be obtained, wherein theta represents an observation inclination angle, and m represents imaging point position information. The dip gather is a gather generated in seismic data processing, the method for acquiring the dip gather of the area to be imaged is not specifically limited in the embodiment of the invention, and a user can select the dip gather according to actual requirements.
In the embodiment of the invention, the abscissa of the dip gather is the observation dip, the value range is [ -90 degrees, 90 degrees ], the ordinate thereof is the imaging point position, and the value of each point in the dip gather is the determined observation dip and the amplitude value of the imaging point position, that is, the dip gather represents the amplitude characteristics of the underground imaging point in different angular directions.
Step S104, determining an amplitude data set of the target imaging point based on the dip gather.
After the dip angle gather is obtained, determining an amplitude data set of a target imaging point according to the dip angle gather; wherein, the target imaging point is any one of all imaging points in the region to be imaged, that is, if the target imaging point is m0Then the amplitude data set for the target imaging point should be represented as X ═ S (θ, m)0)＝(x1,x2,x3…xM-1,xM) And X is a one-dimensional array. x is the number ofiThe ith amplitude data in the amplitude data set X representing the target imaging point, M representing the target imaging point M0The corresponding total number of observed dips.
And S106, determining a diffracted wave imaging result of the target imaging point by combining a preset double-order weight function and the amplitude data set of the target imaging point.
The dip gathers are obtained by conventional seismic migration processes that include both reflected and diffracted wave energy, and therefore, to achieve diffracted wave imaging of the target imaging point, the reflected wave energy must be suppressed and the diffracted wave energy highlighted. In the embodiment of the present invention, after determining the amplitude data set of the target imaging point, combining the amplitude data set with a preset dual-order weight function, where each amplitude data in the amplitude data set includes: diffraction wave amplitude and reflected wave amplitude; the preset double-order weight function value corresponding to the diffracted wave amplitude is larger than a first preset threshold, the preset double-order weight function value corresponding to the reflected wave amplitude is smaller than a second preset threshold, the first preset threshold is larger than the second preset threshold, after the preset double-order weight function is processed, diffracted waves with large preset double-order weight function values (weights) are focused, reflected waves with small preset double-order weight function values are suppressed, and then the diffracted wave imaging result of the target imaging point is obtained.
The embodiment of the invention does not specifically limit the value of the first preset threshold, the value of the second preset threshold and the form of the preset double-order weight function, and a user can set the values according to actual requirements as long as the weight value conditions of the diffracted wave and the reflected wave can be met.
And S108, determining a comprehensive diffracted wave imaging result of the to-be-imaged area based on diffracted wave imaging results of all imaging points in the to-be-imaged area.
Therefore, by using the processing method in steps S104 to S106, the diffracted wave imaging results of all imaging points in the region to be imaged can be obtained, and the set of the diffracted wave imaging results of all imaging points is used as the comprehensive diffracted wave imaging result of the region to be imaged.
According to the diffracted wave imaging method, after an inclination angle gather of a region to be imaged is obtained, an amplitude data set of a target imaging point is determined based on the inclination angle gather; the target imaging point is any one of all imaging points in the region to be imaged; then determining a diffracted wave imaging result of the target imaging point by combining a preset double-order weight function and the amplitude data set of the target imaging point; and finally, determining a comprehensive diffracted wave imaging result of the to-be-imaged area based on diffracted wave imaging results of all imaging points in the to-be-imaged area. In the method of the present invention, each amplitude data in the amplitude data set comprises: diffraction wave amplitude and reflected wave amplitude; the preset double-order weight function value corresponding to the diffracted wave amplitude is larger than a first preset threshold, the preset double-order weight function value corresponding to the reflected wave amplitude is smaller than a second preset threshold, and the first preset threshold is larger than the second preset threshold, so that the method can effectively suppress the energy of reflected waves, the diffracted wave imaging is focused, the diffracted wave high-precision imaging is realized, and the technical problem of low imaging result accuracy in the diffracted wave imaging method in the prior art is effectively solved.
The diffracted wave imaging method provided by the embodiment of the present invention is briefly described above, and an alternative implementation of some of the method steps involved therein is described in detail below.
In an optional embodiment, as shown in fig. 2, in the step S106, determining a diffracted wave imaging result of the target imaging point by combining the preset biquadratic weight function and the amplitude data set of the target imaging point, specifically includes the following steps:
step S1061, processing the amplitude data set by using a preset double-order weight function to obtain a preset double-order weight function value array of the target imaging point.
As can be seen from the above description, the preset dual-order weight function has the function of making the preset dual-order weight function value of the reflected wave be a smaller value and the preset dual-order weight function value of the diffracted wave be a larger value, so as to effectively suppress the reflected wave. Therefore, after the amplitude data set of the target imaging point is obtained, each amplitude data in the amplitude data set is input into the preset double-order weight function, the preset double-order weight function value corresponding to the amplitude data is obtained, and then the preset double-order weight function value array of the target imaging point is obtained.
Step S1062, determining a diffracted wave imaging result of the target imaging point based on the amplitude data set of the target imaging point and the preset double-order weight function value array.
After the preset double-order weight function value array of the target imaging point is obtained, each amplitude data is endowed with a weight, the reflected wave is effectively suppressed, and further, the amplitude data set and the preset double-order weight function value array are integrated, so that the diffracted wave imaging result of the target imaging point can be obtained.
In an optional embodiment, the determining the diffracted wave imaging result of the target imaging point based on the amplitude data set of the target imaging point and the preset double-order weight function value array specifically includes the following steps:
equation of utilization
Determining a target imaging point m0The imaging result of the diffracted wave; it is composed ofIn, I (m)0) Representing a target imaging point m0The diffracted wave imaging result of (1), M represents a target imaging point M0Corresponding total number of observed inclinations, xiThe ith amplitude data, w, in the amplitude data set X representing the target imaging pointd(xi) Represents the ith amplitude data xiAnd the corresponding preset double-order weight function value.
According to the above formula, after the preset dual-order weight function value is obtained, the preset dual-order weight function value is multiplied by the amplitude data in the corresponding amplitude data set, all the amplitude data in the amplitude data set are subjected to the same processing, all the product results are superposed, and finally, the superposed result is used as the imaging result of the target imaging point.
The method for performing diffracted wave imaging on a target imaging point is described in detail above, and an alternative embodiment of how to construct a preset dual-order weight function is described in detail below.
In an alternative embodiment, the method of the present invention further comprises the steps of:
step S201, a preset zeroth-order weight function and a preset first-order weight function are obtained.
Wherein the predetermined zeroth order weight function is expressed as
The preset first-order weight function is expressed as
Wherein x isiI-th amplitude data in the amplitude data set representing the target imaging point, u represents an amplitude data mean value of the amplitude data set, and ε represents a preset minimum value, x'iDenotes xiThe first derivative of (a) is,
u' represents the mean value of the absolute values of the first derivatives of the amplitude data in the amplitude data set, ω represents a predetermined attenuation threshold, optionally 1<ω<2。
Step S202, a preset double-order weight function is constructed based on a preset zero-order weight function and a preset first-order weight function.
Specifically, the invention idea of the embodiment of the present invention is to construct a preset double-order weight function to suppress a reflected wave and to retain a diffracted wave, that is, the preset double-order weight function value of the reflected wave is small, and the preset double-order weight function value of the diffracted wave is large. Ideally, the weight function value for the reflected wave is 0, and the weight function value for the diffracted wave is 1.
The reflected wave is limited by snell's law, the illumination angle is limited, the energy distribution in the angular domain is limited by the first fresnel zone, and the energy is focused in a certain angular range, so that the energy of the reflected wave has a larger difference from the overall (energy at all observation angles) average value, therefore, the preset zeroth-order weight function in the embodiment of the invention can be known from the above expression when the amplitude data x is obtainediWhen the value is close to the mean value u of the amplitude data set, a zeroth order weight function value which is close to 1 can be obtained; on the contrary, when the difference between the two is large, the zeroth order weight function value approaching 0 can be obtained.
In order to avoid the condition that the denominator in the division is 0 (unstable in calculation), the preset zeroth-order weight function in the embodiment of the present invention is provided with a preset minimum value epsilon, and generally, a value of the preset minimum value is much smaller than the amplitude data mean value u. The embodiment of the invention does not specifically limit the form of the preset zeroth-order weight function, and a user can set the zero-order weight function according to actual requirements, so that the processing trend of the amplitude data is met.
Because most of the reflected wave energy has a large difference with the overall mean value and the corresponding zeroth order weight function value is small, the reflected wave energy can be suppressed, but at some angles, the reflected wave energy is close to the mean value, and at this time, the corresponding zeroth order weight function value is large, so that the reflected wave energy is difficult to suppress by presetting the zeroth order weight function, residue is caused, and the residual reflected wave can generate obvious offset arcs in imaging, so that certain artifacts in a diffracted wave imaging section are caused, and geological interpretation is misled.
Considering that when the reflected wave energy is close to the overall average value, the corresponding first derivative or change rate is large, the embodiment of the present invention further needs to use a predetermined first order weight function to further suppress the reflected energy that may remain. According to the expression of the preset first-order weight function, when the reflected wave energy is close to the overall average value, namely, the first-order derivative or the change rate of the amplitude data is large, the first-order weight function value which is close to 0 can be obtained; conversely, when the first derivative or rate of change of the amplitude data is small, a first-order weight function value close to 1 can be obtained. The embodiment of the invention does not specifically limit the form of the preset first-order weight function, and a user can set the function according to actual requirements, so that the processing trend of the first-order derivative of the amplitude data is met.
After the preset zero-order weight function and the preset first-order weight function are obtained, a preset double-order weight function can be constructed, the constructed preset double-order weight function meets the condition that the preset double-order weight function value of the reflected wave is small, and the preset double-order weight function value of the diffracted wave is large.
Optionally, the predetermined dual-order weight function is expressed as
Wherein, wd(xi) Represents the ith amplitude data xiCorresponding predetermined two-order weight function value, s1Representing a preset reflection suppression threshold, s2Indicating a predetermined diffraction damage threshold. As can be seen from the above expression, s1The reflection suppression is more obvious when the value is larger, and the value can be set near 0.2 in reference; s2The diffraction wave damage is smaller as the value is smaller, and the value can be set near 0.6 in reference.
In summary, the embodiments of the present invention provide a diffracted wave imaging method, which utilizes a sample (x)i) The difference between the integral mean value (u) and the preset double-order weight function is constructed, the energy of reflected waves is effectively suppressed, and the distribution of the offset artifacts and the noise in an angle domain is also local, so that the method can eliminate the influence of the offset artifacts and the noise on the diffracted wave imaging, the diffracted wave imaging is focused, the high-precision imaging of the diffracted waves is realized, and the underground small-scale discontinuous geologic body is accurately positioned. In addition, the dependence on the formation is avoided in view of the method of the inventionThe inclination angle information and the filter window information, therefore, the application effect and the application prospect of the diffracted wave imaging are effectively improved.
Example two
The embodiments of the present invention further provide a diffracted wave imaging apparatus, which is mainly used to perform the diffracted wave imaging method provided in the first embodiment, and the diffracted wave imaging apparatus provided in the embodiments of the present invention is specifically described below.
Fig. 3 is a functional block diagram of a diffracted wave imaging apparatus according to an embodiment of the present invention, and as shown in fig. 3, the apparatus mainly includes: the first obtaining module 10, the first determining module 20, the second determining module 30, and the third determining module 40, wherein:
the first acquiring module 10 is configured to acquire a dip gather of a region to be imaged.
A first determining module 20 for determining an amplitude dataset of a target imaging point based on a dip gather; the target imaging point is any one of all imaging points in the region to be imaged.
The second determining module 30 is configured to determine a diffracted wave imaging result of the target imaging point by combining the preset bigram weight function and the amplitude data set of the target imaging point; wherein each amplitude data in the amplitude data set comprises: diffraction wave amplitude and reflected wave amplitude; the preset double-order weight function value corresponding to the diffracted wave amplitude is larger than a first preset threshold, the preset double-order weight function value corresponding to the reflected wave amplitude is smaller than a second preset threshold, and the first preset threshold is larger than the second preset threshold.
And the third determining module 40 is configured to determine a comprehensive diffracted wave imaging result of the to-be-imaged area based on the diffracted wave imaging results of all imaging points in the to-be-imaged area.
The diffracted wave imaging device provided by the embodiment of the invention comprises a first acquisition module 10, a second acquisition module, a third acquisition module and a fourth acquisition module, wherein the first acquisition module is used for acquiring an inclination angle gather of an area to be imaged; a first determining module 20 for determining an amplitude dataset of a target imaging point based on a dip gather; the target imaging point is any one of all imaging points in the region to be imaged; the second determining module 30 is configured to determine a diffracted wave imaging result of the target imaging point by combining the preset bigram weight function and the amplitude data set of the target imaging point; and the third determining module 40 is configured to determine a comprehensive diffracted wave imaging result of the to-be-imaged area based on the diffracted wave imaging results of all imaging points in the to-be-imaged area. Each amplitude data in the amplitude data set comprises: diffraction wave amplitude and reflected wave amplitude; in the embodiment of the invention, the preset double-order weight function value corresponding to the diffracted wave amplitude is greater than the first preset threshold, the preset double-order weight function value corresponding to the reflected wave amplitude is less than the second preset threshold, and the first preset threshold is greater than the second preset threshold, so that the device can effectively suppress the reflected wave energy, and the diffracted wave imaging is more focused, thereby realizing the high-precision diffracted wave imaging and effectively solving the technical problem of low imaging result accuracy in the diffracted wave imaging method in the prior art.
Optionally, the second determining module 30 includes:
and the processing unit is used for processing the amplitude data set by using a preset double-order weight function to obtain a preset double-order weight function value array of the target imaging point.
And the determining unit is used for determining the diffracted wave imaging result of the target imaging point based on the amplitude data set of the target imaging point and the preset double-order weight function value array.
Optionally, the determining unit is specifically configured to:
equation of utilization
Determining a target imaging point m0The imaging result of the diffracted wave; wherein, I (m)0) Representing a target imaging point m0The diffracted wave imaging result of (1), M represents a target imaging point M0Corresponding total number of observed inclinations, xiThe ith amplitude data, w, in the amplitude data set X representing the target imaging pointd(xi) Represents the ith amplitude data xiAnd the corresponding preset double-order weight function value.
Optionally, the apparatus further comprises:
the second acquisition module is used for acquiring a preset zeroth-order weight function and a preset first-order weight function; wherein the zero order weight is presetThe function is expressed as
The preset first-order weight function is expressed as
Wherein x isiI-th amplitude data in the amplitude data set representing the target imaging point, u represents an amplitude data mean value of the amplitude data set, and ε represents a preset minimum value, x'iDenotes xiU' represents the mean of the absolute values of the first derivative of the amplitude data in the amplitude data set, and ω represents the preset attenuation threshold.
And the construction module is used for constructing a preset double-order weight function based on the preset zero-order weight function and the preset first-order weight function.
Optionally, the predetermined dual-order weight function is expressed as
Wherein, wd(xi) Represents the ith amplitude data xiCorresponding predetermined two-order weight function value, s1Representing a preset reflection suppression threshold, s2Indicating a predetermined diffraction damage threshold.
EXAMPLE III
Referring to fig. 4, an embodiment of the present invention provides an electronic device, including: a processor 60, a memory 61, a bus 62 and a communication interface 63, wherein the processor 60, the communication interface 63 and the memory 61 are connected through the bus 62; the processor 60 is arranged to execute executable modules, such as computer programs, stored in the memory 61.
The memory 61 may include a high-speed Random Access Memory (RAM) and may also include a non-volatile memory (non-volatile memory), such as at least one disk memory. The communication connection between the network element of the system and at least one other network element is realized through at least one communication interface 63 (which may be wired or wireless), and the internet, a wide area network, a local network, a metropolitan area network, and the like can be used.
The bus 62 may be an ISA bus, PCI bus, EISA bus, or the like. The bus may be divided into an address bus, a data bus, a control bus, etc. For ease of illustration, only one double-headed arrow is shown in FIG. 4, but that does not indicate only one bus or one type of bus.
The memory 61 is used for storing a program, the processor 60 executes the program after receiving an execution instruction, and the method executed by the apparatus defined by the flow process disclosed in any of the foregoing embodiments of the present invention may be applied to the processor 60, or implemented by the processor 60.
The processor 60 may be an integrated circuit chip having signal processing capabilities. In implementation, the steps of the above method may be performed by integrated logic circuits of hardware or instructions in the form of software in the processor 60. The Processor 60 may be a general-purpose Processor, and includes a Central Processing Unit (CPU), a Network Processor (NP), and the like; the device can also be a Digital Signal Processor (DSP), an Application Specific Integrated Circuit (ASIC), a Field-Programmable Gate Array (FPGA), or other Programmable logic devices, discrete Gate or transistor logic devices, discrete hardware components. The various methods, steps and logic blocks disclosed in the embodiments of the present invention may be implemented or performed. A general purpose processor may be a microprocessor or the processor may be any conventional processor or the like. The steps of the method disclosed in connection with the embodiments of the present invention may be directly implemented by a hardware decoding processor, or implemented by a combination of hardware and software modules in the decoding processor. The software module may be located in ram, flash memory, rom, prom, or eprom, registers, etc. storage media as is well known in the art. The storage medium is located in a memory 61, and the processor 60 reads the information in the memory 61 and, in combination with its hardware, performs the steps of the above method.
The diffracted wave imaging method, apparatus, and computer program product of an electronic device provided in the embodiments of the present invention include a computer-readable storage medium storing processor-executable nonvolatile program code, where instructions included in the program code may be used to execute the method described in the foregoing method embodiments, and specific implementation may refer to the method embodiments, and will not be described herein again.
In addition, functional units in the embodiments of the present invention may be integrated into one processing unit, or each unit may exist alone physically, or two or more units are integrated into one unit.
The functions, if implemented in the form of software functional units and sold or used as a stand-alone product, may be stored in a non-volatile computer-readable storage medium executable by a processor. Based on such understanding, the technical solution of the present invention may be embodied in the form of a software product, which is stored in a storage medium and includes instructions for causing a computer device (which may be a personal computer, a server, or a network device) to execute all or part of the steps of the method according to the embodiments of the present invention. And the aforementioned storage medium includes: a U-disk, a removable hard disk, a Read-Only Memory (ROM), a Random Access Memory (RAM), a magnetic disk or an optical disk, and other various media capable of storing program codes.
It should be noted that: like reference numbers and letters refer to like items in the following figures, and thus, once an item is defined in one figure, it need not be further defined and explained in subsequent figures.
In the description of the present invention, it should be noted that the terms "center", "upper", "lower", "left", "right", "vertical", "horizontal", "inner", "outer", etc. indicate orientations or positional relationships based on the orientations or positional relationships shown in the drawings or the orientations or positional relationships that the products of the present invention are conventionally placed in use, and are only used for convenience in describing the present invention and simplifying the description, but do not indicate or imply that the devices or elements referred to must have a specific orientation, be constructed and operated in a specific orientation, and thus, should not be construed as limiting the present invention. Furthermore, the terms "first," "second," "third," and the like are used solely to distinguish one from another and are not to be construed as indicating or implying relative importance.
Furthermore, the terms "horizontal", "vertical", "overhang" and the like do not imply that the components are required to be absolutely horizontal or overhang, but may be slightly inclined. For example, "horizontal" merely means that the direction is more horizontal than "vertical" and does not mean that the structure must be perfectly horizontal, but may be slightly inclined.
In the description of the present invention, it should also be noted that, unless otherwise explicitly specified or limited, the terms "disposed," "mounted," "connected," and "connected" are to be construed broadly and may, for example, be fixedly connected, detachably connected, or integrally connected; can be mechanically or electrically connected; they may be connected directly or indirectly through intervening media, or they may be interconnected between two elements. The specific meanings of the above terms in the present invention can be understood in specific cases to those skilled in the art.
Finally, it should be noted that: the above embodiments are only used to illustrate the technical solution of the present invention, and not to limit the same; while the invention has been described in detail and with reference to the foregoing embodiments, it will be understood by those skilled in the art that: the technical solutions described in the foregoing embodiments may still be modified, or some or all of the technical features may be equivalently replaced; and the modifications or the substitutions do not make the essence of the corresponding technical solutions depart from the scope of the technical solutions of the embodiments of the present invention.
Claims (10)
1. A diffracted wave imaging method, comprising:
acquiring a dip angle gather of a region to be imaged;
determining an amplitude dataset of a target imaging point based on the dip gather; the target imaging point is any one of all imaging points in the region to be imaged;
determining a diffracted wave imaging result of the target imaging point by combining a preset double-order weight function and the amplitude data set of the target imaging point; wherein each amplitude data in the set of amplitude data comprises: diffraction wave amplitude and reflected wave amplitude; a preset double-order weight function value corresponding to the amplitude of the diffracted wave is larger than a first preset threshold value, a preset double-order weight function value corresponding to the amplitude of the reflected wave is smaller than a second preset threshold value, and the first preset threshold value is larger than the second preset threshold value;
and determining a comprehensive diffracted wave imaging result of the to-be-imaged area based on the diffracted wave imaging results of all imaging points in the to-be-imaged area.
2. The method of claim 1, wherein determining the diffracted wave imaging result of the target imaging point in combination with a preset biquadratic weight function and the amplitude dataset of the target imaging point comprises:
processing the amplitude data set by using the preset double-order weight function to obtain a preset double-order weight function value array of the target imaging point;
and determining a diffracted wave imaging result of the target imaging point based on the amplitude data set of the target imaging point and the preset double-order weight function value array.
3. The method of claim 2, wherein determining a diffracted wave imaging result for the target imaging point based on the amplitude dataset for the target imaging point and the preset array of dual order weight function values comprises:
equation of utilization
Determining a target imaging point m0The imaging result of the diffracted wave; wherein, I (m)0) Representing the target imaging point m0The diffracted wave imaging result of (1), M represents a target imaging point M0Corresponding total number of observed inclinations, xiAn ith amplitude number in an amplitude data set X representing the target imaging pointAccording to wd(xi) Represents the ith amplitude data xiAnd the corresponding preset double-order weight function value.
4. The method of claim 1, further comprising:
acquiring a preset zeroth-order weight function and a preset first-order weight function; wherein the predetermined zeroth order weight function is expressed as
The preset first-order weight function is expressed as
Wherein x isiThe ith amplitude data in the amplitude data set representing the target imaging point, u represents the mean of the amplitude data set, epsilon represents a preset minimum value, and xi' represents xiU' represents the mean value of the absolute values of the first derivatives of the amplitude data in the amplitude data set, and ω represents a preset attenuation threshold;
and constructing the preset double-order weight function based on the preset zero-order weight function and the preset first-order weight function.
5. The method of claim 4,
the preset double-order weight function is expressed as
Wherein, wd(xi) Represents the ith amplitude data xiCorresponding predetermined two-order weight function value, s1Representing a preset reflection suppression threshold, s2Indicating a predetermined diffraction damage threshold.
6. A diffracted wave imaging apparatus, comprising:
the first acquisition module is used for acquiring a dip angle gather of a region to be imaged;
a first determination module to determine an amplitude dataset of a target imaging point based on the dip gather; the target imaging point is any one of all imaging points in the region to be imaged;
the second determining module is used for determining a diffracted wave imaging result of the target imaging point by combining a preset double-order weight function and the amplitude data set of the target imaging point; wherein each amplitude data in the set of amplitude data comprises: diffraction wave amplitude and reflected wave amplitude; a preset double-order weight function value corresponding to the amplitude of the diffracted wave is larger than a first preset threshold value, a preset double-order weight function value corresponding to the amplitude of the reflected wave is smaller than a second preset threshold value, and the first preset threshold value is larger than the second preset threshold value;
and the third determining module is used for determining a comprehensive diffracted wave imaging result of the to-be-imaged area based on the diffracted wave imaging results of all imaging points in the to-be-imaged area.
7. The apparatus of claim 6, wherein the second determining module comprises:
the processing unit is used for processing the amplitude data set by using the preset double-order weight function to obtain a preset double-order weight function value array of the target imaging point;
and the determining unit is used for determining the diffracted wave imaging result of the target imaging point based on the amplitude data set of the target imaging point and the preset double-order weight function value array.
8. The apparatus according to claim 7, wherein the determining unit is specifically configured to:
equation of utilization
Determining a target imaging point m0The imaging result of the diffracted wave; wherein, I (m)0) Representing the target imaging point m0The diffracted wave imaging result of (1), M represents a target imaging point M0Corresponding total number of observed inclinations, xiThe ith amplitude data, w, in the amplitude data set X representing the target imaging pointd(xi) Represents the ith amplitude data xiAnd the corresponding preset double-order weight function value.
9. An electronic device comprising a memory, a processor, and a computer program stored on the memory and executable on the processor, wherein the processor implements the steps of the method of any of claims 1 to 5 when executing the computer program.
10. A computer-readable medium having non-volatile program code executable by a processor, the program code causing the processor to perform the method of any of claims 1 to 5.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110878253.9A CN113608261B (en) | 2021-07-30 | 2021-07-30 | Diffracted wave imaging method and device and electronic equipment |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110878253.9A CN113608261B (en) | 2021-07-30 | 2021-07-30 | Diffracted wave imaging method and device and electronic equipment |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113608261A true CN113608261A (en) | 2021-11-05 |
| CN113608261B CN113608261B (en) | 2021-12-28 |
Family
ID=78306381
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110878253.9A Active CN113608261B (en) | 2021-07-30 | 2021-07-30 | Diffracted wave imaging method and device and electronic equipment |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113608261B (en) |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2008111037A2 (en) * | 2007-03-12 | 2008-09-18 | Geomage 2003 Ltd | A method for identifying and analyzing faults/fractures using reflected and diffracted waves |
| CN105607121A (en) * | 2016-02-02 | 2016-05-25 | 中国矿业大学(北京) | Coal collapse column identification method and apparatus |
| CN106772583A (en) * | 2017-01-10 | 2017-05-31 | 中国科学院地质与地球物理研究所 | A kind of earthquake diffracted wave separation method and device |
| CN107255832A (en) * | 2017-07-10 | 2017-10-17 | 中国石油大学(华东) | A kind of inversion method of subsurface structure |
| CN109143343A (en) * | 2018-11-14 | 2019-01-04 | 中国矿业大学(北京) | Diffraction wave imaging method, device and electronic equipment |
| CN110118993A (en) * | 2019-05-22 | 2019-08-13 | 中国矿业大学(北京) | Diffraction wave imaging method and device |
| CN111880219A (en) * | 2020-08-07 | 2020-11-03 | 中国科学院地质与地球物理研究所 | Diffracted wave imaging method and device based on azimuth-dip angle gather |
| CN111929729A (en) * | 2020-08-20 | 2020-11-13 | 中国矿业大学(北京) | Diffracted wave imaging method and device and electronic equipment |
-
2021
- 2021-07-30 CN CN202110878253.9A patent/CN113608261B/en active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2008111037A2 (en) * | 2007-03-12 | 2008-09-18 | Geomage 2003 Ltd | A method for identifying and analyzing faults/fractures using reflected and diffracted waves |
| CN105607121A (en) * | 2016-02-02 | 2016-05-25 | 中国矿业大学(北京) | Coal collapse column identification method and apparatus |
| CN106772583A (en) * | 2017-01-10 | 2017-05-31 | 中国科学院地质与地球物理研究所 | A kind of earthquake diffracted wave separation method and device |
| CN107255832A (en) * | 2017-07-10 | 2017-10-17 | 中国石油大学(华东) | A kind of inversion method of subsurface structure |
| CN109143343A (en) * | 2018-11-14 | 2019-01-04 | 中国矿业大学(北京) | Diffraction wave imaging method, device and electronic equipment |
| CN110118993A (en) * | 2019-05-22 | 2019-08-13 | 中国矿业大学(北京) | Diffraction wave imaging method and device |
| CN111880219A (en) * | 2020-08-07 | 2020-11-03 | 中国科学院地质与地球物理研究所 | Diffracted wave imaging method and device based on azimuth-dip angle gather |
| CN111929729A (en) * | 2020-08-20 | 2020-11-13 | 中国矿业大学(北京) | Diffracted wave imaging method and device and electronic equipment |
Non-Patent Citations (1)
| Title |
|---|
| 李正伟 等: "基于倾角-偏移距域道集的绕射波成像", 《地球物理学报》 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113608261B (en) | 2021-12-28 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US11409020B2 (en) | Method and apparatus for obtaining residual gravity anomaly | |
| US20050222774A1 (en) | Method, device and software package for smoothing a subsurface property | |
| CN111929729B (en) | Diffracted wave imaging method and device and electronic equipment | |
| EP1555632A2 (en) | Signal noise filtering | |
| CN113640872A (en) | Diffracted wave separation method and device and electronic equipment | |
| CN103927737A (en) | SAR image change detecting method based on non-local mean | |
| CN106772592A (en) | The analysis method and device of diffracted wave focus energy | |
| CN112526470A (en) | Method and device for calibrating radar parameters, electronic equipment and storage medium | |
| CN111337981B (en) | Diffracted wave imaging method and device and electronic equipment | |
| CN110297273B (en) | Pre-stack inversion method and device for reservoir | |
| CN112381940A (en) | Processing method and device for generating digital elevation model from point cloud data and terminal equipment | |
| CN113608261B (en) | Diffracted wave imaging method and device and electronic equipment | |
| CN113945987B (en) | Method and device for detecting disease geologic body and electronic equipment | |
| CN112559559A (en) | List similarity calculation method and device, computer equipment and storage medium | |
| CN111929730B (en) | Small-scale geological anomaly detection method and device | |
| Corsaro et al. | Bayesian peak bagging analysis of 19 low-mass low-luminosity red giants observed with Kepler (Corrigendum) | |
| CN115903043B (en) | Diffracted wave separation method and device | |
| WO2017107695A1 (en) | Method and device for sorting news | |
| CN115685320B (en) | Method and device for denoising seismic shot gather data | |
| CN116381786A (en) | Diffraction wave imaging method and device | |
| CN112578449A (en) | Deep reflection seismic profile random medium parameter determination method and device and electronic equipment | |
| CN111640149A (en) | Analysis method for evaluating point cloud modeling error | |
| CN113031073B (en) | Pre-stack depth migration speed determination method, device and equipment | |
| KR101749834B1 (en) | Apparatus and method for calculating design load | |
| CN115808867A (en) | Comprehensive pulsar time determination method and system based on timing residual error |
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 |