Method, device and electronic equipment for determining leakage channels of mine wastewater

Abstract

The present invention provides a method, device, and electronic device for determining mine wastewater leakage channels, relating to the technical field of geological exploration. The method constructs a reflection wave objective function based on the instantaneous frequency smoothing constraint term, the spectral energy distribution constraint term, and the nuclear norm constraint term of all reflected wave sub-signals. By solving the reflection wave objective function, all reflected wave sub-signals can be accurately extracted, avoiding the signal aliasing problem. After obtaining the reflected wave sub-signals at each frequency, the reflected wave field data in the seismic shot data is then constructed. The reflected 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 in the target mine wastewater treatment area can be obtained. This method can extract effective information about wastewater leakage channels from complex mine seismic wave fields, providing an efficient and accurate solution for the precise identification of mine wastewater leakage channels.

Description

Technical Field

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

[0002] Wastewater generated during mining operations often contains large amounts of heavy metals, toxic chemicals, and other pollutants. This wastewater seeps into groundwater, damaging the mining area's ecological environment and threatening surrounding water resources. Effectively identifying wastewater leakage channels and underground pollution diffusion pathways is a crucial step in mine ecological restoration.

[0003] Mine wastewater leakage channels are characterized by their small size and complex structure, and their wavefield response is scattered waves. Therefore, by accurately identifying the scattered waves in the detection area, high-resolution imaging of the wastewater leakage channels and underground pollution diffusion paths can be achieved. However, since the energy of scattered waves is generally much lower than the energy of interfering wavefields such as reflected waves, if the wavefield of the entire detection area is directly offset and imaged, the scattered waves are likely to be overwhelmed by the more energetic interference signals, resulting in low imaging resolution. Therefore, before imaging the scattered wave signals in the detection area, they must be separated to achieve high-precision imaging of the mine wastewater leakage channels. However, when using traditional signal decomposition methods, the coupling between the reflected and scattered waves leads to energy aliasing between the decomposed subspaces, making it impossible to accurately separate the reflected and scattered waves, thus affecting the accuracy of the separation results. Summary of the Invention

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

[0005] In a first aspect, the present invention provides a method for determining a mine wastewater leakage channel, comprising: obtaining seismic shot gather 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 reflected wave sub-signals; solving the reflection wave target function to obtain reflection wave field data; removing the reflection wave field data from the seismic shot gather data to obtain scattered wave field data; and 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 kThe instantaneous frequency smoothing constraint term of (x, t), δ(x, t) represents the Dirac distribution, x represents the position 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 gather 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 relative updated amounts of energy 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, Represents the set of reflected wave sub-signals with sequence numbers less than k after the n+1th round of 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 the 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 follows: in, represents the square of the second norm of the k-th reflected wave sub-signal after the n-th iteration, It represents the square of the energy difference of the k-th reflected wave 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 target 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 target function to obtain reflection wave field data; a removal module for removing reflection wave field data from the 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 comprising a memory and a processor, wherein the memory stores a computer program that can be run on the processor, and when the processor executes the computer program, it implements the method for determining the mine wastewater leakage channel of any one of the aforementioned embodiments.

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

[0016] The present invention provides a method for determining mine wastewater leakage channels. This method constructs a reflection wave objective function based on the instantaneous frequency smoothing constraint term, spectral energy distribution constraint term, and 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-signal 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 in the target mine wastewater treatment area can be obtained. This method can extract effective information about wastewater leakage channels from complex mine seismic wave fields, providing an efficient and accurate solution for the precise identification of mine wastewater leakage channels. BRIEF DESCRIPTION OF THE DRAWINGS

[0017] In order to more clearly illustrate the specific embodiments of the present invention or the technical solutions in the prior art, the following briefly introduces the drawings required for use in the specific embodiments or the description of the prior art. Obviously, the drawings described below are some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying any 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] To make the objectives, technical solutions, and advantages of the embodiments of the present invention more clear, the technical solutions of the embodiments of the present invention will be clearly and completely described below in conjunction with the accompanying drawings of the embodiments of the present invention. Obviously, the described embodiments are only 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 herein 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 as claimed, but rather merely represents selected embodiments of the present invention. All other embodiments derived by persons of ordinary skill in the art based on the embodiments of the present invention without creative effort shall fall within the scope of protection of the present invention.

