Method and device for determining mine wastewater leakage channel and electronic equipment

Abstract

The invention provides a method and device for determining a mine waste water leakage channel and electronic equipment, and relates to the technical field of geological prospecting, and the method comprises the steps: constructing a reflected wave target function based on the instantaneous frequency smooth constraint term, the spectrum energy distribution constraint term and the nuclear norm constraint term of all reflected wave sub-signals, and obtaining a reflected wave leakage channel; all the reflected wave sub-signals can be accurately extracted by solving the reflected wave target function, the signal aliasing problem is avoided, and after the reflected wave sub-signals under each frequency are obtained, the reflected wave field data in the seismic shot gather data are constructed. And removing the reflected wave field data from the seismic shot gather data to obtain scattered wave field data, and performing migration imaging on the scattered wave field data to obtain the mine wastewater leakage channel of the target mine wastewater treatment area. According to the method, effective information of the wastewater leakage channel can be extracted from a complex mine seismic wave field, and an efficient and accurate solution is provided for accurate identification of the mine wastewater leakage channel.

Description

Technical Field

[0001] The invention relates to the technical field of geological exploration, and in particular to a method, a device and an electronic device for determining a leakage channel of mine wastewater. Background Art

[0002] Wastewater generated during mining usually contains a large amount of heavy metals, toxic chemicals and other pollutants. These wastewaters seep into groundwater, which not only damages the ecological environment of the mining area, but also threatens the safety of surrounding water resources. Effectively identifying wastewater leakage channels and underground pollution diffusion paths is an important part of mine ecological environment restoration.

[0003] The leakage channel of mine wastewater has the characteristics of small scale and complex structure, and its wave field response is scattered wave. Therefore, the scattered wave in the detection area can be identified with high precision, so as to achieve the purpose of high-resolution imaging of the wastewater leakage channel in the detection area and the underground pollution diffusion path. However, since the energy of scattered waves is generally much lower than the energy of interference wave fields such as reflected waves, if the wave field of the entire detection area is directly offset and imaged, it is very likely that the scattered waves will be submerged by interference signals with stronger energy, resulting in low imaging resolution. Therefore, before imaging the scattered wave signal in the detection area, it is necessary to separate it first to achieve the purpose of high-precision imaging of the leakage channel of mine wastewater. However, when using the traditional signal decomposition method, due to the coupling characteristics between the reflected wave and the scattered wave, the energy aliasing phenomenon occurs between the decomposed subspaces, and the reflected wave and the scattered wave cannot be accurately separated, which affects the accuracy of the separation result. Summary of the invention

[0004] The purpose of the present invention is to provide a method, device and electronic equipment for determining a mine wastewater leakage channel, so as to provide an efficient and accurate solution for the accurate identification of the mine wastewater leakage channel.

[0005] In a first aspect, the present invention provides a method for determining a mine wastewater leakage channel, comprising: acquiring seismic shot data of a target mine wastewater treatment area; constructing a reflection wave target function based on instantaneous frequency smoothing constraints, spectrum energy distribution constraints and nuclear norm constraints of all reflection wave sub-signals; solving the reflection wave target function to obtain reflection wave field data; removing reflection wave field data from the seismic shot data to obtain scattered wave field data; performing offset imaging on the scattered wave field data to obtain a mine wastewater leakage channel in the target mine wastewater treatment area.

[0006] Optionally, the reflection wave objective function is expressed as: in, represents the kth reflected wave sub-signal r k(x, t), δ(x, t) represents the Dirac distribution, x represents the location of the seismic trace, t represents the earthquake recording time, j represents the imaginary unit, K represents the number of reflected wave sub-signals, r(x,t) represents the reflected wave field data, Represents r k The analytical signal corresponding to (x, t), ω k Represents r k The center frequency of (x,t), represents the partial derivative of the earthquake record time, s(x,t) represents the seismic shot data; Represents r k The spectral energy distribution constraint term of (x, t), |ω-ω k | 2 Indicates that the frequency ω deviates from the center frequency ω k degree, F(·) represents the Fourier transform, F[r k (x,t)] 2 =|r k (x,w) 2 , represents the spectrum energy at frequency ω; Represents r k The nuclear norm constraint term of (x,t), σ i Represents r k The i-th singular value of (x,t).

[0007] Optionally, solving the reflected wave objective function to obtain the reflected wave wavefield data includes: constructing a Lagrangian function based on the reflected wave objective function; iteratively updating each reflected wave sub-signal and its center frequency based on the Lagrangian function until the sum of the relative update amounts of the energies of all reflected wave sub-signals is less than a preset threshold, thereby obtaining all reflected wave sub-signals; and combining all reflected wave sub-signals to obtain the reflected wavefield data.

[0008] Optionally, the Lagrangian function is expressed as:

[0009]

