Method and device for determining velocity model, electronic equipment and storage medium

Abstract

The application provides a speed model determination method and device, electronic equipment and a storage medium. The application determines the intersection point coordinate value between the straight line equation of the control point vertically downward and the plane coordinate equation corresponding to each three-dimensional horizon in the stratum through the coordinate information of the control point, then determines the intersection point speed corresponding to each intersection point coordinate value based on a pre-established speed function, and finally performs interpolation processing based on the control point speed and the intersection point speed to determine the speed model of the stratum, thereby improving the precision of three-dimensional horizon speed interpolation.

Description

Technical Field

[0001] This application relates to the field of seismic exploration technology, and in particular to a method, apparatus, electronic device and storage medium for determining a velocity model. Background Technology

[0002] The parameter validation of a seismic acquisition and observation system can be divided into three aspects: theoretical calculation, analysis of actual data, and model-based design. These three aspects corroborate each other. Theoretical calculation is the foundation of the observation system parameter validation, a crucial and indispensable step. Currently, mainstream acquisition software in China considers parameter validation based on 3D horizon data as a very important module. However, in practice, it is difficult to obtain dense control point information; sparseness leads to inaccurate 3D horizon velocity interpolation. Summary of the Invention

[0003] To address the aforementioned problems, this application provides a method, apparatus, electronic device, and storage medium for determining a velocity model, which can improve the accuracy of velocity interpolation.

[0004] This application provides a method for determining a velocity model, including:

[0005] Acquire control point information of control points in the formation, wherein the control point information includes: coordinate information and the control point velocity corresponding to the coordinate information;

[0006] Based on the coordinate information, determine the intersection point coordinates of the line equation vertically downward from the control point and the plane coordinate equations corresponding to each three-dimensional layer in the stratum.

[0007] Based on a pre-established velocity function, determine the velocity at each intersection point corresponding to its coordinate values.

[0008] Interpolation is performed based on the control point velocity and the intersection velocity to determine the velocity model of the formation.

[0009] In some embodiments, determining the intersection point coordinates of the intersection point between the equation of the vertically downward line from the control point and the plane coordinate equations corresponding to each three-dimensional layer in the stratum based on the coordinate information includes:

[0010] Based on the coordinate information, construct the equation of a straight line perpendicular to the control point downwards;

[0011] Determine the coordinates of the intersection point between the linear equation and the plane coordinate equations corresponding to each three-dimensional layer in the stratum.

[0012] In some embodiments, determining the velocity model of the formation based on the control point velocity and intersection data includes:

[0013] The velocity model of the formation is determined using the Kriging interpolation algorithm based on the control point velocity and intersection data.

[0014] In some embodiments, the method further includes:

[0015] Obtain sample control point information of sample control points, wherein the sample control point information includes: sample coordinate information and the actual velocity value corresponding to the sample coordinate information, and the actual velocity value corresponding to the sample coordinate information of multiple layers vertically downward from the sample coordinate information.

[0016] A velocity function is constructed based on the sample control point information.

[0017] In some embodiments, constructing the velocity function based on the sample control point information includes:

[0018] Construct a velocity function to be determined, wherein the velocity function to be determined includes: coefficients;

[0019] The sample coordinate information corresponding to the sample control point is input into the velocity function to be determined to determine the calculated velocity value corresponding to each coordinate information.

[0020] The sum of squares of the errors is determined based on each calculated speed value and the corresponding actual speed value;

[0021] Determine the minimum value function of the sum of squares of the errors;

[0022] Based on the coefficients, the partial derivatives of the minimum function are obtained to obtain the equations for calculating the coefficients;

[0023] Based on the calculation equation, the calculation results of the coefficients are determined;

[0024] The velocity function is determined based on the calculation results of the aforementioned coefficients.

[0025] This application provides a velocity model determination device, comprising:

[0026] The acquisition module is used to acquire control point information, coordinate information, and control point velocity corresponding to the coordinate information of control points in the stratum.

[0027] The first determining module is used to determine the intersection point coordinates of the intersection point between the vertically downward straight line equation of the control point and the plane coordinate equation corresponding to each three-dimensional layer in the stratum, based on the coordinate information.

[0028] The second determining module is used to determine the intersection velocity corresponding to the coordinate values ​​of each intersection point based on a pre-established velocity function;

[0029] The third determining module is used to determine the velocity model of the formation based on the control point velocity and the intersection velocity.

