Seismic data matching method and device, storage medium and electronic equipment

Abstract

The application relates to the technical field of physical exploration, in particular to a seismic data matching method and device, a storage medium and an electronic device, which comprises the following steps: determining a plurality of controlled source seismic traces and a plurality of explosive source seismic traces with the same number according to the construction design of a target area; performing phase rotation on each controlled source seismic trace according to a preset angle interval; respectively performing structure similarity evaluation on the signal of each explosive source seismic trace and the signal of the controlled source seismic trace corresponding to the explosive source seismic trace, so as to obtain the best phase rotation angle of each controlled source seismic trace; obtaining the best average phase rotation angle of the target area by averaging; and performing phase rotation on all controlled sources in the target area, so as to obtain matched controlled source data. The application can obtain a more accurate phase rotation angle through quantitative phase rotation and similarity evaluation, can better realize in-phase superposition of the controlled source and the explosive source, and effectively improves the seismic imaging quality.

Description

Technical Field

[0001] This application relates to the field of geophysical exploration technology, and in particular to a seismic data matching method, apparatus, storage medium, and electronic device. Background Technology

[0002] In densely populated areas such as beaches, deserts, lakes, and towns, where surface conditions limit the effectiveness of seismic activity, it is often necessary to alternate between controlled seismic sources and explosive seismic sources. Due to the different source types, the records from explosive and controlled seismic sources acquired by the geophones differ in phase, frequency, and energy, making it impossible to superimpose the two sets of data in phase. To solve the data matching problem between controlled and explosive seismic sources, commonly used methods can be divided into two main categories: one is to perform matching operations on the data from controlled and explosive seismic sources to obtain matching operators; the other is to perform polarity reversal or phase adjustment on the controlled seismic sources.

[0003] The common practice for obtaining matching operators is to first stack the repeating segments of data from two different sources, and then use the stacked data to calculate the matching operator. However, the stacked data is often affected by preprocessing steps such as static correction, velocity picking, and deconvolution, resulting in a generally unsatisfactory matching operator. Furthermore, when the signal-to-noise ratio of the seismic data is too low, the obtained matching operator is unlikely to meet processing requirements, significantly impacting subsequent seismic imaging. Methods using controlled-source polarity reversal or phase adjustment also have limitations. The phase difference between a controlled-source and an explosive source is not always exactly 180 degrees, making direct polarity reversal unreasonable. Additionally, the determination of the phase adjustment angle in phase adjustment methods lacks quantitative analytical tools, thus limiting the effectiveness of this approach in current applications. Summary of the Invention

[0004] To address the aforementioned issues, this application provides a seismic data matching method, apparatus, storage medium, and electronic device, which solves the technical problem of poor imaging quality caused by inaccurate seismic data matching in related technologies.

[0005] In a first aspect, this application provides a seismic data matching method, the method comprising:

[0006] Based on the construction design of the target area, determine the same number of multiple controllable source seismic channels and multiple explosive source seismic channels.

[0007] For each controllable source seismic trace, phase rotation is performed within the range of 0 to 180 degrees according to a preset angle interval to obtain the rotation results of each controllable source seismic trace at different phase rotation angles;

[0008] Structural similarity evaluation was performed on the seismic traces of each explosive source and the corresponding controllable source seismic traces with different phase rotation angles to obtain the optimal phase rotation angle for each controllable source seismic trace.

[0009] The average of the optimal phase rotation angles for each controllable source seismic trace is calculated to obtain the optimal average phase rotation angle for the target area.

[0010] Based on the optimal average phase rotation angle, all controllable seismic sources within the target area are phase rotated to obtain matched controllable seismic source data.

[0011] In some embodiments, determining the controllable source seismic path and the explosive source seismic path based on the construction design of the target area includes:

[0012] The overlapping or adjacent areas of controllable seismic sources and explosive seismic sources are determined according to the construction design of the work area;

[0013] Based on the overlapping or adjacent areas of the controllable seismic source and the explosive seismic source, the same number of multiple controllable seismic source seismic channels and multiple explosive seismic channels are determined.

[0014] In some embodiments, the preset angle interval is 5 degrees.

[0015] In some embodiments, the step of performing phase rotation on each controllable source seismic trace within a range of 0 to 180 degrees at preset angle intervals to obtain rotation results for each controllable source seismic trace at different phase rotation angles includes:

