Reverse time focusing operator updating calculation method and device

Abstract

The invention provides a reverse time focusing operator updating calculation method and device, and relates to the technical field of geophysical exploration. The method includes: periodically performing calculation of a differential time shift panel; calculating a correlation matrix between adjacent channel signals in the differential time shifting panel and a signal amplitude statistical distribution matrix in the differential time shifting panel; and according to the differential time shifting panel, the correlation matrix and the signal amplitude statistical distribution matrix, generating visual content on an interactive interface, so that a user can continuously perform time shifting selection operation according to the visual content until the signal response in the differential time shifting panel is at zero time. And taking the inverse time focusing operator in the process of calculating the differential time shifting panel at the moment as the updated inverse time focusing operator. The device executes the method. According to the inverse time focusing operator updating calculation method and device provided by the embodiment of the invention, the inverse time focusing operator can be accurately and quickly obtained.

Description

Technical Field

[0001] The present invention relates to the technical field of geophysical exploration, and particularly to a method and device for updating and calculating an inverse-time focusing operator. Background Art

[0002] In the common focus point technology, calculating the inverse-time focusing operator and constructing the common focus point gather are the core parts. The common focus point gather is a new seismic data set obtained through partial migration operations. The diffracted energy in the gather has been migrated, and random noise has been suppressed. After obtaining the common focus point gather, the corresponding inverse-time focusing operator can be calculated based on the equal-time principle. The inverse-time focusing operator represents a set of one-way travel times from the common focus point to the ground, including information such as velocity inhomogeneity and anisotropy of the underground structure. Therefore, how to quickly and accurately obtain the inverse-time focusing operator is the most basic and core component in the common focus point technology.

[0003] With the deepening of exploration, the underground space structure has become increasingly complex. Some theoretical assumptions in traditional seismic data processing, such as small source-receiver offset and small dip angle assumptions, are no longer applicable. At the same time, with the widespread application of 3D seismic exploration, the amount of calculation and data has increased rapidly. The low calculation efficiency of traditional inversion imaging methods hinders their application in actual data. Moreover, because the common focus point technology is not mature enough in terms of theory and algorithm for obtaining the inverse-time focusing operator, the domestic research on the common focus point technology has gradually faded. Summary of the Invention

[0004] In view of the problems in the prior art, embodiments of the present invention provide a method and device for updating and calculating an inverse-time focusing operator, which can at least partially solve the problems existing in the prior art.

[0005] On the one hand, the present invention proposes a method for updating and calculating an inverse-time focusing operator, including:

[0006] Periodically performing the calculation of the differential time-shifted panel;

[0007] Calculating the correlation matrix between adjacent trace signals in the differential time-shifted panel, and the signal amplitude statistical distribution matrix in the differential time-shifted panel;

[0008] Generating visual content on an interactive interface according to the differential time-shifted panel, the correlation matrix, and the signal amplitude statistical distribution matrix, for the user to continuously perform time-shift selection operations according to the visual content until the signal response in the differential time-shifted panel is at the zero moment, and taking the inverse-time focusing operator in the process of calculating the differential time-shifted panel at this time as the updated inverse-time focusing operator.

[0009] Wherein, the periodically performing the calculation of the differential time-shifted panel includes:

[0010] Calculate the common focus point gather and the focusing operator;

[0011] Update and calculate the reverse time focusing operator according to the focusing operator;

[0012] Calculate the differential time shift panel according to the updated reverse time focusing operator and the common focus point gather.

[0013] Among them, the updating and calculating the reverse time focusing operator according to the focusing operator includes:

[0014] Take the focusing operator with the reflection depth on the surface as the initial focusing operator, and perform a transpose operation on the initial focusing operator;

[0015] Take the initial focusing operator after the transpose operation as the initial reverse time focusing operator for the first calculation, and replace the original reverse time focusing operator with the latest calculated reverse time focusing operator in each differential time shift panel calculation period;

[0016] Among them, the first original reverse time focusing operator is the initial reverse time focusing operator.

[0017] Among them, the calculating the differential time shift panel according to the updated reverse time focusing operator and the common focus point gather includes:

[0018] Perform a cross-correlation calculation on the updated reverse time focusing operator and the common focus point gather to obtain the differential time shift panel.

[0019] Among them, calculating the correlation matrix between adjacent channel signals in the differential time shift panel includes:

[0020] Perform a cross-correlation calculation between each pair of adjacent channel signals, and take the calculation results as the correlation matrix elements in the correlation matrix in turn.

[0021] Among them, calculating the signal amplitude statistical distribution matrix in the differential time shift panel includes:

[0022] Take the ratio of each signal amplitude to the number of signal channels in the differential time shift panel as the signal amplitude statistical distribution matrix elements in the signal amplitude statistical distribution matrix in turn.

[0023] Among them, the generating the visualization content on the interactive interface according to the differential time shift panel, the correlation matrix and the signal amplitude statistical distribution matrix includes:

[0024] Obtain the correlation matrix elements in the correlation matrix, and select the relatively large correlation matrix elements equal to the preset proportion as the target correlation matrix elements;

[0025] Obtain the signal amplitude statistical distribution matrix elements in the signal amplitude statistical distribution matrix, remove the signal amplitude statistical distribution matrix elements less than zero, calculate the average value of the remaining signal amplitude statistical distribution matrix elements, and use the remaining signal amplitude statistical distribution matrix elements greater than the average value as the target signal amplitude statistical distribution matrix elements;

[0026] Load the differential time shift panel, the target correlation matrix elements, and the target signal amplitude statistical distribution matrix elements in the generated initial interactive interface to obtain the visualization content on the interactive interface.

[0027] On the one hand, the present invention proposes an inverse time focusing operator update calculation device, including:

[0028] A first calculation unit for periodically performing the calculation of the differential time shift panel;

[0029] A second calculation unit for calculating the correlation matrix between adjacent channel signals in the differential time shift panel and the signal amplitude statistical distribution matrix in the differential time shift panel;

[0030] A third calculation unit for generating the visualization content on the interactive interface according to the differential time shift panel, the correlation matrix, and the signal amplitude statistical distribution matrix, so that the user can continuously perform time shift selection operations according to the visualization content until the signal response in the differential time shift panel is at the zero moment, and use the inverse time focusing operator in the process of calculating the differential time shift panel at this time as the updated inverse time focusing operator.

[0031] On the other hand, an embodiment of the present invention provides a computer device, including a memory, a processor, and a computer program stored on the memory and executable on the processor. When the processor executes the computer program, the following method is implemented:

[0032] Periodically perform the calculation of the differential time shift panel;

[0033] Calculate the correlation matrix between adjacent channel signals in the differential time shift panel and the signal amplitude statistical distribution matrix in the differential time shift panel;

[0034] Generate the visualization content on the interactive interface according to the differential time shift panel, the correlation matrix, and the signal amplitude statistical distribution matrix, so that the user can continuously perform time shift selection operations according to the visualization content until the signal response in the differential time shift panel is at the zero moment, and use the inverse time focusing operator in the process of calculating the differential time shift panel at this time as the updated inverse time focusing operator.

[0035] An embodiment of the present invention provides a computer-readable storage medium, including:

[0036] The computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, the following method is implemented:

[0037] Periodically perform the calculation of the differential time-shifted panel;

[0038] Calculate the correlation matrix between adjacent trace signals in the differential time-shifted panel, and the signal amplitude statistical distribution matrix in the differential time-shifted panel;

[0039] Generate visual content on an interactive interface according to the differential time-shifted panel, the correlation matrix, and the signal amplitude statistical distribution matrix, so that the user can continuously perform time-shift selection operations according to the visual content until the signal response in the differential time-shifted panel is at the zero moment. At this time, the inverse time focusing operator in the process of calculating the differential time-shifted panel is used as the updated inverse time focusing operator.

[0040] The inverse time focusing operator update calculation method and device provided by the embodiments of the present invention periodically perform the calculation of the differential time-shifted panel; calculate the correlation matrix between adjacent trace signals in the differential time-shifted panel, and the signal amplitude statistical distribution matrix in the differential time-shifted panel; generate visual content on an interactive interface according to the differential time-shifted panel, the correlation matrix, and the signal amplitude statistical distribution matrix, so that the user can continuously perform time-shift selection operations according to the visual content until the signal response in the differential time-shifted panel is at the zero moment. At this time, the inverse time focusing operator in the process of calculating the differential time-shifted panel is used as the updated inverse time focusing operator, and can accurately and quickly obtain the inverse time focusing operator. Description of the Drawings

[0041] In order to more clearly illustrate the technical solutions in the embodiments of the present invention or the prior art, the following will briefly introduce the drawings required for use in the description of the embodiments or the prior art. Obviously, the following drawings are only some embodiments of the present invention. For those of ordinary skill in the art, without creative efforts, other drawings can be obtained according to these drawings. In the drawings:

[0042] Figure 1 is a schematic flowchart of an inverse time focusing operator update calculation method provided by an embodiment of the present invention.

[0043] Figure 2 is a schematic diagram of an interactive interface description provided by an embodiment of the present invention.

[0044] Figure 3 is a schematic structural diagram of an inverse time focusing operator update calculation device provided by an embodiment of the present invention.

[0045] Figure 4Schematic diagram of the physical structure of the computer device provided by the embodiment of the present invention. Detailed implementation manners

[0046] To make the objectives, technical solutions and advantages of the embodiments of the present invention clearer and more understandable, the embodiments of the present invention will be further described in detail below with reference to the accompanying drawings. Herein, the illustrative embodiments of the present invention and their descriptions are used to explain the present invention, but not to limit the present invention. It should be noted that, without conflict, the embodiments in the present application and the features in the embodiments can be combined arbitrarily with each other.

[0047] Figure 1 is a schematic flowchart of the method for updating and calculating the inverse-time focusing operator provided by an embodiment of the present invention. As Figure 1 shown, the method for updating and calculating the inverse-time focusing operator provided by the embodiment of the present invention includes:

[0048] Step S1: Periodically execute the calculation of the differential time-shifted panel.

[0049] Step S2: Calculate the correlation matrix between adjacent trace signals in the differential time-shifted panel, and the signal amplitude statistical distribution matrix in the differential time-shifted panel.

[0050] Step S3: Generate the visualization content on the interactive interface according to the differential time-shifted panel, the correlation matrix and the signal amplitude statistical distribution matrix, so that the user can continuously perform the time-shift selection operation according to the visualization content until the signal response in the differential time-shifted panel is at the zero moment. At this time, the inverse-time focusing operator in the process of calculating the differential time-shifted panel is used as the updated inverse-time focusing operator.

[0051] In the above step S1, the device periodically executes the calculation of the differential time-shifted panel. The device can be a computer device that executes this method. In the technical solution of the present application, the acquisition, storage, use, processing, etc. of data all comply with relevant regulations. Whether to execute the method for updating and calculating the inverse-time focusing operator again can be triggered by whether the signal response in the differential time-shifted panel is at the zero moment, that is, if the signal response in the differential time-shifted panel is at the zero moment, the method for updating and calculating the inverse-time focusing operator is no longer executed.

[0052] If the signal response in the differential time-shifted panel is not at the zero moment, the method for updating and calculating the inverse-time focusing operator is executed again.

[0053] The differential time-shifted panel can be calculated by using the confocal point technology. The periodic execution of the calculation of the differential time-shifted panel includes:

[0054] Calculate the confocal point gather and the focusing operator;

[0055] Update and calculate the inverse-time focusing operator according to the focusing operator;

[0056] Calculate the differential time-lapse panel based on the updated reverse-time focusing operator and the common focus point gather.

[0057] The updating and calculating the reverse-time focusing operator according to the focusing operator includes:

[0058] Take the focusing operator with the reflection depth at the surface as the initial focusing operator, and perform a transpose operation on the initial focusing operator;

[0059] Take the initially transposed focusing operator as the initially calculated initial reverse-time focusing operator, and replace the original reverse-time focusing operator with the latest calculated reverse-time focusing operator in each differential time-lapse panel calculation period;

[0060] Wherein, the first original reverse-time focusing operator is the initial reverse-time focusing operator.

[0061] The calculating the differential time-lapse panel based on the updated reverse-time focusing operator and the common focus point gather includes:

[0062] Perform cross-correlation calculation on the updated reverse-time focusing operator and the common focus point gather to obtain the differential time-lapse panel.

[0063] The common focus point gather can be calculated according to the following formula:

[0064]

[0065] Wherein, the elements in the vector are the geophone positions at the same frequency, z m represents the reflection layer depth, z s represents the source depth; R + is the matrix representation of the reflection function; W + is the matrix representation of the propagation function, + represents downward propagation; D + is the downward wave field; S represents the source function. - represents upward propagation, is the unit row vector, and there is a 1 at the i-th position at the reflection layer depth z m (the other positions of this row vector are all zero), R + (z m ) is the matrix representation of the reflection function with the reflection layer depth as the independent variable, W + (z m , z s ) is the matrix representation of the propagation function with the reflection layer depth and the source depth as the independent variables, D + (z s ) is the downward wave field with the source depth as the independent variable, S(z s ) is the source function with the source depth as the independent variable.

