Multi-static correction joint optimization correction determination method, device, equipment and medium

Abstract

This invention relates to the field of seismic exploration or surveying technology, and discloses a method, apparatus, equipment, and medium for determining optimized calibration quantities through multi-static calibration combined with other methods. The method includes: determining a first static calibration quantity to fill to the reference surface based on a pre-established static calibration equivalent model using a first static calibration method; determining a second static calibration quantity to fill to the surface elevation based on the first static calibration quantity and the static calibration equivalent model; determining shot-receiver point information and first arrival based on the second static calibration quantity; determining a third static calibration quantity from the surface elevation to the reference surface based on the shot-receiver point information, first arrival, and the static calibration equivalent model using the second static calibration method; and determining optimized calibration quantities based on the second and third static calibration quantities. The solution of this invention adopts a step-by-step filling method, decomposing the filling to the surface at the shot-receiver point into filling to the reference surface, and can employ different static calibration methods to achieve multi-static calibration combined with other methods for comprehensive analysis.

Description

Technical Field

[0001] This invention relates to the field of seismic exploration or surveying technology, specifically to a method, apparatus, equipment, and medium for determining optimized calibration quantities in multi-static calibration combined with other methods. Background Technology

[0002] The existing definition of static correction is to eliminate the influence of factors such as changes in surface elevation, variations in velocity and thickness of low-velocity layers, and differences in shot point inlet depth. Static correction is a qualitative definition, and there is currently no unified calculation formula for solving the static correction problem. Various methods have emerged to address this issue. Currently, the application of each static correction method must first meet the pre-set assumptions of the method. In actual production, under complex near-surface conditions, the static correction values ​​calculated by different static correction methods, when applied to seismic data overlay, result in varying imaging accuracy in the overlay profile.

[0003] In related technologies, high- and low-frequency separation is typically employed, with the low-frequency component of a certain static correction method selected to control the orientation. Imaging accuracy is achieved by fusing the high-frequency components of different static correction methods in segments. However, this approach still suffers from the limitation of using a single static correction method. Due to the incompatibility of the software interfaces of various static correction methods, it is impossible to perform comprehensive analysis using different static correction software. Summary of the Invention

[0004] In view of this, the present invention provides a method, apparatus, equipment, and medium for determining optimized calibration quantities for multiple static calibrations, in order to solve the technical problem that multiple static calibration methods cannot be combined for comprehensive analysis.

[0005] In a first aspect, the present invention provides a method for determining optimized calibration quantities for multiple static calibrations, the method comprising: determining a first static calibration quantity for filling to a reference surface based on a pre-established static calibration equivalent model and using a first static calibration method; determining a second static calibration quantity for filling to the ground surface elevation based on the first static calibration quantity and the static calibration equivalent model; determining shot-receiver point information and initial elevation based on the second static calibration quantity; determining a third static calibration quantity for filling from the ground surface elevation to the reference surface based on the shot-receiver point information, initial elevation, and the static calibration equivalent model and using a second static calibration method; and determining optimized calibration quantities based on the second and third static calibration quantities.

[0006] In conjunction with the first aspect, in one possible implementation of the first aspect, a second static correction amount is determined based on the first static correction amount and the static correction equivalent model, including: determining the filling thickness between the reference surface and the ground surface elevation based on the model data of the static correction equivalent model; determining the filling amount corresponding to the filling thickness based on the model data of the static correction equivalent model; and determining the second static correction amount to the ground surface elevation based on the difference between the first static correction amount and the filling amount.

[0007] In conjunction with the first aspect, in one possible implementation of the first aspect, the shot-receiver information and first arrival are determined based on the second static correction, including: adding the second static correction to the original single-shot first arrival to form a new first arrival; and, based on the second static correction, determining the corresponding shot coordinates, receiver coordinates, shot elevation, receiver elevation, and well depth; setting the well depth to preset data, and using the set well depth, shot coordinates, receiver coordinates, shot elevation, and receiver elevation as the shot-receiver information.

[0008] In conjunction with the first aspect, in one possible implementation of the first aspect, setting the well depth to preset data includes: setting the well depth to zero.

[0009] In conjunction with the first aspect, in one possible implementation of the first aspect, based on shot-receiver information, first arrival, and static correction equivalent model, a third static correction amount is determined by using a second static correction method to fill the surface elevation to the reference surface. This includes: determining the corresponding replacement velocity between the surface elevation and the reference surface based on model data from the static correction equivalent model; and performing inversion calculations using the second static correction method based on shot-receiver information, first arrival, and replacement velocity to determine the third static correction amount filling the surface elevation to the reference surface.

