Method, device and equipment for detecting water guide channel of mine

Abstract

The embodiment of the invention provides a method and device for detecting a mine water guide channel, electronic equipment, a storage medium and a program product. The method comprises the following steps: in response to acquired seismic shot gather data associated with a mine target area, separating a time domain diffracted wave signal from the seismic shot gather data; performing dip angle domain conversion on the time domain diffracted wave signal by using a preset conversion rule to obtain a dip angle domain diffracted wave field; wherein the dip angle domain diffracted wave travel time equation is involved in the dip angle domain diffracted wave field; based on the travel time equation, utilizing Taylor approximation to construct a diffracted wave target function; on the basis of the target function, diffracted wave event inversion operation is executed, and an optimal superposition path is obtained; and based on the dip angle domain diffracted wave field and the optimal superposition path, performing superposition operation, and generating a water guide channel imaging result of the mine target area. By utilizing the method disclosed by the embodiment of the invention, the speed sensitivity limitation of a traditional method can be broken through, and high-resolution imaging of a small-scale water guide structure is realized.

Description

Technical Field

[0001] This disclosure relates to the fields of mine safety production and mining engineering technology, specifically to a method, apparatus, electronic device, computer-readable storage medium, and computer program product for detecting mine water diversion channels. Background Technology

[0002] Mining water-conducting channels (such as small-scale discontinuous structures like hidden fissures, faults, and collapse column boundaries) are key geological factors inducing mine water inrush, shaft flooding, and surface subsidence. If such structures form hydraulic connections with aquifers, they can easily and rapidly conduct water flow under mining disturbances, leading to major safety accidents. Therefore, high-precision identification of water-conducting channels is crucial for achieving safe mine production, precise water hazard prevention, and the construction of green mines.

[0003] Diffraction waves, with their unique physical properties—disregarding Snell's law and propagating without being limited by the angle of incidence—exhibit advantages such as wide-area illumination capabilities, high sensitivity to small-scale structures, and rich boundary detail information, making them a core information source for detecting concealed water-conducting channels. However, the high sensitivity of the diffraction wave phase axis to velocity errors in the dip domain and the difficulty in obtaining accurate offset velocities pose a dual challenge. This causes traditional horizontal stacking methods to suffer severe image defocusing even with minor velocity errors, resulting in the failure to locate underground water-conducting channels. Summary of the Invention

[0004] To address the aforementioned technical problems, the present disclosure provides a solution. Embodiments of this disclosure provide a method, apparatus, electronic device, computer-readable storage medium, and computer program product for detecting water diversion channels in mines.

[0005] According to a first aspect of the present disclosure, a method for detecting a mine water diversion channel is provided, wherein the method includes: In response to acquiring seismic shot gather data associated with the target mining area, a time-domain diffraction wave signal is separated from the seismic shot gather data; Using a preset conversion rule, the time-domain diffracted wave signal is converted to the tilt domain to obtain the tilt domain diffracted wave field; wherein, the tilt domain diffracted wave field involves the tilt domain diffracted wave travel time equation; Based on the travel time equation, a diffraction wave objective function is constructed using the Taylor approximation; wherein, the objective function characterizes the relationship between the phase axis shape of the diffraction wave and the offset velocity error; Based on the objective function, a diffraction wave phase axis inversion operation is performed to obtain the optimal superposition path; Based on the diffraction wave field in the tilt domain and the optimal stacking path, a stacking operation is performed to generate an imaging result of the water conduit channel in the target area of ​​the mine.

[0006] According to a second aspect of the present disclosure, an apparatus for detecting water diversion channels in mines is provided, wherein the apparatus includes: The data preprocessing unit is configured to: in response to acquiring seismic shot gather data associated with the target mining area, separate the time-domain diffraction wave signal from the seismic shot gather data; The data conversion processing unit is configured to: perform tilt domain conversion on the time-domain diffraction wave signal using a preset conversion rule to obtain a tilt domain diffraction wave field; wherein, the tilt domain diffraction wave field includes the tilt domain diffraction wave travel time equation; The objective function construction unit is configured to: construct a diffracted wave objective function using Taylor approximation based on the travel time equation; wherein the objective function characterizes the relationship between the phase axis shape of the diffracted wave and the offset velocity error; The inversion unit is configured to perform a diffraction wave phase axis inversion operation based on the objective function to obtain the optimal superposition path; The overlay unit is configured to perform an overlay operation based on the diffraction wave field in the tilt domain and the optimal overlay path to generate an imaging result of the water conduit channel in the target area of ​​the mine.