[0066] The focusing operator can be calculated according to the following formula:

[0067]

[0068] Wherein, reference can be made to the above description; z r represents the reflection depth; is the focusing operator acting on the geophone position; [D - (z r )] -1 is the inverse matrix form of the up-going wavefield function from the reflection depth to the receiving depth. W + (z m ,z r ) is the matrix expression form of the propagation function with the reflection layer depth and the reflection depth as independent variables.

[0069] The initial reverse-time focusing operator can be calculated according to the following formula

[0070]

[0071] Wherein, z0 represents the reflection depth at the surface of the earth, and correspondingly represents the initial focusing operator, represents the transpose operation on the initial focusing operator.

[0072] The calculation period of each differential time shift panel can be equivalent to the execution period of the reverse-time focusing operator update calculation method. That is, the reverse-time focusing operator calculated in the first differential time shift panel calculation period is The reverse-time focusing operator calculated in the second differential time shift panel calculation period is Use to replace to complete the subsequent steps, and so on.

[0073] The differential time shift panel can be calculated according to the following formula:

[0074]

[0075] DTS is the differential time shift panel, is the cross-correlation calculation. The calculation period of each differential time shift panel can use the corresponding to replace j is the number of differential time shift panel calculation periods. Therefore, the differential time shift panel can be dynamically updated according to the reverse-time focusing operator.

[0076] In the above step S2, the device calculates the correlation matrix between adjacent channel signals in the differential time shift panel and the signal amplitude statistical distribution matrix in the differential time shift panel. Calculating the correlation matrix between adjacent channel signals in the differential time shift panel includes:

[0077] Perform cross-correlation calculation between each pair of adjacent channel signals, and take the calculation results as the correlation matrix elements in the correlation matrix in sequence. The correlation matrix can be calculated according to the following formula:

[0078]

[0079] where R is the correlation matrix, dts i is a signal in the differential time shift panel, dts i+1 is an adjacent signal in the differential time shift panel, and n is the number of signal channels in the differential time shift panel.

[0080] Calculating the signal amplitude statistical distribution matrix in the differential time shift panel includes:

[0081] Take the ratio of each signal amplitude to the number of signal channels in the differential time shift panel as the signal amplitude statistical distribution matrix elements in the signal amplitude statistical distribution matrix in sequence. The signal amplitude statistical distribution matrix can be calculated according to the following formula:

[0082]

[0083] where A is the signal amplitude statistical distribution matrix, a i is the i-th signal amplitude, and n is the number of signal channels in the differential time shift panel.

[0084] In the above step S3, the device generates visual content on the interactive interface according to the differential time shift panel, the correlation matrix, and the signal amplitude statistical distribution matrix, so that the user can continuously perform time shift selection operations according to the visual content until the signal response in the differential time shift panel is at the zero moment. At this time, the inverse time focusing operator in the process of calculating the differential time shift panel is used as the updated inverse time focusing operator. The visual content on the interactive interface, as Figure 2 shown, Figure 2 where the horizontal axis in represents the number of signal channels in the differential time shift panel, and the vertical axis represents the time difference, specifically the difference from the inverse time focusing operator. Figure 2 The two rectangular boxes (the shape is not specifically limited) in can respectively represent two constraint conditions for the time shift selection operation, where one constraint condition corresponds to the correlation matrix, and the other constraint condition corresponds to the signal amplitude statistical distribution matrix.

[0085] The above two constraint conditions can assist the user in continuously performing time-shift selection operations based on the visualization content. Each time-shift selection operation can make the signal response continuously approach the point at 0 on the vertical axis as Figure 2 shown in the vertical axis.

[0086] The rectangular box is the area with strong correlation in the correlation matrix, and the area where the signal amplitude statistical distribution matrix is dense respectively. The curve is the predicted time-shift curve. With the assistance of the correlation matrix, the signal amplitude statistical distribution matrix, and the predicted time-shift curve, the user picks up the time shift. The interactive interface receives the user's picking action to complete the picking of the time shift and completes an update of the reverse-time focusing operator.

[0087] Generating the visualization content on the interactive interface according to the differential time-shift panel, the correlation matrix, and the signal amplitude statistical distribution matrix includes:

[0088] Obtain the correlation matrix elements in the correlation matrix, and select the correlation matrix elements with larger values equal to the preset proportion as the target correlation matrix elements; the preset proportion can be set independently according to the actual situation, and can be 40%. The correlation matrix elements can be sorted in descending order of values, and the correlation matrix elements equal to 40% of the total number of correlation matrix elements can be selected from the front to the back as the target correlation matrix elements.

[0089] Obtain the signal amplitude statistical distribution matrix elements in the signal amplitude statistical distribution matrix, remove the signal amplitude statistical distribution matrix elements less than zero, calculate the average value of the remaining signal amplitude statistical distribution matrix elements, and use the remaining signal amplitude statistical distribution matrix elements greater than the average value as the target signal amplitude statistical distribution matrix elements; the remaining signal amplitude statistical distribution matrix elements retained are all signal amplitude statistical distribution matrix elements greater than zero, and then the above signal amplitude statistical distribution matrix elements greater than zero are further screened by the average value to obtain the target signal amplitude statistical distribution matrix elements.

[0090] Load the differential time-shift panel, the target correlation matrix elements, and the target signal amplitude statistical distribution matrix elements into the generated initial interactive interface to obtain the visualization content on the interactive interface.

[0091] The beneficial effect of the reverse-time focusing operator update calculation method provided by the embodiment of the present invention is that, compared with the prior art, this method uses the interactive interface, as well as the correlation matrix and the amplitude statistical distribution matrix to constrain the picking range, so as to predict the time-shift curve, thereby assisting the user in picking up the time shift. Since this method comprehensively and systematically considers the distribution characteristics of the differential time-shift panel signal and the correlation characteristics between signals, therefore, a more accurate preliminary analysis of the differential time-shift panel can be performed, providing a reference for the user and solving the problem of difficult time-shift picking in the calculation of the reverse-time focusing operator.

[0092] By calculating the inverse time focusing operator, further work such as confocal point tomography velocity modeling can be carried out.

[0093] The inverse time focusing operator update calculation method provided by the embodiments of the present invention periodically performs calculating a differential time shift panel; calculating a correlation matrix between adjacent trace signals in the differential time shift panel, and a signal amplitude statistical distribution matrix in the differential time shift panel; generating visual content on an interactive interface according to the differential time shift panel, the correlation matrix, and the signal amplitude statistical distribution matrix, so that a user can continuously perform time shift selection operations according to the visual content until the signal response in the differential time shift panel is at the zero moment, and taking the inverse time focusing operator in the process of calculating the differential time shift panel at this time as the updated inverse time focusing operator, which can accurately and quickly obtain the inverse time focusing operator.

[0094] Further, the periodically performing calculating the differential time shift panel includes:

[0095] Calculating a confocal point gather and a focusing operator; which can be described with reference to the above embodiments and will not be elaborated here.

[0096] Updating and calculating the inverse time focusing operator according to the focusing operator; which can be described with reference to the above embodiments and will not be elaborated here.

[0097] Calculating the differential time shift panel according to the updated inverse time focusing operator and the confocal point gather; which can be described with reference to the above embodiments and will not be elaborated here.

[0098] Further, the updating and calculating the inverse time focusing operator according to the focusing operator includes:

[0099] Taking the focusing operator with the reflection depth at the surface as the initial focusing operator, and performing a transpose operation on the initial focusing operator; which can be described with reference to the above embodiments and will not be elaborated here.

[0100] Taking the initial focusing operator after the transpose operation as the initial inverse time focusing operator for the first calculation, and replacing the original inverse time focusing operator with the latest calculated inverse time focusing operator in each differential time shift panel calculation period; which can be described with reference to the above embodiments and will not be elaborated here.

[0101] Wherein, the first original inverse time focusing operator is the initial inverse time focusing operator; which can be described with reference to the above embodiments and will not be elaborated here.

[0102] Further, the calculating the differential time shift panel according to the updated inverse time focusing operator and the confocal point gather includes:

[0103] Perform cross - correlation calculation on the updated reverse - time focusing operator and the confocal point gather to obtain the differential time - shift panel. It can be described with reference to the above - mentioned embodiments and will not be elaborated here.

[0104] Further, calculating the correlation matrix between adjacent - trace signals in the differential time - shift panel includes:

[0105] Perform cross - correlation calculation between each pair of adjacent - trace signals, and use the calculation results as the correlation matrix elements in the correlation matrix in sequence. It can be described with reference to the above - mentioned embodiments and will not be elaborated here.