[0016] According to the formula s(t)=x(t)e jφ Phase rotation is performed on each controllable source seismic trace to obtain the rotation result s(t) for each controllable source seismic trace at different phase rotation angles;

[0017] Where s(t) is the controlled-source seismic trace after phase rotation; x(t) is the controlled-source seismic trace before phase rotation; e is the natural index; j is the imaginary unit; and φ is the phase rotation angle.

[0018] In some embodiments, the structural similarity evaluation is performed on the signals of the seismic traces of each explosive source and the signals of the corresponding controlled-source seismic traces with different phase rotation angles to obtain the optimal phase rotation angle for each controlled-source seismic trace, including:

[0019] Select the calculation time window at the first arrival position of the controlled-source seismic trace and the explosive-source seismic trace;

[0020] Within the calculation window, the structural similarity of the seismic traces of each explosive source and the corresponding controllable source seismic traces with different phase rotation angles is evaluated to obtain the optimal phase rotation angle for each controllable source seismic trace.

[0021] In some embodiments, the structural similarity evaluation of the seismic traces of each explosive source and the corresponding controlled-source seismic traces with different phase rotation angles within the calculation window, to obtain the optimal phase rotation angle for each controlled-source seismic trace, includes:

[0022] Based on the average energy of the explosive source seismic trace signals and the controlled source seismic trace signals used in the calculation within the calculation time window, the energy variance of the explosive source seismic trace signals and the controlled source seismic trace signals used in the calculation within the calculation time window is obtained.

[0023] Based on the energy variance of the data of the explosive source seismic trace signal and the controlled source seismic trace signal within the calculation time window, the correlation coefficient of the data of the explosive source seismic trace signal and the controlled source seismic trace signal within the calculation time window is obtained.

[0024] Based on the correlation coefficient and energy variance of the seismic trace signals from explosive sources and controlled sources within the calculation window, the signal structure similarity is calculated to obtain the similarity results.

[0025] The phase rotation angle of the controllable source seismic trace signal with the highest similarity to the explosive source seismic trace signal is taken as the optimal phase rotation angle of the controllable source seismic trace.

[0026] In some embodiments, the step of calculating signal structure similarity based on the correlation coefficient and energy variance of the data between the explosive source seismic trace signal and the controlled source seismic trace signal within the calculation time window, and obtaining the similarity result, includes:

[0027] According to the formula The signal structure similarity (MSS) is calculated.

[0028] Where C is a very small positive number to avoid a denominator of 0, x is the explosive source seismic trace signal, y is the controlled source seismic trace signal, θx is the energy variance of the explosive source seismic trace signal within the calculation window, θy is the energy variance of the controlled source seismic trace signal within the calculation window, and θxy is the correlation coefficient between the explosive source seismic trace signal and the controlled source seismic trace signal within the calculation window.

[0029] Secondly, a seismic data matching device, the device comprising:

[0030] The determination unit is used to determine the same number of multiple controllable source seismic channels and multiple explosive source seismic channels according to the construction design of the target area.

[0031] The phase rotation unit is used to rotate the phase of each controllable source seismic trace within a range of 0 to 180 degrees according to a preset angle interval, so as to obtain the rotation results of each controllable source seismic trace at different phase rotation angles.

[0032] The evaluation unit is used to evaluate the structural similarity of the seismic traces of each explosive source and the corresponding controllable source seismic traces with different phase rotation angles, and to obtain the optimal phase rotation angle for each controllable source seismic trace.

[0033] The calculation unit is used to calculate the average of the optimal phase rotation angle of each controllable source seismic trace, and finally obtain the optimal average phase rotation angle of the target area.

[0034] The matching unit is used to perform phase rotation on all controllable seismic sources within the target area according to the optimal average phase rotation angle, so as to obtain matched controllable seismic source data.

[0035] Thirdly, a storage medium storing a computer program that can be executed by one or more processors to implement the seismic data matching method as described in the first aspect above.

[0036] Fourthly, an electronic device includes a memory and a processor, wherein a computer program is stored on the memory and the processor is communicatively connected to each other, and the computer program, when executed by the processor, performs the seismic data matching method as described in the first aspect above.