[0024] The following embodiments of the present invention are described in detail with reference to the accompanying drawings. In the absence of conflict, the following embodiments and features in the embodiments may be combined with each other.

[0025] The method for identifying mine wastewater leakage channels provided by the embodiments of the present invention is a sophisticated detection method for mine wastewater leakage channels and pollution monitoring. It aims to improve the accuracy of identifying and monitoring mine wastewater leakage channels, effectively preventing the impact of wastewater pollution on groundwater and the surrounding ecological environment. The following describes the method provided by the embodiments of the present invention in detail.

[0026] Example 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 is shown in FIG. Figure 1 As shown, the method specifically includes the following steps:

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

[0029] Specifically, the seismic shot gather data s(x, t) for the target mine wastewater treatment area consists of the reflected wave field data r(x, t) and the scattered wave field data d(x, t), i.e., s(x, t) = r(x, t) + d(x, t). Here, x represents the location of the seismic trace, i.e., the location of the geophone, and t represents the time of the earthquake recording.

[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. To accurately separate scattered wave wavefield data from seismic shot gather data, embodiments of the present invention propose constructing an objective function that can accurately extract all reflection wave sub-signals from the seismic shot gather data. This objective function combines all reflection wave sub-signals to obtain reflection wavefield data, which are then removed from the seismic shot gather data to obtain 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 spectrum change characteristics of the signal at each moment. Smoothness can reflect the continuity and gradual change of the instantaneous frequency curve in time. 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 small-scale discontinuous geological bodies underground, carrying geological information about the complex geological structure and lithological differences of the underground. Therefore, the smooth instantaneous frequency corresponds to the regular structure of the signal (reflected wave), while the instantaneous frequency with relatively large fluctuations is usually related to the scattered wave. Therefore, when constructing the reflected wave objective function, the embodiment of the present invention adopts the instantaneous frequency smoothness constraint term of the reflected wave sub-signal to optimize the instantaneous frequency smoothness of the reflected wave sub-signal.

[0033] Furthermore, to avoid spectral overlap between reflected and scattered waves and improve scattered wave separation, an embodiment of the present invention adds a spectral energy distribution constraint term for the reflected wave sub-signal to the reflected wave objective function. This spectral energy distribution constraint term is used to control the center frequency and bandwidth of each reflected wave sub-signal, constraining the spectral energy of each reflected wave sub-signal to be as concentrated as possible near its own center spectrum, thereby ensuring the independence of the spectral distribution, avoiding interference from irrelevant spectra, and enhancing the accuracy of the separation results.

[0034] Furthermore, given that reflected wave signals exhibit low rank in mathematical representation, this embodiment of the present invention also incorporates a nuclear norm constraint on the reflected wave sub-signals within the reflected wave objective function to achieve low rank constraints on the reflected wave sub-signals. In this 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, the spectral energy distribution constraint, and the nuclear norm constraint for all reflected wave sub-signals.

[0035] Step S106: solving the reflected wave objective function to obtain reflected 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 combining all the reflection wave sub-signals can obtain 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] This embodiment of the present invention does not specifically limit the migration imaging method; users can select one based on their actual needs. Alternatively, using the Kirchhoff migration algorithm, the scattered wave field data can be migrated and imaged to obtain precise imaging results of the mine wastewater leakage path, that is, to obtain the mine wastewater leakage channel in the target mine wastewater treatment area.

[0041] The present invention provides a method for determining mine wastewater leakage channels. This method constructs a reflection wave objective function based on the instantaneous frequency smoothing constraint term, spectral energy distribution constraint term, and 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-signal 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 in the target mine wastewater treatment area can be obtained. This method can extract effective information about wastewater leakage channels from complex mine seismic wave fields, providing an efficient and accurate solution for the precise identification of mine wastewater leakage channels.

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

[0043]

[0044] in, represents the kth reflected wave sub-signal r k The instantaneous frequency smoothing constraint term of (x, t), δ(x, t) represents the Dirac distribution, x represents the position 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 recording time, and s(x,t) represents the seismic shot gather data.