[0106] Further, calculating the signal amplitude statistical distribution matrix in the differential time - shift panel includes:

[0107] Take the ratio of each signal amplitude to the number of signal traces in the differential time - shift panel as the signal amplitude statistical distribution matrix elements in the signal amplitude statistical distribution matrix in sequence. It can be described with reference to the above - mentioned embodiments and will not be elaborated here.

[0108] Further, generating the visualization content on the interactive interface according to the differential time - shift panel, the correlation matrix, and the signal amplitude statistical distribution matrix includes:

[0109] Obtain the correlation matrix elements in the correlation matrix, and select the numerically larger correlation matrix elements equal to the preset proportion as the target correlation matrix elements; it can be described with reference to the above - mentioned embodiments and will not be elaborated here.

[0110] Obtain the signal amplitude statistical distribution matrix elements in the signal amplitude statistical distribution matrix, remove the signal amplitude statistical distribution matrix elements less than zero, calculate the average value of the remaining signal amplitude statistical distribution matrix elements, and take the remaining signal amplitude statistical distribution matrix elements greater than the average value as the target signal amplitude statistical distribution matrix elements; it can be described with reference to the above - mentioned embodiments and will not be elaborated here.

[0111] Load the differential time - shift panel, the target correlation matrix elements, and the target signal amplitude statistical distribution matrix elements into the generated initial interactive interface to obtain the visualization content on the interactive interface. It can be described with reference to the above - mentioned embodiments and will not be elaborated here.

[0112] Figure 3 It is a schematic structural diagram of the reverse - time focusing operator update calculation device provided by an embodiment of the present invention. As Figure 3 shown, the reverse - time focusing operator update calculation device provided by the embodiment of the present invention includes a first calculation unit 301, a second calculation unit 302, and a third calculation unit 303, where:

[0113] The first calculation unit 301 is used to periodically perform the calculation of the differential time-shifted panel; the second calculation unit 302 is used to calculate the correlation matrix between adjacent channel signals in the differential time-shifted panel and the signal amplitude statistical distribution matrix in the differential time-shifted panel; the third calculation unit 303 is used to generate visual content on the interactive interface according to the differential time-shifted panel, the correlation matrix and the signal amplitude statistical distribution matrix, so that the user can continuously perform time-shift selection operations according to the visual content until the signal response in the differential time-shifted panel is at the zero moment, and the inverse time focusing operator in the process of calculating the differential time-shifted panel at this time is used as the updated inverse time focusing operator.

[0114] Specifically, the first calculation unit 301 in the device is used to periodically perform the calculation of the differential time-shifted panel; the second calculation unit 302 is used to calculate the correlation matrix between adjacent channel signals in the differential time-shifted panel and the signal amplitude statistical distribution matrix in the differential time-shifted panel; the third calculation unit 303 is used to generate visual content on the interactive interface according to the differential time-shifted panel, the correlation matrix and the signal amplitude statistical distribution matrix, so that the user can continuously perform time-shift selection operations according to the visual content until the signal response in the differential time-shifted panel is at the zero moment, and the inverse time focusing operator in the process of calculating the differential time-shifted panel at this time is used as the updated inverse time focusing operator.

[0115] The inverse time focusing operator update calculation device provided by the embodiment of the present invention periodically performs the calculation of the differential time-shifted panel; calculates the correlation matrix between adjacent channel signals in the differential time-shifted panel and the signal amplitude statistical distribution matrix in the differential time-shifted panel; generates visual content on the interactive interface according to the differential time-shifted panel, the correlation matrix and the signal amplitude statistical distribution matrix, so that the user can continuously perform time-shift selection operations according to the visual content until the signal response in the differential time-shifted panel is at the zero moment, and the inverse time focusing operator in the process of calculating the differential time-shifted panel at this time is used as the updated inverse time focusing operator, and can accurately and quickly obtain the inverse time focusing operator.

[0116] Figure 4 It is a schematic diagram of the entity structure of the computer device provided by the embodiment of the present invention, as Figure 4 shown, the computer device includes: a memory 401, a processor 402, and a computer program stored on the memory 401 and executable on the processor 402. When the processor 402 executes the computer program, the following method is implemented:

[0117] Periodically perform the calculation of the differential time-shifted panel;

[0118] Calculate the correlation matrix between adjacent trace signals in the differential time - shift panel, and the signal amplitude statistical distribution matrix in the differential time - shift panel;

