A two-dimensional offset seismic profile generation method, device, equipment and medium

Abstract

The present application relates to the field of oil and gas exploration seismic data processing, and discloses a two-dimensional migration seismic profile generation method, device, equipment and medium, which can select a target underground area including a volcanic channel according to the volcanic channel, perform numerical simulation based on the two-dimensional vertical profile of the target underground area, generate a volcanic deformation matrix, and then generate a two-dimensional migration seismic profile corresponding to the two-dimensional vertical profile based on the volcanic deformation matrix. The present application can effectively realize numerical simulation of a two-dimensional migration seismic profile containing a volcanic channel, and can improve the accuracy of the generated two-dimensional migration seismic profile containing a volcanic channel.

Description

Technical Field

[0001] The present invention relates to the field of oil and gas exploration seismic data processing, and in particular to a two-dimensional offset seismic profile generation method, device, equipment and medium. Background Art

[0002] A two-dimensional migration seismic profile refers to a profile obtained by performing migration processing on seismic data obtained through two-dimensional seismic exploration technology.

[0003] 2D migration seismic profiles are widely used in various geological exploration scenarios, especially oil and gas exploration. They enable exploration workers to more accurately understand and interpret subsurface structures, identify stratigraphic interfaces and structural features, and optimize oil and gas extraction plans.

[0004] However, the two-dimensional offset seismic profiles generated by related technologies for underground areas containing volcanic channels have low accuracy. Summary of the Invention

[0005] The present invention provides a two-dimensional migration seismic profile generation method, device, equipment and medium, which are used to solve the defect of low accuracy of two-dimensional migration seismic profiles containing volcanic channels generated in related technologies, and improve the accuracy of two-dimensional migration seismic profiles containing volcanic channels.

[0006] In a first aspect, the present invention provides a method for generating a two-dimensional migration seismic profile, comprising:

[0007] Meshing the target underground area in a two-dimensional vertical section in the horizontal and vertical directions to obtain a corresponding meshing interface; wherein the two-dimensional vertical section includes a volcanic channel;

[0008] Determining the number of horizontal grid rows and the number of vertical grid columns in the grid division interface, and constructing a zero matrix of corresponding dimensions according to the number of horizontal grid rows and the number of vertical grid columns as a matrix to be assigned;

[0009] According to the location distribution of the volcanic channel, determining a vertical starting grid corresponding to the volcanic channel on the grid division interface, and generating a left boundary curve and a right boundary curve corresponding to the volcanic channel;

[0010] Traversing all the grids from the vertical starting grid to the last grid on the vertical line in sequence along the vertical line where the vertical starting grid is located;

[0011] Assigning values ​​to the elements in the matrix to be assigned according to each of the traversed grids, obtaining a matrix after the assignment is completed and serving as a volcano deformation matrix;

[0012] generate a two-dimensional offset seismic profile corresponding to the two-dimensional vertical profile according to the volcanic deformation matrix.

[0013] Optionally, the assigning values to elements in the matrix to be valued according to each grid traversed to obtain a matrix after the values are assigned, comprises:

[0014] For the first grid traversed, it is determined whether the first grid reaches the bottom of the volcanic conduit;

[0015] If it is determined that the first grid does not reach the bottom of the volcanic conduit, a volcanic deformation value corresponding to the first grid is calculated, the elements in the matrix to be valued are assigned using the volcanic deformation value, a new matrix to be valued is obtained, and a second grid is traversed.

[0016] For the second grid traversed, it is determined whether the second grid reaches the bottom of the volcanic conduit until it is determined that the grid newly traversed reaches the bottom of the volcanic conduit, and the matrix newly to be valued is taken as the matrix after the values are assigned.

[0017] Optionally, the determining whether the first grid reaches the bottom of the volcanic conduit comprises:

[0018] According to the vertical coordinates of the first grid, corresponding left boundary points and right boundary points are determined on the left boundary curve and the right boundary curve respectively;

[0019] The left boundary points and the right boundary points are vertically displaced by perturbing the vertical coordinates of the left boundary points and the right boundary points respectively to obtain perturbed left boundary points and perturbed right boundary points;

[0020] An intermediate position point is randomly created between the perturbed left boundary points and the perturbed right boundary points;

[0021] The first distance and the second distance between the intermediate position point and the perturbed left boundary points and the perturbed right boundary points respectively are calculated;

[0022] If the first distance and the second distance are both greater than 0, it is determined that the first grid does not reach the bottom of the volcanic conduit;

[0023] If at least one of the first distance and the second distance is not greater than 0, it is determined that the first grid reaches the bottom of the volcanic conduit.

[0024] Optionally, the calculating the volcanic deformation value corresponding to the first grid comprises:

[0025] determining a projection point of the intermediate position point to a target line segment, the target line segment being from the left boundary point to the right boundary point;

[0026] constructing a first Bezier function according to the left boundary point, the projection point and the intermediate position point, and constructing a second Bezier function according to the right boundary point, the projection point and the intermediate position point;

[0027] obtaining horizontal coordinates of the left boundary point and the right boundary point, inputting the horizontal coordinate of the left boundary point into the first Bezier function to calculate a first function value, and inputting the horizontal coordinate of the right boundary point into the second Bezier function to calculate a second function value;

[0028] respectively determining the first function value and the second function value as first and second volcano-shaped deformation values corresponding to the left boundary point and the right boundary point, and taking the first and second volcano-shaped deformation values as a whole as a volcano-shaped deformation value corresponding to the first grid.

[0029] Optionally, the using the volcano-shaped deformation value to assign values to elements in the to-be-assigned matrix to obtain a new to-be-assigned matrix, comprises:

[0030] In the grid division interface, respectively determining a first grid and a second grid corresponding to the position coordinates of the left boundary point and the position coordinates of the right boundary point;

[0031] In the to-be-assigned matrix, respectively determining a first element position and a second element position corresponding to the first grid and the second grid;

[0032] respectively replacing zero values at the first element position and the second element position with the first and second volcano-shaped deformation values to obtain a new to-be-assigned matrix.

[0033] Optionally, the generating a two-dimensional offset seismic profile corresponding to the two-dimensional vertical profile according to the volcano-shaped deformation matrix, comprises:

[0034] obtaining position data of a plurality of key horizons in the target underground region and elastic parameter data of the underground medium;

[0035] generating a reflectivity matrix corresponding to the two-dimensional vertical profile according to the position data of each key horizon, the elastic parameter data, the number of horizontal grid rows and the number of vertical grid columns;

[0036] introducing local fold features and fault features in the reflectivity matrix to obtain a fault reflectivity matrix;

[0037] mapping the fault reflectivity matrix based on the volcano deformation matrix to obtain a volcano reflectivity matrix;

[0038] creating a seismic wavelet function corresponding to the two-dimensional vertical profile, and performing convolution between the volcano reflectivity matrix and the seismic wavelet function to obtain the two-dimensional migration seismic profile.

[0039] Optionally, the position data of each key horizon includes position coordinates of a plurality of horizon points, and the elastic parameter data includes an elastic parameter value corresponding to each key horizon;

[0040] The generating of the reflectivity matrix corresponding to the two-dimensional vertical profile according to the position data of each key horizon, the elastic parameter data, the number of horizontal grid rows and the number of vertical grid columns includes:

[0041] According to the position coordinates of the horizon points in the position data of each key horizon and the number of horizontal grid rows, the horizon points are expanded to generate a horizon point array of each key horizon, and the number of elements in each horizon point array is equal to the number of horizontal grid rows;

[0042] According to the elastic parameter data and the number of vertical grid columns, the elastic parameter value is expanded to obtain an elastic parameter array, and the number of elements in the elastic parameter array is equal to the number of vertical grid columns;

[0043] Global fold features are introduced into each horizon array to obtain a plurality of global fold horizon arrays corresponding thereto;

[0044] A two-dimensional horizon matrix is generated according to the plurality of global fold horizon arrays;

[0045] A two-dimensional elastic parameter matrix is generated according to the two-dimensional horizon matrix and the elastic parameter array;

[0046] The reflectivity matrix is generated based on the two-dimensional elastic parameter matrix.

[0047] In a second aspect, the present application provides a two-dimensional migration seismic profile generation device, which comprises:

[0048] A division unit is configured to perform grid division along a horizontal direction and a vertical direction on a two-dimensional vertical profile of a target underground area to obtain a corresponding grid division interface, wherein the two-dimensional vertical profile includes a volcano channel;

[0049] A determination unit is configured to determine a number of horizontal grid rows and a number of vertical grid columns in the grid division interface;

[0050] A construction unit is configured to construct a zero matrix of a corresponding dimension according to the number of horizontal grid rows and the number of vertical grid columns, and use the zero matrix as a to-be-assigned matrix.

[0051] a processing unit, configured to determine, on the grid division interface, a vertical starting grid corresponding to the volcanic channel according to the position distribution of the volcanic channel, and to generate a left boundary curve and a right boundary curve corresponding to the volcanic channel;

[0052] a traversal unit, configured to traverse all the grids from the vertical starting grid to the last grid on the vertical line in sequence along the vertical line where the vertical starting grid is located;

[0053] An assignment unit, configured to assign values ​​to elements in the matrix to be assigned according to each of the traversed grids, to obtain a matrix after the assignment is completed;

[0054] As a unit, it is used to use the matrix after the assignment as the volcano deformation matrix;

[0055] A generating unit is used to generate a two-dimensional offset seismic profile corresponding to the two-dimensional vertical profile according to the volcanic deformation matrix.

[0056] In a third aspect, the present invention provides a computer device comprising: a memory and a processor, the memory and the processor being communicatively connected to each other, the memory storing computer instructions, and the processor executing the two-dimensional offset seismic profile generation method of the first aspect or any corresponding embodiment thereof by executing the computer instructions.

[0057] In a fourth aspect, the present invention provides a computer-readable storage medium having computer instructions stored thereon, the computer instructions being used to enable a computer to execute the two-dimensional offset seismic profile generation method of the first aspect or any corresponding embodiment thereof.

[0058] The method, device, equipment, and medium for generating 2D migration seismic profiles provided by the present invention can select a target underground area containing a volcanic channel based on the volcanic channel, perform numerical simulation based on the 2D vertical profile of the target underground area, generate a volcanic deformation matrix, and then generate a 2D migration seismic profile corresponding to the 2D vertical profile based on the volcanic deformation matrix. This invention can effectively perform numerical simulation of 2D migration seismic profiles containing volcanic channels and improve the accuracy of the generated 2D migration seismic profiles containing volcanic channels. BRIEF DESCRIPTION OF THE DRAWINGS

[0059] In order to more clearly illustrate the technical solutions in the present invention or the prior art, a brief introduction is given below to the drawings required for use in the embodiments or the description of the prior art. Obviously, the drawings described below are some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying any creative work.

[0060] Figure 1 A flow chart of a method for generating a two-dimensional migration seismic profile provided by an embodiment of the present invention;

[0061] Figure 2 A schematic cross-sectional view corresponding to a fault reflectivity matrix provided in an embodiment of the present invention;

[0062] Figure 3 A schematic cross-sectional view corresponding to a volcanic deformation matrix provided by an embodiment of the present invention;

[0063] Figure 4 A time series graph of a seismic wavelet provided by an embodiment of the present invention;

[0064] Figure 5 A two-dimensional migration seismic profile provided by an embodiment of the present invention;

[0065] Figure 6 A schematic structural diagram of a two-dimensional migration seismic profile generating device provided by an embodiment of the present invention;

[0066] Figure 7 A schematic structural diagram of a computer device provided in an embodiment of the present invention. DETAILED DESCRIPTION

[0067] To make the objectives, technical solutions, and advantages of the present invention more clear, the technical solutions of the present invention will be clearly and completely described below in conjunction with the accompanying drawings. Obviously, the embodiments described are only some of the embodiments of the present invention, not all of them. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without making creative efforts shall fall within the scope of protection of the present invention.

[0068] The following combination Figure 1-Figure 5 The method for generating a two-dimensional migration seismic profile according to the present invention is described.

[0069] like Figure 1 As shown, this embodiment proposes a first method for generating a two-dimensional migration seismic profile, which may include the following steps:

[0070] S101. Perform grid division in the horizontal and vertical directions on a two-dimensional vertical section of the target underground area to obtain a corresponding grid division interface; wherein the two-dimensional vertical section includes a volcanic channel.

[0071] The target underground area may be an underground area including the volcanic channel to be studied. In this embodiment, the target underground area may be selected by technicians based on the area where the volcanic channel to be studied is located.

[0072] Specifically, the two-dimensional vertical profile may be an actual or inferred section along the vertical direction of the target underground area, showing the geological structure within a certain depth of the target underground area.

[0073] Alternatively, the 2D vertical section may be a depth-domain section, in which case the horizontal axis of the 2D vertical section is the length axis and the vertical axis is the depth axis. Alternatively, the 2D vertical section may be a time-domain section, in which case the horizontal axis of the 2D vertical section is the length axis and the vertical axis is the time axis.

[0074] Specifically, this embodiment can evenly divide the two-dimensional vertical section along the horizontal and vertical directions to obtain a corresponding meshing interface. It is understood that the meshing interface can include multiple meshes, each of which has the same size. The horizontal and vertical dimensions of a single mesh can be the same or different.

[0075] S102: Determine the number of horizontal grid rows and vertical grid columns in the grid division interface.

[0076] The number of horizontal grid rows may be the total number of grid rows in the grid division interface, and the number of vertical grid columns may be the total number of grid columns in the grid division interface.

[0077] Specifically, this embodiment can count the number of horizontal grid rows and the number of vertical grid columns in the grid division interface.

[0078] S103. Construct a zero matrix of corresponding dimensions according to the number of horizontal grid rows and the number of vertical grid columns and use it as the matrix to be assigned.

[0079] Specifically, this embodiment can construct a zero matrix with the number of element rows equal to the number of horizontal grid rows and the number of element columns equal to the number of vertical grid columns based on the number of horizontal grid rows and the number of vertical grid columns. It can be understood that the value of each element in the zero matrix is ​​0.

[0080] Specifically, in this embodiment, the constructed zero matrix can be used as the matrix to be assigned.

[0081] S104. According to the location distribution of the volcanic channel, a vertical starting grid corresponding to the volcanic channel is determined on the grid division interface, and a left boundary curve and a right boundary curve corresponding to the volcanic channel are generated.

