High-precision seismic signal enhancement method, device and electronic equipment
By constructing quadratic equations and using filtering techniques, the problem of insufficient accuracy in seismic signal enhancement with curved wavefronts in existing technologies has been solved, and high-precision seismic signal enhancement has been achieved.
Patent Information
- Authority / Receiving Office
- CN · China
- Patent Type
- Patents(China)
- Current Assignee / Owner
- BEIJING INFORMATION SCI & TECH UNIV
- Filing Date
- 2025-05-28
- Publication Date
- 2026-05-19
AI Technical Summary
Existing seismic signal enhancement methods assume that seismic waves are plane waves, which makes it difficult to effectively fit seismic signals with curved wavefronts, resulting in the loss of effective information.
A quadratic equation characterizing the travel time of the wavefront of a seismic signal is constructed. The accuracy of the seismic signal is improved through solving and filtering techniques. This method is suitable for enhancing seismic signals with a certain curvature on the wavefront.
It improves the accuracy of earthquake signal prediction, reduces the loss of effective information, and achieves high-precision earthquake signal enhancement.
Smart Images

Figure CN120468940B_ABST
Abstract
Description
Technical Field
[0001] This invention relates to the field of earthquake data processing technology, and in particular to a high-precision earthquake signal enhancement method, apparatus, and electronic equipment. Background Technology
[0002] Seismic data is frequently affected by noise, which can disrupt the continuity of the data and obscure valuable information. Furthermore, many subsequent processing tasks, such as seismic attribute analysis, amplitude versus migration (AVO) analysis, and automated interpretation, are negatively impacted by noise. Therefore, improving the signal-to-noise ratio of seismic data to enhance its valuable information is crucial for improving the efficiency of seismic exploration.
[0003] Denoising methods based on underground geological information and the physical laws of seismic waves can enhance the effective signal by predicting the effective signal or interfering noise. Existing seismic signal enhancement methods generally assume that seismic waves are plane waves. However, seismic signal enhancement methods based on the plane wave assumption are only suitable for seismic signals with obvious local linear characteristics of the phase axis. They are difficult to effectively fit seismic signals with curved wavefronts. In the case of curved wavefronts, the inability to accurately predict the effective seismic information leads to the loss of effective information during the seismic signal enhancement process. Summary of the Invention
[0004] In view of this, the purpose of the present invention is to provide a high-precision seismic signal enhancement method, apparatus and electronic device to alleviate the above-mentioned problems existing in the existing seismic signal enhancement technology.
[0005] In a first aspect, embodiments of the present invention provide a high-precision seismic signal enhancement method, comprising: acquiring raw seismic data; constructing a quadratic equation characterizing the wavefront travel time of the seismic signal; solving the equation based on the raw seismic data and the quadratic equation to obtain the wavefront travel time of the seismic signal in the raw seismic data; and filtering the raw seismic data based on the wavefront travel time to obtain the enhanced seismic signal.
[0006] Secondly, embodiments of the present invention also provide a high-precision seismic signal enhancement device, comprising: an acquisition module for acquiring raw seismic data; a construction module for constructing a quadratic equation characterizing the wavefront travel time of a seismic signal; a solution module for solving the equation based on the raw seismic data and the quadratic equation to obtain the wavefront travel time of the seismic signal from the raw seismic data; and a filtering module for filtering the raw seismic data based on the wavefront travel time of the seismic signal to obtain an enhanced seismic signal.
[0007] Thirdly, embodiments of the present invention also provide an electronic device, including a processor and a memory, wherein the memory stores computer-executable instructions that can be executed by the processor, and the processor executes the computer-executable instructions to implement the high-precision seismic signal enhancement method described in the first aspect above.
[0008] This invention provides a high-precision seismic signal enhancement method, apparatus, and electronic device. The method involves acquiring raw seismic data, constructing a quadratic equation characterizing the wavefront travel time of the seismic signal, solving the equation based on the raw seismic data and the quadratic equation to obtain the wavefront travel time of the original seismic signal, and filtering the raw seismic data based on the wavefront travel time to obtain the enhanced seismic signal. This technique is suitable for enhancing seismic signals with a certain curvature on the wavefront. Compared to existing seismic signal enhancement techniques, it improves the accuracy of seismic signal prediction and effectively reduces the loss of effective information during the seismic signal enhancement process, thereby achieving high-precision seismic signal enhancement.
[0009] Other features and advantages of the invention will be set forth in the description which follows, and will be apparent in part from the description, or may be learned by practicing the invention. The objects and other advantages of the invention are realized and obtained in accordance with the structures particularly pointed out in the description, claims and drawings.
[0010] To make the above-mentioned objects, features and advantages of the present invention more apparent and understandable, preferred embodiments are described below in detail with reference to the accompanying drawings. Attached Figure Description
[0011] To more clearly illustrate the specific embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the specific embodiments or the prior art will be briefly introduced below. Obviously, the drawings described below are some embodiments of the present invention. For those skilled in the art, other drawings can be obtained from these drawings without creative effort.
[0012] Figure 1 This is a flowchart illustrating a high-precision seismic signal enhancement method according to an embodiment of the present invention;
[0013] Figure 2 This is a flowchart illustrating the implementation of high-precision seismic signal enhancement in an embodiment of the present invention;
[0014] Figure 3 This is a schematic diagram of the structure of a high-precision seismic signal enhancement device according to an embodiment of the present invention;
[0015] Figure 4 This is a schematic diagram of the structure of an electronic device according to an embodiment of the present invention. Detailed Implementation
[0016] To make the objectives, technical solutions, and advantages of the embodiments of the present invention clearer, the technical solutions of the present invention will be clearly and completely described below in conjunction with the embodiments. Obviously, the described embodiments are only some embodiments of the present invention, not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of the present invention.
[0017] Current seismic signal enhancement methods generally assume that seismic waves are plane waves. However, these methods are only suitable for seismic signals with obvious local linear characteristics along the same phase axis, and are difficult to effectively fit seismic signals with curved wavefronts. Furthermore, in cases with curved wavefronts, the inability to accurately predict effective seismic information leads to the loss of this information during the enhancement process. Therefore, this invention provides a high-precision seismic signal enhancement method, apparatus, and electronic equipment that can alleviate the aforementioned problems in existing seismic signal enhancement technologies.
[0018] To facilitate understanding of this embodiment, a high-precision seismic signal enhancement method disclosed in this invention will first be described in detail. (See [link to relevant documentation]). Figure 1 As shown, the method may include the following steps:
[0019] Step S102: Obtain raw seismic data.
[0020] The original seismic data can be two-dimensional data with time as the vertical axis and trace number as the horizontal axis, such as two-dimensional common offset data (i.e., data obtained through common offset processing technology). There are no restrictions on the acquisition method and data type of the original seismic data.
[0021] Step S104: Construct a quadratic equation characterizing the travel time pattern of the seismic signal wavefront.
[0022] When constructing the quadratic equation, the coefficients of the quadratic equation are unknown. It is necessary to solve for the coefficients of the quadratic equation to more accurately characterize the travel time pattern of the seismic signal wavefront.
[0023] Step S106: Solve the original seismic data and quadratic equation to obtain the travel time of the seismic signal wavefront from the original seismic data.
[0024] The original seismic data can be directly substituted into the quadratic equation for solution, or the quadratic equation can be processed before the original seismic data is substituted into the processed equation for solution, thereby obtaining the travel time of the seismic signal wavefront from the original seismic data.
[0025] Step S108: Based on the wavefront travel time of the seismic signal, the original seismic data is filtered to obtain the enhanced seismic signal.
[0026] The filtering can be done using moving average filtering, median filtering, etc., and there is no limitation on the filtering method here.
[0027] After obtaining the wavefront travel time of the original seismic signal, the original seismic data can be filtered along the obtained wavefront travel time to enhance the seismic signal and obtain the enhanced seismic signal.
[0028] This invention provides a high-precision seismic signal enhancement method. The method involves acquiring raw seismic data, constructing a quadratic equation characterizing the wavefront travel time of the seismic signal, solving the equation based on the raw seismic data to obtain the wavefront travel time of the original seismic signal, and then filtering the raw seismic data based on the wavefront travel time to obtain the enhanced seismic signal. This method is suitable for enhancing seismic signals with a certain curvature on the wavefront. Compared with existing seismic signal enhancement techniques, it improves the accuracy of seismic signal prediction and effectively reduces the loss of effective information during the enhancement process, thereby achieving high-precision seismic signal enhancement.
[0029] As one possible implementation, the raw seismic data can be two-dimensional data with time as the vertical axis and trace number as the horizontal axis; for example, two-dimensional seismic data... Extracted as raw seismic data x The track number indicating the earthquake record. t Indicates the time sample number recorded. x and t The sampling interval is 1. Based on this, we can assume that the travel time of the seismic signal wavefront follows a quadratic equation, which can be expressed as:
[0030]
[0031] in, t For time, x As a Taoist name, a , b and c The coefficients of the quadratic equation;
[0032] For the same wavefront, if its amplitude is constant, then the quadratic equation needs to satisfy the following constraint equations:
[0033]
[0034] in, P For amplitude, C It is a constant.
[0035] As one possible implementation, step S106 (i.e., solving the original seismic data and the quadratic equation to obtain the seismic signal wavefront travel time of the original seismic data) may include:
[0036] Step A1: Find the first and second derivatives of the constraint equations respectively.
[0037] Continuing from the previous example, we can apply the constraint equations By taking the first derivative and the second derivative respectively, the result of taking the first derivative of the constraint equation can be:
[0038]
[0039] The result of taking the second derivative of the constraint equations can be:
[0040] .
[0041] Step A2: Based on the original seismic data and the results of obtaining the first and second derivatives of the constraint equations, solve the coefficients of the quadratic equations to obtain the travel time of the seismic signal wavefront from the original seismic data.
[0042] For example, step A2 above may include steps A21 to A24:
[0043] Step A21: Combine the results of finding the first and second derivatives of the constraint equations to obtain a system of simultaneous equations.
[0044] Continuing from the previous example, we can obtain the first derivative result of the constraint equation. The result of finding the second derivative The equations are combined into a system of equations.
[0045] Step A22: Solve the simultaneous equations to obtain a system of partial differential equations.
[0046] Continuing from the previous example, the system of simultaneous equations obtained in step A21 can be solved to obtain... and These two partial differential equations (that is, the expressions for the coefficients of the quadratic equations) together form a system of partial differential equations.
[0047] Step A23: Based on the original seismic data, the partial differential parts of the partial differential equation system are solved using the finite difference method to obtain the coefficient values of the quadratic equation.
[0048] Following the previous example, raw seismic data can be used. The partial differential equations obtained in step A22 are solved using the finite difference method, yielding the partial differential terms of the system:
[0049]
[0050]
[0051]
[0052]
[0053]
[0054] In this formula, , , , , , , , , Equivalent to , , , , , , , , Substituting the obtained partial differential terms into the system of partial differential equations obtained in step A22 will yield the coefficients of the quadratic equation. a and b Their respective values.
[0055] Step A24: Based on the quadratic equation and the obtained coefficient values, determine the travel time of the seismic signal wavefront from the original seismic data.
[0056] Continuing from the previous example, we can use a quadratic equation. and its coefficients a and b The individual values are then processed to obtain the travel time of the seismic signal wavefront from the original seismic data. :
[0057]
[0058] in, i The change value of the channel number, i Take the integer part.
[0059] As one possible implementation, the filtering in step S108 (i.e., filtering the original seismic data based on the wavefront travel time of the seismic signal to obtain the enhanced seismic signal) can be moving average filtering or median filtering.
[0060] Continuing from the previous example, the filter window radius can be set. k When traveling along the front of the seismic signal wave Raw seismic data Moving average filtering is performed to enhance the seismic signal; the enhanced seismic signal... It can be represented as: .
[0061] To facilitate understanding, the process of achieving high-precision seismic signal enhancement using the above-mentioned high-precision seismic signal enhancement method is described in the following exemplary description using a specific application example.
[0062] See Figure 2 As shown, the high-precision seismic signal enhancement process mainly includes the following steps:
[0063] The first step is to read the earthquake data.
[0064] Reading earthquake data , This is two-dimensional stacked data (such as two-dimensional common offset data). x The track number indicating the earthquake record. t Indicates the time sample number recorded. x and t The sampling interval is 1.
[0065] The second step is to establish the travel time pattern of the seismic signal wavefront.
[0066] The plane wave assumption is often used for seismic signal enhancement with good local linearity. To address seismic signal enhancement with poor local linearity, the travel time of the seismic signal wavefront can be assumed to follow a quadratic equation:
[0067] (1)
[0068] In equation (1), a , b and c These are the three coefficients of the quadratic equation.
[0069] For the same wavefront, its amplitude is constant, that is:
[0070] (2)
[0071] In equation (2), P For amplitude, C It is a constant.
[0072] The third step is to calculate the travel time of the seismic signal wavefront.
[0073] Taking the first and second derivatives of equation (2) respectively, we get:
[0074] (3)
[0075] Solving equation (3), we get:
[0076] (4)
[0077] The partial differential in equation (4) can be obtained using the finite difference method:
[0078] (5)
[0079] In equation (5), , , , , , , , , Equivalent to , , , , , , , , .
[0080] Substituting equation (5) into equation (4) yields the coefficients of the quadratic equation. a and b .
[0081] After finding the coefficients of the quadratic equation a and b Then, from equation (1), we can obtain:
[0082] (6)
[0083] Equation (6) can represent the travel time of the seismic signal wavefront. .
[0084] The fourth step is to enhance the seismic signal by utilizing the wavefront travel time.
[0085] When walking along the obtained wave front earthquake data Seismic signal enhancement can be achieved by performing moving average filtering; the enhanced seismic signal... It can be represented as:
[0086] (7)
[0087] In equation (7), k This represents the radius of the filtering window.
[0088] The advantages of the above-mentioned high-precision seismic signal enhancement method are:
[0089] Traditional seismic signal enhancement methods based on the plane wave assumption are difficult to effectively fit seismic signals with curved wavefronts. The high-precision seismic signal enhancement method described above constructs a quadratic equation form for wavefront travel time, which is suitable for seismic signal enhancement with a certain curvature of the wavefront. It has higher accuracy in seismic signal prediction, can effectively reduce the loss of effective information during the signal enhancement process, and achieves high-precision seismic signal enhancement.
[0090] Based on the above-described high-precision seismic signal enhancement method, this invention also provides a high-precision seismic signal enhancement device, see [link to relevant documentation]. Figure 3 As shown, the device may include:
[0091] The acquisition module 302 is used to acquire raw seismic data.
[0092] Module 304 is used to construct a quadratic equation characterizing the travel time of the wavefront of a seismic signal.
[0093] The solver module 306 is used to solve the original seismic data and the quadratic equation to obtain the travel time of the seismic signal wavefront of the original seismic data.
[0094] The filtering module 308 is used to filter the original seismic data based on the wavefront travel time of the seismic signal to obtain an enhanced seismic signal.
[0095] This invention provides a high-precision seismic signal enhancement device. It acquires raw seismic data, constructs a quadratic equation characterizing the wavefront travel time of the seismic signal, solves the equation based on the raw seismic data and the quadratic equation to obtain the wavefront travel time of the original seismic signal, and filters the raw seismic data based on the wavefront travel time to obtain the enhanced seismic signal. This device is suitable for enhancing seismic signals with a certain curvature on the wavefront. Compared with existing seismic signal enhancement technologies, it improves the accuracy of seismic signal prediction and effectively reduces the loss of effective information during the seismic signal enhancement process, thereby achieving high-precision seismic signal enhancement.
[0096] The aforementioned raw seismic data can be considered as two-dimensional data with time as the ordinate and trace number as the abscissa; based on this, the expression for the aforementioned quadratic equation can be: ;in, t For time, x As a Taoist name, a , b and c Let be the coefficients of the quadratic equation; the quadratic equation can satisfy the following constraint equations: ;in, P For amplitude, C It is a constant.
[0097] The aforementioned solution module 306 can also be used to: calculate the first and second derivatives of the constraint equations respectively; and, based on the original seismic data and the results of calculating the first and second derivatives of the constraint equations, solve the coefficients of the quadratic equations to obtain the travel time of the seismic signal wavefront of the original seismic data.
[0098] The solution module 306 described above can also be used to: combine the results of finding the first and second derivatives of the constraint equations to obtain a system of simultaneous equations; solve the system of simultaneous equations to obtain a system of partial differential equations; based on the original seismic data, use the finite difference method to solve the partial differential parts of the system of partial differential equations to obtain the coefficient values of the quadratic equations; and based on the quadratic equations and the obtained coefficient values, determine the travel time of the seismic signal wavefront of the original seismic data.
[0099] The first derivative of the above constraint equations can be obtained as follows:
[0100]
[0101] The result of taking the second derivative of the above constraint equations can be:
[0102] .
[0103] The expression for the partial differential part of the above system of partial differential equations can be:
[0104]
[0105]
[0106]
[0107]
[0108] .
[0109] The expression for the travel time of the seismic signal wavefront in the above raw seismic data can be: ;in, When the seismic signal wave travels ahead, i The change value of the channel number, i Take the integer part.
[0110] The above filtering can be a moving average filter; based on this, the enhanced seismic signal expression can be: ;in, For the enhanced seismic signal, k is the radius of the filtering window.
[0111] The high-precision seismic signal enhancement device provided in this embodiment of the invention has the same implementation principle and technical effect as the aforementioned high-precision seismic signal enhancement method embodiment. For the sake of brevity, any parts not mentioned in the high-precision seismic signal enhancement device embodiment can be referred to the corresponding content in the aforementioned high-precision seismic signal enhancement method embodiment.
[0112] This invention also provides an electronic device, such as... Figure 4 The diagram shows the structure of the electronic device, which includes a processor 41 and a memory 40. The memory 40 stores computer-executable instructions that can be executed by the processor 41. The processor 41 executes the computer-executable instructions to implement the above-mentioned high-precision seismic signal enhancement method.
[0113] exist Figure 4 In the illustrated embodiment, the electronic device further includes a bus 42 and a communication interface 43, wherein the processor 41, the communication interface 43, and the memory 40 are connected via the bus 42.
[0114] The memory 40 may include high-speed random access memory (RAM) and may also include non-volatile memory, such as at least one disk storage device. Communication between this system network element and at least one other network element is achieved through at least one communication interface 43 (which can be wired or wireless), such as the Internet, wide area network, local area network, metropolitan area network, etc. The bus 42 may be an ISA (Industry Standard Architecture) bus, a PCI (Peripheral Component Interconnect) bus, or an EISA (Extended Industry Standard Architecture) bus, etc. The bus 42 can be divided into an address bus, a data bus, a control bus, etc. For ease of representation, Figure 4 The symbol is represented by a single double-headed arrow, but this does not mean that there is only one bus or one type of bus.
[0115] Processor 41 may be an integrated circuit chip with signal processing capabilities. In implementation, each step of the above method can be completed by the integrated logic circuitry in the hardware of processor 41 or by software instructions. Processor 41 can be a general-purpose processor, including a Central Processing Unit (CPU), a Network Processor (NP), etc.; it can also be a Digital Signal Processor (DSP), an Application Specific Integrated Circuit (ASIC), a Field-Programmable Gate Array (FPGA), or other programmable logic devices, discrete gate or transistor logic devices, or discrete hardware components. The general-purpose processor can be a microprocessor or any conventional processor. The steps of the method disclosed in the embodiments of this invention can be directly implemented by a hardware decoding processor, or implemented by a combination of hardware and software modules in the decoding processor. The software modules can reside in random access memory, flash memory, read-only memory, programmable read-only memory, electrically erasable programmable memory, registers, or other mature storage media in the art. The storage medium is located in the memory. The processor 41 reads the information in the memory and, in conjunction with its hardware, completes the steps of the high-precision seismic signal enhancement method described in the aforementioned embodiment.
[0116] Unless otherwise specifically stated, the relative steps, numerical expressions, and values of the components and steps described in these embodiments do not limit the scope of the invention.
[0117] If the aforementioned functions are implemented as software functional units and sold or used as independent products, they can be stored in a processor-executable, non-volatile, computer-readable storage medium. Based on this understanding, the technical solution of this invention, essentially, or the part that contributes to the prior art, or a portion of the technical solution, can be embodied in the form of a software product. This computer software product is stored in a storage medium and includes several instructions to cause a computer device (which may be a personal computer, server, or network device, etc.) to execute all or part of the steps of the methods described in the various embodiments of this invention. The aforementioned storage medium includes various media capable of storing program code, such as USB flash drives, portable hard drives, read-only memory (ROM), random access memory (RAM), magnetic disks, or optical disks.
[0118] In the description of this invention, it should be noted that the terms "center," "upper," "lower," "left," "right," "vertical," "horizontal," "inner," and "outer," etc., indicate the orientation or positional relationship based on the orientation or positional relationship shown in the accompanying drawings. They are used only for the convenience of describing the invention and for simplifying the description, and do not indicate or imply that the device or element referred to must have a specific orientation, or be constructed and operated in a specific orientation. Therefore, they should not be construed as limitations on the invention. Furthermore, the terms "first," "second," and "third" are used for descriptive purposes only and should not be construed as indicating or implying relative importance.
[0119] Finally, it should be noted that the above-described embodiments are merely specific implementations of the present invention, used to illustrate the technical solutions of the present invention, and not to limit it. The scope of protection of the present invention is not limited thereto. Although the present invention has been described in detail with reference to the foregoing embodiments, those skilled in the art should understand that any person skilled in the art can still modify or easily conceive of changes to the technical solutions described in the foregoing embodiments within the technical scope disclosed in the present invention, or make equivalent substitutions for some of the technical features; and these modifications, changes, or substitutions do not cause the essence of the corresponding technical solutions to deviate from the spirit and scope of the technical solutions of the embodiments of the present invention, and should all be covered within the scope of protection of the present invention. Therefore, the scope of protection of the present invention should be determined by the scope of the claims.
Claims
1. A high-precision seismic signal enhancement method, characterized in that, include: Obtain raw earthquake data; Construct a quadratic equation characterizing the travel time of the wavefront of seismic signals; The travel time of the seismic signal wavefront from the original seismic data is obtained by solving the original seismic equation based on the original seismic data. Based on the wavefront travel time of the earthquake signal, the original earthquake data is filtered to obtain an enhanced earthquake signal. The original seismic data is two-dimensional data with time as the vertical axis and track number as the horizontal axis; The expression for the quadratic equation is: in, t For time, x As a Taoist name, a , b and c The coefficients of the quadratic equation are denoted as . The quadratic equation satisfies the following constraint equations: in, P For amplitude, C It is a constant; The process of obtaining the wavefront travel time of the seismic signal from the original seismic data by solving the quadratic equation includes: calculating the first and second derivatives of the constraint equations; simultaneously solving the first and second derivatives of the constraint equations to obtain a system of simultaneous equations; solving the system of simultaneous equations to obtain a system of partial differential equations; using the finite difference method to solve the partial differential parts of the system of partial differential equations based on the original seismic data to obtain the coefficients of the quadratic equation; and determining the wavefront travel time of the seismic signal from the original seismic data based on the quadratic equation and the obtained coefficients. The expression for the wavefront travel time of the seismic signal in the original seismic data is: in, When the seismic signal wave travels ahead, i The change value of the channel number, i Take the integer part.
2. The high-precision seismic signal enhancement method according to claim 1, characterized in that, The result of taking the first derivative of the constraint equation is: The result of taking the second derivative of the constraint equation is: 。 3. The high-precision seismic signal enhancement method according to claim 2, characterized in that, The expression for the partial differential part of the system of partial differential equations is as follows: 。 4. The high-precision seismic signal enhancement method according to claim 3, characterized in that, The filtering method is a moving average filtering method; the enhanced seismic signal expression is: in, For the enhanced seismic signal, k Let be the radius of the filtering window.
5. A high-precision seismic signal enhancement device, characterized in that, include: The acquisition module is used to acquire raw seismic data; The module is used to construct a quadratic equation characterizing the travel time of the wavefront of a seismic signal. The solution module is used to solve the original seismic data and the quadratic equation to obtain the wavefront travel time of the seismic signal from the original seismic data. The filtering module is used to filter the original seismic data based on the wavefront travel time of the seismic signal to obtain an enhanced seismic signal. The original seismic data is two-dimensional data with time as the vertical axis and track number as the horizontal axis; The expression for the quadratic equation is: in, t For time, x As a Taoist name, a , b and c The coefficients of the quadratic equation are denoted as . The quadratic equation satisfies the following constraint equations: in, P For amplitude, C It is a constant; The solution module is further configured to: calculate the first and second derivatives of the constraint equations; combine the results of calculating the first and second derivatives of the constraint equations to obtain a system of simultaneous equations; solve the system of simultaneous equations to obtain a system of partial differential equations; based on the original seismic data, solve the partial differential parts of the system of partial differential equations using the finite difference method to obtain the coefficient values of the quadratic equations; and determine the travel time of the seismic signal wavefront from the original seismic data based on the quadratic equations and the obtained coefficient values. The expression for the wavefront travel time of the seismic signal in the original seismic data is: in, When the seismic signal wave travels ahead, i The change value of the channel number, i Take the integer part.
6. An electronic device, characterized in that, It includes a processor and a memory, the memory storing computer-executable instructions that can be executed by the processor, the processor executing the computer-executable instructions to implement the high-precision seismic signal enhancement method according to any one of claims 1 to 4.