[0030] In some embodiments, determining the intersection point coordinates of the intersection point between the equation of the vertically downward line from the control point and the plane coordinate equations corresponding to each three-dimensional layer in the stratum based on the coordinate information includes:

[0031] Based on the coordinate information, construct the equation of a straight line perpendicular to the control point downwards;

[0032] Determine the coordinates of the intersection point between the linear equation and the plane coordinate equations corresponding to each three-dimensional layer in the stratum.

[0033] In some embodiments, determining the velocity model of the formation based on the control point velocity and intersection data includes:

[0034] The velocity model of the formation is determined using the Kriging interpolation algorithm based on the control point velocity and intersection data.

[0035] This application provides an electronic device, including a memory and a processor. The memory stores a computer program, which, when executed by the processor, performs the speed model determination method described in any of the above-mentioned embodiments.

[0036] This application provides a storage medium storing a computer program that can be executed by one or more processors and can be used to implement the speed model determination method described in any of the above claims.

[0037] This application provides a method, apparatus, electronic device, and storage medium for determining a velocity model. It determines the intersection coordinates between the vertically downward linear equation of the control point and the planar coordinate equations corresponding to each three-dimensional layer in the formation by using the coordinate information of the control point. Then, based on a pre-established velocity function, it determines the intersection velocity corresponding to each intersection coordinate value. Finally, it performs interpolation processing based on the control point velocity and the intersection velocity to determine the velocity model of the formation, thereby improving the accuracy of three-dimensional layer velocity interpolation. Attached Figure Description

[0038] The present application will be described in more detail below based on embodiments and with reference to the accompanying drawings.

[0039] Figure 1 A schematic diagram illustrating the implementation process of a method for determining a velocity model provided in an embodiment of this application;

[0040] Figure 2 A schematic diagram illustrating the implementation process of constructing a velocity function based on the sample control point information, provided in an embodiment of this application;

[0041] Figure 3A schematic diagram illustrating the implementation process of a method for determining a velocity model provided in an embodiment of this application;

[0042] Figure 4 This is a schematic diagram of the composition structure of the electronic device provided in the embodiments of this application.

[0043] In the accompanying drawings, the same parts are referred to by the same reference numerals, and the drawings are not drawn to scale. Detailed Implementation

[0044] To make the objectives, technical solutions, and advantages of this application clearer, the application will be further described in detail below with reference to the accompanying drawings. The described embodiments should not be regarded as limitations on this application. All other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of this application.

[0045] In the following description, references are made to “some embodiments,” which describe a subset of all possible embodiments. However, it is understood that “some embodiments” may be the same subset or different subsets of all possible embodiments and may be combined with each other without conflict.

[0046] If the application documents contain similar descriptions such as "first, second, third", the following explanation shall be added: In the following description, the terms "first, second, third" are used only to distinguish similar objects and do not represent a specific order of objects. It is understood that "first, second, third" may be interchanged in a specific order or sequence where permitted, so that the embodiments of this application described herein can be implemented in an order other than that illustrated or described herein.

[0047] Unless otherwise defined, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this application belongs. The terminology used herein is for the purpose of describing embodiments of this application only and is not intended to limit this application.

[0048] To address the problems existing in related technologies, this application provides a method for determining a speed model. This method is applied to electronic devices, such as computers and mobile terminals. The function implemented by the speed model determination method provided in this application can be achieved by the processor of the electronic device calling program code, wherein the program code can be stored in a computer storage medium.

[0049] Example 1

[0050] This application provides a method for determining a velocity model. Figure 1 This is a schematic diagram illustrating the implementation process of a method for determining a velocity model provided in an embodiment of this application, as shown below. Figure 1As shown, it includes the following steps:

[0051] Step S101: Obtain control point information of control points in the stratum. The control point information includes: coordinate information and the control point velocity corresponding to the coordinate information.

[0052] In this embodiment of the application, the control point refers to the measurement control point. The measurement control point refers to a series of points set up in the area to be measured before the measurement operation is carried out to complete the measurement operation of the entire area. When selecting the points, the existing topographic maps and control point data of the survey area are first investigated and collected.

[0053] In this embodiment of the application, measurement operations can be performed for each control point to obtain the coordinate information and speed of the control point.

[0054] In this embodiment of the application, the electronic device can communicate with the measuring device to obtain control point information of the control points in the formation from the measuring device.

[0055] In this embodiment of the application, there may be multiple control points, and the formation may be a three-dimensional formation.

