Two-dimensional scattered wave imaging method, system and device based on nonlinear bunching filtering
By using a nonlinear clustering filtering method, setting parameter windows and parabolic parameters, calculating similar energy spectra, and generating optimal reflected wave seismic profiles, the problem of difficult scattered wave imaging is solved, more accurate scattered wave imaging is achieved, and a new tool for understanding geological structures is provided.
Patent Information
- Authority / Receiving Office
- CN · China
- Patent Type
- Patents(China)
- Current Assignee / Owner
- CHINA UNIV OF PETROLEUM (EAST CHINA)
- Filing Date
- 2024-01-18
- Publication Date
- 2026-07-21
AI Technical Summary
Existing technologies struggle to effectively image scattered waves, especially since their energy is weak and their propagation direction does not satisfy Snell's law, which increases the difficulty of imaging.
A nonlinear clustering filter-based method is adopted. By setting the parameter window and the parabolic parameter range, the similar energy spectrum is calculated, the optimal parabolic parameters are determined, the optimal reflected wave seismic profile is generated, and finally the reflected wave information is subtracted from the original seismic profile to obtain the scattered wave seismic profile.
It improves the accuracy of scattered wave imaging, ensures the precision of scattered wave imaging, provides a more in-depth tool for understanding underground geological structures, and offers a new perspective for seismology and geological exploration.
Smart Images
Figure CN117890974B_ABST
Abstract
Description
Technical Field
[0001] This invention relates to the field of seismic exploration and development technology, and in particular to a two-dimensional scattered wave imaging method, system and equipment based on nonlinear beam filtering. Background Technology
[0002] Seismic images contain valuable information about geological structures, including faults, fractures, reflection coefficients, and oil and gas conduits. Conventional seismic exploration techniques use artificial sources to generate seismic waves at the surface. As these waves propagate underground, they encounter wave impedance reflection interfaces, producing reflected waves. By deploying geophones on the surface to collect these reflected wave information, and then processing the data in the laboratory, subsurface structural information can be obtained. This information is used to locate oil and gas resources and prepare for subsequent oil and gas development.
[0003] Besides reflected wave information, seismic waves generate scattered waves when they encounter scattering bodies such as faults, cracks, and cavities during propagation. Compared with reflected waves, scattered waves have the following characteristics: 1) Scattered waves have very weak energy, usually only one-tenth of the energy of reflected waves, making it difficult to extract information from them; 2) Scattered waves do not obey Snell's law when propagating at the scattering point, so their propagation direction is across all space, which further increases the difficulty of imaging scattered waves. Summary of the Invention
[0004] The purpose of this invention is to provide a two-dimensional scattered wave imaging method, system, and device based on nonlinear beam filtering, which improves the accuracy of scattered wave imaging.
[0005] To achieve the above objectives, the present invention provides the following solution:
[0006] A two-dimensional scattered wave imaging method based on nonlinear beamforming filtering includes:
[0007] Obtain the raw two-dimensional seismic imaging profile; the raw two-dimensional seismic imaging profile includes multiple actual seismic traces; each actual seismic trace includes multiple points; each point corresponds to one seismic data; the horizontal coordinate of the point corresponds to the location of the seismic data, the vertical coordinate of the point corresponds to the travel time of the seismic data, and the seismic data is energy;
[0008] Set multiple parameter windows and the range of parabolic parameters;
[0009] Any parameter window is designated as the current parameter window, the parameter channel corresponding to the current parameter window is designated as the current channel, any point on the current channel is designated as the scanning target point, and all points within the current parameter window are designated as energy spectrum calculation points; any parabolic parameter within the range of the parabolic parameters is designated as the current parameter;
[0010] Based on the current parameters, the seismic data corresponding to the scanned target point, and the abscissa of all the energy spectrum calculation points, calculate the similar energy spectrum corresponding to the scanned target point under the current parameters;
[0011] Based on all similar energy spectra corresponding to the scanned target point, determine the optimal parabolic parameters;
[0012] Based on the optimal parabolic parameters, the seismic data corresponding to the scanning target point, and the abscissa of all the energy spectrum calculation points, the optimal reflected wave seismic profile is calculated.
[0013] Based on the original two-dimensional seismic imaging profile and the optimal reflected wave seismic profile, the scattered wave seismic profile is determined.
[0014] Optionally, the parameter path corresponding to the current parameter window is determined based on the midpoint of the current parameter window.
[0015] Optionally, the formula for calculating the similarity energy spectrum is:
[0016]
[0017] Wherein, S(t) ip x ip () represents the similarity energy spectrum of the target point being scanned; win est The parameter window size is denoted by ; x is the x-coordinate of a point within the parameter window; u(t) ip +Δt(x,x ip ), x) represents the seismic data corresponding to the point with abscissa x; Δt(x, x) ip Let Δt(x, x) be the travel time difference between the point with x-coordinate and the target point being scanned. ip ) = A(xx ip )+B(xx ip ) 2 A and B are parabolic parameters; t ip x is the ordinate of the target point being scanned; ip The x-coordinate of the target point being scanned.
[0018] Optionally, based on all similar energy spectra corresponding to the scanned target point, the optimal parabolic parameters are determined, including:
[0019] The parabolic parameters that maximize the similarity energy spectrum corresponding to the scanned target point are determined as the optimal parabolic parameters.
[0020] Optionally, based on the original two-dimensional seismic imaging profile and the optimal reflected wave seismic profile, a scattered wave seismic profile is determined, including:
[0021] The difference is obtained by subtracting the original two-dimensional seismic imaging profile from the optimal reflected wave seismic profile.
[0022] The difference is determined as the scattered wave seismic profile.
[0023] A two-dimensional scattered wave imaging system based on nonlinear beam filtering includes:
[0024] The raw profile acquisition module is used to acquire raw two-dimensional seismic imaging profiles; the raw two-dimensional seismic imaging profiles include multiple actual seismic traces; each actual seismic trace includes multiple points; each point corresponds to one seismic data; the horizontal coordinate of the point corresponds to the location of the seismic data, and the vertical coordinate of the point corresponds to the travel time of the seismic data, wherein the seismic data is energy.
[0025] The settings module is used to set multiple parameter windows and the range of parabolic parameters;
[0026] The current parameter module is used to determine any parameter window as the current parameter window, the parameter channel corresponding to the current parameter window as the current channel, any point on the current channel as the scanning target point, and all points within the current parameter window as energy spectrum calculation points; and to determine any parabolic parameter within the range of the parabolic parameters as the current parameter.
[0027] The similar energy spectrum calculation module is used to calculate the similar energy spectrum of the scanning target point under the current parameters based on the current parameters, the seismic data corresponding to the scanning target point and the abscissa of all the energy spectrum calculation points.
[0028] The optimal parabola parameter determination module is used to determine the optimal parabola parameters based on all similar energy spectra corresponding to the scanned target point.
[0029] The optimal reflection wave seismic profile determination module is used to calculate the optimal reflection wave seismic profile based on the optimal parabolic parameters, the seismic data corresponding to the scanning target point, and the abscissa of all the energy spectrum calculation points.
[0030] The scattered wave seismic profile determination module is used to determine the scattered wave seismic profile based on the original two-dimensional seismic imaging profile and the optimal reflected wave seismic profile.
[0031] An apparatus includes a memory and a processor, the memory storing a computer program, and the processor running the computer program to enable the apparatus to perform the two-dimensional scattering wave imaging method based on nonlinear beam filtering as described above.
[0032] Optionally, the memory is a readable storage medium.
[0033] According to specific embodiments provided by the present invention, the present invention discloses the following technical effects:
[0034] This invention discloses a two-dimensional scattered wave imaging method, system, and device based on nonlinear clustering filtering. First, an original two-dimensional seismic imaging profile is acquired; multiple parameter windows and a range of parabolic parameters are set; any parameter window is designated as the current parameter window, the parameter trace corresponding to the current parameter window is designated as the current trace, any point on the current trace is designated as the scanning target point, and all points within the current parameter window are designated as energy spectrum calculation points; any parabolic parameter within the range of parabolic parameters is designated as the current parameter; second, based on the current parameter, the seismic data corresponding to the scanning target point, and the abscissa of all energy spectrum calculation points, the similar energy spectrum corresponding to the scanning target point at the current parameter is calculated; based on all similar energy spectra corresponding to the scanning target point, the optimal parabolic parameters are determined; third, based on the optimal parabolic parameters, the seismic data corresponding to the scanning target point, and the abscissa of all energy spectrum calculation points, the optimal reflected wave seismic profile is calculated; finally, based on the original two-dimensional seismic imaging profile and the optimal reflected wave seismic profile, the scattered wave seismic profile is determined. This invention first images the reflected wave, and then uses the original two-dimensional seismic imaging profile and the optimal reflected wave seismic profile to determine the scattered wave seismic profile. Since the energy of the reflected wave is much stronger than that of the scattered wave, the final scattered wave seismic profile is more accurate, thus improving the accuracy of the scattered wave imaging. Attached Figure Description
[0035] To more clearly illustrate the technical solutions in the embodiments of the present invention or the prior art, the drawings used in the embodiments will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present invention. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.
[0036] Figure 1 This is a schematic diagram of the two-dimensional scattered wave imaging method based on nonlinear beam filtering provided in Embodiment 1 of the present invention;
[0037] Figure 2 This is a schematic diagram illustrating the principle of nonlinear beam filtering.
[0038] Figure 3 This is a schematic diagram of the two-dimensional scattered wave imaging signal extraction technology based on nonlinear beam filtering.
[0039] Figure 4 A schematic diagram for determining the parameters of the channel;
[0040] Figure 5 A schematic diagram for generating the optimal reflected wave superposition seismic profile;
[0041] Figure 6 This is a schematic diagram of a two-dimensional imaging profile;
[0042] Figure 7 This is a schematic diagram of the superimposed gathers obtained by considering the inclination angle and curvature.
[0043] Figure 8 This is a schematic diagram of similar energy spectra obtained by considering tilt angle and curvature.
[0044] Figure 9 This is the original input two-dimensional imaging profile.
[0045] Figure 10 To extract the scattering imaging information map using the nonlinear beamforming filtering method. Detailed Implementation
[0046] The technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of the present invention, and not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of the present invention.
[0047] The purpose of this invention is to provide a two-dimensional scattered wave imaging method, system, and device based on nonlinear beam filtering, which aims to improve the accuracy of scattered wave imaging.
[0048] To make the above-mentioned objects, features and advantages of the present invention more apparent and understandable, the present invention will be further described in detail below with reference to the accompanying drawings and specific embodiments.
[0049] Example 1
[0050] Figure 1 This is a schematic diagram of the two-dimensional scattered wave imaging method based on nonlinear beamforming filtering provided in Embodiment 1 of the present invention. Figure 1 As shown, the two-dimensional scattered wave imaging method based on nonlinear beam filtering in this embodiment includes:
[0051] Step 101: Obtain the original two-dimensional seismic imaging profile.
[0052] The original 2D seismic imaging profile includes multiple actual seismic traces; each actual seismic trace includes multiple points; each point corresponds to one seismic data point. The x-coordinate of the point corresponds to the location of the seismic data, and the y-coordinate corresponds to the travel time of the seismic data. The seismic data is in the form of energy.
[0053] Step 102: Set multiple parameter windows and the range of parabolic parameters.
[0054] Step 103: Determine the target scanning point, energy spectrum calculation point, and current parameters. Specifically:
[0055] Define any parameter window as the current parameter window, define the parameter channel corresponding to the current parameter window as the current channel, define any point on the current channel as the scanning target point, and define all points within the current parameter window as energy spectrum calculation points; define any parabolic parameter within the range of the parabolic parameters as the current parameter.
[0056] As an optional implementation, the parameter path corresponding to the current parameter window is determined based on the midpoint of the current parameter window.
[0057] Step 104: Based on the current parameters, the seismic data corresponding to the target point, and the abscissa of all energy spectrum calculation points, calculate the similar energy spectrum corresponding to the target point under the current parameters.
[0058] As an optional implementation method, the formula for calculating the similarity energy spectrum is:
[0059]
[0060] Wherein, S(t) ip x ip () represents the similarity energy spectrum of the target point being scanned; win ext The parameter window size is denoted by ; x is the x-coordinate of a point within the parameter window; u(t) ip +Δt(x,x ip ), x) represents the seismic data corresponding to the point with abscissa x; Δt(x, x) ip Let Δt(x, x) be the travel time difference between the point with x-coordinate and the target point being scanned. ip ) = A(xx ip )+B(xx ip ) 2 A and B are parabolic parameters; t ip x is the ordinate of the target point being scanned; ip The x-coordinate of the target point being scanned.
[0061] Step 105: Determine the optimal parabolic parameters based on all similar energy spectra corresponding to the target point.
[0062] As an optional implementation, step 105 includes:
[0063] The optimal parabola parameter is determined based on the parabola parameter that maximizes the similarity energy spectrum corresponding to the scanned target point.
[0064] Step 106: Calculate the optimal reflected wave seismic profile based on the optimal parabolic parameters, the seismic data corresponding to the target scanning point, and the abscissa of all energy spectrum calculation points.
[0065] Step 107: Determine the scattered wave seismic profile based on the original two-dimensional seismic imaging profile and the optimal reflected wave seismic profile.
[0066] As an optional implementation, step 107 includes:
[0067] The difference is obtained by subtracting the original two-dimensional seismic imaging profile from the optimal reflection wave seismic profile.
[0068] The difference was determined as the scattered wave seismic profile.
[0069] In fact, the principle of the nonlinear beam filtering method used in Example 1 is as follows: Figure 2 As shown, actual seismic traces are represented by vertical solid lines, parametric scan trajectories by dashed lines, and optimal scan parameter trajectories by solid lines with squares at both ends. By performing nonlinear scans on multiple sets of different parameter scan trajectories, the corresponding optimal scan parameter trajectory can be determined. Superimposing these optimal superimposed trajectories helps extract the energy information of reflected waves, ultimately generating a reflected wave imaging profile. The essence of the nonlinear clustering filter method lies in subtracting the energy of reflected waves from the original seismic record, thereby separating reflected and scattered waves to obtain a scattered wave imaging profile. This profile contains information from the original seismic record that was not occupied by reflected waves, which plays a crucial role in more accurate imaging of subsurface structures. This method provides a deeper understanding of subsurface geological structures, offering new tools and perspectives for research in seismology and geological exploration.
[0070] like Figure 3 As shown, the present invention also provides a two-dimensional scattered wave imaging signal extraction technique based on a nonlinear beamforming filtering method, including:
[0071] 1. It begins with a raw 2D seismic imaging profile obtained through seismic exploration. The raw 2D seismic imaging profile includes multiple actual seismic traces; each actual seismic trace includes multiple points; each point corresponds to a seismic data u(t, x) representing the seismic response of the subsurface structure; the x-coordinate of the point corresponds to the location of the seismic data, and the t-coordinate of the point corresponds to the travel time of the seismic data. The seismic data is in energy form. u(t, x) is the data, representing the seismic response of the subsurface structure.
[0072] 2. To perform parametric scanning, it is necessary to define the window size and the range of scanning tilt and curvature parameters. The selection of these parameters is crucial for subsequent imaging and data analysis.
[0073] win est ∈[win sum ,2win sum (1).
[0074] win op ∈[win sum 1.5win sum (2).
[0075] Among them, win est Indicates the size of the parameter window; win op Indicates the size of the operator window; win sum This indicates the size of the summation window. Given a maximum scan tilt angle θ, the maximum curvature parameter calculated is... The window size for parameter scanning.
[0076] 3. Scan the target point (t) on the parameter track. ip x ip ) Calculate the travel time difference Δt from the start.
[0077] Δt(x 0 x 0 +Δx)=t(x 0 +Δx)-t(x 0 )=AΔx+BΔx 2 (3).
[0078]
[0079]
[0080] Among them, t ip x is the ordinate (during travel) of the i-th point (i.e., the target point) on the p-th parameter path; ip Let t(x) be the x-coordinate (position) of the i-th point in the p-th parameter path; 0 +Δx) is related to the position x 0 The ordinate of the point at a distance Δx (during travel); t(x 0 ) is x 0 The corresponding travel time; Δx is the distance between the scan point and a point on the parameter track; x 0 Let A be the x-coordinate of a point on the parabola; A and B are parabola parameters, where A is the inclination angle and B is the curvature.
[0081] Scan target point (t) ip x ip () represents a point on the parameter path. The vertical line passing through the midpoint of the parameter window is the parameter path. The parameter path is determined as follows: Figure 4 As shown. Figure 4 In the diagram, the dashed lines represent the parameter paths.
[0082] 4. Next, similar energy spectra are calculated along different scan trajectories. Similar energy spectra help to understand the similarity between data at different locations, while stacked gathers are a method for synthesizing seismic trace data.
[0083]
[0084]
[0085]
[0086] Wherein, S(t) ip x ip () represents the similarity energy spectrum of the target point being scanned; win est The parameter window size is denoted by ; x is the x-coordinate of a point within the parameter window; u(t) ip +Δt(x,x ip ), x) represents the seismic data corresponding to the point with abscissa x; Δt(x, x) ip Let Δt(x, x) be the travel time difference between the point with x-coordinate and the target point being scanned. ip ) = A(xx ip )+B(xx ip ) 2 ; This represents the square of the total energy in the parameter window. w(x, x) represents the sum of squares of the energies in the parameter window. ip The weights are based on the parameters and the seismic traces within the parameter calculation range, used to suppress noise and enhance seismic data.
[0087] The scanning trajectory of formula (6) is determined by the range of formulas (3), (4), and (5). When S approaches its maximum value, it indicates that the energy distribution among different components is more uniform, and the energy will not be too concentrated at a certain point, satisfying the energy propagation law of waves; the superimposed gather is all the channels belonging to the parameter window. From the similar energy spectrum S(t ip x ip The optimal tilt angle A and curvature B are selected from the data to obtain the optimal parabolic parameters. The selection of these parameters will affect subsequent data processing and imaging results.
[0088] 5. Using the optimal parabolic parameters, from the parameter trace superposition trace set u(t) p x p (u(t) of all scan target points on the p-th parameter track) ip x ipExtracting the optimal stacking seismic traces from the (structure) process helps generate the best reflection wave stacked seismic profile. The location with the strongest energy on the similarity energy spectrum is the optimal stacking seismic trace. The A and B values corresponding to this optimal point are the picked parameter values. These parameter values can be used to pick the corresponding stacking traces on the stacking trace set. Putting all the optimal stacking seismic traces together gives the final stacked seismic profile. A schematic diagram of generating the best reflection wave stacked seismic profile is shown below. Figure 5 As shown. Figure 5 In the diagram, the dotted line curve represents the curve determined by the optimal parabolic parameters.
[0089] The points in the summation window are calculated according to formula (9), where, It depends on the positional relationship between the reflected wave superimposed seismic trace and the actual seismic trace within the operator window.
[0090]
[0091] Where u(t, x) represents the optimal two-dimensional seismic data; ∏win op For a certain operator window; w(x, x p ) represents the weight; Adding the time difference between the summation channel and the parameter channel to a point on a given channel, and subtracting the time difference between that channel and the parameter channel, yields the time difference between that channel and the summation channel; Δt(x, x p ) for x and x p The time difference between them; for With x p The time difference between them.
[0092] 6. Finally, the final scattered wave imaging profile is obtained by subtracting the optimal stacked profile from the original seismic profile. This step allows for the separation of reflected and scattered waves, leading to a better understanding of subsurface geological structures. The final scattered wave imaging information is the profile. The imaging is based on the energy of the scattered waves. The energy obtained is actually the energy of the reflected waves, as the energy of reflected waves is several times greater than that of scattered waves and is continuous, making it easier to obtain. After obtaining this energy, the scattered wave energy is obtained by removing the reflected wave energy from the entire profile. This research provides new tools and insights for the fields of seismology and geological exploration.
[0093] u sca (t, x) = u(t, x) - u opt (t p x p (10).
[0094] Among them, u sca (t, x) represents the two-dimensional seismic data of the scattered waves; u opt (t p x p() represents two-dimensional seismic data of reflected waves.
[0095] The invention also tested actual land-based two-dimensional data.
[0096] like Figure 6 As shown, a local depth-domain imaging profile is presented, where the vertical solid lines represent the target trace locations in the parametric scan. It is clearly visible that strong scattering imaging information is located at a depth of approximately 6.8 kilometers.
[0097] on the other hand, Figure 7 and Figure 8 (Dip represents tilt angle, Curvature represents curvature, and depth represents depth) This presents superimposed gathers and similar energy spectra obtained using different scanning parameters. By observing these superimposed gathers and similar energy spectra, it can be clearly seen that the location of the reflected wave imaging has relatively concentrated energy, while the location of the scattered wave imaging exhibits energy dispersion characteristics. Figure 7 This represents the superimposed trace set obtained by scanning the target trace under different dip angles and curvatures, that is, the trace set obtained by superimposing all traces in its corresponding parameter region under different dip angles and curvatures. Figure 8 The displayed energy spectrum is similar in range to Figure 7 A consistent and optimal energy spectrum profile. This invention utilizes this characteristic to completely remove reflected wave information while retaining scattered wave information.
[0098] Figure 9 This demonstrates the use of the technique of the present invention to extract data from raw seismic imaging profiles. Figure 10 The results of scattering imaging information extracted using the method of the present invention are shown. (Comparison) Figure 9 and Figure 10 It was found that after applying the method of the present invention, the reflected wave imaging information was effectively removed, leaving only the scattered wave imaging information.
[0099] Example 2
[0100] The two-dimensional scattering wave imaging system based on nonlinear beam filtering in this embodiment includes:
[0101] The raw profile acquisition module is used to acquire raw 2D seismic imaging profiles. The raw 2D seismic imaging profile includes multiple actual seismic traces. Each actual seismic trace includes multiple points. Each point corresponds to one seismic data point. The x-coordinate of the point corresponds to the location of the seismic data point, and the y-coordinate of the point corresponds to the travel time of the seismic data point. The seismic data is in the form of energy.
[0102] The settings module is used to set multiple parameter windows and the range of parabolic parameters.
[0103] The current parameter module is used to define any parameter window as the current parameter window, define the parameter channel corresponding to the current parameter window as the current channel, define any point on the current channel as the scanning target point, define all points within the current parameter window as energy spectrum calculation points, and define any parabolic parameter within the range of the parabolic parameters as the current parameter.
[0104] The similar energy spectrum calculation module is used to calculate the similar energy spectrum of the scanned target point under the current parameters, based on the current parameters, the seismic data corresponding to the scanned target point, and the abscissa of all energy spectrum calculation points.
[0105] The optimal parabola parameter determination module is used to determine the optimal parabola parameters based on all similar energy spectra corresponding to the scanned target point.
[0106] The optimal reflection wave seismic profile determination module is used to calculate the optimal reflection wave seismic profile based on the optimal parabolic parameters, the seismic data corresponding to the scanned target point, and the abscissa of all energy spectrum calculation points.
[0107] The scattered wave seismic profile determination module is used to determine the scattered wave seismic profile based on the original two-dimensional seismic imaging profile and the optimal reflected wave seismic profile.
[0108] Example 3
[0109] An apparatus includes a memory and a processor, the memory storing a computer program and the processor running the computer program to cause the apparatus to perform the two-dimensional scattering wave imaging method based on nonlinear beam filtering as described in Embodiment 1.
[0110] As an optional implementation, the memory is a readable storage medium.
[0111] Advantages of this invention:
[0112] 1. The method of this invention has wide applicability in the field of post-stack seismic imaging, whether for time-domain or depth-domain profiles, and can even be extended to process 3D seismic data. This multi-dimensional application flexibility makes the method a powerful tool in seismology and geological exploration.
[0113] 2. Compared with traditional pre-stack algorithms, the method of this invention has a significant computational efficiency advantage in extracting scattered wave imaging information. This means that this invention can process large-scale seismic data more quickly, thereby improving the speed and efficiency of data processing.
[0114] 3. The algorithm of this invention employs a local scanning strategy, which enables it to handle data volumes of various sizes and facilitates parallel computation. This means that large seismic datasets can be easily processed without concern for computational complexity, thus better meeting the needs of modern earthquake research.
[0115] The various embodiments in this specification are described in a progressive manner, with each embodiment focusing on its differences from other embodiments. Similar or identical parts between embodiments can be referred to interchangeably. For the systems disclosed in the embodiments, since they correspond to the methods disclosed in the embodiments, the descriptions are relatively simple; relevant parts can be referred to the method section.
[0116] This document uses specific examples to illustrate the principles and implementation methods of the present invention. The descriptions of the above embodiments are only for the purpose of helping to understand the method and core ideas of the present invention. Furthermore, those skilled in the art will recognize that, based on the ideas of the present invention, there will be changes in the specific implementation methods and application scope. Therefore, the content of this specification should not be construed as a limitation of the present invention.
Claims
1. A two-dimensional scattered wave imaging method based on nonlinear beamforming filtering, characterized in that, The method includes: Obtain the raw two-dimensional seismic imaging profile; the raw two-dimensional seismic imaging profile includes multiple actual seismic traces; each actual seismic trace includes multiple points; each point corresponds to one seismic data; the horizontal coordinate of the point corresponds to the location of the seismic data, the vertical coordinate of the point corresponds to the travel time of the seismic data, and the seismic data is energy; Set multiple parameter windows and the range of parabolic parameters; Any parameter window is designated as the current parameter window, the parameter channel corresponding to the current parameter window is designated as the current channel, any point on the current channel is designated as the scanning target point, and all points within the current parameter window are designated as energy spectrum calculation points; any parabolic parameter within the range of the parabolic parameters is designated as the current parameter; Based on the current parameters, the seismic data corresponding to the scanned target point, and the abscissa of all the energy spectrum calculation points, calculate the similar energy spectrum corresponding to the scanned target point under the current parameters; Based on all similar energy spectra corresponding to the scanned target point, determine the optimal parabolic parameters; Based on the optimal parabolic parameters, the seismic data corresponding to the scanning target point, and the abscissa of all the energy spectrum calculation points, the optimal reflected wave seismic profile is calculated. Based on the original two-dimensional seismic imaging profile and the optimal reflected wave seismic profile, the scattered wave seismic profile is determined.
2. The two-dimensional scattered wave imaging method based on nonlinear beamforming filtering according to claim 1, characterized in that, The parameter path corresponding to the current parameter window is determined based on the midpoint of the current parameter window.
3. The two-dimensional scattered wave imaging method based on nonlinear beamforming filtering according to claim 1, characterized in that, The formula for calculating similar energy spectra is: Wherein, S(t) ip ,x ip () represents the similarity energy spectrum of the target point being scanned; win est The parameter window size is denoted by ; x is the x-coordinate of a point within the parameter window; u(t) ip +Δt(x,x ip ),x) represents the seismic data corresponding to the point with x-coordinate; Δt(x,x) ip Let Δt(x, x) be the travel time difference between the point with x-coordinate and the target point being scanned. ip ) = A(xx ip )+B(xx ip ) 2 A and B are parabolic parameters; t ip x is the ordinate of the target point being scanned; ip The x-coordinate of the target point being scanned.
4. The two-dimensional scattered wave imaging method based on nonlinear beamforming filtering according to claim 1, characterized in that, Based on all similar energy spectra corresponding to the scanned target point, the optimal parabolic parameters are determined, including: The parabolic parameters that maximize the similarity energy spectrum corresponding to the scanned target point are determined as the optimal parabolic parameters.
5. The two-dimensional scattered wave imaging method based on nonlinear beamforming filtering according to claim 1, characterized in that, Based on the original two-dimensional seismic imaging profile and the optimal reflected wave seismic profile, the scattered wave seismic profile is determined, including: The difference is obtained by subtracting the original two-dimensional seismic imaging profile from the optimal reflected wave seismic profile. The difference is determined as the scattered wave seismic profile.
6. A two-dimensional scattered wave imaging system based on nonlinear beamforming filtering, characterized in that, The system includes: The raw profile acquisition module is used to acquire raw two-dimensional seismic imaging profiles; the raw two-dimensional seismic imaging profiles include multiple actual seismic traces; each actual seismic trace includes multiple points; each point corresponds to one seismic data; the horizontal coordinate of the point corresponds to the location of the seismic data, and the vertical coordinate of the point corresponds to the travel time of the seismic data, wherein the seismic data is energy. The settings module is used to set multiple parameter windows and the range of parabolic parameters; The current parameter module is used to determine any parameter window as the current parameter window, the parameter channel corresponding to the current parameter window as the current channel, any point on the current channel as the scanning target point, and all points within the current parameter window as energy spectrum calculation points; and to determine any parabolic parameter within the range of the parabolic parameters as the current parameter. The similar energy spectrum calculation module is used to calculate the similar energy spectrum of the scanning target point under the current parameters based on the current parameters, the seismic data corresponding to the scanning target point and the abscissa of all the energy spectrum calculation points. The optimal parabola parameter determination module is used to determine the optimal parabola parameters based on all similar energy spectra corresponding to the scanned target point. The optimal reflection wave seismic profile determination module is used to calculate the optimal reflection wave seismic profile based on the optimal parabolic parameters, the seismic data corresponding to the scanning target point, and the abscissa of all the energy spectrum calculation points. The scattered wave seismic profile determination module is used to determine the scattered wave seismic profile based on the original two-dimensional seismic imaging profile and the optimal reflected wave seismic profile.
7. A device, characterized in that, The device includes a memory and a processor, the memory being used to store a computer program, and the processor running the computer program to enable the device to perform the two-dimensional scattering wave imaging method based on nonlinear beam filtering as described in any one of claims 1 to 5.
8. The device according to claim 7, characterized in that, The memory is a readable storage medium.