[0037] This application provides a seismic data matching method, apparatus, storage medium, and electronic device, comprising: determining the same number of multiple controlled-source seismic traces and multiple explosive-source seismic traces based on the construction design of the target area; performing phase rotation on each controlled-source seismic trace within a range of 0 to 180 degrees at preset angle intervals to obtain rotation results for each controlled-source seismic trace at different phase rotation angles; performing structural similarity evaluation on the signals of each explosive-source seismic trace and the signals of its corresponding controlled-source seismic traces at different phase rotation angles to obtain the optimal phase rotation angle for each controlled-source seismic trace; averaging the optimal phase rotation angles of each controlled-source seismic trace to obtain the optimal average phase rotation angle for the target area; and performing phase rotation on all controlled sources within the target area based on the optimal average phase rotation angle to obtain matched controlled-source data. This application, through quantitative phase rotation and similarity evaluation, obtains more accurate phase rotation angles, enabling better in-phase superposition of controlled-source and explosive-source seismic data and effectively improving seismic imaging quality. Attached Figure Description

[0038] To more clearly illustrate the technical solutions in the embodiments of the present invention or the prior art, the drawings used in the description of the embodiments or the prior art will be briefly introduced below. Obviously, the drawings described below are only embodiments of the present invention. For those skilled in the art, other drawings can be obtained based on the provided drawings without creative effort.

[0039] Figure 1 A flowchart illustrating a seismic data matching method provided in an embodiment of this application;

[0040] Figure 2 A schematic diagram of seismic records of an explosive source in a target area provided in an embodiment of this application;

[0041] Figure 3 A schematic diagram of seismic records of a controllable source in the target area provided for an embodiment of this application;

[0042] Figure 4 Comparison diagram of controllable source seismic traces with partial phase rotation angles and corresponding explosive source seismic traces provided for embodiments of this application;

[0043] Figure 5 A statistical graph showing the similarity of different phase rotation angles and signal structures of a controllable source seismic trace within a selected time window, provided for embodiments of this application.

[0044] Figure 6 Superimposed profile images of controllable source data before and after matching, provided for embodiments of this application;

[0045] Figure 7 This is a schematic diagram of the structure of a seismic data matching device provided in an embodiment of this application;

[0046] Figure 8 This is a connection block diagram of an electronic device provided in an embodiment of this application. Detailed Implementation

[0047] The following detailed description of the embodiments of this application, in conjunction with the accompanying drawings, will provide a thorough understanding of how this application uses technical means to solve technical problems and achieve corresponding technical effects, enabling its implementation. The embodiments of this application and the various features within them can be combined with each other without conflict, and all resulting technical solutions are within the protection scope of this application.

[0048] As the background technology indicates, the common practice for obtaining matching operators is to first stack the repeating segments of data from two different sources, and then use the stacked data to calculate the matching operator. However, the stacked data is usually affected by preprocessing steps such as static correction, velocity picking, and deconvolution, resulting in a generally unsatisfactory matching operator. Furthermore, when the signal-to-noise ratio of the seismic data is too low, the obtained matching operator is unlikely to meet processing requirements, significantly impacting subsequent seismic imaging. Methods using controlled-source polarity reversal or phase adjustment also have limitations. The phase difference between a controlled-source and an explosive source is generally not exactly 180 degrees, making direct polarity reversal unreasonable. Moreover, the determination of the phase adjustment angle in phase adjustment methods lacks quantitative analytical tools, thus this method has not yet achieved satisfactory application results.

[0049] In view of this, this application provides a seismic data matching method, apparatus, storage medium, and electronic device, which solves the technical problem of poor imaging quality caused by inaccurate seismic data matching in related technologies.

[0050] Example 1

[0051] Figure 1 This is a flowchart illustrating a seismic data matching method provided in an embodiment of this application, as shown below. Figure 1 As shown, this method includes:

[0052] S101. Determine the same number of multiple controllable source seismic channels and multiple explosive source seismic channels based on the construction design of the target area.

[0053] S102. For each controllable source seismic trace, perform phase rotation within the range of 0 to 180 degrees according to a preset angle interval to obtain the rotation results of each controllable source seismic trace at different phase rotation angles.

[0054] S103. For the signal of the seismic trace of each explosive source and the signal of the corresponding controllable source seismic trace with different phase rotation angles, perform structural similarity evaluation to obtain the optimal phase rotation angle of each controllable source seismic trace.

[0055] S104. Calculate the average of the optimal phase rotation angles for each controllable source seismic trace to obtain the optimal average phase rotation angle for the target area.

[0056] S105. Based on the optimal average phase rotation angle, perform phase rotation on all controllable seismic sources within the target area to obtain matched controllable seismic source data.