[0007] According to a third aspect of the present disclosure, an electronic device is provided, the electronic device comprising: a processor; a memory for storing executable instructions of the processor; the processor being configured to read the executable instructions from the memory and execute the instructions to implement the method for detecting mine water diversion channels as described in the present disclosure.

[0008] According to a fourth aspect of the present disclosure, a computer-readable storage medium is provided, the storage medium storing a computer program for executing the method for detecting mine water diversion channels as described in the present disclosure.

[0009] According to a fifth aspect of the present disclosure, a computer program product is provided, including a computer program, wherein when executed by a processor, the computer program implements the method for detecting a mine water diversion channel as described in the present disclosure.

[0010] In summary, the method for detecting water-conducting channels in mines provided in this disclosure is based on the inventive concept of "travel time approximation-inversion superposition." First, it constructs an approximate travel time equation for diffracted waves in the dip domain through Taylor expansion, quantitatively revealing the influence mechanism of velocity errors on the morphology of the in-phase axis. Then, it uses this equation as the objective function to invert and solve for the optimal diffracted wave propagation path. Finally, it performs focused superposition of diffracted wave energy along this path, overcoming the velocity sensitivity limitations of traditional methods and achieving high-resolution imaging of small-scale water-conducting structures. This method provides a geological exploration solution with universal engineering applicability for mine disaster prevention and control. Attached Figure Description

[0011] The above and other objects, features, and advantages of this disclosure will become more apparent from the more detailed description of the embodiments thereof in conjunction with the accompanying drawings. The drawings are provided to further illustrate the embodiments of this disclosure and form part of the specification. They are used together with the embodiments of this disclosure to explain the disclosure and do not constitute a limitation thereof. In the drawings, the same reference numerals generally represent the same components or steps.

[0012] Figure 1 This is a flowchart illustrating a method for detecting a mine water diversion channel provided in an exemplary embodiment of this disclosure; Figure 2 This is a public announcement Figure 1 An exemplary flowchart of a method for detecting water diversion channels in a mine, provided in this embodiment; Figure 3 This is a public announcement Figure 1 Another exemplary flowchart of the method for detecting mine water diversion channels provided in the embodiments; Figure 4 This is a schematic diagram of the structure of an apparatus for detecting water diversion channels in a mine, provided in an exemplary embodiment of this disclosure; Figure 5 This is a schematic diagram of the structure of an application embodiment of the electronic device disclosed herein. Detailed Implementation

[0013] The present disclosure will be further described below with reference to the embodiments shown in the accompanying drawings. Obviously, the described embodiments are merely some embodiments of the present disclosure, and not all embodiments of the present disclosure. It should be understood that the present disclosure is not limited to the exemplary embodiments described herein.

[0014] It should be noted that, unless otherwise specifically stated, the relative arrangement, numerical expressions, and values ​​of the components and steps set forth in these embodiments do not limit the scope of this disclosure.

[0015] Those skilled in the art will understand that the terms "first," "second," etc., in the embodiments of this disclosure are only used to distinguish different steps, devices, or modules, and do not represent any specific technical meaning, nor do they indicate a necessary logical order between them.

[0016] It should also be understood that in the embodiments disclosed herein, "a plurality of" may refer to two or more, and "at least one" may refer to one, two or more.

[0017] It should also be understood that any component, data or structure mentioned in the embodiments of this disclosure can generally be understood as one or more unless expressly defined or given to the contrary in the context.

[0018] Furthermore, the term "and / or" in this disclosure is merely a description of the relationship between related objects, indicating that three relationships can exist. For example, A and / or B can represent: A existing alone, A and B existing simultaneously, or B existing alone. Additionally, the character " / " in this disclosure generally indicates that the preceding and following related objects have an "or" relationship.