[0010] Optionally, the formula for iterative update of the reflected wave sub-signal is: in, represents the kth reflected wave sub-signal after the n+1th iteration, It represents the set of reflected wave sub-signals with sequence numbers less than k after the n+1th iteration. Represents the set of reflected wave sub-signals with sequence numbers greater than or equal to k after the nth round of iteration, Represents the set of central frequencies of all reflected wave sub-signals after the nth iteration.

[0011] Optionally, the formula for iteratively updating the center frequency of the reflected wave sub-signal is: in, represents the center frequency of the kth reflected wave sub-signal after the n+1th iteration, represents the set of all reflected wave sub-signals after the n+1th iteration, It represents the set of center frequencies of reflected wave sub-signals with sequence numbers less than k after the n+1th iteration. Represents the set of center frequencies of reflected wave sub-signals with sequence numbers greater than or equal to k after the nth iteration.

[0012] Optionally, the relative update amount of all reflected wave sub-signal energies is calculated as: in, represents the square of the second norm of the k-th reflected wave sub-signal after the n-th iteration, Represents the square of the energy difference of the k-th reflected wave sub-signal before and after the n+1th iteration.

[0013] In a second aspect, the present invention provides a device for determining a mine wastewater leakage channel, comprising: an acquisition module for acquiring seismic shot gather data of a target mine wastewater treatment area; a construction module for constructing a reflection wave objective function based on instantaneous frequency smoothing constraints, spectrum energy distribution constraints and nuclear norm constraints of all reflection wave sub-signals; a solution module for solving the reflection wave objective function to obtain reflection wave field data; a removal module for removing reflection wave field data from seismic shot gather data to obtain scattered wave field data; and an imaging module for performing offset imaging on the scattered wave field data to obtain a mine wastewater leakage channel in the target mine wastewater treatment area.

[0014] In a third aspect, the present invention provides an electronic device including a memory and a processor, wherein the memory stores a computer program executable on the processor, and when the processor executes the computer program, the method for determining a mine wastewater leakage channel of any one of the aforementioned embodiments is implemented.

[0015] In a fourth aspect, the present invention provides a computer-readable storage medium, which stores computer instructions. When the computer instructions are executed by a processor, the method for determining the mine wastewater leakage channel of any one of the aforementioned embodiments is implemented.

[0016] The present invention provides a method for determining a mine wastewater leakage channel. The method constructs a reflection wave objective function based on the instantaneous frequency smoothing constraint term, the spectrum energy distribution constraint term and the nuclear norm constraint term of all reflection wave sub-signals. By solving the reflection wave objective function, all reflection wave sub-signals can be accurately extracted, avoiding the signal aliasing problem. After obtaining the reflection wave sub-signals at each frequency, the reflection wave field data in the seismic shot data is then constructed. The reflection wave field data is removed from the seismic shot data to obtain the scattered wave field data. By performing offset imaging on the scattered wave field data, the mine wastewater leakage channel of the target mine wastewater treatment area can be obtained. The method can extract effective information of the wastewater leakage channel from the complex mine seismic wave field, and provides an efficient and accurate solution for the accurate identification of the mine wastewater leakage channel. BRIEF DESCRIPTION OF THE DRAWINGS

[0017] In order to more clearly illustrate the specific implementation methods of the present invention or the technical solutions in the prior art, the drawings required for use in the specific implementation methods or the description of the prior art will be briefly introduced below. Obviously, the drawings described below are some implementation methods of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative work.

[0018] Figure 1 A flow chart of a method for determining a mine wastewater leakage channel provided by an embodiment of the present invention;

[0019] Figure 2 A flow chart of solving a reflection wave objective function and obtaining reflection wave field data provided by an embodiment of the present invention;

[0020] Figure 3 A functional module diagram of a device for determining a mine wastewater leakage channel provided by an embodiment of the present invention;

[0021] Figure 4 A schematic diagram of an electronic device provided by an embodiment of the present invention. DETAILED DESCRIPTION

[0022] In order to make the purpose, 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 in conjunction with the drawings in the embodiments of the present invention. Obviously, the described embodiments are part of the embodiments of the present invention, not all of the embodiments. Generally, the components of the embodiments of the present invention described and shown in the drawings here can be arranged and designed in various different configurations.

[0023] Therefore, the following detailed description of the embodiments of the present invention provided in the accompanying drawings is not intended to limit the scope of the invention claimed for protection, but merely represents selected embodiments of the present invention. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0024] Some embodiments of the present invention are described in detail below in conjunction with the accompanying drawings. In the absence of conflict, the following embodiments and features in the embodiments can be combined with each other.