[0057] In some embodiments, determining the controllable source seismic path and the explosive source seismic path based on the construction design of the target area includes:

[0058] The overlapping or adjacent areas of controllable seismic sources and explosive seismic sources are determined according to the construction design of the work area;

[0059] Based on the overlapping or adjacent areas of the controllable seismic source and the explosive seismic source, the same number of multiple controllable seismic source seismic channels and multiple explosive seismic channels are determined.

[0060] It should be noted that, as Figure 2 The image shows a schematic diagram of the seismic record of an explosive source in a target area. Figure 3 The image shows a schematic diagram of seismic records from controllable sources in the target area. Figure 2 and Figure 3 It can be seen that there is a significant phase difference between the two records at the initial arrival position. Assuming that there are approximately 500 controlled source excitation points in this region, an average of 10 controlled source seismic traces were selected from the region as matching calculation data. At the same time, 10 corresponding seismic traces were also selected from the explosive source data of adjacent or repeated segments. The number of seismic traces selected can be adjusted according to needs.

[0061] In some embodiments, the preset angle interval is 5 degrees.

[0062] It should be noted that the preset angle interval is preferably 5 degrees, but it can be adjusted as needed.

[0063] In some embodiments, the step of performing phase rotation on each controllable source seismic trace within a range of 0 to 180 degrees at preset angle intervals to obtain rotation results for each controllable source seismic trace at different phase rotation angles includes:

[0064] According to the formula s(t)=x(t)e jφPhase rotation is performed on each controllable source seismic trace to obtain the rotation result s(t) for each controllable source seismic trace at different phase rotation angles;

[0065] Where s(t) is the controlled-source seismic trace after phase rotation; x(t) is the controlled-source seismic trace before phase rotation; e is the natural index; j is the imaginary unit; and φ is the phase rotation angle.

[0066] It should be noted that, as Figure 4 The image shows a comparison of the seismic traces of a controlled source and the corresponding explosive source with partial phase rotation angles. It can be seen that the morphological similarity between the controlled source signal and the explosive source signal varies with different phase angles. The rectangular boxes in the image represent the selected matching calculation time windows.

[0067] In some embodiments, the structural similarity evaluation is performed on the signals of the seismic traces of each explosive source and the signals of the corresponding controlled-source seismic traces with different phase rotation angles to obtain the optimal phase rotation angle for each controlled-source seismic trace, including:

[0068] Select the calculation time window at the first arrival position of the controlled-source seismic trace and the explosive-source seismic trace;

[0069] Within the calculation window, the structural similarity of the seismic traces of each explosive source and the corresponding controllable source seismic traces with different phase rotation angles is evaluated to obtain the optimal phase rotation angle for each controllable source seismic trace.

[0070] It should be noted that by utilizing the first arrival information (i.e., first arrival position) of seismic data, the problem of inaccurate data matching due to low signal-to-noise ratio is avoided to the greatest extent. Compared with the conventional method of obtaining matching operators using stacked profiles, the method of this invention has stronger anti-interference ability.

[0071] It should be further explained that, such as Figure 5 The figure shown is a statistical graph of the similarity between different phase rotation angles and signal structure of a controllable source seismic trace within a selected time window. Figure 5 The results show that the structural similarity between the two seismic traces is highest at approximately a phase rotation of 130 degrees. Therefore, the optimal phase rotation angle for the seismic traces calculated in this study is 130 degrees. Similarly, signal structure similarity calculations were performed on other selected seismic traces, and the average of all the obtained optimal phase rotation angles was calculated. The final calculated optimal average phase rotation angle for this region is 135 degrees.

[0072] In some embodiments, the structural similarity evaluation of the seismic traces of each explosive source and the corresponding controlled-source seismic traces with different phase rotation angles within the calculation window, to obtain the optimal phase rotation angle for each controlled-source seismic trace, includes:

[0073] Based on the average energy of the explosive source seismic trace signals and the controlled source seismic trace signals used in the calculation within the calculation time window, the energy variance of the explosive source seismic trace signals and the controlled source seismic trace signals used in the calculation within the calculation time window is obtained.

[0074] Based on the energy variance of the data of the explosive source seismic trace signal and the controlled source seismic trace signal within the calculation time window, the correlation coefficient of the data of the explosive source seismic trace signal and the controlled source seismic trace signal within the calculation time window is obtained.