[0019] It should also be understood that the description of the various embodiments in this disclosure emphasizes the differences between the various embodiments, and the similarities or similarities can be referred to each other. For the sake of brevity, they will not be described in detail.

[0020] At the same time, it should be understood that, for ease of description, the dimensions of the various parts shown in the accompanying drawings are not drawn according to actual scale.

[0021] The following description of at least one exemplary embodiment is merely illustrative and is in no way intended to limit this disclosure or its application or use.

[0022] Techniques, methods, and equipment known to those skilled in the art may not be discussed in detail, but where appropriate, such techniques, methods, and equipment should be considered part of the specification.

[0023] It should be noted that similar labels and letters in the following figures indicate similar items; therefore, once an item is defined in one figure, it does not need to be discussed further in subsequent figures.

[0024] Example 1 Figure 1 This is a schematic flowchart of a method for detecting water diversion channels in a mine, provided by an exemplary embodiment of this disclosure. The method can be executed on a server or an edge computing terminal, wherein the server may include, but is not limited to, a server or a cloud computing platform.

[0025] Specifically, refer to Figure 1 The method for detecting mine water diversion channels includes: S110. In response to acquiring seismic shot gather data associated with the target mining area, the time-domain diffraction wave signal is separated from the seismic shot gather data.

[0026] Among them, seismic shot gather data refers to the reflection and diffraction signals recorded during seismic exploration.

[0027] Since diffracted waves are sensitive to small-scale structures but are usually mixed with reflected waves, it is necessary to process the seismic shot gather data to separate the diffracted wave field (target signal).

[0028] Optionally, the method for separating the time-domain diffraction wave signal from the seismic shot gather data may include, but is not limited to, any one or any combination of wave field vector decomposition, tilt angle filtering, and adaptive plane wave destruction filtering.

[0029] Optionally, the server or edge computing terminal, acting as the execution entity, can obtain the seismic shot gather data by communicating with the storage service that stores the "seismic shot gather data".

[0030] S120. Using a preset conversion rule, perform tilt domain conversion on the time domain diffraction wave signal to obtain the tilt domain diffraction wave field.

[0031] The diffracted wave field in the tilt domain involves the travel time equation of the diffracted wave in the tilt domain.

[0032] Optionally, step S120 can be implemented in the following way: Using the following calculation formula, the time-domain diffracted wave signal is subjected to tilt domain transformation to obtain the tilt domain diffracted wave field; (1) (2)

[0033]

[0034]

[0035] in, This represents the diffracted wave field in the tilt region; This represents the time-domain diffracted wave signal; Represents the Dirac function; Represent the travel time equation for diffracted waves in the tilt domain; Indicates the offset angle; Indicates the depth of the diffraction point; Indicates the actual velocity of the formation; Indicates the offset speed; Spatial location of the imaging point; Spatial location of the diffraction point; Indicates spatial location; Indicates time; This indicates the vertical two-way travel time of the reflection point.

[0036] S130. Based on the travel time equation, construct the diffraction wave objective function using the Taylor approximation.

[0037] The objective function characterizes the relationship between the phase axis shape of the diffracted wave and the offset velocity error.

[0038] Optionally, step S130 can be implemented as follows: performing a Taylor expansion of the travel time equation at a reference point to obtain the diffracted wave travel time formula, which is used as the objective function.

[0039] The reference point satisfies the following condition: the offset angle The value is 0; the ratio of the offset velocity to the true velocity of the formation is 1; the ratio of the lateral distance between the imaging point and the diffraction point to the depth of the diffraction point is 0. That is, it can be expressed using the calculation formula: .

[0040] Optionally, refer to Figure 2 The phrase "performing a Taylor expansion of the travel time equation at a reference point to obtain the diffracted wave travel time formula" may include: S1310. Perform Taylor expansion on the cosine, numerator, and denominator terms in the travel time equation at the reference point to obtain the corresponding cosine Taylor expansion, numerator Taylor expansion, and denominator Taylor expansion.