[0010] In conjunction with the first aspect, in one possible implementation of the first aspect, after determining the third static correction amount from the surface elevation to the reference surface using the second static correction method based on the shot-receiver information, the initial arrival and the static correction equivalent model, the method further includes: using the sum of the second static correction and the third static correction as the new first static correction amount; determining the new second static correction amount from the surface elevation to the reference surface based on the new first static correction amount and the static correction equivalent model; determining new shot-receiver information and a new initial arrival based on the new second static correction amount; determining the new third static correction amount from the surface elevation to the reference surface using the third static correction method based on the new shot-receiver information, the new initial arrival and the static correction equivalent model, and using the new third static correction amount as the third static correction amount.

[0011] In conjunction with the first aspect, in one possible implementation of the first aspect, the method further includes: using multiple different first static correction methods and second static correction methods to iteratively determine multiple optimized correction quantities; superimposing the multiple optimized correction quantities, comparing the profile imaging accuracy with the attitude, and determining the optimal optimized correction quantity.

[0012] Secondly, the present invention provides an optimized calibration quantity determination device for multiple static calibrations, the device comprising: a first static calibration quantity determination module, used to determine a first static calibration quantity to fill to a reference surface based on a pre-established static calibration equivalent model and using a first static calibration method; a second static calibration quantity determination module, used to determine a second static calibration quantity to fill to the ground surface elevation based on the first static calibration quantity and the static calibration equivalent model; an information determination module, used to determine shot-receiver point information and initial arrival based on the second static calibration quantity; a third static calibration quantity determination module, used to determine a third static calibration quantity to fill from the ground surface elevation to the reference surface based on the shot-receiver point information, initial arrival, and the static calibration equivalent model and using a second static calibration method; and an optimized calibration quantity determination module, used to determine an optimized calibration quantity based on the second and third static calibration quantities.

[0013] Thirdly, the present invention provides a computer device, comprising: a memory and a processor, wherein the memory and the processor are communicatively connected to each other, the memory stores computer instructions, and the processor executes the computer instructions to perform the optimized calibration determination method for multi-static calibration linkage described in the first aspect or any corresponding embodiment.

[0014] Fourthly, the present invention provides a computer-readable storage medium storing computer instructions for causing a computer to perform the optimized calibration determination method for multi-static calibration linkage described in the first aspect or any corresponding embodiment thereof.

[0015] The technical solution of this invention has the following advantages: This invention provides a method, apparatus, equipment, and medium for determining optimized calibration quantities through multiple static corrections. The method uses a pre-established static correction equivalent model to determine a first static correction quantity filled to a reference surface using a first static correction method. Based on this, a second static correction quantity filled to the ground surface elevation is determined. Furthermore, using the second static correction method, a third static correction quantity filled from the ground surface elevation to the reference surface is determined. The optimized calibration quantities are then determined using the second and third static correction quantities. This process employs a step-by-step filling method, decomposing the filling calculation calibration process into two inversion processes: determining the second and third static correction quantities. The determination of the second static correction quantity corresponds to filling to the ground surface at the shot receiver point, while the determination of the third static correction quantity corresponds to filling to a fixed horizontal reference surface or a floating reference surface. Because the filling calculation calibration process is decomposed, the two inversion processes can use the same or different first and second static correction methods, thus achieving the combined application of multiple static correction methods. This allows for comprehensive analysis through multiple combined static correction methods to determine the optimized calibration quantities. Attached Figure Description

[0016] To more clearly illustrate the specific embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the specific embodiments or the prior art will be briefly introduced below. Obviously, the drawings described below are some embodiments of the present invention. For those skilled in the art, other drawings can be obtained from these drawings without creative effort.

[0017] Figure 1 This is a flowchart illustrating an optimized calibration determination method for multiple static calibrations provided in an embodiment of the present invention. Figure 2 This is a schematic diagram of the calibration calculation of a static correction equivalent model provided in an embodiment of the present invention; Figure 3 This is a schematic diagram of the residuals of the first arrival wave participating in the inversion model provided by an embodiment of the present invention; Figure 4 This is a schematic diagram comparing the accuracy of CMP overlay profile imaging according to an embodiment of the present invention; Figure 5 This is another CMP overlay profile imaging accuracy comparison diagram provided by an embodiment of the present invention; Figure 6 This is a structural block diagram of an optimized calibration determination device for multi-static calibration combined according to an embodiment of the present invention; Figure 7 This is a schematic diagram of the hardware structure of a computer device according to an embodiment of the present invention. Detailed Implementation

[0018] To make the objectives, technical solutions, and advantages of the embodiments of the present invention clearer, the technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of the present invention, not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of the present invention.

[0019] According to an embodiment of the present invention, an embodiment of an optimized calibration determination method for multiple static calibrations is provided. It should be noted that the steps shown in the flowchart in the accompanying drawings can be executed in a computer system such as a set of computer-executable instructions. Furthermore, although a logical order is shown in the flowchart, in some cases, the steps shown or described may be executed in a different order than that shown here.