[0075] Based on the correlation coefficient and energy variance of the seismic trace signals from explosive sources and controlled sources within the calculation window, the signal structure similarity is calculated to obtain the similarity results.

[0076] The phase rotation angle of the controllable source seismic trace signal with the highest similarity to the explosive source seismic trace signal is taken as the optimal phase rotation angle of the controllable source seismic trace.

[0077] It should be noted that performing phase rotation on a signal is essentially filtering it, and the specific formula is as follows:

[0078] s(t)=x(t)e jφ (1)

[0079] In the formula: s(t) is the phase-rotated seismic record; x(t) is the original seismic record; e is the natural index; j is the imaginary unit; φ is the phase rotation angle.

[0080] Previously, phase adjustment of controlled seismic sources typically involved quality control based on the stacked profiles after phase adjustment, selecting the optimal phase adjustment angle. This method lacks quantitative analysis, and human subjectivity significantly influences the results. To address this issue, this invention proposes a structural similarity evaluation method that can more accurately determine the optimal phase adjustment angle for controlled seismic sources. Based on the attribute definition used by Chopra et al. for describing seismic texture attributes, μ can be defined. x and μ y The expressions for the mean energy values ​​of data x and y within the corresponding time window are as follows:

[0081]

[0082]

[0083] In the formula: x k and y k Let θx and θy be the amplitude values ​​at the k-th sample point within the time window for data x and data y, respectively; N is the total number of samples within the time window. Define θx and θy as the energy variances of the signals in data x and data y within the corresponding time window, respectively, with the following expressions:

[0084]

[0085]

[0086] The physical meaning of the above formula is a measure of the local signal change within a time window. Define θxy as the correlation coefficient between data x and data y within the corresponding time window, expressed as:

[0087]

[0088] By analyzing the correlation of data, the similarity of signals can be reflected. Therefore, the mean structural similarity (MSS) of seismic data x and y within a time window can be further defined:

[0089] In the formula, C is a very small positive number to avoid the denominator being zero. MSS can be used to compare the similarity of two data sets within the same time window. The parameter used for similarity evaluation must also satisfy that MSS takes a value between 0 and 1, that is, MSS = 1 when the two data sets are completely similar, and MSS = 0 when the two data sets are completely dissimilar. Usually, the boundedness of MSS can be satisfied by adjusting the parameter C.

[0090] In some embodiments, the step of calculating signal structure similarity based on the correlation coefficient and energy variance of the data between the explosive source seismic trace signal and the controlled source seismic trace signal within the calculation time window, and obtaining the similarity result, includes:

[0091] According to the formula The signal structure similarity (MSS) is calculated.

[0092] Where C is a very small positive number to avoid a denominator of 0, x is the explosive source seismic trace signal, y is the controlled source seismic trace signal, θx is the energy variance of the explosive source seismic trace signal within the calculation window, θy is the energy variance of the controlled source seismic trace signal within the calculation window, and θxy is the correlation coefficient between the explosive source seismic trace signal and the controlled source seismic trace signal within the calculation window.

[0093] Furthermore, to verify the effectiveness of the invention, data from adjacent segments of the controllable seismic source and the explosive seismic source were superimposed. For example... Figure 6The image shows the stacked profiles before and after matching controllable source data. The three sub-images, from left to right, represent: the stacked profile before matching, the stacked profile after matching using conventional matched filtering, and the stacked profile after matching using the method of this invention. Figure 6 As can be seen from the area indicated by the arrows and rectangles, compared with conventional matched filtering methods, the method provided by this invention has better continuity of the in-phase axis in the superimposed profile after matching, achieving in-phase superposition, thus proving the effectiveness of the method of this invention.

[0094] In summary, this application provides a seismic data matching method, comprising: determining the same number of controlled-source seismic traces and multiple explosive-source seismic traces based on the construction design of the target area; performing phase rotation on each controlled-source seismic trace within a range of 0 to 180 degrees at preset angle intervals to obtain rotation results for each controlled-source seismic trace at different phase rotation angles; performing structural similarity evaluation on the signals of each explosive-source seismic trace and the signals of its corresponding controlled-source seismic traces at different phase rotation angles to obtain the optimal phase rotation angle for each controlled-source seismic trace; averaging the optimal phase rotation angles for each controlled-source seismic trace to obtain the optimal average phase rotation angle for the target area; and performing phase rotation on all controlled sources within the target area based on the optimal average phase rotation angle to obtain matched controlled-source data. This application, through quantitative phase rotation and similarity evaluation, obtains a more accurate phase rotation angle, enabling better in-phase superposition of controlled-source and explosive-source seismic data and effectively improving seismic imaging quality.