[0041] Among them, with Using the reference point, the standard expansion of the multi-parameter Taylor expansion of the travel time equation in the diffracted wave offset tilt domain can be calculated as follows (3): (3)

[0042]

[0043] in, Indicates the offset angle; same , representing the ratio of the offset velocity to the actual velocity of the formation; same , which represents the ratio of the lateral distance between the imaging point and the diffraction point to the depth of the diffraction point; This represents a higher-order remainder operator.

[0044] As a specific example, based on equation (3), the numerator term in the time-travel equation (i.e., equation (2)) is analyzed. Performing a Taylor expansion at the reference point yields the corresponding second-order approximation. Specifically, this can be expressed using the following formula:

[0045] Due to the offset tilt angle Satisfying the small angle assumption, therefore, it can be... Furthermore, the numerator of the diffraction travel time equation (i.e., equation (2)) in the tilt domain can be... Approximately:

[0046] As a specific example, the numerator and denominator of equation (2) are first multiplied together. Then, by simplifying the trigonometric identities, we obtain the following recombined calculation formula (4): (4) As a specific example, based on equation (3), a geometric series expansion is performed on the denominator term of the travel time equation (i.e., equation (4)) at the reference point to obtain the corresponding higher-order approximation. Specifically, it can be expressed using the following formula:

[0047]

[0048] As a specific example, based on equation (3), a Taylor expansion is performed on the cosine term in the travel time equation (i.e., equation (2) or (4)) at the reference point to obtain the corresponding higher-order approximation. Specifically, it can be expressed using the following formula:

[0049] S1320. Substitute the cosine Taylor expansion, the numerator Taylor expansion, and the denominator Taylor expansion into the travel time equation to obtain the following calculation formula for the diffraction wave travel time. (5)

[0050]

[0051] in, Indicates the offset angle; Indicates the depth of the diffraction point; Indicates the actual velocity of the formation; This indicates the offset speed.

[0052] Specifically, substituting the expansions of the numerator, denominator, and cosine terms obtained in step S1310 into equation (4) yields: (6)

[0053] in, This represents the error between the image point of the diffracted wave and the actual diffracted point. Indicates the depth of the diffraction point; same .

[0054] Since the offset velocity is the result of multiple rounds of iterative calculations, it can be considered that... Much smaller than the depth of the diffraction point Based on this, it can be Since it is approximately 0, the travel time formula of diffracted wave in the tilt domain can be approximated from equation (6) to equation (5).

[0055] It should be noted that the above approximate result (i.e., equation (5)) conforms to the parabolic expression, where, , The above equation clearly shows that, in the offset tilt domain, the shape of the diffracted wave's phase axis is closely related to the offset velocity. When the offset velocity is accurate (i.e., When the value is 1), the phase axis of the concentrated diffracted wave in the tilt channel appears as a straight line, that is... When the migration velocity is inaccurate, the phase axis of the diffracted wave is no longer a straight line, but exhibits an upward or downward parabolic shape. Furthermore, when dealing with actual seismic data with low signal-to-noise ratios, the overlap of multiple wavefields complicates the velocity spectrum. No matter how much the velocity analysis method is improved, the accuracy of the migration velocity can only be increased to a certain extent, and absolute accuracy cannot be guaranteed. The phase axis shape of the diffracted wave in the dip domain is highly sensitive to the migration velocity. Therefore, when using traditional horizontal stacking methods for imaging, even small velocity errors can lead to severe defocusing of the diffracted imaging results, resulting in low resolution.

[0056] S140. Based on the objective function, perform the diffraction wave phase axis inversion operation to obtain the optimal superposition path.

[0057] Optionally, refer to Figure 3 Step S140 may include: S1410. Determine the shape parameter A of the phase axis of the diffracted wave based on the objective function.

[0058] S1420. Determine the parameter range of the shape parameter A based on the preset error range.

[0059] S1430. Within the parameter range, search for the shape parameter A that maximizes the superposition energy of the diffracted waves, and determine the parabolic path based on the shape parameter A, denoted as the optimal superposition path.