[0056] Step S102: Based on the coordinate information, determine the coordinate values ​​of the intersection point between the vertical downward straight line equation of the control point and the plane coordinate equation corresponding to each three-dimensional layer in the stratum.

[0057] In this embodiment, a vertically downward straight line equation for the control point can be constructed based on the coordinate information, and the coordinate values ​​of the intersection points between the straight line equation and the planar coordinate equations corresponding to each three-dimensional layer in the stratum can be determined. In this embodiment, the stratum is a three-dimensional stratum, and the three-dimensional stratum can be subdivided into different depths to obtain each three-dimensional layer.

[0058] In this embodiment of the application, the strata can be divided into multiple three-dimensional layers. Each three-dimensional layer can be represented by a plane coordinate equation.

[0059] In this embodiment of the application, the intersection point coordinates of the intersection point between the linear equation and the plane coordinate equation corresponding to each three-dimensional layer in the stratum can be calculated.

[0060] Step S103: Based on the pre-established velocity function, determine the intersection velocity corresponding to the coordinate values ​​of each intersection point.

[0061] In this embodiment, the velocity function includes the calculated relationship between velocity and coordinates. In this embodiment, the velocity function can be established based on sample control points.

[0062] In this embodiment of the application, the coordinate values ​​of each intersection point can be input into the velocity function, thereby calculating the intersection velocity corresponding to each intersection point coordinate value.

[0063] In this embodiment of the application, the intersection velocity represents the velocity of the control point perpendicularly downward to the next three-dimensional layer plane.

[0064] Step S104: Perform interpolation based on the control point velocity and intersection velocity to determine the velocity model of the formation.

[0065] In this embodiment, the velocity model of the formation is determined by the Kriging interpolation algorithm based on the control point velocity and the intersection velocity.

[0066] In this embodiment of the application, all control points can be processed sequentially using the Kriging interpolation algorithm to determine the velocity corresponding to each control point, and based on the velocity corresponding to each control point, a velocity model of the formation can be obtained.

[0067] In this embodiment, the Kriging interpolation algorithm first considers the variation distribution of spatial attributes in spatial location to determine the distance range that affects the value of a point to be interpolated, and then uses sampling points within this range to estimate the attribute value of the point to be interpolated. Based on the different spatial locations of the samples and the different degrees of correlation between samples, different weights are assigned to each sample grade, and a moving weighted average is performed to estimate the average velocity of the central block segment.

[0068] This application provides a method for determining a velocity model. By using the coordinate information of control points, the method determines the intersection coordinates between the equation of the vertically downward line from the control points and the plane coordinate equations corresponding to each three-dimensional layer in the stratum. Then, based on a pre-established velocity function, the method determines the intersection velocity corresponding to each intersection coordinate value. Finally, the method performs interpolation processing based on the control point velocity and the intersection velocity to determine the velocity model of the stratum, which can improve the accuracy of three-dimensional layer velocity interpolation.

[0069] Example 2

[0070] Based on the foregoing embodiments, this application further provides a method for determining a velocity model, the method comprising:

[0071] Step S201: Obtain sample control point information of sample control points, wherein the sample control point information includes: sample coordinate information and the actual velocity value corresponding to the sample coordinate information, and the actual velocity value corresponding to the sample coordinate information of multiple layers vertically downward from the sample coordinate information.

[0072] In this embodiment of the application, the sample control information of the sample control point can be measured by a measuring device.

[0073] In this embodiment of the application, the sample control point information can be represented as: control point array v i[j], where i represents the control point number, i = {1, 2, 3, ..., n}, and J is the velocity array of each layer below the control point, n = {1, 2, 3, ..., m}.

[0074] In this embodiment of the application, there are n control points and m layers.

[0075] Step S202: Construct a velocity function based on the sample control point information.

[0076] In this embodiment of the application, step S202 can be implemented through the following steps. Figure 2 This application provides a schematic diagram of an implementation process for constructing a velocity function based on the sample control point information, as shown in the embodiment of the present application. Figure 2 As shown:

[0077] Step S2021: Input the sample coordinate information corresponding to the sample control point into the velocity function to be determined to determine the calculated velocity value corresponding to each coordinate information.

[0078] In this embodiment, the velocity function to be determined can be pre-established as a quadratic function equation. For example, the velocity function to be determined is expressed as:

[0079] h θ (x)=θ0+θ1x+θ2x 2 ;