[0095] Example 2

[0096] Based on the earthquake data matching method disclosed in the above embodiments of the present invention, Figure 7 The seismic data matching device that applies this seismic data matching method is specifically disclosed.

[0097] like Figure 7 As shown in the figure, an embodiment of the present invention discloses a seismic data matching device, which includes:

[0098] Unit 701 is used to determine the same number of multiple controllable source seismic channels and multiple explosive source seismic channels according to the construction design of the target area.

[0099] The phase rotation unit 702 is used to rotate the phase of each controllable source seismic trace within the range of 0 to 180 degrees according to a preset angle interval, so as to obtain the rotation results of each controllable source seismic trace at different phase rotation angles.

[0100] Evaluation unit 703 is used to evaluate the structural similarity of the seismic traces of each explosive source and the corresponding controllable source seismic traces with different phase rotation angles, and to obtain the optimal phase rotation angle of each controllable source seismic trace.

[0101] The calculation unit 704 is used to calculate the average of the optimal phase rotation angle of each controllable source seismic trace, and finally obtain the optimal average phase rotation angle of the target area.

[0102] The matching unit 705 is used to perform phase rotation on all controllable seismic sources within the target area according to the optimal average phase rotation angle, so as to obtain matched controllable seismic source data.

[0103] The specific working process of each unit in the seismic data matching device disclosed in the above embodiments of the present invention can be found in the corresponding content of the seismic data matching method disclosed in the above embodiments of the present invention, and will not be repeated here.

[0104] In summary, this application provides a seismic data matching device, comprising: determining the same number of multiple controllable source seismic traces and multiple explosive source seismic traces according to the construction design of the target area; performing phase rotation on each controllable source seismic trace within a range of 0 to 180 degrees at preset angle intervals to obtain rotation results for each controllable source seismic trace at different phase rotation angles; performing structural similarity evaluation on the signals of each explosive source seismic trace and the signals of the corresponding controllable source seismic traces at different phase rotation angles to obtain the optimal phase rotation angle for each controllable source seismic trace; averaging the optimal phase rotation angles of each controllable source seismic trace to obtain the optimal average phase rotation angle for the target area; and performing phase rotation on all controllable sources within the target area based on the optimal average phase rotation angle to obtain matched controllable source data. This application, through quantitative phase rotation and similarity evaluation, obtains a more accurate phase rotation angle, enabling better in-phase superposition of controllable and explosive sources, and effectively improving seismic imaging quality.

[0105] Example 3

[0106] This embodiment also provides a computer-readable storage medium, such as flash memory, hard disk, multimedia card, card-type memory (e.g., SD or DX memory), random access memory (RAM), static random access memory (SRAM), read-only memory (ROM), electrically erasable programmable read-only memory (EEPROM), programmable read-only memory (PROM), magnetic memory, disk, optical disk, server, App application store, etc., which stores a computer program. When the computer program is executed by a processor, it can implement the method steps as described in Embodiment 1. This embodiment will not repeat the description here.

[0107] Example 4

[0108] Figure 8 A connection block diagram of an electronic device 800 provided in an embodiment of this application is shown below. Figure 8 As shown, the electronic device 800 may include: a processor 801, a memory 802, a multimedia component 803, an input / output (I / O) interface 804, and a communication component 805.

[0109] The processor 801 is used to execute all or part of the steps in the seismic data matching method as described in Embodiment 1. The memory 802 is used to store various types of data, which may include, for example, instructions for any application or method in an electronic device, as well as application-related data.

[0110] The processor 801 may be implemented as an application-specific integrated circuit (ASIC), a digital signal processor (DSP), a digital signal processing device (DSPD), a programmable logic device (PLD), a field-programmable gate array (FPGA), a controller, a microcontroller, a microprocessor, or other electronic components, and is used to execute the seismic data matching method in Embodiment 1 above.

[0111] The memory 802 can be implemented by any type of volatile or non-volatile storage device or a combination thereof, such as static random access memory (SRAM), electrically erasable programmable read-only memory (EEPROM), erasable programmable read-only memory (EPROM), programmable read-only memory (PROM), read-only memory (ROM), magnetic storage, flash memory, magnetic disk, or optical disk.