[0060] Wherein, the shape parameter A is the aforementioned .

[0061] Since the velocity analysis is the result of multiple iterations, its error will not be too large. As an optional example, compared to the accurate velocity, the preset error range can be ±15%, that is, Therefore, we can know .

[0062] In the separated diffracted wave field, the diffracted wave, as the target wave field, has the highest energy. Therefore, finding the parabola with the highest superposition energy is sufficient to determine the phase axis of the diffracted wave. As mentioned above, the phase axis of the diffracted wave in the tilt domain exhibits a parabolic shape. However, due to the unknowns in the actual data, the size and direction of the parabola's opening cannot be determined. Unable to determine. To accurately identify the shape of each diffracted wave's phase axis, we need to further investigate... The system performs a traversal search within the path, calculating the cumulative energy value of the corresponding path. This searches for the parabolic path that maximizes the superposition energy of the diffracted waves, which is then considered the optimal superposition path.

[0063] S150. Based on the diffraction wave field in the tilt domain and the optimal superposition path, perform a superposition operation to generate the water channel imaging result of the target area of ​​the mine.

[0064] Optionally, step S150 can be implemented by performing nonlinear superposition of the diffracted wave field in the tilt domain along the optimal superposition path to generate the water channel imaging result of the target area of ​​the mine.

[0065] The nonlinear superposition satisfies the following conditions: when the offset velocity error is absent, the phase axis of the diffracted wave is a straight line; when the offset velocity error exists, the phase axis of the diffracted wave is a parabola.

[0066] In summary, the method for detecting water-conducting channels in mines provided in this disclosure is based on the inventive concept of "travel time approximation-inversion superposition." First, it constructs an approximate travel time equation for diffracted waves in the dip domain through Taylor expansion, quantitatively revealing the influence mechanism of velocity errors on the morphology of the in-phase axis. Then, it uses this equation as the objective function to invert and solve for the optimal diffracted wave propagation path. Finally, it performs focused superposition of diffracted wave energy along this path, overcoming the velocity sensitivity limitations of traditional methods and achieving high-resolution imaging of small-scale water-conducting structures. This method provides a geological exploration solution with universal engineering applicability for mine disaster prevention and control.

[0067] Example 2 It should be understood that the methods described in the foregoing embodiments for detecting mine water diversion channels can also be similarly applied to the following apparatus for detecting mine water diversion channels, with similar extensions. For simplicity, they are not described in detail.

[0068] Figure 4 This is a schematic diagram of an apparatus for detecting water diversion channels in a mine, provided in an exemplary embodiment of this disclosure. (Refer to...) Figure 4 The device includes: The data preprocessing unit 410 is configured to: in response to acquiring seismic shot gather data associated with the target mining area, separate the time-domain diffraction wave signal from the seismic shot gather data; The data conversion processing unit 420 is configured to: perform tilt domain conversion on the time domain diffraction wave signal using a preset conversion rule to obtain a tilt domain diffraction wave field; wherein, the tilt domain diffraction wave field includes the tilt domain diffraction wave travel time equation; The objective function construction unit 430 is configured to: construct a diffracted wave objective function using Taylor approximation based on the travel time equation; wherein the objective function characterizes the relationship between the phase axis shape of the diffracted wave and the offset velocity error; Inversion unit 440 is configured to perform diffraction wave phase axis inversion operation based on the objective function to obtain the optimal superposition path; The overlay unit 450 is configured to perform an overlay operation based on the diffraction wave field in the tilt domain and the optimal overlay path to generate an imaging result of the water channel in the target area of ​​the mine.

[0069] In summary, the device for detecting water-conducting channels in mines provided in this disclosure is based on the inventive concept of "travel time approximation-inversion superposition." First, it constructs an approximate travel time equation for diffracted waves in the dip domain through Taylor expansion, quantitatively revealing the influence mechanism of velocity errors on the morphology of the in-phase axis. Then, it uses this equation as the objective function to invert and solve for the optimal diffracted wave propagation path. Finally, it performs focused superposition of diffracted wave energy along this path, overcoming the velocity sensitivity limitations of traditional methods and achieving high-resolution imaging of small-scale water-conducting structures. This method provides a geological exploration solution with universal engineering applicability for mine disaster prevention and control.