[0025] The method for determining the mine wastewater leakage channel provided by the embodiment of the present invention is a fine detection method for mine wastewater leakage channel and pollution monitoring, which aims to improve the identification and monitoring accuracy of mine wastewater leakage channel and effectively prevent the impact of wastewater pollution on groundwater and surrounding ecological environment. The method provided by the embodiment of the present invention is specifically introduced below.

[0026] Embodiment 1

[0027] Figure 1 A flow chart of a method for determining a mine wastewater leakage channel provided by an embodiment of the present invention, such as Figure 1 As shown, the method specifically comprises the following steps:

[0028] Step S102, obtaining seismic shot gather data of the target mine wastewater treatment area.

[0029] Specifically, the seismic shot data s(x, t) of the target mine wastewater treatment area is composed of the reflected wave field data r(x, t) and the scattered wave field data d(x, t), that is, s(x, t) = r(x, t) + d(x, t). Where x represents the position of the seismic trace, that is, the position of the detector, and t represents the earthquake recording time.

[0030] Step S104: constructing a reflected wave objective function based on the instantaneous frequency smoothing constraint term, the spectrum energy distribution constraint term and the nuclear norm constraint term of all reflected wave sub-signals.

[0031] In the embodiment of the present invention, the reflected wave field data r(x, t) is composed of a plurality of reflected wave sub-signals, that is, K represents the number of reflected wave sub-signals, r k (x, t) represents the kth reflection wave sub-signal, and each reflection wave sub-signal corresponds to a specific center frequency. In order to accurately separate the scattered wave wavefield data from the seismic shot data, the embodiment of the present invention proposes to construct an objective function that can accurately extract all reflection wave sub-signals from the seismic shot data, so as to obtain the reflection wave field data by combining all the reflection wave sub-signals, and finally remove them from the seismic shot data to obtain the scattered wave wavefield data.

[0032] It is known that the instantaneous frequency represents the local frequency of the signal at a certain moment, reflecting the spectral variation characteristics of the signal at each moment. The smoothness can reflect the continuity and gradual change of the instantaneous frequency curve in time, and the smoothness of the instantaneous frequency reflects the continuity of the underground medium. The reflected wave has strong energy, good continuity, and linear characteristics. The scattered wave is the seismic response generated by the underground small-scale discontinuous geological body, carrying the geological information of the underground complex geological structure and lithology differences. Therefore, the smooth instantaneous frequency corresponds to the regular structure of the signal (reflected wave), and the instantaneous frequency with relatively large fluctuations is usually related to the scattered wave. Therefore, when constructing the reflected wave objective function, the instantaneous frequency smoothing constraint term of the reflected wave sub-signal is adopted to optimize the instantaneous frequency smoothness of the reflected wave sub-signal.

[0033] Furthermore, in order to avoid the spectrum overlap between the reflected wave and the scattered wave and improve the scattered wave separation effect, the embodiment of the present invention adds a spectrum energy distribution constraint term of the reflected wave sub-signal in the reflected wave objective function, so as to control the center frequency and bandwidth of each reflected wave sub-signal through the spectrum energy distribution constraint term, constrain the spectrum energy of each reflected wave sub-signal to be concentrated as much as possible near its own center spectrum, thereby ensuring the independence of the spectrum distribution, avoiding the interference of irrelevant spectrum, and enhancing the accuracy of the separation result.

[0034] In addition, it is known that the reflected wave signal has a low-rank characteristic in mathematical representation. Therefore, the embodiment of the present invention also adopts the nuclear norm constraint term of the reflected wave sub-signal in the reflected wave objective function to realize the low-rank constraint on the reflected wave sub-signal. In the embodiment of the present invention, the reflected wave objective function is constructed with the goal of minimizing the sum of the instantaneous frequency smoothing constraint term, the spectrum energy distribution constraint term and the nuclear norm constraint term of all reflected wave sub-signals.

[0035] Step S106, solving the reflected wave objective function to obtain reflected wave wave field data.

[0036] From the introduction above, we can know that by solving the reflection wave objective function, we can obtain all the reflection wave sub-signals, and by combining all the reflection wave sub-signals, we can get the reflection wave field data, that is,

[0037] Step S108, removing the reflected wave field data from the seismic shot gather data to obtain the scattered wave field data.

[0038] Specifically, the scattered wave field data is determined by the following formula: d(x, t) = s(x, t) - r(x, t).

[0039] Step S110 , performing migration imaging on the scattered wave field data to obtain the mine wastewater leakage channel of the target mine wastewater treatment area.

[0040] The embodiment of the present invention does not specifically limit the migration imaging method, and the user can select it according to actual needs. Optionally, by using the Kirchhoff migration algorithm, the scattered wave field data is migrated and imaged to obtain an accurate imaging result of the mine wastewater leakage path, that is, to obtain the mine wastewater leakage channel of the target mine wastewater treatment area.