[0112] Multimedia component 803 may include a screen, which may be a touchscreen, and an audio component for outputting and / or inputting audio signals. For example, the audio component may include a microphone for receiving external audio signals. The received audio signals may be further stored in memory or transmitted via a communication component. The audio component also includes at least one speaker for outputting audio signals.

[0113] I / O interface 804 provides an interface between processor 801 and other interface modules, such as keyboards, mice, and buttons. These buttons can be virtual or physical buttons.

[0114] The communication component 805 is used for wired or wireless communication between the electronic device 800 and other devices. Wireless communication includes, for example, Wi-Fi, Bluetooth, Near Field Communication (NFC), 2G, 3G, or 4G, or a combination thereof. Therefore, the corresponding communication component 805 may include: a Wi-Fi module, a Bluetooth module, or an NFC module.

[0115] In summary, this application provides a seismic data matching method, apparatus, storage medium, and electronic device. The method includes: determining the same number of multiple controlled-source seismic traces and multiple explosive-source seismic traces based on the construction design of the target area; performing phase rotation on each controlled-source seismic trace within a range of 0 to 180 degrees at preset angle intervals to obtain rotation results for each controlled-source seismic trace at different phase rotation angles; performing structural similarity evaluation on the signals of each explosive-source seismic trace and the signals of its corresponding controlled-source seismic traces at different phase rotation angles to obtain the optimal phase rotation angle for each controlled-source seismic trace; averaging the optimal phase rotation angles for each controlled-source seismic trace to obtain the optimal average phase rotation angle for the target area; and performing phase rotation on all controlled sources within the target area based on the optimal average phase rotation angle to obtain matched controlled-source data. This application, through quantitative phase rotation and similarity evaluation, obtains a more accurate phase rotation angle, enabling better in-phase superposition of controlled-source and explosive-source seismic data and effectively improving seismic imaging quality.

[0116] In the several embodiments provided in this application, it should be understood that the disclosed methods can also be implemented in other ways. The method embodiments described above are merely illustrative.

[0117] It should be noted that, in this document, the terms "comprising," "including," or any other variations thereof are intended to cover non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements includes not only those elements but also other elements not expressly listed, or elements inherent to such a process, method, article, or apparatus. Unless otherwise specified, an element defined by the phrase "comprising one..." does not exclude the presence of other identical elements in the process, method, article, or apparatus that includes that element.

[0118] Although the embodiments disclosed in this application are as described above, the above content is merely for the purpose of facilitating understanding of this application and is not intended to limit this application. Any person skilled in the art to which this application pertains may make any modifications and changes in the form and details of the implementation without departing from the spirit and scope disclosed in this application; however, the scope of patent protection of this application shall still be determined by the scope defined in the appended claims.

Claims

1. A method of seismic data matching, characterized by, The method comprises: determining a plurality of controllable source seismic traces and a plurality of explosive source seismic traces according to a construction design of a target region; performing phase rotation on each controllable source seismic trace at a preset angle interval within a range of 0 to 180 degrees to obtain a rotation result of each controllable source seismic trace at different phase rotation angles; respectively performing structural similarity evaluation on a signal of each explosive source seismic trace and a signal of a controllable source seismic trace corresponding to the explosive source seismic trace to obtain an optimal phase rotation angle of each controllable source seismic trace; averaging the optimal phase rotation angles of each controllable source seismic trace to obtain an optimal average phase rotation angle of the target region; performing phase rotation on all controllable sources in the target region according to the optimal average phase rotation angle to obtain matched controllable source data; wherein the structural similarity evaluation on the signal of each explosive source seismic trace and the signal of the controllable source seismic trace corresponding to the explosive source seismic trace comprises:

2. The method of claim 1, wherein, selecting a calculation time window at a first arrival position of the controllable source seismic trace and the explosive source seismic trace; and respectively performing structural similarity evaluation on a signal of each explosive source seismic trace and a signal of a controllable source seismic trace corresponding to the explosive source seismic trace within the calculation time window to obtain an optimal phase rotation angle of each controllable source seismic trace. The method comprises:

3. The method of claim 1, wherein, determining a plurality of controllable source seismic traces and a plurality of explosive source seismic traces according to a construction design of a target region; 4. The method of claim 1, wherein, determining a repeating or adjacent region of the controllable source and the explosive source according to a construction design of a work area; According to the formula Phase rotation is performed on each controllable source seismic trace to obtain each controllable source seismic trace, and the rotation results of different phase rotation angles s ( t ) Wherein, s(t) is the controllable seismic source seismic trace after phase rotation; x(t) is the controllable seismic source seismic trace before phase rotation; e is natural index; j is imaginary unit; is the phase rotation angle.

5. The method of claim 1, wherein, determining a plurality of controllable source seismic traces and a plurality of explosive source seismic traces according to the repeating or adjacent region of the controllable source and the explosive source. The preset angle interval is 5 degrees. The method comprises: determining a plurality of controllable source seismic traces and a plurality of explosive source seismic traces according to a construction design of a target region; determining a repeating or adjacent region of the controllable source and the explosive source according to a construction design of a work area; determining a plurality of controllable source seismic traces and a plurality of explosive source seismic traces according to the repeating or adjacent region of the controllable source and the explosive source. The preset angle interval is 5 degrees. The method comprises: determining a plurality of controllable source seismic traces and a plurality of explosive source seismic traces according to a construction design of a target region; determining a repeating or adjacent region of the controllable source and the explosive source according to a construction design of a work area; determining a plurality of controllable source seismic traces and a plurality of explosive source seismic traces according to the repeating or adjacent region of the controllable source and the explosive source. The preset angle interval is 5 degrees. The method comprises: determining a plurality of controllable source seismic traces and a plurality of explosive source seismic traces according to a construction design of a target region; determining a repeating or adjacent region of the controllable source and the explosive source according to a construction design of a work area; determining a plurality of controllable source seismic traces and a plurality of explosive source seismic traces according to the repeating or adjacent region of the controllable source and the explosive source. The preset angle interval is 5 degrees. The method comprises: determining a plurality of controllable source seismic traces and a plurality of explosive source seismic traces according to a construction design of a target region; determining a repeating or adjacent region of the controllable source and the explosive source according to a construction design of a work area; determining a plurality of controllable source seismic traces and a plurality of explosive source seismic traces according to the repeating or adjacent region of the controllable source and the explosive source. The preset angle interval is 5 degrees. According to the correlation coefficient and energy variance of explosive seismic trace signals and controllable seismic trace signals in a calculation time window, signal structure similarity is calculated to obtain a similarity result; A phase rotation angle of a controllable seismic trace signal with the highest similarity to the explosive seismic trace signal in the similarity result is taken as an optimal phase rotation angle of the controllable seismic trace.

6. The method of claim 5, wherein, The similarity result is obtained by calculating the signal structure similarity according to the correlation coefficient and energy variance of the explosive seismic trace signals and the controllable seismic trace signals in the calculation time window, and the similarity result includes: According to the formula The signal structure similarity MSS is calculated. Wherein, C is a minimum positive number to avoid a denominator of 0, x is the explosive seismic trace signal, y is the controllable seismic trace signal, θx is the energy variance of the explosive seismic trace signal in the calculation time window, θy is the energy variance of the controllable seismic trace signal in the calculation time window, and θxy is the correlation coefficient of the explosive seismic trace signal and the controllable seismic trace signal in the calculation time window.

7. A seismic data matching apparatus, characterized by, The device includes: A determination unit is configured to determine a plurality of controllable seismic traces and a plurality of explosive seismic traces according to a construction design of a target area; A phase rotation unit is configured to rotate each controllable seismic trace by a preset angle interval within a range of 0 to 180 degrees to obtain a rotation result of each controllable seismic trace at different phase rotation angles. An evaluation unit is configured to evaluate the signal of each explosive seismic trace and the signal of the controllable seismic trace corresponding to each explosive seismic trace at different phase rotation angles to obtain an optimal phase rotation angle of each controllable seismic trace. A calculation unit is configured to average the optimal phase rotation angles of the controllable seismic traces to obtain an optimal average phase rotation angle of the target area. A matching unit is configured to rotate all controllable seismic sources in the target area according to the optimal average phase rotation angle to obtain matched controllable seismic source data.

8. A storage medium, characterized by The computer program stored in the storage medium can be executed by one or more processors and can be used to implement the seismic data matching method of any one of claims 1-6.

9. An electronic device, comprising: The computer program stored in the storage medium can be executed by one or more processors and can be used to implement the seismic data matching method of any one of claims 1-6.

Images