[0070] Example 3 In addition, this disclosure also provides an electronic device, including: a memory for storing a computer program; and a processor for executing the computer program stored in the memory, wherein when the computer program is executed, it implements the method for detecting a mine water diversion channel as described in any of the above embodiments of this disclosure.

[0071] Figure 5 This is a schematic diagram of the structure of an application embodiment of the electronic device disclosed herein. Below, reference is made to… Figure 5 This describes an electronic device according to embodiments of the present disclosure. The electronic device may be either or both of a first device and a second device, or a standalone device independent of them, which may communicate with the first device and the second device to receive acquired input signals from them.

[0072] like Figure 5As shown, the electronic device includes one or more processors and memory. The processor may be a central processing unit (CPU) or other processing unit with data processing capabilities and / or instruction execution capabilities, and may control other components in the electronic device to perform desired functions. The memory may include one or more computer program products, which may include various forms of computer-readable storage media, such as volatile memory and / or non-volatile memory. The volatile memory may, for example, include random access memory (RAM) and / or cache memory. The non-volatile memory may, for example, include read-only memory (ROM), hard disk, flash memory, etc. One or more computer program instructions may be stored on the computer-readable storage medium, and the processor may execute the program instructions to implement the methods for detecting mine water diversion channels and / or other desired functions described in the various embodiments of this disclosure above.

[0073] In one example, the electronic device may further include input and output devices, which are interconnected via a bus system and / or other forms of connection mechanisms (not shown). Furthermore, the input device may include, for example, a keyboard, a mouse, etc. The output device can output various information to the outside, including determined distance information, direction information, etc. The output device may include, for example, a display, a speaker, a printer, and a communication network and its connected remote output devices, etc.

[0074] Of course, for the sake of simplicity, Figure 5 Only some of the components of the electronic device relevant to this disclosure are shown, omitting components such as buses, input / output interfaces, etc. In addition, the electronic device may include any other suitable components depending on the specific application.

[0075] In addition to the methods and devices described above, embodiments of this disclosure may also be computer program products comprising computer program instructions that, when executed by a processor, cause the processor to perform the steps of the methods for detecting mine water diversion channels according to various embodiments of this disclosure as described in the foregoing portion of this specification.

[0076] The computer program product can be written in any combination of one or more programming languages ​​to perform the operations of the embodiments of this disclosure. The programming languages ​​include object-oriented programming languages ​​such as Java and C++, as well as conventional procedural programming languages ​​such as C or similar languages. The program code can be executed entirely on a user's computing device, partially on a user's computing device, as a standalone software package, partially on a user's computing device and partially on a remote computing device, or entirely on a remote computing device or server.

[0077] Furthermore, embodiments of this disclosure may also be computer-readable storage media storing computer program instructions that, when executed by a processor, cause the processor to perform the steps of the method for detecting a mine water diversion channel according to various embodiments of this disclosure as described in the foregoing portion of this specification.

[0078] The computer-readable storage medium may be any combination of one or more readable media. A readable medium may be a readable signal medium or a readable storage medium. A readable storage medium may be, for example, an electrical, magnetic, optical, electromagnetic, infrared, or semiconductor system, apparatus, or device, or any combination thereof. More specific examples of readable storage media (a non-exhaustive list) include: an electrical connection having one or more wires, a portable disk, a hard disk, random access memory (RAM), read-only memory (ROM), erasable programmable read-only memory (EPROM or flash memory), optical fiber, portable compact disk read-only memory (CD-ROM), optical storage device, magnetic storage device, or any suitable combination thereof.

[0079] Those skilled in the art will understand that all or part of the steps of the above method embodiments can be implemented by hardware related to program instructions. The aforementioned program can be stored in a computer-readable storage medium. When the program is executed, it performs the steps of the above method embodiments. The aforementioned storage medium includes various media that can store program code, such as ROM, RAM, magnetic disk, or optical disk.