[0041] The present invention provides a method for determining a mine wastewater leakage channel. The method constructs a reflection wave objective function based on the instantaneous frequency smoothing constraint term, the spectrum energy distribution constraint term and the nuclear norm constraint term of all reflection wave sub-signals. By solving the reflection wave objective function, all reflection wave sub-signals can be accurately extracted, avoiding the signal aliasing problem. After obtaining the reflection wave sub-signals at each frequency, the reflection wave field data in the seismic shot data is then constructed. The reflection wave field data is removed from the seismic shot data to obtain the scattered wave field data. By performing offset imaging on the scattered wave field data, the mine wastewater leakage channel of the target mine wastewater treatment area can be obtained. The method can extract effective information of the wastewater leakage channel from the complex mine seismic wave field, and provides an efficient and accurate solution for the accurate identification of the mine wastewater leakage channel.

[0042] In an optional implementation, the reflected wave objective function is expressed as:

[0043]

[0044] in, represents the kth reflected wave sub-signal r k (x, t), δ(x, t) represents the Dirac distribution, x represents the location of the seismic trace, t represents the earthquake recording time, j represents the imaginary unit, K represents the number of reflected wave sub-signals, r(x,t) represents the reflected wave field data, Represents r k The analytical signal corresponding to (x, t), ω k Represents r k The center frequency of (x,t), It represents the partial derivative of earthquake record time, and s(x,t) represents the seismic shot data.

[0045] In the instantaneous frequency smoothing constraint, Represents r k The analytical signal corresponding to (x, t) is a complex-valued signal, which represents the original signal through the real part and provides phase information related to the original signal through the imaginary part. is the kernel function for constructing the analytical signal, The role of the reflected wave signal r kThe spectrum corresponding to (x, t) is moved to zero frequency (baseband). Using the kernel function and center frequency adjustment The complex-valued signal of the reflected wave sub-signal is constructed to better extract and describe the instantaneous frequency characteristics of the reflected wave sub-signal.

[0046] represents the partial derivative of earthquake record time, Represents r k The derivative of the analytical signal corresponding to (x, t) with respect to the seismic recording time is obtained by Calculate r k The instantaneous frequency of (x, t) to emphasize the smoothness of the instantaneous frequency; It represents the L2 norm, which is used to penalize the drastic change of instantaneous frequency during the iteration process, encourage the instantaneous frequency to change in a smooth direction, and constrain the smoothness of the instantaneous frequency of the reflected wave signal.

[0047] Represents r k The spectral energy distribution constraint term of (x, t), |ω-ω k | 2 Indicates that the frequency ω deviates from the center frequency ω k degree, F(·) represents the Fourier transform, F[r k (x,t)] 2 =|r k (x,w) 2 , represents the spectrum energy at frequency ω.

[0048] Based on the expression of the spectral energy distribution constraint term, it can be seen that by controlling the center frequency and bandwidth of the reflected wave sub-signal component and minimizing the value of the constraint term, the spectral energy of each reflected wave sub-signal can be constrained to be concentrated as close to the center spectrum as possible, avoiding spectral overlap between different sub-signals, ensuring the independence of spectral distribution, and enhancing the accuracy of the separation result.

[0049] Represents r k The nuclear norm constraint term of (x,t), σ i Represents r k The i-th singular value of (x,t). k The nuclear norm of (x,t) is r k The sum of all singular values ​​of (x, t) can achieve the low rank constraint on the reflected wave sub-signal by minimizing the nuclear norm.

[0050] In an alternative embodiment, Figure 2 As shown, the above step S106, solving the reflected wave objective function to obtain the reflected wave field data, specifically includes the following steps:

[0051] Step S1061, constructing a Lagrangian function based on the reflected wave objective function.

[0052] Constructing the Lagrangian function according to the reflected wave objective function is equivalent to converting the constrained optimization problem into an unconstrained optimization problem. In the embodiment of the present invention, the Lagrangian function is expressed as:

[0053]

[0054] Step S1062, iteratively updating each reflected wave sub-signal and its central frequency based on the Lagrangian function until the sum of the relative update amounts of the energies of all reflected wave sub-signals is less than a preset threshold, thereby obtaining all reflected wave sub-signals.

[0055] Step S1063: combine all reflected wave sub-signals to obtain reflected wave field data.

[0056] The embodiment of the present invention needs to control the reflected wave sub-signals and their center frequencies to perform multiple rounds of iterative updates. After all reflected wave sub-signals and their center frequencies are initialized, in each round of iteration, all reflected wave sub-signals are first controlled to be updated in ascending order of serial numbers, and then the center frequencies of all reflected wave sub-signals are controlled to be updated in ascending order of serial numbers.