[0020] This embodiment provides an optimized calibration determination method for multiple static calibrations, such as... Figure 1 As shown, the method includes the following steps: S101. Based on the pre-established static correction equivalent model, the first static correction amount to be filled to the reference surface is determined using the first static correction method.

[0021] Specifically, the establishment of a static correction equivalent model requires the following elements: near-surface shot point and receiver elevation, elevation of the high-velocity top interface and average velocity from the receiver to the high-velocity top interface, shot point firing depth, and average velocity from the actual shot point firing position to the high-velocity top interface. The current model obtained from tomographic static correction inversion is the average velocity within each small grid cell downwards from the near-surface plumb line. It can also be converted into a near-surface equivalent model through simple calculations. The static correction equivalent model allows for rapid analysis of the static correction calculation process. For example... Figure 2 As shown, where, T S T represents the shot point correction amount filled to the datum elevation. S填 This indicates the amount of fill from the shot point at the top elevation of the high-speed highway to the horizontal reference surface. H S地表高程 This indicates the elevation of the ground surface at the point where the plumb bob points upwards. H S高速顶 Indicates the elevation of the high-speed peak of the artillery position. H S井深 Indicates the depth of the firing point. V S平均 This represents the average velocity of the shot point from the firing point to the high-speed peak elevation of the shot point. V 替换速度 This indicates the rate of constant substitution in the calculation of the calibration. T R This represents the correction amount for the receiver point filled to the elevation of the reference plane. T R剥 T represents the amount of stripping from the receiver at the surface to the corresponding elevation of the high-speed peak at the receiver. R填 This indicates the amount of fill from the detector point at the top elevation of the high-speed highway to the horizontal reference surface. H R地表高程 This indicates the surface elevation at the receiver point. H R高速顶 This indicates the elevation of the high-speed peak corresponding to the detector point. V R平均 This represents the average velocity of the receiver point from the ground surface to the corresponding high-speed peak elevation.

[0022] Specifically, the first static correction method includes: refraction wave method, tomography method, generalized linear method, Cronowski inversion method, etc. Typically, multiple static correction methods and software are selected to perform static correction inversion calculations on the two-dimensional or three-dimensional survey lines of the work area, calculating the model quantities and first-arrival residual quantities. The obtained static correction quantities are then superimposed on the CMP (Common Middle Point) gather to perform comparative analysis of the accuracy and attitude of the superimposed profile imaging. Several static corrections suitable for this work area are selected for combined use of multiple static corrections, i.e., as the first, second, and third static correction methods.

[0023] Specifically, static corrections include shot point corrections and receiver point corrections. In the inverted near-surface model, these can be approximated as the geometric process of the self-excited and self-received vertical propagation path from the excitation point to the receiver point in the corresponding geodetic coordinates. This process is completed through a stripping-filling calculation. Specifically: In the static correction equivalent model, the receiver point correction calculation first involves stripping from the receiver point surface to the high-velocity top interface from the receiver point's vertical direction (negative value). Then, it involves filling from the high-velocity top interface from the receiver point's vertical direction to the fixed or floating reference surface using a high-velocity replacement rate (positive value). The shot point correction calculation first involves stripping from the shot point's well depth to the high-velocity top interface from the shot point's vertical direction (negative value). Then, it involves filling from the high-velocity top interface from the shot point's vertical direction to the fixed or floating reference surface using a high-velocity replacement rate (positive value). Therefore, the first static correction amount filling to the reference surface refers to the static correction amount that is stripped and then filled directly to the reference surface.

[0024] S102. Based on the first static correction amount and the static correction equivalent model, determine the second static correction amount to fill to the ground surface elevation.

[0025] Specifically, the second static correction amount for filling to the ground surface elevation refers to the static correction amount for filling to the ground surface elevation by decomposing the filling process in the relevant technology.

[0026] Specifically, determining the second static correction amount based on the first static correction amount and the static correction equivalent model involves using the model data from the static correction equivalent model to determine the filling amount between the reference surface and the ground surface elevation, and then determining the second static correction amount by the difference between the first static correction amount and the filling amount.

[0027] S103. Based on the second static correction, determine the shot receiver information and initial arrival.

[0028] Specifically, determining the shot-receiver information and initial arrival based on the second static correction means adding the shot point correction and receiver correction corresponding to the second static correction to the original single shot initial arrival to form a new initial arrival. The shot-receiver information, including shot point coordinates, receiver coordinates, shot point elevation, receiver elevation, and well depth, is determined through the second static correction, where the well depth is set to a preset value.

[0029] S104. Based on shot receiver information and the equivalent model of initial arrival and static correction, the third static correction amount from the ground elevation to the reference surface is determined using the second static correction method.