[0080] Among them, h θ (x) Calculation speed. θ0, θ1, and θ2 are the parameters to be calculated, and x is the coordinate information (i.e., depth value).

[0081] In this embodiment of the application, the sample coordinate information corresponding to each control point can be input into the velocity function to be determined to determine the calculated velocity value corresponding to each coordinate information.

[0082] In this embodiment of the application, the calculated speed value is a calculated value including coefficients θ0, θ1 and θ2 to be calculated.

[0083] Step S2022: Determine the sum of squares of the errors based on each calculated speed value and the corresponding actual speed value.

[0084] In this embodiment, the difference between each calculated speed value and the corresponding actual speed value is obtained. Then, the square of the difference is calculated to obtain the square of each error. Finally, the squares are summed to obtain the sum of squares of the errors.

[0085] In this embodiment of the application, the sum of squares of the errors can be expressed as:

[0086]

[0087] Among them, h θ (x i ) represents the calculated speed value, and yi represents the actual speed value.

[0088] Step S2023: Determine the minimum value function of the sum of squares of the error.

[0089] In this embodiment of the application, the minimum value function can be expressed as:

[0090]

[0091] Step S2024: Based on the coefficients, take the partial derivative of the minimum value function to obtain the calculation equation for the coefficients.

[0092] In this embodiment of the application, the partial derivative formula can be expressed as:

[0093] In this embodiment, the minimum function can be transformed into a coefficient calculation equation using the partial derivative formula.

[0094] Step S2025: Based on the calculation equation, determine the calculation result of the coefficient.

[0095] In this embodiment of the application, the calculation equation can be solved to determine the calculation results of the coefficients.

[0096] Step S2026: Determine the velocity function based on the calculation results of the coefficients.

[0097] In this embodiment of the application, the obtained coefficients are substituted into the velocity function to be determined, and the optimal velocity function can be determined.

[0098] Step S203: Obtain control point information of control points in the stratum. The control point information includes: coordinate information and the control point velocity corresponding to the coordinate information.

[0099] In this embodiment of the application, the control point refers to the measurement control point. The measurement control point refers to a series of points set up in the area to be measured before the measurement operation is carried out to complete the measurement operation of the entire area. When selecting the points, the existing topographic maps and control point data of the survey area are first investigated and collected.

[0100] In this embodiment of the application, measurement operations can be performed for each control point to obtain the coordinate information and speed of the control point.

[0101] In this embodiment of the application, the electronic device can communicate with the measuring device to obtain control point information of the control points in the formation from the measuring device.

[0102] In this embodiment, there may be multiple control points. The formation is a three-dimensional formation.

[0103] Step S204: Based on the coordinate information, determine the coordinate values ​​of the intersection point between the equation of the vertically downward line of the control point and the plane coordinate equation corresponding to each three-dimensional layer in the stratum.

[0104] In this embodiment of the application, a straight line equation perpendicular to the control point can be constructed based on the coordinate information, and the coordinate values ​​of the intersection point between the straight line equation and the plane coordinate equation corresponding to each three-dimensional layer in the stratum can be determined.

[0105] In this embodiment of the application, the strata can be divided into multiple three-dimensional layers. Each three-dimensional layer can be represented by a plane coordinate equation.

[0106] In this embodiment of the application, the intersection point coordinates of the intersection point between the linear equation and the plane coordinate equation corresponding to each three-dimensional layer in the stratum can be calculated.

[0107] Step S205: Based on the pre-established velocity function, determine the intersection velocity corresponding to the coordinate values ​​of each intersection point.

[0108] In this embodiment, the velocity function includes the calculated relationship between velocity and coordinates. In this embodiment, the velocity function can be established based on sample control points.

[0109] In this embodiment of the application, the coordinate values ​​of each intersection point can be input into the velocity function, thereby calculating the intersection velocity corresponding to each intersection point coordinate value.

[0110] In this embodiment of the application, the velocity of the intersection point represents the velocity of the control point perpendicularly downward to the next three-dimensional layer plane.

[0111] Step S206: Perform interpolation based on the control point velocity and intersection velocity to determine the velocity model of the formation.

[0112] In this embodiment, the velocity model of the formation is determined by the Kriging interpolation algorithm based on the control point velocity and the intersection velocity.

[0113] In this embodiment of the application, all control points can be processed sequentially using the Kriging interpolation algorithm to determine the velocity corresponding to each control point, and based on the velocity corresponding to each control point, a velocity model of the formation can be obtained.