[0082] Specifically, this embodiment can determine the vertical starting grid of the volcanic channel on the grid division interface according to the position distribution of the volcanic channel in the two-dimensional vertical section.

[0083] It can be understood that the vertical starting grid can be a grid in the grid division interface that corresponds to the vertical starting position of the volcanic channel in the two-dimensional vertical section.

[0084] The left boundary curve can be a curve corresponding to the left side boundary of the volcanic conduit in the two-dimensional vertical profile, and the right boundary curve can be a curve corresponding to the right side boundary of the volcanic conduit in the two-dimensional vertical profile.

[0085] S105, along the vertical line where the vertical starting grid is located, sequentially traversing all grids between the vertical starting grid and the last grid on the vertical line.

[0086] Specifically, the embodiment can start from the vertical starting grid in the grid division interface, and traverse the grids along the vertical line where the vertical starting grid is located.

[0087] S106, accordinging to each grid traversed, assigning values to elements in the to-be-assigned matrix to obtain a matrix after assignment.

[0088] Specifically, the embodiment can assign values to elements at different positions in the to-be-assigned matrix according to each grid traversed.

[0089] Optionally, step S106 can include:

[0090] For the first grid traversed, it is determined whether the first grid reaches the bottom of the volcanic conduit.

[0091] If it is determined that the first grid does not reach the bottom of the volcanic conduit, a volcanic deformation value corresponding to the first grid is calculated, the elements in the to-be-assigned matrix are assigned using the volcanic deformation value, a new to-be-assigned matrix is obtained, and the second grid is traversed.

[0092] For the second grid traversed, it is determined whether the second grid reaches the bottom of the volcanic conduit until it is determined that the latest traversed grid reaches the bottom of the volcanic conduit, and the latest to-be-assigned matrix is taken as the matrix after assignment.

[0093] It can be understood that the embodiment can determine that the assignment is completed when it is determined that a certain grid reaches the bottom of the volcanic conduit, and determine the current latest to-be-assigned matrix as the matrix after assignment.

[0094] Optionally, the above determination of whether the first grid reaches the bottom of the volcanic conduit includes:

[0095] According to the vertical coordinates of the first grid, corresponding left boundary points and right boundary points are determined on the left boundary curve and the right boundary curve respectively.

[0096] The vertical coordinates of the left boundary points and the right boundary points are disturbed to displace the left boundary points and the right boundary points vertically respectively to obtain disturbed left boundary points and disturbed right boundary points.

[0097] Randomly create intermediate position points between the left boundary point after perturbation and the right boundary point after perturbation;

[0098] Calculate the first distance and the second distance between the middle position point and the left boundary point after disturbance and the right boundary point after disturbance respectively;

[0099] If both the first distance and the second distance are greater than 0, it is determined that the first grid has not reached the bottom of the volcanic channel;

[0100] If at least one of the first distance and the second distance is not greater than 0, it is determined that the first one reaches the bottom of the volcanic channel.

[0101] Specifically, this embodiment can perform a small-scale perturbation on the vertical coordinates of the left boundary point and the boundary point based on the random displacement function, thereby performing a small-scale random displacement on the left boundary point and the boundary point to obtain the perturbed left boundary point and the perturbed right boundary point.

[0102] It should be noted that the horizontal coordinate of the middle point is between the horizontal coordinate of the left boundary point after disturbance and the horizontal coordinate of the right boundary point after disturbance, and the vertical coordinate of the middle point is between the vertical coordinate of the left boundary point after disturbance and the vertical coordinate of the right boundary point after disturbance.

[0103] Specifically, the first distance is the distance between the middle position point and the left boundary point after the disturbance, and the second distance is the distance between the middle position point and the right boundary point after the disturbance.

[0104] Optionally, the above calculation of the volcano deformation value corresponding to the first grid includes:

[0105] Determine the projection point from the middle position point to the target line segment, where the target line segment is from the left boundary point to the right boundary point;

[0106] Constructing a corresponding first Bessel function according to the left boundary point, the projection point and the middle position point, and constructing a corresponding second Bessel function according to the right boundary point, the projection point and the middle position point;

[0107] Obtaining the horizontal coordinates of the left boundary point and the right boundary point, inputting the horizontal coordinate of the left boundary point into the first Bessel function to calculate the corresponding first function value, and inputting the horizontal coordinate of the right boundary point into the second Bessel function to calculate the corresponding second function value;

[0108] The first function value and the second function value are respectively determined as the first volcanic deformation value and the second volcanic deformation value corresponding to the left boundary point and the right boundary point, and the first volcanic deformation value and the second volcanic deformation value are taken as the volcanic deformation value corresponding to the first grid as a whole.

[0109] Specifically, the embodiment can draw a perpendicular line of the intermediate position point to the target line segment, determine an intersection point of the perpendicular line and the target line segment, and determine the intersection point as a projection point of the intermediate position point to the target line segment.

[0110] Specifically, the embodiment can construct a first Bezier curve according to the left boundary point, the projection point and the intermediate position point, and construct a second Bezier curve according to the right boundary point, the projection point and the intermediate position point.

[0111] Specifically, the first Bezier function corresponds to the first Bezier curve, and the second Bezier function corresponds to the second Bezier curve.

[0112] Optionally, the above assigning values to elements in the to-be-assigned matrix by using the volcano deformation values to obtain a new to-be-assigned matrix comprises:

[0113] In the grid division interface, a first grid and a second grid corresponding to the position coordinates of the left boundary point and the position coordinates of the right boundary point are determined respectively.

[0114] In the to-be-assigned matrix, a first element position and a second element position corresponding to the first grid and the second grid are determined respectively.

[0115] The zero values at the first element position and the second element position are replaced by the first volcano deformation value and the second volcano deformation value respectively to obtain a new to-be-assigned matrix.

[0116] Specifically, the embodiment can determine a grid containing the left boundary point as the first grid and a grid containing the right boundary point as the second grid in the grid division interface.

[0117] It should be noted that the to-be-assigned matrix is constructed by the embodiment according to the number of horizontal grid rows and the number of vertical grid columns in the grid division interface, and the number of element rows and the number of element columns in the to-be-assigned matrix are equal to the number of horizontal grid rows and the number of vertical grid columns in the grid division interface. The position of any element in the to-be-assigned matrix can correspond to the position of a certain grid in the grid division interface, and all elements in the to-be-assigned matrix are in one-to-one correspondence with all grids in the grid division interface.

[0118] Specifically, the embodiment can determine the first element position and the second element position in the to-be-assigned matrix according to the first grid and the second grid. Then, the embodiment can replace 0 at the first element position with the first volcano deformation value and replace 0 at the second element position with the second volcano deformation value.

[0119] S107, the matrix after the assignment is completed is taken as a volcano deformation matrix.

[0120] Specifically, this embodiment can stop traversal and determine that the assignment is completed when it is determined that a certain grid reaches the bottom of the volcanic channel, and determine the current latest matrix to be assigned as the matrix after the assignment is completed and used as the volcano deformation matrix.

[0121] S108. Generate a two-dimensional migration seismic profile corresponding to the two-dimensional vertical profile according to the volcanic deformation matrix.

[0122] Specifically, this embodiment can generate a two-dimensional remote seismic profile corresponding to the two-dimensional vertical profile based on the volcanic deformation matrix and the migration seismic data of the target underground area.