[0030] Specifically, the third static correction amount, which is the amount of static correction from the ground surface elevation to the reference surface, refers to the amount of static correction obtained by decomposing the filling process in the relevant technology and filling it from the ground surface elevation to the reference surface. The second static correction method can be the same as or different from the correction method used in the first static correction method, and the specific selection depends on the actual working conditions. This embodiment does not impose any specific limitations on this.

[0031] Specifically, based on shot-receiver point information, first arrival and static correction equivalent model, the third static correction amount from the surface elevation to the reference surface is determined using the second static correction method. This refers to determining the corresponding replacement rate between the surface elevation and the reference surface based on the model data of the static correction equivalent model, and combining the determined shot-receiver point information and first arrival, the third static correction amount from the surface elevation to the reference surface is determined by inversion using the second static correction method.

[0032] S105. Based on the second static correction value and the third static correction value, determine the optimized correction value.

[0033] Specifically, determining the optimal calibration value based on the second static calibration value and the third static calibration value means adding the second static calibration value and the third static calibration value as the optimal calibration value.

[0034] This invention provides a method, apparatus, equipment, and medium for determining optimized calibration quantities through multiple static corrections. The method uses a pre-established static correction equivalent model to determine a first static correction quantity filled to a reference surface using a first static correction method. Based on this, a second static correction quantity filled to the ground surface elevation is determined. Then, using the second static correction method, a third static correction quantity filled from the ground surface elevation to the reference surface is determined. The optimized calibration quantities are then determined using the second and third static correction quantities. This process employs a step-by-step filling method, decomposing the filling calculation calibration process into two inversion processes: determining the second and third static correction quantities. The determination of the second static correction quantity corresponds to filling to the ground surface elevation position of the shot receiver point, while the determination of the third static correction quantity corresponds to filling to a fixed horizontal reference surface or a floating reference surface. Furthermore, because the filling calculation calibration process is decomposed, the two inversion processes can use the same or different first and second static correction methods, thus achieving the joint application of multiple static correction methods. This allows for comprehensive analysis through multiple combined static correction methods to determine the optimized calibration quantities.

[0035] In one optional implementation, a second static correction amount is determined based on a first static correction amount and a static correction equivalent model, including: Based on the model data of the static correction equivalent model, the filling thickness between the reference surface and the ground surface elevation is determined; based on the model data of the static correction equivalent model, the filling amount corresponding to the filling thickness is determined; based on the difference between the first static correction amount and the filling amount, the second static correction amount to fill to the ground surface elevation is determined.

[0036] Specifically, based on the model data of the static correction equivalent model, determining the fill thickness between the datum and the ground surface elevation means that, since establishing the static correction equivalent model requires determining the elevation of the fixed or floating datum and the ground surface elevation of the shot-receiver point, the fill thickness between the datum and the ground surface elevation is the difference between the elevation of the fixed or floating datum and the ground surface elevation of the shot-receiver point. Figure 2 As shown, that is H 基准面 and H 地表高程 The difference between them.

[0037] Specifically, determining the fill amount corresponding to the fill thickness based on the model data of the static correction equivalent model means that, since establishing the static correction equivalent model requires determining the corresponding replacement rate, determining the fill amount corresponding to the fill thickness is done by dividing the fill thickness by the corresponding replacement rate to determine the fill amount within the fill thickness range, i.e., the fill amount corresponding to the fill thickness.

[0038] Specifically, the determined second static correction amount is as follows: Figure 3 As shown in the first round of inversion by the medium-to-high-level analysis, based on the difference between the first static correction and the filling amount, the second static correction amount for filling to the surface elevation is determined to be the static correction amount for filling directly to the reference surface after the first static correction amount is stripped and then filled. After subtracting the filling amount corresponding to the filling thickness from this, the remaining filling amount corresponds to the filling amount from the top interface of the high-speed road to the surface elevation, which is the second static correction amount for filling to the surface elevation.

[0039] Specifically, the static correction equivalent model includes: a first static correction equivalent model, a second static correction equivalent model, and a third static correction equivalent model. The first static correction equivalent model refers to the firing point being activated within a low-velocity zone, such as... Figure 2 As shown, S represents the firing position of the shot point. The second static correction equivalent model refers to the shot point firing at the ground surface, and the third static correction equivalent model refers to the shot point firing below the high-velocity top. Therefore, based on the model data of the static correction equivalent model, the steps to determine the second static correction amount to fill to the ground surface elevation can be expressed by a formula.