[0045] In the instantaneous frequency smoothing constraint, Represents r k The analytical signal corresponding to (x, t) is a complex-valued signal that represents the original signal through its real part and provides phase information related to the original signal through its imaginary part. is the kernel function for constructing the analytical signal, The function is to reflect the sub-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 earthquake 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 sub-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 components and minimizing the value of this 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 results.

[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 low-rank constraints 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 a Lagrangian function based on the reflected wave objective function is equivalent to converting a 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 update each reflected wave sub-signal and its center frequency based on the Lagrangian function until the sum of the relative update amounts of all reflected wave sub-signal energies 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] In the embodiment of the present invention, it is necessary to control the reflected wave sub-signals and their center frequencies to perform multiple rounds of iterative updates. After initializing all reflected wave sub-signals and their center frequencies, 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, the formula for iteratively updating the reflected wave sub-signal is: in, represents the kth reflected wave sub-signal after the n+1th iteration, Represents the set of reflected wave sub-signals with sequence numbers less than k after the n+1th round of 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 sequentially in ascending order of serial number. During 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. To facilitate understanding, the sequential updating of reflected wave sub-signals is illustrated by an example. Starting from the n+1th round, the first reflected wave sub-signal is to be updated (i.e., when k=1). The reflected wave signals used are all the results of the nth round of iteration: When updating the second reflected wave sub-signal (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, 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 the 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 their updates, the center frequencies of the reflected wave sub-signals are sequentially updated. Since all reflected wave sub-signals have completed their updates in this round, when the center frequencies are updated, all reflected wave sub-signals use the updated results. To facilitate understanding, the sequential update of the center frequencies of the reflected wave sub-signals is illustrated as follows. In the (n+1)th round, the center frequency of the first reflected wave sub-signal is updated (that is, when k=1). The center frequencies of the reflected wave signals used are all the results of the nth round of iteration: 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. Among them, 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 squared diagonal of the energy difference between the kth reflected wave sub-signal before and after the (n+1) iteration reflects the magnitude of the energy change in the kth reflected wave sub-signal before and after the update. Therefore, the relative update E of the reflected wave sub-signal energy can reflect the magnitude of the energy change in the reflected wave sub-signal before and after the update.

[0062] As can be seen from the above description, in the embodiment of the present invention, the iteration termination condition of the reflected wave sub-signals and their center frequencies 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 is completed, all reflected wave sub-signals and their corresponding center frequencies are obtained. All reflected wave sub-signals are combined to obtain the reflected wave field data. This reflected wave field data is then removed from the seismic shot gather data to separate the scattered wave field data. The separated scattered wave field data is then migrated and imaged using the Kirchhoff migration algorithm to determine the mine wastewater leakage channel.

[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, spectral energy distribution constraint term, and nuclear norm constraint term of all reflected wave sub-signals, achieving accurate extraction of all reflected wave sub-signals and avoiding signal aliasing problems. When solving the Lagrangian function, the convergence accuracy is improved through an efficient iterative algorithm, making the embodiment of the present invention highly robust and suitable for processing complex non-stationary seismic signals. It provides an efficient and accurate solution for the precise identification and monitoring of mine wastewater leakage channels and contaminated areas.

[0065] Example 2

[0066] An embodiment of the present invention further 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 detailed 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 solving 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] An embodiment of the present invention provides a device for determining mine wastewater leakage channels. The device 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 in the target mine wastewater treatment area can be obtained. The device can extract effective information about the wastewater leakage channel from the complex mine seismic wave field, providing an efficient and accurate solution for the precise identification of mine wastewater leakage channels.

[0074] Optionally, the reflected wave objective function is expressed as: in, represents the kth reflected wave sub-signal r k The instantaneous frequency smoothing constraint term of (x, t), δ(x, t) represents the Dirac distribution, x represents the position 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 gather 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 configured to:

[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 updated 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 the 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, Represents the set of reflected wave sub-signals with sequence numbers less than k after the n+1th round of 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 the 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 follows: in, represents the square of the second norm of the k-th reflected wave sub-signal after the n-th iteration, It represents the square of the energy difference of the k-th reflected wave signal before and after the (n+1)th iteration.

[0084] Example 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, wherein the processor 60, the communication interface 63 and the memory 61 are connected via 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 system network element and at least one other network element are connected via at least one communication interface 63 (which may be wired or wireless), and the Internet, a wide area network, a local area network, a 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 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. During implementation, the steps of the above method may be performed by hardware integrated logic circuits or software instructions within the processor 60. The processor 60 may be a general-purpose processor, including a central processing unit (CPU), a network processor (NP), etc.; it may 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, or discrete hardware components. It may implement or execute the methods, steps, and logic block diagrams disclosed in the embodiments of the present invention. The general-purpose processor may be a microprocessor or any conventional processor. The steps of the methods disclosed in conjunction with the embodiments of the present invention may be directly implemented and executed by a hardware decoding processor, or by a combination of hardware and software modules within the decoding processor. The software modules may be located in storage media well-known in the art, such as random access memory, flash memory, read-only memory, programmable read-only memory, electrically erasable programmable memory, registers, or the like. 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] The embodiments of the present invention provide a computer program product for a method, device, and electronic device for determining a mine wastewater leakage channel, including a computer-readable storage medium storing 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 and 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 enabling a computer device (which can be a personal computer, server, or network device, etc.) to execute 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, and other media that can store program code.

[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, it does not need to be further defined or explained 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," "inner," and "outer" and the like indicate positions or locations based on the positions shown in the accompanying drawings, or the positions or locations in which the inventive product is typically placed when in use. These terms are intended solely to facilitate the description of the present invention and to simplify the description, and are not intended to indicate or imply that the devices or components referred to must have a specific orientation, be constructed, or operate in a specific orientation. Therefore, they should not be construed as limitations on the present invention. Furthermore, the terms "first," "second," and "third," etc., are used solely to distinguish descriptions and should not be construed as indicating or implying relative importance.

[0095] Furthermore, terms such as "horizontal," "vertical," and "overhanging" do not necessarily imply that a component must be absolutely horizontal or overhanging, but rather that it can be slightly tilted. For example, "horizontal" simply means that its direction is more horizontal than "vertical," and does not mean that the structure must be completely horizontal, but rather that it can be slightly tilted.

[0096] In the description of the present invention, it should also be noted that, unless otherwise expressly specified or limited, the terms "disposed," "installed," "connected," and "connected" should be understood in a broad sense. For example, they may refer to fixed connections, detachable connections, or integral connections; they may refer to mechanical connections or electrical connections; they may refer to direct connections or indirect connections through an intermediate medium; and they may refer to internal communication between two components. Those skilled in the art will understand the specific meanings of the above terms in the present invention based on the 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 above embodiments, those skilled in the art should understand that they can still modify the technical solutions described in the above 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 for 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 field data; Eliminating the reflected wave field data from the seismic shot gather data to obtain scattered wave field data; Migration imaging is performed on the scattered wave field data 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 signal The instantaneous frequency smoothing constraint term, represents the Dirac distribution, represents the location of the seismic trace, represents the earthquake recording time, represents the imaginary unit, , K represents the number of reflected wave sub-signals, Represents the reflected wave field data, express The corresponding analytical signal, express The center frequency, represents the partial derivative of earthquake record time, represents seismic shot gather data; ,express The spectrum energy distribution constraint term is: Indicates frequency Deviation from center frequency degree, represents the Fourier transform, , indicating the frequency The spectrum energy under ,express The nuclear norm constraint term, express The i-th singular value of .

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 updated 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, Represents the set of reflected wave sub-signals with sequence numbers less than k after the n+1th round of 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, represents the Lagrangian function.

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 the reflected wave sub-signals with sequence numbers less than k after the n+1th iteration. Represents the set of center frequencies of the reflected wave sub-signals with sequence numbers greater than or equal to k after the nth round of iteration, represents the Lagrangian function.

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, It represents the square of the energy difference of the k-th reflected wave 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 the target mine wastewater treatment area; A construction module is used to construct 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, configured to remove the reflected wave field data from the seismic shot gather data to obtain scattered wave field data; An imaging module is used to perform migration 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