[0123] The 2D migration seismic profile generation method proposed in this embodiment can select a target underground region containing a volcanic conduit based on the volcanic conduit, perform numerical simulation based on the 2D vertical profile of the target underground region, generate a volcanic deformation matrix, and then generate a 2D migration seismic profile corresponding to the 2D vertical profile based on the volcanic deformation matrix. This embodiment can effectively perform numerical simulation of 2D migration seismic profiles containing a volcanic conduit and improve the accuracy of the generated 2D migration seismic profiles containing the volcanic conduit.

[0124] based on Figure 1 This embodiment proposes a second method for generating a two-dimensional migration seismic profile. In this method, step S108 may include S1081, S1082, S1083, S1084, S1085, and S1086.

[0125] S1081. Acquire location data of multiple key layers in the target underground area and elastic parameter data of the underground medium.

[0126] Optionally, the position data of each key layer includes the position coordinates of multiple layer points, and the elastic parameter data includes the elastic parameter value corresponding to each key layer.

[0127] Specifically, the position data of each key layer includes the position coordinates of multiple layer points.

[0128] The elastic parameter data of the underground medium may include the elastic parameter value of each key layer.

[0129] Optionally, elastic parameter values ​​may include longitudinal wave velocity and density.

[0130] Specifically, this embodiment determines the key layer information based on the actual offset seismic data of the target underground area, and sets the number of key layers to NL. Each layer i, i∈NL information should contain at least the horizontal position coordinates and vertical position coordinates of two layer points to form a layer dictionary Z(i), i∈NL. It should be noted that two layer points determine a layer straight line. For layer lines of other shapes, multiple layer point controls can be added in this step, or a straight layer can be formed first, and the layer morphology can be enriched in the subsequent steps of adding global and local folds.

[0131] Specifically, this embodiment can obtain the one-dimensional elastic parameter information of the underground medium in the target underground area to form the longitudinal wave velocity dictionary k p (i),i∈NL and density dictionary ρ(i),i∈NL.

[0132] S1082. Generate a reflectivity matrix corresponding to a two-dimensional vertical profile based on the position data, elastic parameter data, number of horizontal grid rows, and number of vertical grid columns of each key layer.

[0133] Optionally, step S1082 may include:

[0134] According to the position coordinates of the layer sites and the number of horizontal grid rows in the position data of each key layer, the layer site array of each key layer is expanded to generate the layer site array. The number of elements in each layer site array is equal to the number of horizontal grid rows.

[0135] Expanding the elastic parameter value according to the elastic parameter data and the number of vertical grid columns to obtain an elastic parameter array, wherein the number of elements in the elastic parameter array is equal to the number of vertical grid columns;

[0136] Introducing global wrinkle features into each layer array to obtain corresponding multiple global wrinkle layer arrays;

[0137] Generate a two-dimensional horizon matrix based on multiple global fold horizon arrays;

[0138] Generate a two-dimensional elastic parameter matrix according to the two-dimensional layer matrix and the elastic parameter array;

[0139] Generates a reflectivity matrix based on a 2D elastic parameter matrix.

[0140] Specifically, in the process of expanding the layer site to obtain the layer site array of each key layer, this embodiment can perform linear spline fitting on the position coordinates of all layer sites in each layer dictionary Z(i). In the i-th layer, there are K i layer sites, each layer site k, k∈K i Contains the horizontal position x k and vertical position z k , each layer fits a polynomial Si (x). Where:

[0141]

[0142] Where x is the horizontal coordinate. The polynomial after fitting each layer is interpolated separately, that is, when x is taken from [1, NX] in ascending order, the corresponding S is obtained. i The vertical position of the layer (x) is calculated and placed into the layer array in sequence, converting NL unequal-length layer dictionaries Z(i), i∈NL into NL equal-length layer arrays Z(x), x∈NX. After interpolation, the layer array length NX is consistent with the number of horizontal grids after the actual data is segmented, ensuring the scale and accuracy of the simulation, making the simulation results more accurate and reliable.

[0143] Specifically, in the process of expanding the elastic parameter values ​​to obtain the elastic parameter array, the elastic parameter dictionary obtained, i.e., the longitudinal wave velocity dictionary, may be used. p (i),i∈NL and density dictionary ρ(i),i∈NL. According to the vertical position coordinates of the key layer array Z(x), a random elastic parameter increment is added to the random position to supplement the intermediate elastic parameter values ​​between the key layers. The random elastic parameter increment is implemented as follows:

[0144] A(z)=(U(0,1)-0.5)*s*V0(z).

[0145] Among them, U(0,1) represents the random number randomly selected in the uniform distribution interval [0,1], s is the elastic parameter increment coefficient, which is generally set to 0.6 and can be determined according to the waveform distribution of seismic data in the actual work area. The value range of the independent variable z is [1,NZ]. After random increment processing at random positions, a one-dimensional longitudinal wave velocity array V is formed. p (z),z∈NZ, and the density array ρ(z),z∈NZ, which contains the corresponding P-wave velocity and density values ​​and is arranged in ascending order according to z. The one-dimensional elastic parameter array, after random position processing, remains consistent with the actual data on the vertical scale, ensuring the scale of the simulation. At the same time, the random increments introduced between layers simulate the uneven variations in actual subsurface lithology, enhancing the elastic parameters' ability to characterize complex subsurface structures. In the subsequent wavefield convolution simulation, the random size variations of the elastic parameters are equivalent to adding scattering points and reasonable background noise to the simulation process, making the simulation results more consistent with actual conditions.

[0146] Specifically, in this embodiment, a global fold feature is introduced into each layer array. In the process of obtaining corresponding multiple global fold layer arrays, a global stratum fold deformation can be randomly generated and introduced into the layer array based on a three-point Bessel function. The deformation equation is:

[0147] B(x) = (1+x) 2 P0+2(1-x)xP1+x 2 P2.

[0148] wherein x is an independent variable, the number of values of x corresponds to the range of wrinkle deformation, for the global lateral wrinkle deformation, the number of x is NX. p0 is the starting point coordinate corresponding to the Bezier base function, p1 is the control point coordinate corresponding to the Bezier base function, and P2 is the terminal point coordinate corresponding to the Bezier base function.

[0149] After the horizon array is transformed by the random Bezier function, the horizon array Z g (x) introducing global wrinkle is obtained, x∈NX. The horizon array Z g (x) introducing global wrinkle is:

[0150] Z g (x) = Z(x) + B(x).

[0151] Subsequently, the two-dimensional profile interpolation is performed on the NL horizon arrays Z g (x), x∈NX by using the gradient linear interpolation method, so as to realize the continuity of the vertical gradient change of the horizon, thereby avoiding the step-type abnormality of the horizon.

[0152] Specifically, in the process of generating a two-dimensional horizon matrix according to a plurality of global wrinkle horizon arrays, a midpoint can be selected as a reference position Z g (x0) in a single global wrinkle horizon array, the difference Z g (x)-Z g (x0) between the horizontal horizon and the reference position is calculated. The reference position Z g (x0) of each key horizon is taken as an independent variable, and the difference is taken as a dependent variable, and a polynomial fitting is performed to obtain a target polynomial:

[0153] f x (t) = a x +b x t+c x t 2 +d x .

[0154] wherein t is an independent variable, f x (t) is a dependent variable, a x ,b x ,c x ,d x are coefficients of the polynomial, and the coefficients are obtained by bringing the values Z g (x0) of the NL independent variables into t and the dependent variable Z g(x)-Z g (x0) Substitute f x (t) is obtained. After fitting the polynomial, the vertical position coordinate z, z∈NZ is substituted into the polynomial to obtain the layer gradient Z g ′ (x,z),x∈NX,z∈NZ. Each point in the layer gradient matrix is ​​added to its corresponding depth value to obtain a two-dimensional layer matrix Z g (x,z),x∈NX,z∈NZ.

[0155] Afterwards, in the process of generating the two-dimensional elastic parameter matrix, the two-dimensional layer matrix Z g (x,z)The midpoint of the horizontal axis Z g (x0,z) as the independent variable, interpolation Z g (x,z)-Z g (x0,z) is used as the dependent variable to fit a new polynomial, and then the longitudinal wave velocity array V p (z),z∈NZ and the density array ρ(z),x∈NZ are substituted into the new polynomial to form the longitudinal wave velocity matrix V p (x,z),x∈NX,z∈NZ and density matrix ρ(x,z),x∈NX,z∈NZ.

[0156] Specifically, in the process of generating a reflectivity matrix based on a two-dimensional elastic parameter matrix, according to the Zoeppritz equation, the simplified model is the reflection coefficient of vertically incident longitudinal waves, and the two-dimensional reflectivity matrix R0(x,z) can be calculated from the two-dimensional elastic parameter matrix, x∈NX, z∈NZ.

[0157]

[0158] The vertical axis subscript z ranges from [1, NZ-1], the horizontal axis subscript x ranges from [1, NX], and the reflectivity coefficient of the grid points not involved in the calculation is 0.

[0159] S1083. Introduce local fold characteristics and fault characteristics into the reflectivity matrix to obtain the fault reflectivity matrix.

[0160] Specifically, this embodiment can first introduce local wrinkle features into the reflectivity matrix. This embodiment can add local small-scale wrinkles into the reflectivity matrix to obtain an updated reflectivity matrix R1(x,z), x∈NX, z∈NZ. For the introduction of local wrinkles, the deformation function is the same as the above deformation equation. On this basis, it is also necessary to limit the scope of the local wrinkles. Assuming that the boundary of a local wrinkle is fully represented by the left boundary and right boundary functions, for a wedge-shaped deformation, its boundary representation function is:

[0161]

[0162] Among them, z is the vertical coordinate, the range is [1, NZ], a l ,b l , and a r ,b r They are the left boundary parameter and the bounded parameter respectively.

[0163] Specifically, this embodiment can introduce fault features into the reflectivity matrix R1(x,z), x∈NX, z∈NZ, to form a reflectivity model that includes fault structures, namely the fault reflectivity matrix R2(x,z), x∈NX, z∈NZ. This embodiment can introduce fault features by executing Algorithm 1. The input of Algorithm 1 includes:

[0164] reflectivity: A two-dimensional array representing the reflectivity matrix of the formation.

[0165] fault_position: One-dimensional array, representing the coordinates [z,x] of the fault starting position.

[0166] fault_angle: A numeric value indicating the inclination angle of the fault.

[0167] fault_direction: An integer value indicating the direction of fault displacement (positive values ​​indicate upward displacement, negative values ​​indicate downward displacement).

[0168] fault_displacement: A numeric value representing the displacement of the fault.

[0169] Among them, the output of Algorithm 1 is:

[0170] new_reflectivity: A two-dimensional array representing the new reflectivity matrix after considering the impact of the fault.

[0171] Specifically, the execution process of Algorithm 1 includes:

[0172] Get the height NZ and width NX of the input reflectivity matrix reflectivity.

[0173] Create a new reflectivity matrix new_reflectivity that is a copy of reflectivity.

[0174] Calculate the vertical displacement dz of the fault using fault_displacement and fault_angle.

[0175] For each x from fault_position[1] to NX:

[0176] a. Calculate the new vertical position of the fault at its current x position, new_fault_pos_z.

[0177] b. If new_fault_pos_z is in the valid vertical range [0, NZ]:

[0178] For each z from new_fault_pos_z to NZ:

[0179] Calculate the new vertical position new_z, taking into account the direction of the fault displacement.

[0180] If new_z is in the valid vertical range [0,NZ-1], copy the value of reflectivity[z,x] to new_reflectivity[new_z,x].

[0181] Returns the new reflectivity matrix new_reflectivity.

[0182] The cross-sectional schematic diagram corresponding to the fault reflectivity matrix obtained in this embodiment can be shown as follows Figure 2 shown.

[0183] Specifically, after generating the volcanic deformation matrix, this embodiment can use the coordinate mapping method to apply the volcanic deformation matrix to the reflection coefficient R2(x,z),x∈NX,z∈NZ matrix to form the reflectivity matrix R3(x,z),x∈NX,z∈NZ containing the volcanic channel structure.

[0184] In this embodiment, the volcano deformation matrix can be generated based on Algorithm 2. The input of Algorithm 2 is:

[0185] indata: 2D reflectivity matrix.

[0186] The output of Algorithm 2 is:

[0187] oudata: two-dimensional deformation matrix Q.

[0188] The execution process of Algorithm 2 includes:

[0189] Step 1: Initialize constants.

[0190] NZ and NX: Get the vertical and horizontal grid numbers of the two-dimensional reflectivity matrix.

[0191] z0: The vertical starting position coordinate of the volcanic channel.

[0192] oudata: deformation matrix used to represent the deformation of the strata in the volcanic channel area.

[0193] Step 2: Generate the volcanic channel side boundary curve.

[0194] left_bond and right_bond: Calculate the left and right boundary curves of the volcanic channel to be simulated, indicating the left and right boundary positions of the volcanic channel on the two-dimensional section.

[0195] left_shift and right_shift: Generate random displacement functions, which are used to subsequently generate small-scale perturbations to the vertical coordinates of the boundary points P0 and P2.

[0196] Step 3: Generate internal stratigraphic control points of the volcanic channel.

[0197] Traverse the vertical grid points [z0,NZ]:

[0198] Calculate the position coordinates P0 and P2 of the left and right boundaries left_bond and right_bond under the current grid;

[0199] Add random shifts left_shift and right_shift to the vertical coordinates of P0 and P2 to obtain P0′ and P2′.

[0200] A random middle position P1 is created as a reference point between the left and right boundaries.

[0201] Calculate the displacements X0 and X2 between the middle position P1 and the two boundary positions, and determine whether the bottom of the volcanic channel has been reached based on the size of X0 and X2.

[0202] Step 4: Generate the internal stratigraphic control curve of the volcanic channel.

[0203] Traverse the vertical grid points [z0,NZ]:

[0204] Determine whether the displacements X0 and X2 between the left and right boundaries are both greater than zero:

[0205] If so, calculate two Bezier curves. Take the projection of P1 on the line connecting P0 and P2 as the middle point P k , respectively (P0′,P k ,P1) and (P1,P k ,P2′) are the endpoints to calculate two Bessel functions, which are used to describe the displacement of the stratum from the left boundary to the middle position and from the middle position to the right boundary.