[0057] In an optional implementation manner, the formula for iteratively updating the reflected wave sub-signal is: in, represents the kth reflected wave sub-signal after the n+1th iteration, It represents the set of reflected wave sub-signals with sequence numbers less than k after the n+1th iteration. Represents the set of reflected wave sub-signals with sequence numbers greater than or equal to k after the nth round of iteration, Represents the set of central frequencies of all reflected wave sub-signals after the nth iteration.

[0058] In the embodiment of the present invention, all reflected wave sub-signals need to be updated in order of the serial numbers from small to large. In any round of updating of all reflected wave sub-signals, the center frequency of each reflected wave sub-signal uses the result of the previous round of iteration. For ease of understanding, an example is given for sequential updating of reflected wave sub-signals. Starting from the n+1th round, the first reflected wave sub-signal is to be updated (that is, when k=1), and the reflected wave signals used are all the results of the nth round of iteration: When the second reflected wave sub-signal is updated (i.e. when k=2), since the first reflected wave sub-signal has been updated, the first reflected wave sub-signal should use the result of the n+1th round of iteration. The remaining reflected wave sub-signals still use the iteration results of the nth round, and so on, to complete the update of all reflected wave sub-signals in the n+1th round. Among them, when all reflected wave sub-signals are updated in the n+1th round, the center frequencies of all reflected wave sub-signals use the iteration results of the nth round.

[0059] In an optional implementation manner, the formula for iteratively updating the center frequency of the reflected wave sub-signal is: in, represents the center frequency of the kth reflected wave sub-signal after the n+1th iteration, represents the set of all reflected wave sub-signals after the n+1th iteration, It represents the set of center frequencies of reflected wave sub-signals with sequence numbers less than k after the n+1th iteration. Represents the set of center frequencies of reflected wave sub-signals with sequence numbers greater than or equal to k after the nth iteration.

[0060] In any round of iteration, after all reflected wave sub-signals have completed the update, the center frequencies of the reflected wave sub-signals are sequentially updated. Since all reflected wave sub-signals have completed the update in this round, when the center frequencies are updated, all reflected wave sub-signals use the updated results. For ease of understanding, an example is given of the sequential update of the center frequencies of the reflected wave sub-signals. In the n+1th round, the center frequency of the first reflected wave sub-signal is to be updated (that is, when k=1), and the center frequencies of the reflected wave signals used are all the iteration results of the nth round: When the center frequency of the second reflected wave sub-signal is updated (that is, when k=2), since the center frequency of the first reflected wave sub-signal has been updated, the center frequency of the first reflected wave sub-signal should use the result of the n+1th round of iteration. The center frequencies of the remaining reflected wave sub-signals still use the iteration results of the nth round, and so on, to complete the update of the center frequencies of all reflected wave sub-signals in the n+1th round. When the center frequencies of all reflected wave sub-signals are updated in the n+1th round, all reflected wave sub-signals use the iteration results of the n+1th round.

[0061] In an optional implementation manner, the formula for the relative update amount of all reflected wave sub-signal energies is: in, It represents the square of the second norm of the k-th reflected wave sub-signal after the n-th iteration, reflecting the energy of the k-th reflected wave sub-signal. The square of the second norm of the energy difference of the kth reflected wave sub-signal before and after the n+1th iteration reflects the magnitude of the energy change of the kth reflected wave sub-signal before and after this iteration. Therefore, the relative update amount E of the reflected wave sub-signal energy can reflect the magnitude of the energy change of the reflected wave sub-signal before and after the iteration update.

[0062] It can be seen from the above description that in the embodiment of the present invention, the iteration termination condition of the reflected wave sub-signal and its center frequency is: E<ε, where E represents the relative update amount of the energy of all reflected wave sub-signals, and ε represents the preset threshold.

[0063] After the iterative update, all the reflected wave sub-signals and the corresponding center frequencies can be obtained. All the reflected wave sub-signals are combined to obtain the reflected wave field data, and then the reflected wave field data is removed from the seismic shot gather data to separate the scattered wave field data. The separated scattered wave field data is migrated and imaged through the Kirchhoff migration algorithm to determine the leakage channel of mine wastewater.

[0064] In summary, the embodiment of the present invention constructs and solves the reflection wave objective function based on the instantaneous frequency smoothing constraint term, spectrum energy distribution constraint term and nuclear norm constraint term of all reflection wave sub-signals, thereby achieving accurate extraction of all reflection wave sub-signals and avoiding the signal aliasing problem. When solving the Lagrangian function, the convergence accuracy is improved through an efficient iterative algorithm, so that the embodiment of the present invention has good robustness, is suitable for the processing of complex non-stationary seismic signals, and provides an efficient and accurate solution for the accurate identification and monitoring of mine wastewater leakage channels and contaminated areas.

[0065] Embodiment 2