[0119] Generate visualization content on an interactive interface based on the differential time - shift panel, the correlation matrix, and the signal amplitude statistical distribution matrix, for the user to continuously perform time - shift selection operations according to the visualization content until the signal response in the differential time - shift panel is at the zero moment. Take the inverse - time focusing operator in the process of calculating the differential time - shift panel at this time as the updated inverse - time focusing operator.

[0120] This embodiment discloses a computer program product. The computer program product includes a computer program, and when the computer program is executed by a processor, the following method is implemented:

[0121] Periodically execute the calculation of the differential time - shift panel;

[0122] Calculate the correlation matrix between adjacent trace signals in the differential time - shift panel, and the signal amplitude statistical distribution matrix in the differential time - shift panel;

[0123] Generate visualization content on an interactive interface based on the differential time - shift panel, the correlation matrix, and the signal amplitude statistical distribution matrix, for the user to continuously perform time - shift selection operations according to the visualization content until the signal response in the differential time - shift panel is at the zero moment. Take the inverse - time focusing operator in the process of calculating the differential time - shift panel at this time as the updated inverse - time focusing operator.

[0124] This embodiment provides a computer - readable storage medium. The computer - readable storage medium stores a computer program, and when the computer program is executed by a processor, the following method is implemented:

[0125] Periodically execute the calculation of the differential time - shift panel;

[0126] Calculate the correlation matrix between adjacent trace signals in the differential time - shift panel, and the signal amplitude statistical distribution matrix in the differential time - shift panel;

[0127] Generate visualization content on an interactive interface based on the differential time - shift panel, the correlation matrix, and the signal amplitude statistical distribution matrix, for the user to continuously perform time - shift selection operations according to the visualization content until the signal response in the differential time - shift panel is at the zero moment. Take the inverse - time focusing operator in the process of calculating the differential time - shift panel at this time as the updated inverse - time focusing operator.

[0128] Those skilled in the art should understand that the embodiments of the present invention can be provided as a method, a system, or a computer program product. Therefore, the present invention can take the form of a complete hardware embodiment, a complete software embodiment, or an embodiment combining software and hardware aspects. Moreover, the present invention can take the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) that contain computer-usable program code.

[0129] The present invention is described with reference to the flowcharts and / or block diagrams of methods, apparatuses (systems), and computer program products according to embodiments of the present invention. It should be understood that each flow and / or block in the flowchart and / or block diagram can be implemented by computer program instructions, and the combination of flows and / or blocks in the flowchart and / or block diagram can also be implemented by computer program instructions. These computer program instructions can be provided to the processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing devices to generate a machine, such that the instructions executed by the processor of the computer or other programmable data processing devices generate means for implementing the functions specified in Figure 1 one flow or multiple flows and / or blocks Figure 1 one block or multiple blocks.

[0130] These computer program instructions can also be stored in a computer-readable memory that can direct a computer or other programmable data processing device to work in a specific manner, such that the instructions stored in the computer-readable memory generate a manufactured article including instruction means that implement the functions specified in Figure 1 one flow or multiple flows and / or blocks Figure 1 one block or multiple blocks.

[0131] These computer program instructions can also be loaded onto a computer or other programmable data processing device, such that a series of operation steps are executed on the computer or other programmable device to generate a computer-implemented process, so that the instructions executed on the computer or other programmable device provide steps for implementing the functions specified in Figure 1 one flow or multiple flows and / or blocks Figure 1 one block or multiple blocks.

[0132] In the description of this specification, the description with reference to terms such as "one embodiment", "a specific embodiment", "some embodiments", "for example", "example", "specific example", or "some examples", etc. means that the specific features, structures, materials, or characteristics described in connection with the embodiment or example are included in at least one embodiment or example of the present invention. In this specification, the schematic expressions of the above terms do not necessarily refer to the same embodiment or example. Moreover, the specific features, structures, materials, or characteristics described can be combined in a suitable manner in any one or more embodiments or examples.

[0133] The specific embodiments described above further elaborate on the purpose, technical solutions, and beneficial effects of the present invention. It should be understood that the above are only specific embodiments of the present invention and are not used to limit the protection scope of the present invention. Any modifications, equivalent replacements, improvements, etc. made within the spirit and principles of the present invention shall be included in the protection scope of the present invention.

Claims