[0206] At the current vertical position, within the horizontal positions corresponding to P0 and P2, the values ​​of the two Bezier curves are filled into the matrix to be assigned.

[0207] Step 5. Return the generated volcano deformation matrix.

[0208] The cross-sectional diagram corresponding to the volcano deformation matrix generated in this embodiment can be shown as follows Figure 3shown.

[0209] S1084. Perform coordinate mapping on the fault reflectivity matrix based on the volcanic deformation matrix to obtain the volcanic reflectivity matrix.

[0210] Specifically, after generating the volcanic deformation matrix Q(x,z), this embodiment can obtain the fault reflectivity matrix R2(x,z), x∈NX,z∈NZ, and the position coordinate matrix Q0(x,z), x∈NX,z∈NZ. Then, the new coordinate matrix Q1(x,z) with the volcanic channel is:

[0211] Q1(x,z)=Q0(x,z)+Q(x,z).

[0212] Next, in this embodiment, Q0(x,z) can be substituted into the target polynomial as the independent variable and Q1(x,z) as the dependent variable to perform polynomial fitting. Substituting R2(x,z) into the fitted polynomial yields the reflectivity matrix for the volcanic channel, i.e., the volcanic reflectivity matrix.

[0213] S1085. Create a seismic wavelet function corresponding to the two-dimensional vertical profile.

[0214] S1086. Convolve the volcanic reflectivity matrix with the seismic wavelet function to obtain a two-dimensional migration seismic profile.

[0215] Specifically, in this embodiment, the seismic wavelet function w(z) can be set, and the time series curve of the seismic wavelet is as follows: Figure 4 As shown. Then, the seismic wavelet can be convolved with the volcanic reflectivity matrix R3(x,z),x∈NX,z∈NZ to obtain the final simulated migration seismic data, that is, the two-dimensional migration seismic profile D(x,z),x∈NX,z∈NZ. Where:

[0216] D(x,z)=w(z)*R3(x,z).

[0217] The convolution is performed along the vertical axis (depth or time axis) to obtain a two-dimensional migration seismic profile, such as Figure 5 shown.

[0218] It should be noted that in seismic data from related technologies, volcanic conduits have relatively clear identification features. Relatively narrow volcanic conduit phases exhibit chaotic high-frequency signals, are at a certain angle to the continuity axis of the surrounding strata, and have a relatively clear upward flow trend. Therefore, related technologies can identify volcanic conduits in migrated seismic data through feature recognition methods. However, volcanic conduit phases run through areas with extensive volcanic rock development. Seismic attribute methods have difficulty directly distinguishing volcanic rocks of different microfacies, and manual interpretation of volcanic conduits is labor-intensive. Supervised feature recognition methods require a large amount of volcanic conduit phase data for feature learning. Therefore, obtaining migrated seismic data containing volcanic conduits has become an important research issue.

[0219] The two-dimensional migration seismic profile generation method proposed in this embodiment can achieve characteristic simulation of the internal reflection coefficient of the volcanic channel in the migration seismic data. Through the random design of the stratigraphic structure, faults and volcanic channels, a large amount of complex two-dimensional migration seismic data containing volcanic channels can be generated.

[0220] like Figure 6 As shown, this embodiment provides a two-dimensional migration seismic profile generating device, which may include:

[0221] A division unit 601 is configured to perform grid division in the horizontal and vertical directions on a two-dimensional vertical section of the target underground area to obtain a corresponding grid division interface; wherein the two-dimensional vertical section includes a volcanic channel;

[0222] A determining unit 602 is used to determine the number of horizontal grid rows and vertical grid columns in the grid division interface;

[0223] A construction unit 603 is configured to construct a zero matrix of corresponding dimensions according to the number of horizontal grid rows and the number of vertical grid columns and use the zero matrix as the matrix to be assigned;

[0224] The processing unit 604 is used to determine the vertical starting grid corresponding to the volcanic channel on the grid division interface according to the location distribution of the volcanic channel, and generate the left boundary curve and the right boundary curve corresponding to the volcanic channel;

[0225] A traversal unit 605 is configured to traverse all grids from the vertical starting grid to the last grid on the vertical line in sequence along the vertical line where the vertical starting grid is located;

[0226] An assignment unit 606 is configured to assign values ​​to elements in the to-be-assigned matrix according to each traversed grid, to obtain a matrix after the assignment is completed;

[0227] As unit 607, used to use the matrix after the assignment as the volcano deformation matrix;

[0228] The generating unit 608 is configured to generate a two-dimensional migration seismic profile corresponding to the two-dimensional vertical profile according to the volcanic deformation matrix.

[0229] It should be noted that the processing procedures of the division unit 601, the determination unit 602, the construction unit 603, the processing unit 604, the traversal unit 605, the assignment unit 606, the as unit 607 and the generation unit 608 and the beneficial effects thereof can be referred to in detail. Figure 1 Steps S101 to S108 in the above are not described in detail.

[0230] Optionally, the assignment unit 606 is further configured to:

[0231] For the first grid traversed, determine whether the first grid reaches the bottom of the volcanic channel;

[0232] If it is determined that the first grid has not reached the bottom of the volcanic channel, the volcanic deformation value corresponding to the first grid is calculated, and the volcanic deformation value is used to assign the elements in the to-be-assigned matrix to obtain a new to-be-assigned matrix, and the second grid is traversed;

[0233] For the second grid traversed, determine whether the second grid reaches the bottom of the volcanic channel, until it is determined that the latest traversed grid reaches the bottom of the volcanic channel, and use the latest matrix to be assigned as the matrix after the assignment is completed.

[0234] Optionally, the assignment unit 606 is further configured to:

[0235] According to the vertical coordinates of the first grid, the corresponding left boundary point and right boundary point are determined on the left boundary curve and the right boundary curve respectively;

[0236] By perturbing the vertical coordinates of the left boundary point and the right boundary point, the left boundary point and the right boundary point are vertically displaced respectively to obtain the disturbed left boundary point and the disturbed right boundary point;

[0237] Randomly create intermediate position points between the left boundary point after perturbation and the right boundary point after perturbation;

[0238] Calculate the first distance and the second distance between the middle position point and the left boundary point after disturbance and the right boundary point after disturbance respectively;

[0239] If both the first distance and the second distance are greater than 0, it is determined that the first grid has not reached the bottom of the volcanic channel;

[0240] If at least one of the first distance and the second distance is not greater than 0, it is determined that the first one reaches the bottom of the volcanic channel.

[0241] Optionally, the assignment unit 606 is further configured to:

[0242] determining a projection point of the intermediate position point to a target line segment, the target line segment being from the left boundary point to the right boundary point;

[0243] constructing a first Bezier function according to the left boundary point, the projection point and the intermediate position point, and constructing a second Bezier function according to the right boundary point, the projection point and the intermediate position point;

[0244] obtaining horizontal coordinates of the left boundary point and the right boundary point, inputting the horizontal coordinate of the left boundary point into the first Bezier function to calculate a corresponding first function value, and inputting the horizontal coordinate of the right boundary point into the second Bezier function to calculate a corresponding second function value;

[0245] respectively determining the first function value and the second function value as first and second volcano deformation values corresponding to the left boundary point and the right boundary point, and taking the first and second volcano deformation values as a whole as volcano deformation values corresponding to a first grid.