[0066] An embodiment of the present invention also provides a device for determining a mine wastewater leakage channel, which is mainly used to execute the method for determining a mine wastewater leakage channel provided in the above-mentioned embodiment 1. The following is a specific introduction to the device for determining a mine wastewater leakage channel provided in the embodiment of the present invention.

[0067] Figure 3 A functional module diagram of a device for determining a mine wastewater leakage channel provided by an embodiment of the present invention, such as Figure 3 As shown, the device mainly includes: an acquisition module 10, a construction module 20, a solution module 30, a rejection module 40, and an imaging module 50, wherein:

[0068] The acquisition module 10 is used to acquire seismic shot gather data of the target mine wastewater treatment area.

[0069] The construction module 20 is used to construct a reflected wave objective function based on the instantaneous frequency smoothing constraint term, the spectrum energy distribution constraint term and the nuclear norm constraint term of all reflected wave sub-signals.

[0070] The solution module 30 is used to solve the reflection wave objective function and obtain the reflection wave field data.

[0071] The elimination module 40 is used to eliminate the reflected wave field data from the seismic shot gather data to obtain the scattered wave field data.

[0072] The imaging module 50 is used to perform migration imaging on the scattered wave field data to obtain the mine wastewater leakage channel in the target mine wastewater treatment area.

[0073] The embodiment of the present invention provides a device for determining a mine wastewater leakage channel. The device constructs a reflection wave target function based on the instantaneous frequency smoothing constraint term, the spectrum energy distribution constraint term and the nuclear norm constraint term of all reflection wave sub-signals. By solving the reflection wave target function, all reflection wave sub-signals can be accurately extracted, avoiding the signal aliasing problem. After obtaining the reflection wave sub-signals at each frequency, the reflection wave field data in the seismic shot data is then constructed. The reflection wave field data is removed from the seismic shot data to obtain the scattered wave field data. By performing offset imaging on the scattered wave field data, the mine wastewater leakage channel of the target mine wastewater treatment area can be obtained. The device can extract effective information of the wastewater leakage channel from the complex mine seismic wave field, and provides an efficient and accurate solution for the accurate identification of the mine wastewater leakage channel.

[0074] Optionally, the reflected wave objective function is expressed as: in, represents the kth reflected wave sub-signal r k (x, t), δ(x, t) represents the Dirac distribution, x represents the location of the seismic trace, t represents the earthquake recording time, j represents the imaginary unit, K represents the number of reflected wave sub-signals, r(x,t) represents the reflected wave field data, Represents r k The analytical signal corresponding to (x, t), ω k Represents r k The center frequency of (x,t), represents the partial derivative of the earthquake record time, s(x,t) represents the seismic shot data; Represents r k The spectral energy distribution constraint term of (x, t), rω-ω k | 2 Indicates that the frequency ω deviates from the center frequency ω k degree, F(·) represents the Fourier transform, F[r k (x,t)] 2 =|r k (x,w) 2, represents the spectrum energy at frequency ω; Represents r k The nuclear norm constraint term of (x,t), σ i Represents r k The i-th singular value of (x,t).

[0075] Optionally, the solution module 30 is specifically used for:

[0076] The Lagrangian function is constructed based on the reflected wave objective function.

[0077] Each reflected wave sub-signal and its central frequency are iteratively updated based on the Lagrangian function until the sum of the relative update amounts of the energies of all reflected wave sub-signals is less than a preset threshold, thereby obtaining all reflected wave sub-signals.

[0078] All reflected wave sub-signals are combined to obtain reflected wave field data.

[0079] Alternatively, the Lagrangian function is expressed as:

[0080]

[0081] Optionally, the formula for iteratively updating the reflected wave sub-signal is: in, represents the kth reflected wave sub-signal after the n+1th iteration, It represents the set of reflected wave sub-signals with sequence numbers less than k after the n+1th iteration. Represents the set of reflected wave sub-signals with sequence numbers greater than or equal to k after the nth round of iteration, Represents the set of central frequencies of all reflected wave sub-signals after the nth iteration.

[0082] Optionally, the formula for iteratively updating the center frequency of the reflected wave sub-signal is: in, represents the center frequency of the kth reflected wave sub-signal after the n+1th iteration, represents the set of all reflected wave sub-signals after the n+1th iteration, It represents the set of center frequencies of reflected wave sub-signals with sequence numbers less than k after the n+1th iteration. Represents the set of center frequencies of reflected wave sub-signals with sequence numbers greater than or equal to k after the nth iteration.

[0083] Optionally, the relative update amount of all reflected wave sub-signal energies is calculated as: in, represents the square of the second norm of the k-th reflected wave sub-signal after the n-th iteration, Represents the square of the energy difference of the k-th reflected wave sub-signal before and after the n+1th iteration.

[0084] Embodiment 3