1. A method for updating and calculating an inverse-time focusing operator, characterized in that, Including: Periodically perform the calculation of the differential time-shifted panel; Calculate the correlation matrix between adjacent trace signals in the differential time-shifted panel, and the signal amplitude statistical distribution matrix in the differential time-shifted panel; Generate visual content on the interactive interface according to the differential time-shifted panel, the correlation matrix, and the signal amplitude statistical distribution matrix, so that the user can continuously perform time-shift selection operations according to the visual content until the signal response in the differential time-shifted panel is at the zero moment, and use the reverse-time focusing operator in the process of calculating the differential time-shifted panel at this time as the updated reverse-time focusing operator.

2. The method for updating the calculation of the inverse time focusing operator according to claim 1, wherein The periodically performing the calculation of the differential time-shifted panel includes: Calculate the common focus point gather and the focusing operator; Update and calculate the reverse-time focusing operator according to the focusing operator; Calculate the differential time-shifted panel according to the updated reverse-time focusing operator and the common focus point gather.

3. The calculation method for updating the inverse time focusing operator according to claim 2, characterized in that The updating and calculating the reverse-time focusing operator according to the focusing operator includes: Use the focusing operator with the reflection depth at the surface as the initial focusing operator, and perform a transpose operation on the initial focusing operator; Use the transposed initial focusing operator as the initial reverse-time focusing operator for the first calculation, and replace the original reverse-time focusing operator with the latest calculated reverse-time focusing operator in each differential time-shifted panel calculation period; Wherein, the first original reverse-time focusing operator is the initial reverse-time focusing operator.

4. The calculation method for updating the inverse time focusing operator according to claim 2, characterized in that The calculating the differential time-shifted panel according to the updated reverse-time focusing operator and the common focus point gather includes: Perform cross-correlation calculation on the updated reverse-time focusing operator and the common focus point gather to obtain the differential time-shifted panel.

5. The method for updating and calculating the inverse time focusing operator according to claim 1, wherein Calculating the correlation matrix between adjacent trace signals in the differential time-shifted panel includes: Perform cross-correlation calculation on each pair of adjacent trace signals, and use the calculation results as the correlation matrix elements in the correlation matrix in sequence.

6. The calculation method for updating the inverse-time focusing operator according to claim 1, wherein Calculating the signal amplitude statistical distribution matrix in the differential time-shifted panel includes: Use the ratio of each signal amplitude to the number of signal traces in the differential time-shifted panel as the signal amplitude statistical distribution matrix elements in the signal amplitude statistical distribution matrix in sequence.

7. The method for updating and calculating the inverse time focusing operator according to any one of claims 1 to 6, characterized in that The generating the visual content on the interactive interface according to the differential time-shifted panel, the correlation matrix, and the signal amplitude statistical distribution matrix includes: Obtain the correlation matrix elements in the correlation matrix, and select the relatively large correlation matrix elements equal to the preset proportion as the target correlation matrix elements; Obtain the signal amplitude statistical distribution matrix elements in the signal amplitude statistical distribution matrix, remove the signal amplitude statistical distribution matrix elements less than zero, calculate the average value of the remaining signal amplitude statistical distribution matrix elements, and use the remaining signal amplitude statistical distribution matrix elements greater than the average value as the target signal amplitude statistical distribution matrix elements; Load the differential time-shifted panel, the target correlation matrix elements, and the target signal amplitude statistical distribution matrix elements into the generated initial interactive interface to obtain the visual content on the interactive interface.

8. An inverse time focusing operator update calculation device, characterized in that, Including: A first calculation unit for periodically performing the calculation of the differential time-shifted panel; A second calculation unit, configured to calculate a correlation matrix between adjacent channel signals in the differential time-shifted panel, and a signal amplitude statistical distribution matrix in the differential time-shifted panel; A third calculation unit, configured to generate visualization content on an interactive interface according to the differential time-shifted panel, the correlation matrix, and the signal amplitude statistical distribution matrix, so that a user can continuously perform time-shift selection operations according to the visualization content until the signal response in the differential time-shifted panel is at the zero moment, and use the reverse-time focusing operator in the process of calculating the differential time-shifted panel at this time as the updated reverse-time focusing operator.

9. A computer device, comprising a memory, a processor, and a computer program stored on the memory and executable on the processor, characterized in that, When the processor executes the computer program, the method according to any one of claims 1 to 7 is implemented.

10. A computer-readable storage medium, characterized in that, The computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, the method according to any one of claims 1 to 7 is implemented.

Images