[0246] Optionally, the assignment unit 606 is further configured to:

[0247] respectively determining a first grid and a second grid corresponding to the position coordinate of the left boundary point and the position coordinate of the right boundary point in the grid division interface;

[0248] respectively determining a first element position and a second element position corresponding to the first grid and the second grid in the to-be-assigned matrix;

[0249] respectively replacing zero values at the first element position and the second element position with the first and second volcano deformation values to obtain a new to-be-assigned matrix.

[0250] Optionally, the generation unit 608 is further configured to:

[0251] obtaining position data of a plurality of key horizons and elastic parameter data of a subsurface medium in a target subsurface region;

[0252] generating a reflectivity matrix corresponding to a two-dimensional vertical profile according to the position data of each key horizon, the elastic parameter data, a horizontal grid row number and a vertical grid column number;

[0253] introducing local fold features and fault features into the reflectivity matrix to obtain a fault reflectivity matrix;

[0254] performing coordinate mapping on the fault reflectivity matrix based on the volcano deformation matrix to obtain a volcano reflectivity matrix;

[0255] creating a seismic wavelet function corresponding to the two-dimensional vertical profile, and convolving the volcano reflectivity matrix with the seismic wavelet function to obtain a two-dimensional migration seismic profile.

[0256] Optionally, the position data of each key layer includes the position coordinates of multiple layer points, and the elastic parameter data includes the elastic parameter value corresponding to each key layer;

[0257] The generating unit 608 is further configured to:

[0258] According to the position coordinates of the layer sites and the number of horizontal grid rows in the position data of each key layer, the layer site array of each key layer is expanded to generate the layer site array. The number of elements in each layer site array is equal to the number of horizontal grid rows.

[0259] Expanding the elastic parameter value according to the elastic parameter data and the number of vertical grid columns to obtain an elastic parameter array, wherein the number of elements in the elastic parameter array is equal to the number of vertical grid columns;

[0260] Introducing global wrinkle features into each layer array to obtain corresponding multiple global wrinkle layer arrays;

[0261] Generate a two-dimensional horizon matrix based on multiple global fold horizon arrays;

[0262] Generate a two-dimensional elastic parameter matrix according to the two-dimensional layer matrix and the elastic parameter array;

[0263] Generates a reflectivity matrix based on a 2D elastic parameter matrix.

[0264] The 2D migration seismic profile generation device proposed in this embodiment can select a target underground region containing a volcanic conduit based on the volcanic conduit, perform numerical simulation based on the 2D vertical profile of the target underground region, generate a volcanic deformation matrix, and then generate a 2D migration seismic profile corresponding to the 2D vertical profile based on the volcanic deformation matrix. This embodiment can effectively perform numerical simulation of 2D migration seismic profiles containing a volcanic conduit and improve the accuracy of the generated 2D migration seismic profiles containing the volcanic conduit.

[0265] The two-dimensional offset seismic profile generation device in this embodiment is presented in the form of a functional unit, where the unit refers to an ASIC (Application Specific Integrated Circuit) circuit, a processor and memory that executes one or more software or fixed programs, and / or other devices that can provide the above functions.

[0266] The embodiment of the present invention also provides a computer device having the above Figure 6 The two-dimensional migration seismic profile generating device shown.

[0267] See also Figure 7, a structural diagram of a computer device provided by an optional embodiment of the present invention, the computer device includes: one or more processors 10, a memory 20, and interfaces for connecting various components, including high-speed interfaces and low-speed interfaces. The various components are connected to each other using different buses for communication, and can be installed on a common motherboard or installed in other ways as needed. The processor can process instructions executed in the computer device, including instructions stored in or on the memory to display graphical information of the GUI on an external input / output device (such as a display device coupled to the interface). In some optional embodiments, if necessary, multiple processors and / or multiple buses can be used together with multiple memories and multiple memories. Similarly, multiple computer devices can be connected, and each device provides part of the necessary operations (for example, as a server array, a group of blade servers, or a multi-processor system). Figure 7 A processor 10 is taken as an example.

[0268] The processor 10 may be a central processing unit, a grid processor, or a combination thereof. The processor 10 may further include a hardware chip. The hardware chip may be an application-specific integrated circuit, a programmable logic device, or a combination thereof. The programmable logic device may be a complex programmable logic device, a field programmable gate array, a general purpose array logic, or any combination thereof.

[0269] The memory 20 stores instructions that can be executed by at least one processor 10, so as to enable at least one processor 10 to execute the method shown in the above embodiment.

[0270] The memory 20 may include a program storage area and a data storage area, wherein the program storage area may store an operating system and application programs required for at least one function. The data storage area may store data created based on the use of the computer device, etc. In addition, the memory 20 may include a high-speed random access memory and may also include a non-transient memory, such as at least one disk storage device, a flash memory device, or other non-transient solid-state storage device. In some optional embodiments, the memory 20 may optionally include a memory remotely located relative to the processor 10, and these remote memories may be connected to the computer device via a grid. Examples of the above-mentioned grid include, but are not limited to, the Internet, an intranet, a local area network, a mobile communication network, and combinations thereof.

[0271] The memory 20 may include volatile memory, such as random access memory. The memory may also include non-volatile memory, such as flash memory, a hard disk, or a solid-state drive. The memory 20 may also include a combination of the above types of memory.

[0272] The computer device further comprises a communication interface 30 for communicating with other devices or a communication grid.

[0273] The embodiment of the present invention also provides a computer-readable storage medium. The above-mentioned method according to the embodiment of the present invention can be implemented in hardware, firmware, or implemented as a computer code that can be recorded in a storage medium, or implemented as a computer code that is originally stored in a remote storage medium or a non-temporary machine-readable storage medium and will be stored in a local storage medium through grid download, so that the method described herein can be stored in such software processing on a storage medium using a general-purpose computer, a dedicated processor, or programmable or dedicated hardware. Among them, the storage medium can be a magnetic disk, an optical disk, a read-only storage memory, a random access memory, a flash memory, a hard disk or a solid-state drive, etc.; further, the storage medium can also include a combination of the above-mentioned types of memory. It can be understood that a computer, a processor, a microprocessor controller or programmable hardware includes a storage component that can store or receive software or computer code. When the software or computer code is accessed and executed by a computer, a processor or hardware, the method shown in the above embodiment is implemented.

[0274] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention, rather than to limit it. Although the present invention has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that they can still modify the technical solutions described in the aforementioned embodiments, or make equivalent replacements for some of the technical features therein. However, these modifications or replacements do not deviate the essence of the corresponding technical solutions from the spirit and scope of the technical solutions of the various embodiments of the present invention.

Claims

1. A method for generating a two-dimensional migration seismic profile, characterized in that: include: Meshing the target underground area in a two-dimensional vertical section in the horizontal and vertical directions to obtain a corresponding meshing interface; wherein the two-dimensional vertical section includes a volcanic channel; Determining the number of horizontal grid rows and the number of vertical grid columns in the grid division interface, and constructing a zero matrix of corresponding dimensions according to the number of horizontal grid rows and the number of vertical grid columns as a matrix to be assigned; According to the location distribution of the volcanic channel, determining a vertical starting grid corresponding to the volcanic channel on the grid division interface, and generating a left boundary curve and a right boundary curve corresponding to the volcanic channel; Traversing all the grids from the vertical starting grid to the last grid on the vertical line in sequence along the vertical line where the vertical starting grid is located; Assigning values ​​to the elements in the matrix to be assigned according to each of the traversed grids, obtaining a matrix after the assignment is completed and serving as a volcano deformation matrix; A two-dimensional offset seismic profile corresponding to the two-dimensional vertical profile is generated according to the volcanic deformation matrix.