[0040] When the static correction equivalent model is the first static correction equivalent model: The shot point correction amount for filling to the ground surface elevation is expressed by the following formula: T S=- T S剥 + T S填 =-( H S地表高程 - H s高速顶 - H S井深 ) / V s平均 +( H S地表高程 - H S高速顶 ) / V 替换速度 in, T S This represents the shot point correction amount filled to the reference plane elevation. T S剥 T represents the amount of stripping from the firing point to the high-velocity apex of the shot point. S填 This indicates the amount of fill from the shot point at the top elevation of the high-speed highway to the horizontal reference surface. H S地表高程 This indicates the elevation of the ground surface at the point where the plumb bob points upwards. H S高速顶 Indicates the elevation of the high-speed peak of the artillery position. H S井深 Indicates the depth of the firing point. V S平均 This represents the average velocity of the shot point from the firing point to the high-speed peak elevation of the shot point. V 替换速度 This indicates the rate of constant substitution in the calculation of the calibration. The correction amount for the receiver points filled to the ground surface elevation is expressed by the following formula: T R =- T R剥 + T R填 =-( H R地表高程 - H R高速顶 ) / V R平均 +( H R地表高程 - H R高速顶 ) / V 替换速度 in, T R This represents the correction amount for the receiver point filled to the elevation of the reference plane. T R剥T represents the amount of stripping from the receiver at the surface to the corresponding elevation of the high-speed peak at the receiver. R填 This indicates the amount of fill from the detector point at the top elevation of the high-speed highway to the horizontal reference surface. H R地表高程 This indicates the surface elevation at the receiver point. H R高速顶 This indicates the elevation of the high-speed peak corresponding to the detector point. V R平均 This represents the average velocity of the receiver point from the ground surface to the corresponding high-speed peak elevation.

[0041] When the static correction equivalent model is the second static correction equivalent model: The shot point correction amount for filling to the ground surface elevation is expressed by the following formula: T S =- T S剥 + T S填 =-( H S地表高程 - H s高速顶 ) / V s平均 +( H S基准面 - H S高速顶 ) / V 替换速度 in, H S基准面 Indicates the elevation of the reference surface at the firing point; The correction amount for the receiver points filled to the ground surface elevation is expressed by the following formula: T R =- T R剥 + T R填 =-( H R地表高程 - H R高速顶 ) / V R平均 +( H R基准面 - H R高速顶 ) / V 替换速度 in, H R基准面 This indicates the elevation of the reference surface of the receiver point.

[0042] Specifically, since the blasting point is located on the surface, the well depth does not need to be considered.

[0043] When the static correction equivalent model is the third static correction equivalent model: The shot point correction amount for filling to the ground surface elevation is expressed by the following formula: T S =+ T S填 = ( H 基准面 - H 高速顶 ) / V 替换速度 in, H 基准面 Indicates the elevation of the reference surface. H 高速顶 Indicates the altitude of the highway summit; The correction amount for the receiver points filled to the ground surface elevation is expressed by the following formula: T R =- T R剥 + T R填 =-( H R地表高程 - H R高速顶 ) / V R平均 +( H 基准面 - H 高速顶 ) / V 替换速度 Specifically, since the shot point excitation is below the high-speed peak in the third static correction equivalent model, the shot point correction calculation does not require a stripping process, but directly performs a filling process.

[0044] Specifically, the determined second static correction amount is as follows: Figure 3 As shown in the first round of inversion in the meso-to-surface analysis, since the calculation of the static correction involves a stripping and then filling process, and in this embodiment, the second static correction is the filling to the surface elevation, the process of obtaining the second static correction can also be represented by a formula. T S剥 T S填 , H S地表高程 , H S高速顶 , H S井深 , V S平均 , V 替换速度 , TR剥 T R填 , H R地表高程 , V R平均 All of these can be determined based on the pre-established static correction equivalent model, that is, the model data of the static correction equivalent model.

[0045] In one alternative implementation, determining the shot-receiver point information and first arrival based on the second static correction value includes: The second static correction is added to the original single-shot first arrival to form the first arrival; and, based on the second static correction, the corresponding shot point coordinates, receiver coordinates, shot point elevation, receiver elevation, and well depth are determined; the well depth is set to zero, and the set well depth, shot point coordinates, receiver coordinates, shot point elevation, and receiver elevation are used as shot-receiver information.

[0046] Specifically, adding the second static correction to the original single-shot initial arrival to form the initial arrival means adding the shot point correction and receiver point correction corresponding to the second static correction to the original single-shot initial arrival through data headers or custom software to form a new initial arrival.

[0047] Specifically, determining the corresponding shot point coordinates, receiver point coordinates, shot point elevation, receiver point elevation, and well depth based on the second static correction refers to determining the receiver point information and shot point information for the second static correction. Receiver point information and shot point information include: receiver point coordinates and elevation, shot point coordinates and elevation, and well depth.

[0048] Specifically, setting the well depth to preset data and using the set well depth, shot point coordinates, receiver point coordinates, shot point elevation, and receiver point elevation as shot-receiver point information means adjusting the well depth in the shot point information to zero, while keeping the other shot point and receiver point information unchanged, and using this as shot-receiver point information for subsequent determination of the third static correction. The preset data can be set according to actual working conditions; this embodiment does not impose specific limitations on this.