[0114] In this embodiment, the Kriging interpolation algorithm first considers the variation distribution of spatial attributes in spatial location to determine the distance range that affects the value of a point to be interpolated, and then uses sampling points within this range to estimate the attribute value of the point to be interpolated. Based on the different spatial locations of the samples and the different degrees of correlation between samples, different weights are assigned to each sample grade, and a moving weighted average is performed to estimate the average velocity of the central block segment.

[0115] This application provides a method for determining a velocity model. By using the coordinate information of control points, the method determines the intersection coordinates between the equation of the vertically downward line from the control points and the plane coordinate equations corresponding to each three-dimensional layer in the stratum. Then, based on a pre-established velocity function, the method determines the intersection velocity corresponding to each intersection coordinate value. Finally, the method performs interpolation processing based on the control point velocity and the intersection velocity to determine the velocity model of the stratum, which can improve the accuracy of three-dimensional layer velocity interpolation.

[0116] Example 3

[0117] Based on the foregoing embodiments, this application provides another method for determining a velocity model. This method constructs a velocity model for the entire three-dimensional layer data by using the velocities of different layers below the control point coordinates, thus solving the problem of layer velocity interpolation accuracy when there are too few control points. Figure 3 This is a schematic diagram illustrating the implementation process of a method for determining a velocity model provided in an embodiment of this application, as shown below. Figure 3 As shown, the method includes:

[0118] Step S301: Construct the velocity function below the control point using the least squares method.

[0119] In this embodiment of the application, a control point array can be obtained, the control point array v i [j] includes: coordinate values ​​and velocity values ​​for all downward layers, where i is the control point number, i = {1, 2, 3, ..., n}.

[0120] J is the array of layer velocities below the control point, n = {1, 2, 3, ..., m}.

[0121] In this embodiment of the application, the velocity function of each control point can be constructed using the least squares method.