[0085] See also Figure 4 An embodiment of the present invention provides an electronic device, which includes: a processor 60, a memory 61, a bus 62 and a communication interface 63, the processor 60, the communication interface 63 and the memory 61 are connected through the bus 62; the processor 60 is used to execute an executable module stored in the memory 61, such as a computer program.

[0086] 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 storage. The communication connection between the system network element 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, wide area network, local area network, metropolitan area network, etc. may be used.

[0087] The bus 62 may be an ISA bus, a PCI bus, or an EISA bus. The bus may be divided into an address bus, a data bus, a control bus, and the like. For ease of representation, Figure 4 Only one bidirectional arrow is used in the diagram, but this does not mean that there is only one bus or only one type of bus.

[0088] Among them, the memory 61 is used to store programs, and the processor 60 executes the program after receiving the execution instruction. The method executed by the device defined by the process disclosed in any embodiment of the above-mentioned embodiment of the present invention can be applied to the processor 60 or implemented by the processor 60.

[0089] The processor 60 may be an integrated circuit chip with signal processing capabilities. In the implementation process, each step of the above method can be completed by the hardware integrated logic circuit or software instructions in the processor 60. The above processor 60 can be a general-purpose processor, including a central processing unit (CPU), a network processor (NP), etc.; it 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 gates or transistor logic devices, discrete hardware components. The methods, steps and logic block diagrams disclosed in the embodiments of the present invention can be implemented or executed. The general-purpose processor can be a microprocessor or the processor can also be any conventional processor. The steps of the method disclosed in the embodiment of the present invention can be directly embodied as a hardware decoding processor to execute, or can be executed by a combination of hardware and software modules in the decoding processor. The software module can be located in a mature storage medium in the field such as a random access memory, a flash memory, a read-only memory, a programmable read-only memory, or an electrically erasable programmable memory, a register, etc. The storage medium is located in the memory 61, and the processor 60 reads the information in the memory 61 and completes the steps of the above method in combination with its hardware.

[0090] A computer program product of a method, device and electronic device for determining a mine wastewater leakage channel provided in an embodiment of the present invention includes a computer-readable storage medium storing a non-volatile program code executable by a processor. The instructions included in the program code can be used to execute the method in the previous method embodiment. For specific implementation, please refer to the method embodiment, which will not be repeated here.

[0091] In addition, each functional unit in each embodiment of the present invention may be integrated into one processing unit, or each unit may exist physically separately, or two or more units may be integrated into one unit.

[0092] If the function is implemented in the form of a software functional unit and sold or used as an independent product, it can be stored in a non-volatile computer-readable storage medium that is executable by a processor. Based on this understanding, the technical solution of the present invention, or the part that contributes to the prior art, or the part of the technical solution, can be embodied in the form of a software product. The computer software product is stored in a storage medium and includes several instructions for a computer device (which can be a personal computer, a server, or a network device, etc.) to perform all or part of the steps of the methods of each embodiment of the present invention. The aforementioned storage medium includes: U disk, mobile hard disk, read-only memory (ROM, Read-Only Memory), random access memory (RAM, Random Access Memory), disk or optical disk, etc. Various media that can store program codes.

[0093] It should be noted that similar reference numerals and letters denote similar items in the following drawings, and therefore, once an item is defined in one drawing, further definition and explanation thereof is not required in subsequent drawings.

[0094] In the description of the present invention, it should be noted that the terms "center", "upper", "lower", "left", "right", "vertical", "horizontal", "inside", "outside", etc. indicate positions or positional relationships based on the positions or positional relationships shown in the accompanying drawings, or the positions or positional relationships in which the product of the invention is usually placed when in use. They are only for the convenience of describing the present invention and simplifying the description, and do not indicate or imply that the device or element referred to must have a specific position, be constructed and operated in a specific position, and therefore cannot be understood as limiting the present invention. In addition, the terms "first", "second", "third", etc. are only used to distinguish the description, and cannot be understood as indicating or implying relative importance.

[0095] In addition, the terms "horizontal", "vertical", "overhanging" and the like do not mean that the components are required to be absolutely horizontal or overhanging, but can be slightly tilted. For example, "horizontal" only means that its direction is more horizontal than "vertical", and does not mean that the structure must be completely horizontal, but can be slightly tilted.

[0096] In the description of the present invention, it is also necessary to explain that, unless otherwise clearly specified and limited, the terms "set", "install", "connect", and "connect" should be understood in a broad sense, for example, it can be a fixed connection, a detachable connection, or an integral connection; it can be a mechanical connection or an electrical connection; it can be a direct connection, or it can be indirectly connected through an intermediate medium, or it can be the internal communication of two elements. For ordinary technicians in this field, the specific meanings of the above terms in the present invention can be understood according to specific circumstances.