[0080] The basic principles of this disclosure have been described above with reference to specific embodiments. However, it should be noted that the advantages, benefits, and effects mentioned in this disclosure are merely examples and not limitations, and should not be considered as essential features of each embodiment of this disclosure. Furthermore, the specific details disclosed above are for illustrative and facilitative purposes only, and are not limitations. These details do not limit the scope of this disclosure to the necessity of employing the aforementioned specific details for implementation.

[0081] The various embodiments in this specification are described in a progressive manner, with each embodiment focusing on its differences from other embodiments. Similar or identical parts between embodiments can be referred to interchangeably. For system embodiments, since they largely correspond to method embodiments, the description is relatively simple; relevant parts can be referred to the descriptions in the method embodiments.

[0082] The block diagrams of devices, apparatuses, devices, and systems disclosed herein are merely illustrative examples and are not intended to require or imply that they must be connected, arranged, or configured in the manner shown in the block diagrams. As those skilled in the art will recognize, these devices, apparatuses, devices, and systems can be connected, arranged, and configured in any manner. Words such as “comprising,” “including,” “having,” etc., are open-ended terms meaning “including but not limited to,” and are used interchangeably with them. The terms “or” and “and” as used herein refer to the terms “and / or,” and are used interchangeably with them unless the context clearly indicates otherwise. The term “such as” as used herein refers to the phrase “such as but not limited to,” and is used interchangeably with it.

[0083] The methods and apparatus of this disclosure may be implemented in many ways. For example, they may be implemented by software, hardware, firmware, or any combination of software, hardware, and firmware. The above-described order of steps for the methods is for illustrative purposes only, and the steps of the methods of this disclosure are not limited to the order specifically described above unless otherwise specifically stated. Furthermore, in some embodiments, this disclosure may also be implemented as a program recorded on a recording medium, the program including machine-readable instructions for implementing the methods according to this disclosure. Thus, this disclosure also covers recording media storing programs for performing the methods according to this disclosure.

[0084] It should also be noted that in the apparatus, devices, and methods of this disclosure, the components or steps can be disassembled and / or recombined. These disassemblies and / or recombinations should be considered as equivalent solutions to this disclosure.

[0085] The above description of the disclosed aspects is provided to enable any person skilled in the art to make or use this disclosure. Various modifications to these aspects will be readily apparent to those skilled in the art, and the general principles defined herein may be applied to other aspects without departing from the scope of this disclosure. Therefore, this disclosure is not intended to be limited to the aspects shown herein, but rather to be carried out within the widest scope consistent with the principles and novel features disclosed herein.

[0086] The above description has been given for purposes of illustration and description. Furthermore, this description is not intended to limit the embodiments of this disclosure to the forms disclosed herein. Although numerous exemplary aspects and embodiments have been discussed above, those skilled in the art will recognize certain variations, modifications, alterations, additions, and sub-combinations therein.

Claims

1. A method for detecting water diversion channels in mines, characterized in that, The method includes: In response to acquiring seismic shot gather data associated with the target mining area, a time-domain diffraction wave signal is separated from the seismic shot gather data; Using a preset conversion rule, the time-domain diffracted wave signal is converted to the tilt domain to obtain the tilt domain diffracted wave field; wherein, the tilt domain diffracted wave field involves the tilt domain diffracted wave travel time equation; Based on the travel time equation, a diffraction wave objective function is constructed using the Taylor approximation; wherein, the objective function characterizes the relationship between the phase axis shape of the diffraction wave and the offset velocity error; Based on the objective function, a diffraction wave phase axis inversion operation is performed to obtain the optimal superposition path; Based on the diffraction wave field in the tilt domain and the optimal stacking path, a stacking operation is performed to generate an imaging result of the water conduit channel in the target area of ​​the mine.