2. The method according to claim 1, characterized in that The step of assigning values ​​to the elements in the matrix to be assigned according to each of the traversed grids to obtain a matrix after the assignment is completed includes: For the first grid traversed, determining whether the first grid reaches the bottom of the volcanic channel; If it is determined that the first grid does not reach the bottom of the volcanic channel, then calculating the volcanic deformation value corresponding to the first grid, using the volcanic deformation value to assign values ​​to the elements in the matrix to be assigned, obtaining a new matrix to be assigned, and traversing the second grid; For the second grid traversed, determine whether the second grid reaches the bottom of the volcanic channel, until it is determined that the latest traversed grid reaches the bottom of the volcanic channel, and use the latest matrix to be assigned as the matrix after the assignment is completed.

3. The method according to claim 2, characterized in that Determining whether the first grid reaches the bottom of the volcanic channel includes: Determining corresponding left boundary points and right boundary points on the left boundary curve and the right boundary curve respectively according to the vertical coordinates of the first grid; Performing vertical displacement on the left boundary point and the right boundary point respectively by perturbing the vertical coordinates of the left boundary point and the right boundary point to obtain a perturbed left boundary point and a perturbed right boundary point; Randomly create an intermediate position point between the left boundary point after the disturbance and the right boundary point after the disturbance; Calculating a first distance and a second distance between the intermediate position point and the left boundary point after disturbance and the right boundary point after disturbance, respectively; If both the first distance and the second distance are greater than 0, it is determined that the first grid has not reached the bottom of the volcanic channel; If at least one of the first distance and the second distance is not greater than 0, it is determined that the first one has reached the bottom of the volcanic channel.

4. The method according to claim 3, characterized in that Calculating the volcano deformation value corresponding to the first grid includes: Determine a projection point from the middle position point to a target line segment, where the target line segment is from the left boundary point to the right boundary point; Constructing a corresponding first Bessel function according to the left boundary point, the projection point, and the middle position point, and constructing a corresponding second Bessel function according to the right boundary point, the projection point, and the middle position point; Obtaining horizontal coordinates of the left boundary point and the right boundary point, inputting the horizontal coordinate of the left boundary point into the first Bessel function to calculate a corresponding first function value, and inputting the horizontal coordinate of the right boundary point into the second Bessel function to calculate a corresponding second function value; The first function value and the second function value are respectively determined as the first volcanic deformation value and the second volcanic deformation value corresponding to the left boundary point and the right boundary point, and the first volcanic deformation value and the second volcanic deformation value are taken as the volcanic deformation value corresponding to the first grid as a whole.

5. The method according to claim 4, characterized in that The step of assigning values ​​to the elements in the matrix to be assigned using the volcano deformation value to obtain a new matrix to be assigned includes: In the grid division interface, determining a first grid and a second grid corresponding to the position coordinates of the left boundary point and the position coordinates of the right boundary point respectively; In the matrix to be assigned, determining a first element position and a second element position corresponding to the first grid and the second grid respectively; The zero values ​​at the first element position and the second element position are replaced by the first volcanic deformation value and the second volcanic deformation value respectively to obtain a new matrix to be assigned.

6. The method according to claim 1, characterized in that Generating a two-dimensional migration seismic profile corresponding to the two-dimensional vertical profile according to the volcanic deformation matrix includes: Acquiring location data of a plurality of key layers in the target underground area and elastic parameter data of underground media; Generate a reflectivity matrix corresponding to the two-dimensional vertical profile according to the position data of each key layer, the elastic parameter data, the number of horizontal grid rows and the number of vertical grid columns; Introducing local fold characteristics and fault characteristics into the reflectivity matrix to obtain a fault reflectivity matrix; Performing coordinate mapping on the fault reflectivity matrix based on the volcanic deformation matrix to obtain a volcanic reflectivity matrix; A seismic wavelet function corresponding to the two-dimensional vertical profile is created, and the volcanic reflectivity matrix is ​​convolved with the seismic wavelet function to obtain the two-dimensional migrated seismic profile.

7. The method according to claim 6, characterized in that The position data of each key layer includes the position coordinates of a plurality of layer points, and the elastic parameter data includes the elastic parameter value corresponding to each key layer; Generating a reflectivity matrix corresponding to the two-dimensional vertical profile according to the position data of each key layer, the elastic parameter data, the number of horizontal grid rows, and the number of vertical grid columns includes: Expanding the layer site according to the position coordinates of the layer site in the position data of each key layer site and the number of horizontal grid rows to generate a layer site array for each key layer site, wherein the number of elements in each layer site array is equal to the number of horizontal grid rows; Expanding the elastic parameter value according to the elastic parameter data and the number of vertical grid columns to obtain an elastic parameter array, wherein the number of elements in the elastic parameter array is equal to the number of vertical grid columns; Introducing a global wrinkle feature into each of the layer arrays to obtain a corresponding plurality of global wrinkle layer arrays; generating a two-dimensional horizon matrix according to the plurality of global fold horizon arrays; generating a two-dimensional elastic parameter matrix according to the two-dimensional layer matrix and the elastic parameter array; The reflectivity matrix is ​​generated based on the two-dimensional elastic parameter matrix.

8. A two-dimensional migration seismic profile generating device, characterized in that: include: A partitioning unit is used to perform grid division in the horizontal direction and the vertical direction on a two-dimensional vertical section of the target underground area to obtain a corresponding grid division interface; wherein the two-dimensional vertical section includes a volcanic channel; A determination unit, configured to determine the number of horizontal grid rows and vertical grid columns in the grid division interface; A construction unit, configured to construct a zero matrix of corresponding dimensions according to the number of horizontal grid rows and the number of vertical grid columns and use the zero matrix as a matrix to be assigned; a processing unit, configured to determine, on the grid division interface, a vertical starting grid corresponding to the volcanic channel according to the position distribution of the volcanic channel, and to generate a left boundary curve and a right boundary curve corresponding to the volcanic channel; a traversal unit, configured to traverse all the grids from the vertical starting grid to the last grid on the vertical line in sequence along the vertical line where the vertical starting grid is located; An assignment unit, configured to assign values ​​to elements in the matrix to be assigned according to each of the traversed grids, to obtain a matrix after the assignment is completed; As a unit, it is used to use the matrix after the assignment as the volcano deformation matrix; A generating unit is used to generate a two-dimensional offset seismic profile corresponding to the two-dimensional vertical profile according to the volcanic deformation matrix.

9. A computer device, characterized in that: include: 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 two-dimensional migration seismic profile generation method according to any one of claims 1 to 7 by executing the computer instructions.

10. A computer-readable storage medium, characterized in that The computer-readable storage medium stores computer instructions, and the computer instructions are used to enable a computer to execute the two-dimensional migration seismic profile generation method according to any one of claims 1 to 7.

Images