[0097] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention, rather than to limit it. Although the present invention has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that they can still modify the technical solutions described in the aforementioned embodiments, or replace some or all of the technical features therein with equivalents. However, these modifications or replacements do not cause the essence of the corresponding technical solutions to deviate from the scope of the technical solutions of the embodiments of the present invention.

Claims

1. A method for determining a mine wastewater leakage channel, characterized in that: include: Acquire seismic shot gather data of the target mine wastewater treatment area; Based on the instantaneous frequency smoothing constraint term, spectrum energy distribution constraint term and nuclear norm constraint term of all reflected wave sub-signals, the reflected wave objective function is constructed; Solving the reflected wave objective function to obtain reflected wave wavefield data; Eliminating the reflected wave field data from the seismic shot gather data to obtain scattered wave field data; The scattered wave field data is subjected to migration imaging to obtain the mine wastewater leakage channel of the target mine wastewater treatment area.

2. The method for determining the leakage channel of mine wastewater according to claim 1, characterized in that: The reflected wave objective function is expressed as: in, represents the kth reflected wave sub-signal r k (x, t), δ(x, t) represents the Dirac distribution, x represents the location of the seismic trace, t represents the earthquake recording time, j represents the imaginary unit, K represents the number of reflected wave sub-signals, r(x,t) represents the reflected wave field data, Represents r k The analytical signal corresponding to (x, t), ω k Represents r k The center frequency of (x,t), represents the partial derivative of the earthquake record time, s(x,t) represents the seismic shot data; Represents r k The spectral energy distribution constraint term of (x, t), |ω-ω k 2 Indicates that the frequency ω deviates from the center frequency ω k degree, F(·) represents the Fourier transform, F[r k (x,t)] 2 =|r k (x,w) 2 , represents the spectrum energy at frequency ω; Represents r k The nuclear norm constraint term of (x,t), σ i Represents r k The i-th singular value of (x,t).

3. The method for determining the leakage channel of mine wastewater according to claim 2, characterized in that: Solving the reflected wave objective function to obtain reflected wave field data includes: Constructing a Lagrangian function based on the reflected wave objective function; Iteratively updating each reflected wave sub-signal and its center frequency based on the Lagrangian function until the sum of the relative update amounts of the energies of all reflected wave sub-signals is less than a preset threshold, thereby obtaining all reflected wave sub-signals; All the reflected wave sub-signals are combined to obtain reflected wave field data.

4. The method for determining the leakage channel of mine wastewater according to claim 3, characterized in that: The Lagrangian function is expressed as:

5. The method for determining the leakage channel of mine wastewater according to claim 3, characterized in that: The formula for iterative update of the reflected wave sub-signal is: in, represents the kth reflected wave sub-signal after the n+1th iteration, It represents the set of reflected wave sub-signals with sequence numbers less than k after the n+1th iteration. Represents the set of reflected wave sub-signals with sequence numbers greater than or equal to k after the nth round of iteration, Represents the set of central frequencies of all reflected wave sub-signals after the nth iteration.

6. The method for determining the leakage channel of mine wastewater according to claim 3, characterized in that: The formula for iterative update of the center frequency of the reflected wave sub-signal is: in, represents the center frequency of the kth reflected wave sub-signal after the n+1th iteration, represents the set of all reflected wave sub-signals after the n+1th iteration, It represents the set of center frequencies of reflected wave sub-signals with sequence numbers less than k after the n+1th iteration. Represents the set of center frequencies of reflected wave sub-signals with sequence numbers greater than or equal to k after the nth iteration.

7. The method for determining the leakage channel of mine wastewater according to claim 3, characterized in that: The formula for the relative update of the energy of all reflected wave sub-signals is: in, represents the square of the second norm of the k-th reflected wave sub-signal after the n-th iteration, Represents the square of the energy difference of the k-th reflected wave sub-signal before and after the n+1th iteration.

8. A device for determining a mine wastewater leakage channel, characterized in that: include: An acquisition module, used to acquire seismic shot gather data of a target mine wastewater treatment area; A construction module, used for constructing a reflection wave objective function based on instantaneous frequency smoothing constraints, spectrum energy distribution constraints and nuclear norm constraints of all reflection wave sub-signals; A solution module, used for solving the reflection wave objective function to obtain reflection wave field data; A removal module, used for removing the reflected wave field data from the seismic shot gather data to obtain scattered wave field data; The imaging module is used to perform offset imaging on the scattered wave field data to obtain the mine wastewater leakage channel of the target mine wastewater treatment area.

9. An electronic device comprising a memory and a processor, wherein the memory stores a computer program that can be run on the processor, wherein: When the processor executes the computer program, the method for determining the mine wastewater leakage channel according to any one of claims 1 to 7 is implemented.

10. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores computer instructions, and when the computer instructions are executed by a processor, the method for determining the mine wastewater leakage channel according to any one of claims 1 to 7 is implemented.

Images