[0049] In one alternative implementation, setting the well depth to preset data includes setting the well depth to zero.

[0050] Specifically, the well depth is usually set to zero because the second static correction is the static correction amount for filling to the surface elevation. After filling to the surface elevation, the influence of the well depth has been eliminated, and the shot receiver information provides the data basis for determining the third static correction amount. The third static correction amount corresponds to the process of filling from the surface elevation to the reference surface. In this process, the same or different static correction methods may be adopted as the first static correction method. Therefore, it is necessary to set the well depth to zero to avoid the generation of errors and ensure the accuracy of the data in the step-by-step filling.

[0051] In one optional implementation, based on shot receiver information, the initial arrival and static correction equivalent model, a third static correction amount is determined from the surface elevation fill to the reference surface using a second static correction method, including: Based on the model data of the static correction equivalent model, the corresponding replacement velocity between the surface elevation and the reference surface is determined; based on the shot receiver information, first arrival and replacement velocity, the third static correction amount is determined by inversion calculation using the second static correction method, which fills the surface elevation to the reference surface.

[0052] Specifically, the determined third static correction is as follows: Figure 3 The second round of inversion of the first arrival wave in the meso-to-tomography shows that, based on the model data of the statically corrected equivalent model, the corresponding replacement rate between the surface elevation and the reference surface needs to be determined through the statically corrected equivalent model.

[0053] Specifically, based on shot receiver information, first arrival and replacement velocity, the third static correction amount is determined by inversion calculation using the second static correction method. This means inputting shot receiver information, first arrival and replacement velocity as known information into the software corresponding to the second static correction method, performing inversion calculation, and determining the third static correction amount from the surface elevation to the reference surface.

[0054] In one alternative implementation, after determining the third static correction amount from the surface elevation to the reference surface using a second static correction method based on shot receiver information, an equivalent model of initial arrival and static correction, the method further includes: The sum of the second and third static corrections is used as the new first static correction amount. Based on the new first static correction amount and the static correction equivalent model, a new second static correction amount is determined to fill to the ground surface elevation. Based on the new second static correction amount, new shot-receiver point information and new initial elevation are determined. Based on the new shot-receiver point information, the new initial elevation and the static correction equivalent model, a new third static correction amount is determined to fill from the ground surface elevation to the reference surface using the third static correction method, and the new third static correction amount is used as the third static correction amount.

[0055] Specifically, the process of determining the new second static correction amount based on the sum of the second and third static corrections as the new first static correction amount, and then determining the new second static correction amount to fill the surface elevation based on the new first static correction amount and the static correction equivalent model, is the same as in the above embodiments and will not be repeated here. The process of determining the new shot-receiver point information and the new first arrival based on the new second static correction amount is the same as in the above embodiments and will not be repeated here. The process of determining the new third static correction amount to fill the reference surface from the surface elevation using the third static correction method, based on the new shot-receiver point information, the new first arrival, and the static correction equivalent model, and using the new third static correction amount as the third static correction amount, is the same as in the above embodiments and will not be repeated here. The third static correction method is a static correction method different from the first and second static correction methods.

[0056] In one alternative implementation, the method further includes: Multiple first static correction methods and second static correction methods are used to iteratively determine multiple optimized correction values. These multiple optimized correction values ​​are then superimposed, and the profile imaging accuracy is compared with the attitude to determine the optimal optimized correction value.

[0057] Specifically, the method of using multiple different first static correction methods and second static correction methods to iteratively determine multiple optimized calibration values ​​is the same as in the above embodiments, and will not be repeated here.

[0058] Specifically, such as Figure 4 As shown, an exemplary illustration is presented of a first-order inversion of refracted wave static correction determined by the direct filling method in related technologies, and a CMP gather superimposed profile determined by the step-filling method of this embodiment, which combines refracted wave static correction with tomography-based second-order inversion. It can be seen that, by comparing the profile imaging accuracy and attitude, the static correction amount determined by the step-filling method of this embodiment is more effective.

[0059] like Figure 5 As shown, an exemplary model is presented, demonstrating the CMP gather overlay profiles of the first generalized linear inversion followed by tomographic inversion and the first tomographic inversion followed by generalized linear inversion. It can be seen that, through comparison of profile imaging accuracy and attitude, the static correction effect of the first tomographic inversion followed by generalized linear inversion is better.

[0060] By implementing this embodiment, data support can be provided for the establishment of a platform that comprehensively utilizes various static corrections, thereby enabling the rapid selection of the optimal combination of static correction methods and static correction amounts, and improving the progress and quality of seismic data processing.