2. The method according to claim 1, characterized in that, The step of performing tilt domain transformation on the time-domain diffracted wave signal using a preset transformation rule to obtain the tilt domain diffracted wave field includes: Using the following calculation formula, the time-domain diffracted wave signal is subjected to tilt domain transformation to obtain the tilt domain diffracted wave field; in, This represents the diffracted wave field in the tilt region; This represents the time-domain diffracted wave signal; Represents the Dirac function; Represent the travel time equation for diffracted waves in the tilt domain; Indicates the offset angle; Indicates the depth of the diffraction point; Indicates the actual velocity of the formation; Indicates the offset speed; Spatial location of the imaging point; Spatial location of the diffraction point; Indicates spatial location; Indicates time.

3. The method according to claim 1, characterized in that, The construction of the diffracted wave objective function based on the travel time equation using the Taylor approximation includes: The travel time equation is expanded using Taylor at a reference point to obtain the travel time formula of the diffracted wave, which serves as the objective function.

4. The method according to claim 3, characterized in that, The reference point satisfies the following condition: the offset angle The value is 0; the ratio of the offset velocity to the actual velocity of the stratum is 1; the ratio of the lateral distance between the imaging point and the diffraction point to the depth of the diffraction point is 0.

5. The method according to claim 3, characterized in that, The Taylor expansion of the travel time equation at the reference point to obtain the diffracted wave travel time formula includes: For the cosine term, numerator term, and denominator term in the travel time equation, Taylor expansion is performed at the reference point to obtain the corresponding cosine Taylor expansion, numerator Taylor expansion, and denominator Taylor expansion. Substituting the cosine Taylor expansion, the numerator Taylor expansion, and the denominator Taylor expansion into the travel time equation, we obtain the following formula for calculating the travel time of the diffracted wave. in, Indicates the offset angle; Indicates the depth of the diffraction point; Indicates the actual velocity of the formation; This indicates the offset speed.

6. The method according to claim 1, characterized in that, The step of performing diffracted wave phase axis inversion operation based on the objective function to obtain the optimal superposition path includes: The shape parameter A of the phase axis of the diffracted wave is determined based on the objective function; The parameter range of the shape parameter A is determined based on a preset error range; Within the parameter range, a shape parameter A is searched to maximize the superposition energy of the diffracted waves, and a parabolic path is determined based on the shape parameter A, denoted as the optimal superposition path.

7. The method according to claim 5, characterized in that, The process of performing a superposition operation based on the diffracted wave field in the tilt domain and the optimal superposition path to generate an imaging result of the water-conducting channel in the target area of ​​the mine includes: Nonlinear superposition of the diffraction wave field in the tilt domain is performed along the optimal superposition path to generate the water-conducting channel imaging result of the target area of ​​the mine. Wherein, the nonlinear superposition satisfies: When the offset velocity error is absent, the phase axis of the diffracted wave is straight. When the offset velocity error exists, the phase axis of the diffracted wave is parabolic.

8. The method according to claim 1, characterized in that, Methods for separating time-domain diffracted wave signals from the seismic shot gather data include: wave field vector decomposition, tilt angle filtering, and adaptive plane wave destruction filtering.

9. A device for detecting water diversion channels in mines, characterized in that, The device includes: The data preprocessing unit is configured to: in response to acquiring seismic shot gather data associated with the target mining area, separate the time-domain diffraction wave signal from the seismic shot gather data; The data conversion processing unit is configured to: perform tilt domain conversion on the time-domain diffraction wave signal using a preset conversion rule to obtain a tilt domain diffraction wave field; wherein, the tilt domain diffraction wave field includes the tilt domain diffraction wave travel time equation; The objective function construction unit is configured to: construct a diffracted wave objective function using Taylor approximation based on the travel time equation; wherein the objective function characterizes the relationship between the phase axis shape of the diffracted wave and the offset velocity error; The inversion unit is configured to perform a diffraction wave phase axis inversion operation based on the objective function to obtain the optimal superposition path; The overlay unit is configured to perform an overlay operation based on the diffraction wave field in the tilt domain and the optimal overlay path to generate an imaging result of the water conduit channel in the target area of ​​the mine.

10. An electronic device, the electronic device comprising: processor; Memory used to store the processor's executable instructions; The processor is configured to read the executable instructions from the memory and execute the instructions to implement the method for detecting mine water diversion channels as described in claims 1 to 8.

Images