[0122] In this embodiment of the application, h θ (x)=θ0+θ1x+θ2x 2 , where h θ (x i Let be the velocity function to be constructed. Substituting each depth value at the control point into the function yields the velocity value at the current depth. Taking the partial derivative of the accumulated result gives the optimal quadratic function. The coefficients of the constructed velocity function can be calculated using the following formula, thus determining the velocity function. The formula includes:

[0123]

[0124]

[0125] Step S302: Calculate the coordinates of the intersection point between the control point and the three-dimensional layer vertically downwards.

[0126] In this embodiment, all meshes of the 3D stratigraphic data are traversed. If the projected coordinates of a control point are within the current mesh, the planar coordinate equation of the current mesh is returned. Then, the equation of the line perpendicularly downwards from the control point is constructed, and the intersection point (z-value) with the planar coordinate equation of the mesh is calculated. The above steps are repeated to traverse and calculate the z-values ​​of the intersection points of all control points and the stratigraphic data.

[0127] Step S303: Calculate the velocity value at the intersection point based on the velocity function.

[0128] In this embodiment of the application, the z-value of the intersection point is substituted into the velocity function constructed in step S301 to obtain the velocity value at the corresponding depth.

[0129] It can iterate through all control points and calculate the corresponding speed value.

[0130] Step S304: Use the Kriging interpolation algorithm to interpolate the three-dimensional layer velocity.

[0131] In this embodiment of the application, the velocity value can be substituted into the Kriging interpolation algorithm to generate the velocity value of the three-dimensional stratigraphic grid, thereby obtaining the velocity model of the strata.

[0132] This application provides a method for determining a velocity model, which constructs a velocity model for the entire three-dimensional layer data by using different layer velocities below the control point coordinates, thus solving the problem of layer velocity interpolation accuracy when there are too few control points.

[0133] Example 4

[0134] Based on the foregoing embodiments, this application provides a speed model determination device. The various modules and units included in the device can be implemented by a processor in a computer device; of course, they can also be implemented by specific logic circuits. In the implementation process, the processor can be a central processing unit (CPU), a microprocessor (MPU), a digital signal processor (DSP), or a field programmable gate array (FPGA), etc.

[0135] This application provides a velocity model determination device, which includes:

[0136] The acquisition module is used to acquire control point information of control points in the stratum, wherein the control point information includes: coordinate information and the control point velocity corresponding to the coordinate information;

[0137] The first determining module is used to determine the intersection point coordinates of the intersection point between the vertically downward straight line equation of the control point and the plane coordinate equation corresponding to each three-dimensional layer in the stratum, based on the coordinate information.

[0138] The second determining module is used to determine the intersection velocity corresponding to the coordinate values ​​of each intersection point based on a pre-established velocity function;

[0139] The third determining module is used to perform interpolation processing based on the control point velocity and the intersection velocity to determine the velocity model of the formation.

[0140] In this embodiment of the application, the control point refers to the measurement control point. The measurement control point refers to a series of points set up in the area to be measured before the measurement operation is carried out to complete the measurement operation of the entire area. When selecting the points, the existing topographic maps and control point data of the survey area are first investigated and collected.

[0141] In this embodiment of the application, measurement operations can be performed for each control point to obtain the coordinate information and speed of the control point.

[0142] In this embodiment of the application, the electronic device can communicate with the measuring device to obtain control point information of the control points in the formation from the measuring device.

[0143] In this embodiment of the application, there may be multiple control points, and the formation may be a three-dimensional formation.

[0144] In this embodiment of the application, the coordinates of the intersection point between the vertically downward line equation of the control point and the plane coordinate equation corresponding to each three-dimensional layer in the stratum are determined based on the coordinate information.

[0145] In this embodiment of the application, the strata can be divided into multiple three-dimensional layers. Each three-dimensional layer can be represented by a plane coordinate equation.

[0146] In this embodiment of the application, the intersection point coordinates of the intersection point between the linear equation and the plane coordinate equation corresponding to each three-dimensional layer in the stratum can be calculated.

[0147] In this embodiment, the velocity function includes the calculated relationship between velocity and coordinates. In this embodiment, the velocity function can be established based on sample control points.

[0148] In this embodiment of the application, the coordinate values ​​of each intersection point can be input into the velocity function, thereby calculating the intersection velocity corresponding to each intersection point coordinate value.

[0149] In this embodiment of the application, the intersection velocity represents the velocity of the control point perpendicularly downward to the next three-dimensional layer plane.

[0150] In this embodiment, the velocity model of the formation is determined by the Kriging interpolation algorithm based on the control point velocity and the intersection velocity.

[0151] In this embodiment of the application, all control points can be processed sequentially using the Kriging interpolation algorithm to determine the velocity corresponding to each control point, and based on the velocity corresponding to each control point, a velocity model of the formation can be obtained.

[0152] In this embodiment, the Kriging interpolation algorithm first considers the variation distribution of spatial attributes in spatial location to determine the distance range that affects the value of a point to be interpolated, and then uses sampling points within this range to estimate the attribute value of the point to be interpolated. Based on the different spatial locations of the samples and the different degrees of correlation between samples, different weights are assigned to each sample grade, and a moving weighted average is performed to estimate the average velocity of the central block segment.

[0153] This application provides a velocity model determination device that determines the intersection coordinates between the vertically downward linear equation of the control point and the planar coordinate equations corresponding to each three-dimensional layer in the stratum using the coordinate information of the control point. Then, based on a pre-established velocity function, the intersection velocity corresponding to each intersection coordinate value is determined. Finally, interpolation processing is performed based on the control point velocity and the intersection velocity to determine the velocity model of the stratum, which can improve the accuracy of three-dimensional layer velocity interpolation.

[0154] In some embodiments, determining the intersection point coordinates of the intersection point between the equation of the vertically downward line from the control point and the plane coordinate equations corresponding to each three-dimensional layer in the stratum based on the coordinate information includes:

[0155] Based on the coordinate information, construct the equation of a straight line perpendicular to the control point downwards;

[0156] Determine the coordinates of the intersection point between the linear equation and the plane coordinate equations corresponding to each three-dimensional layer in the stratum.

[0157] In some embodiments, the interpolation process based on the control point velocity and the intersection velocity to determine the velocity model of the formation includes:

[0158] The velocity model of the formation is determined using the Kriging interpolation algorithm based on the control point velocity and intersection data.

[0159] In some embodiments, the velocity model determining device is further used for:

[0160] Obtain sample control point information of sample control points, wherein the sample control point information includes: sample coordinate information and the actual velocity value corresponding to the sample coordinate information, and the actual velocity value corresponding to the sample coordinate information of multiple layers vertically downward from the sample coordinate information.

[0161] A velocity function is constructed based on the sample control point information.

[0162] In some embodiments, constructing the velocity function based on the sample control point information includes:

[0163] Construct a velocity function to be determined, wherein the velocity function to be determined includes: coefficients;

[0164] The sample coordinate information corresponding to the sample control point is input into the velocity function to be determined to determine the calculated velocity value corresponding to each coordinate information.

[0165] The sum of squares of the errors is determined based on each calculated speed value and the corresponding actual speed value;

[0166] Determine the minimum value function of the sum of squares of the errors;

[0167] Based on the coefficients, the partial derivatives of the minimum function are obtained to obtain the equations for calculating the coefficients;

[0168] Based on the calculation equation, the calculation results of the coefficients are determined;

[0169] The velocity function is determined based on the calculation results of the aforementioned coefficients.

[0170] It should be noted that, in the embodiments of this application, if the above-described speed model determination method is implemented as a software functional module and sold or used as an independent product, it can also be stored in a computer-readable storage medium. Based on this understanding, the technical solution of the embodiments of this application, or the part that contributes to the prior art, can be embodied in the form of a software product. This computer software product is stored in a storage medium and includes several instructions to cause a computer device (which may be a personal computer, server, or network device, etc.) to execute all or part of the methods described in the various embodiments of this application. The aforementioned storage medium includes various media capable of storing program code, such as USB flash drives, portable hard drives, read-only memory (ROM), magnetic disks, or optical disks. Thus, the embodiments of this application are not limited to any specific hardware and software combination.

[0171] Accordingly, this application provides a storage medium storing a computer program thereon, characterized in that the computer program, when executed by a processor, implements the steps in the speed model determination method provided in the above embodiments.

[0172] Example 4

[0173] This application provides an electronic device; Figure 4 This is a schematic diagram of the composition structure of the electronic device provided in the embodiments of this application, such as... Figure 4 As shown, the electronic device 700 includes: a processor 701, at least one communication bus 702, a user interface 703, at least one external communication interface 704, and a memory 705. The communication bus 702 is configured to enable communication between these components. The user interface 703 may include a display screen, and the external communication interface 704 may include standard wired and wireless interfaces. The processor 701 is configured to execute a program for determining a speed model stored in the memory, to implement the steps in the speed model determination method provided in the above embodiment.

[0174] The descriptions of the display device and storage medium embodiments above are similar to those of the method embodiments above, and have similar beneficial effects. For technical details not disclosed in the computer device and storage medium embodiments of this application, please refer to the descriptions of the method embodiments of this application for understanding.

[0175] It should be noted that the descriptions of the storage medium and device embodiments above are similar to the descriptions of the method embodiments above, and have similar beneficial effects. For technical details not disclosed in the storage medium and device embodiments of this application, please refer to the descriptions of the method embodiments of this application for understanding.

[0176] It should be understood that the phrase "one embodiment" or "an embodiment" throughout the specification means that a specific feature, structure, or characteristic related to the embodiment is included in at least one embodiment of this application. Therefore, "in one embodiment" or "in an embodiment" appearing throughout the specification does not necessarily refer to the same embodiment. Furthermore, these specific features, structures, or characteristics can be combined in any suitable manner in one or more embodiments. It should be understood that in the various embodiments of this application, the sequence numbers of the above-described processes do not imply a sequential order of execution; the execution order of each process should be determined by its function and internal logic, and should not constitute any limitation on the implementation process of the embodiments of this application. The sequence numbers of the above-described embodiments are merely descriptive and do not represent the superiority or inferiority of the embodiments.

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

[0178] In the several embodiments provided in this application, it should be understood that the disclosed devices and methods can be implemented in other ways. The device embodiments described above are merely illustrative. For example, the division of units is only a logical functional division, and in actual implementation, there may be other division methods, such as: multiple units or components can be combined, or integrated into another system, or some features can be ignored or not executed. In addition, the coupling, direct coupling, or communication connection between the various components shown or discussed can be through some interfaces, and the indirect coupling or communication connection between devices or units can be electrical, mechanical, or other forms.

[0179] The units described above as separate components may or may not be physically separate. The components shown as units may or may not be physical units. They may be located in one place or distributed across multiple network units. Some or all of the units may be selected to achieve the purpose of this embodiment according to actual needs.

[0180] In addition, each functional unit in the various embodiments of this application can be integrated into one processing unit, or each unit can be a separate unit, or two or more units can be integrated into one unit; the integrated unit can be implemented in hardware or in the form of hardware plus software functional units.

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

[0182] Alternatively, if the integrated units described above are implemented as software functional modules and sold or used as independent products, they can also be stored in a computer-readable storage medium. Based on this understanding, the technical solutions of the embodiments of this application, or the parts that contribute to the prior art, can be embodied in the form of a software product. This computer software product is stored in a storage medium and includes several instructions to cause a controller to execute all or part of the methods described in the various embodiments of this application. The aforementioned storage medium includes various media capable of storing program code, such as mobile storage devices, ROMs, magnetic disks, or optical disks.

[0183] The above description is merely an embodiment of this application, but the scope of protection of this application is not limited thereto. Any variations or substitutions that can be easily conceived by those skilled in the art within the scope of the technology disclosed in this application should be included within the scope of protection of this application. Therefore, the scope of protection of this application should be determined by the scope of the claims.

Claims

1. A method for determining a velocity model, characterized in that, include: Acquire control point information of control points in the formation, wherein the control point information includes: coordinate information and the control point velocity corresponding to the coordinate information; Based on the coordinate information, determine the intersection point coordinates of the line equation vertically downward from the control point and the plane coordinate equations corresponding to each three-dimensional layer in the stratum. Based on a pre-established velocity function, determine the velocity at each intersection point corresponding to its coordinate values. Interpolation is performed based on the control point velocity and the intersection velocity to determine the velocity model of the formation; The step of determining the intersection point coordinates between the equation of the vertically downward line from the control point and the plane coordinate equations corresponding to each three-dimensional layer in the stratum based on the coordinate information includes: Based on the coordinate information, construct the equation of a straight line perpendicular to the control point downwards; Determine the intersection point coordinates of the linear equation and the plane coordinate equations corresponding to each three-dimensional layer in the stratum.

2. The method according to claim 1, characterized in that, The process of interpolating the control point velocity and intersection velocity to determine the velocity model of the formation includes: The velocity model of the formation is determined using the Kriging interpolation algorithm based on the control point velocity and intersection data.

3. The method according to claim 1, characterized in that, The method further includes: Obtain sample control point information of sample control points, wherein the sample control point information includes: sample coordinate information and the actual velocity value corresponding to the sample coordinate information, and the actual velocity value corresponding to the sample coordinate information of multiple layers vertically downward from the sample coordinate information. A velocity function is constructed based on the sample control point information.

4. The method according to claim 3, characterized in that, The construction of the velocity function based on the sample control point information includes: Construct a velocity function to be determined, wherein the velocity function to be determined includes: coefficients; The sample coordinate information corresponding to the sample control point is input into the velocity function to be determined to determine the calculated velocity value corresponding to each coordinate information. The sum of squares of the errors is determined based on each calculated speed value and its corresponding actual speed value; Determine the minimum value function of the sum of squares of the errors; Based on the coefficients, the partial derivatives of the minimum function are obtained to obtain the equations for calculating the coefficients; Based on the calculation equation, the calculation results of the coefficients are determined; The velocity function is determined based on the calculation results of the aforementioned coefficients.

5. A device for determining a velocity model, characterized in that, include: The acquisition module is used to acquire control point information, coordinate information, and control point velocity corresponding to the coordinate information of control points in the stratum. The first determining module is used to determine the intersection point coordinates of the intersection point between the vertically downward straight line equation of the control point and the plane coordinate equation corresponding to each three-dimensional layer in the stratum, based on the coordinate information. The second determining module is used to determine the intersection velocity corresponding to the coordinate values ​​of each intersection point based on a pre-established velocity function; The third determining module is used to perform interpolation processing based on the control point velocity and the intersection velocity to determine the velocity model of the formation. The step of determining the intersection point coordinates between the equation of the vertically downward line from the control point and the plane coordinate equations corresponding to each three-dimensional layer in the stratum based on the coordinate information includes: Based on the coordinate information, construct the equation of a straight line perpendicular to the control point downwards; Determine the coordinates of the intersection point between the linear equation and the plane coordinate equations corresponding to each three-dimensional layer in the stratum.

6. The velocity model determining device according to claim 5, characterized in that, The process of interpolating the control point velocity and intersection velocity to determine the velocity model of the formation includes: The velocity model of the formation is determined using the Kriging interpolation algorithm based on the control point velocity and intersection data.

7. An electronic device, characterized in that, It includes a memory and a processor, wherein the memory stores a computer program that, when executed by the processor, performs the method for determining the speed model as described in any one of claims 1 to 4.

8. A storage medium, characterized in that, The computer program stored in the storage medium can be executed by one or more processors and can be used to implement the method for determining the speed model as described in any one of claims 1 to 4.

Images