[0061] This embodiment also provides an optimized calibration determination device for multi-static calibration linkage. This device is used to implement the above embodiments and preferred embodiments, and will not be repeated for details already described. As used below, the term "module" can be a combination of software and / or hardware that implements a predetermined function. Although the device described in the following embodiments is preferably implemented in software, hardware implementation, or a combination of software and hardware, is also possible and contemplated.

[0062] This embodiment provides an optimized calibration determination device for multiple static calibrations, such as... Figure 6 As shown, it includes: The first static correction amount determination module 201 is used to determine the first static correction amount to fill the reference surface based on a pre-established static correction equivalent model and using a first static correction method. For details, please refer to the relevant description of step S101 in the above embodiments, which will not be repeated here.

[0063] The second static correction determination module 202 is used to determine the second static correction amount to fill the ground surface elevation based on the first static correction amount and the static correction equivalent model. For details, please refer to the relevant description of step S102 in the above embodiments, which will not be repeated here.

[0064] The information determination module 203 is used to determine the shot-receiver point information and first arrival based on the second static correction value. For details, please refer to the description of step S103 in the above embodiments, which will not be repeated here.

[0065] The third static correction determination module 204 is used to determine the third static correction amount from the ground surface elevation to the reference surface based on the shot receiver information, the initial arrival and static correction equivalent model, and using the second static correction method. For details, please refer to the relevant description of step S104 in the above embodiments, which will not be repeated here.

[0066] The optimization calibration determination module 205 is used to determine the optimization calibration based on the second static calibration and the third static calibration. For details, please refer to the description of step S105 in the above embodiments, which will not be repeated here.

[0067] In this embodiment, the multi-static calibration linkage optimization calibration determination device is presented in the form of a functional unit. Here, a unit refers to an ASIC (Application Specific Integrated Circuit) circuit, a processor and memory that execute one or more software or fixed programs, and / or other devices that can provide the above functions.

[0068] This invention also provides a computer device having the above-described features. Figure 6 The device shown is an optimized calibration determination device for multiple static calibrations.

[0069] Please see Figure 7 , Figure 7 This is a schematic diagram of the structure of a computer device provided in an optional embodiment of the present invention, such as... Figure 7 As shown, the computer device includes one or more processors 301, memory 302, and interfaces for connecting the components, including high-speed interfaces and low-speed interfaces. The components communicate with each other via different buses and can be mounted on a common motherboard or otherwise installed as needed. The processors can process instructions executed within the computer device, including instructions stored in or on memory to display graphical information of a GUI on external input / output devices (such as display devices coupled to the interfaces). In some alternative implementations, multiple processors and / or multiple buses can be used with multiple memories and multiple memory modules, if desired. Similarly, multiple computer devices can be connected, each providing some of the necessary operations (e.g., as a server array, a group of blade servers, or a multiprocessor system). Figure 7 Take processor 301 as an example.

[0070] Processor 301 may be a central processing unit, a network processor, or a combination thereof. Processor 301 may further include a hardware chip. The hardware chip may be an application-specific integrated circuit (ASIC), a programmable logic device (PLD), or a combination thereof. The programmable logic device may be a complex programmable logic device (CAMP), a field-programmable gate array (FPGA), a general-purpose array logic (GDA), or any combination thereof.

[0071] The memory 302 stores instructions executable by at least one processor 301 to cause the at least one processor 301 to perform the method shown in the above embodiments.

[0072] Memory 302 may include a program storage area and a data storage area. The program storage area may store the operating system and applications required for at least one function; the data storage area may store data created based on the use of the computer device. Furthermore, memory 302 may include high-speed random access memory and may also include non-transitory memory, such as at least one disk storage device, flash memory device, or other non-transitory solid-state storage device. In some alternative embodiments, memory 302 may optionally include memory remotely located relative to processor 301, and this remote memory may be connected to the computer device via a network. Examples of such networks include, but are not limited to, the Internet, intranets, local area networks, mobile communication networks, and combinations thereof.

[0073] The memory 302 may include volatile memory, such as random access memory; the memory may also include non-volatile memory, such as flash memory, hard disk or solid-state drive; the memory 302 may also include combinations of the above types of memory. The computer device also includes a communication interface 303 for communicating with other devices or communication networks.

[0074] This invention also provides a computer-readable storage medium. The methods described above according to embodiments of the invention can be implemented in hardware or firmware, or implemented as computer code that can be recorded on a storage medium, or implemented as computer code downloaded via a network and originally stored on a remote storage medium or a non-transitory machine-readable storage medium and then stored on a local storage medium. Thus, the methods described herein can be processed by software stored on a storage medium using a general-purpose computer, a dedicated processor, or programmable or dedicated hardware. The storage medium can be a magnetic disk, optical disk, read-only memory, random access memory, flash memory, hard disk, or solid-state drive, etc.; further, the storage medium can also include combinations of the above types of memory. It is understood that computers, processors, microprocessor controllers, or programmable hardware include storage components capable of storing or receiving software or computer code, which, when accessed and executed by the computer, processor, or hardware, implements the methods shown in the above embodiments.

[0075] Although embodiments of the invention have been described in conjunction with the accompanying drawings, those skilled in the art can make various modifications and variations without departing from the spirit and scope of the invention, and such modifications and variations all fall within the scope defined by the appended claims.

Claims

1. A method for determining optimized calibration quantities in a multi-static calibration linkage, characterized in that, The method includes: Based on the pre-established static correction equivalent model, the first static correction amount to fill the reference surface is determined using the first static correction method. Based on the first static correction amount and the static correction equivalent model, a second static correction amount is determined to fill to the ground surface elevation. Based on the second static correction, the shot receiver information and initial arrival are determined; Based on the shot receiver information, the initial arrival and the static correction equivalent model, the third static correction amount from the ground elevation to the reference surface is determined using the second static correction method. Based on the second static correction value and the third static correction value, the optimized correction value is determined.

2. The method according to claim 1, characterized in that, The step of determining the second static correction amount based on the first static correction amount and the static correction equivalent model includes: Based on the model data of the static correction equivalent model, the filling thickness between the reference surface and the ground surface elevation is determined; Based on the model data of the static correction equivalent model, the filling amount corresponding to the filling thickness is determined; Based on the difference between the first static correction amount and the filling amount, a second static correction amount is determined to fill to the ground surface elevation.

3. The method according to claim 1, characterized in that, The determination of shot receiver information and first arrival based on the second static correction value includes: The second static correction is added to the original single-shot first arrival to form the first arrival; and, Based on the second static correction, the corresponding shot point coordinates, receiver point coordinates, shot point elevation, receiver point elevation, and well depth are determined. The well depth is set to preset data, and the set well depth, the shot point coordinates, the receiver point coordinates, the shot point elevation, and the receiver point elevation are used as shot-receiver point information.

4. The method according to claim 3, characterized in that, Setting the well depth to preset data includes setting the well depth to zero.

5. The method according to claim 1, characterized in that, The method of determining the third static correction amount from the surface elevation to the reference surface based on the shot receiver information, the initial arrival, and the static correction equivalent model, using a second static correction method, includes: Based on the model data of the static correction equivalent model, the corresponding replacement rate between the ground elevation and the reference surface is determined; Based on the shot receiver information, the initial arrival and the replacement velocity, an inversion calculation is performed using the second static correction method to determine the third static correction amount that fills the reference surface from the ground elevation.

6. The method according to claim 1, characterized in that, After determining the third static correction amount from the surface elevation to the reference surface using the second static correction method based on the shot receiver information, the initial arrival, and the static correction equivalent model, the method further includes: The sum of the second static correction and the third static correction is used as the new first static correction amount. Based on the new first static correction amount and the static correction equivalent model, a new second static correction amount is determined to fill the elevation to the ground surface. Based on the new second static correction, new shot receiver information and new first arrival are determined; Based on the new shot-receiver information, the new initial arrival and the static correction equivalent model, a new third static correction amount is determined by using the third static correction method to fill the reference surface from the ground elevation, and the new third static correction amount is used as the third static correction amount.

7. The method according to claim 1, characterized in that, The method further includes: Multiple different first static correction methods and second static correction methods are used to iteratively determine multiple optimized correction quantities; By superimposing multiple optimized calibrations and comparing the profile imaging accuracy with the attitude, the optimal optimized calibration is determined.

8. A device for determining optimized calibration quantities in multi-static calibration linkage, characterized in that, The device includes: The first static correction amount determination module is used to determine the first static correction amount to be filled to the reference surface based on the pre-established static correction equivalent model and using the first static correction method. The second static correction determination module is used to determine the second static correction amount to fill the ground surface elevation based on the first static correction amount and the static correction equivalent model. The information determination module is used to determine the shot receiver information and initial arrival based on the second static correction value; The third static correction determination module is used to determine the third static correction amount from the surface elevation to the reference surface based on the shot receiver information, the initial arrival and the static correction equivalent model, using the second static correction method. The optimization calibration determination module is used to determine the optimization calibration based on the second static calibration and the third static calibration.

9. A computer device, characterized in that, include: The system includes a memory and a processor, which are interconnected. The memory stores computer instructions, and the processor executes the computer instructions to perform the optimized calibration determination method for multi-static calibration linkage as described in any one of claims 1 to 7.

10. A computer-readable storage medium, characterized in that, The computer-readable storage medium stores computer instructions for causing the computer to execute the optimized calibration determination method for multiple static calibrations as described in any one of claims 1 to 7.

Images