Two-dimensional migration seismic section generation method, device, equipment and medium

Abstract

The invention relates to the field of oil and gas exploration seismic data processing, and discloses a two-dimensional offset seismic section generation method, device, equipment and medium, which can select a target underground area comprising a volcano channel according to the volcano channel, perform numerical simulation based on a two-dimensional vertical section of the target underground area, generate a volcano deformation matrix, and generate a seismic data matrix. And then a two-dimensional offset seismic section corresponding to the two-dimensional vertical section is generated based on the volcano deformation matrix. According to the method, the numerical simulation of the two-dimensional offset seismic section containing the volcano channel can be effectively realized, and the accuracy of the generated two-dimensional offset seismic section containing the volcano channel can be improved.

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 after migration processing of seismic data obtained through two-dimensional seismic exploration technology.

[0003] Two-dimensional migration seismic profiles can be widely used in various geological exploration scenarios, especially oil and gas exploration scenarios. Exploration workers can use two-dimensional migration seismic profiles to more accurately understand and interpret underground structures, identify stratum interfaces and structural features, and optimize oil and gas production 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 offset seismic profile generation method, device, equipment and medium, which are used to solve the defect of low accuracy of two-dimensional offset seismic profiles containing volcanic channels generated in related technologies, and improve the accuracy of two-dimensional offset 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] On a two-dimensional vertical section of the target underground area, meshing is performed in the horizontal direction and the vertical direction to obtain a corresponding meshing interface; wherein the two-dimensional vertical section includes a volcanic channel;

[0008] Determine the number of horizontal grid rows and the number of vertical grid columns in the grid division interface, and construct 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 position distribution of the volcanic channel, determining the vertical starting grid corresponding to the volcanic channel on the grid division interface, and generating the left boundary curve and the right boundary curve corresponding to the volcanic channel;

[0010] Along the vertical line where the vertical starting grid is located, all the grids between the vertical starting grid and the last grid on the vertical line are traversed in sequence;

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

[0012] A two-dimensional offset seismic profile corresponding to the two-dimensional vertical profile is generated according to the volcanic deformation matrix.

[0013] Optionally, assigning values ​​to elements in the matrix to be assigned values ​​according to each of the traversed grids to obtain a matrix after the assignment is completed includes:

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

[0015] 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, the volcanic deformation value is used to assign values ​​to the elements in the matrix to be assigned, a new matrix to be assigned is obtained, and the second grid is traversed;

[0016] For the second traversed grid, 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.

[0017] Optionally, determining whether the first grid reaches the bottom of the volcanic passage includes:

[0018] Determine 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;

[0019] 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 respectively vertically displaced to obtain a perturbed left boundary point and a perturbed right boundary point;

[0020] Randomly create an intermediate position point between the left boundary point after the disturbance and the right boundary point after the disturbance;

[0021] Calculating a first distance and a second distance between the middle position point and the left boundary point after disturbance and the right boundary point after disturbance respectively;

[0022] 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;

[0023] 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 passage.

[0024] Optionally, calculating the volcano deformation value corresponding to the first grid includes:

[0025] 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;

[0026] 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;

[0027] Acquire the horizontal coordinates of the left boundary point and the right boundary point, input the horizontal coordinate of the left boundary point into the first Bessel function to calculate the corresponding first function value, and input the horizontal coordinate of the right boundary point into the second Bessel function to calculate the corresponding second function value;

[0028] 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 a whole as the volcanic deformation value corresponding to the first grid.

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

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

[0031] 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;

[0032] The zero values ​​at the first element position and the second element position are replaced with the first volcanic deformation value and the second volcanic deformation value respectively to obtain a new matrix to be assigned.

[0033] Optionally, generating the two-dimensional migration seismic profile corresponding to the two-dimensional vertical profile according to the volcanic deformation matrix includes:

[0034] Acquire location data of multiple key layers in the target underground area and elastic parameter data of underground media;

[0035] Generate a reflectivity matrix corresponding to the two-dimensional vertical section 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;

[0036] Introducing local fold features and fault features into the reflectivity matrix to obtain a fault reflectivity matrix;

[0037] Performing coordinate mapping on the fault reflectivity matrix based on the volcanic deformation matrix to obtain a volcanic reflectivity matrix;

[0038] 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 migration seismic profile.

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

[0040] 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:

[0041] Expanding the layer site according to the position coordinates of the layer site in the position data of each key layer and the number of horizontal grid rows to generate a layer site array for each key layer, wherein the number of elements in each layer site array is equal to the number of horizontal grid rows;

[0042] 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;

[0043] Introducing a global wrinkle feature into each of the layer arrays to obtain a corresponding plurality of global wrinkle layer arrays;

[0044] generating a two-dimensional horizon matrix according to the plurality of global fold horizon arrays;

[0045] Generate a two-dimensional elastic parameter matrix according to the two-dimensional layer 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 invention provides a two-dimensional migration seismic profile generating device, comprising:

[0048] A division 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;

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

[0050] A construction unit, used to construct a zero matrix of corresponding dimension according to the number of horizontal grid rows and the number of vertical grid columns as a matrix to be assigned;

[0051] A processing unit, configured to determine a vertical starting grid corresponding to the volcanic channel on the grid division interface 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, used for 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;

[0053] An assignment unit, used for assigning values ​​to the elements in the matrix to be assigned values ​​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 two-dimensional migration seismic profile generation method, device, equipment and medium provided by the present invention 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 invention can effectively realize the numerical simulation of the two-dimensional migration seismic profile containing the volcanic channel, and can improve the accuracy of the generated two-dimensional migration seismic profile containing the volcanic channel. BRIEF DESCRIPTION OF THE DRAWINGS

[0059] In order to more clearly illustrate the technical solutions in the present invention or the prior art, the following briefly introduces 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 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 cross-sectional schematic diagram corresponding to a fault reflectivity matrix provided in an embodiment of the present invention;

[0062] Figure 3 A cross-sectional schematic diagram corresponding to a volcano deformation matrix provided in an embodiment of the present invention;

[0063] Figure 4 A time series curve diagram 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 diagram of the structure of a computer device provided in an embodiment of the present invention. DETAILED DESCRIPTION

[0067] In order to make the purpose, technical solution and advantages of the present invention clearer, the technical solution of the present invention will be clearly and completely described below in conjunction with the drawings of the present invention. Obviously, the described embodiments are part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0068] Combine the following Figure 1-Figure 5 The two-dimensional migration seismic profile generation method of the present invention is described.

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

[0070] S101. On a two-dimensional vertical section of the target underground area, grid division is performed in the horizontal direction and the vertical direction 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 according to 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] Optionally, the two-dimensional vertical section can be a depth domain section, in which case the horizontal axis of the two-dimensional vertical section is the length axis and the vertical axis is the depth axis. The two-dimensional vertical section can also be a time domain section, in which case the horizontal axis of the two-dimensional 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 direction and the vertical direction to obtain a corresponding grid division interface. It can be understood that the grid division interface can include multiple grids, each grid has the same size, and the lateral size and vertical size of a single grid 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, constructing a zero matrix of corresponding dimension according to the number of horizontal grid rows and the number of vertical grid columns and using 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 according to 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 a 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] Among them, the left boundary curve can be a curve corresponding to the left boundary of the volcanic channel in the two-dimensional vertical section, and the right boundary curve can be a curve corresponding to the right boundary of the volcanic channel in the two-dimensional vertical section.

[0085] S105 , traversing 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.

[0086] Specifically, in this embodiment, in the mesh division interface, starting from the vertical starting mesh, the mesh traversal can be performed downward along the vertical line where the vertical starting mesh is located.

[0087] S106 , assigning values ​​to the elements in the matrix to be assigned according to each traversed grid, to obtain a matrix after the assignment is completed.

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

[0089] Optionally, step S106 may include:

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

[0091] 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 values ​​to the elements in the matrix to be assigned, a new matrix to be assigned is obtained, and the second grid is traversed;

[0092] 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.

[0093] It is understandable that, in this embodiment, when it is determined that a certain grid reaches the bottom of the volcanic channel, it is determined that the assignment is completed, and the latest matrix to be assigned is determined as the matrix after the assignment is completed.

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

[0095] 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;

[0096] 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;

[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 disturbed left boundary point and the disturbed 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, this embodiment may draw a vertical line from the middle position point to the target line segment, determine the intersection of the vertical line and the target line segment, and determine the intersection as the projection point from the middle position point to the target line segment.

[0110] Specifically, in this embodiment, a corresponding first Bezier curve may be constructed according to the left boundary point, the projection point and the middle position point, and a corresponding second Bezier curve may be constructed according to the right boundary point, the projection point and the middle 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 volcano deformation values ​​are used to assign values ​​to the elements in the matrix to be assigned, so as to obtain a new matrix to be assigned, including:

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

[0114] 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;

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

[0116] Specifically, in this embodiment, in the mesh division interface, a mesh including left boundary points may be determined as a first mesh, and a mesh including right boundary points may be determined as a second mesh.

[0117] It should be noted that the matrix to be assigned is constructed according to the number of horizontal grid rows and vertical grid columns in the grid division interface in this embodiment, and the number of element rows and the number of element columns in the matrix to be assigned are respectively 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 matrix to be assigned can correspond to the position of a grid in the grid division interface, and all elements in the matrix to be assigned are in a one-to-one correspondence with all grids in the grid division interface.

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

[0119] S107, using the matrix after the assignment as the volcano deformation matrix.

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

[0121] S108. Generate a two-dimensional offset 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 offset seismic data of the target underground area.

[0123] The two-dimensional migration seismic profile generation method proposed in this embodiment 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. This embodiment can effectively realize the numerical simulation of the two-dimensional migration seismic profile containing the volcanic channel, and can improve the accuracy of the generated two-dimensional migration seismic profile containing the volcanic channel.

[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. Among them:

[0125] S1081. Obtain 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, horizontal grid row number, and vertical grid column number 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 sites are expanded to generate a layer site array of each key layer, and the number of elements in each layer site array is equal to the number of horizontal grid rows;

[0135] The elastic parameter value is expanded 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, obtaining 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] Generate a reflectivity matrix based on a 2D elastic parameter matrix.

[0140] Specifically, in this embodiment, when expanding the layer site to obtain the layer site array of each key layer, a linear spline fitting can be performed on the position coordinates of all layer sites in each layer dictionary Z(i) obtained. 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 polynomials after fitting in each layer are 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 put into the layer array in turn, so that NL unequal-length layer dictionaries Z(i), i∈NL are converted 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 to ensure 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, that is, the longitudinal wave velocity dictionary 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 value between the key layers. The implementation method of the random elastic parameter increment is:

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

[0145] Among them, U(0,1) represents the random selection of random numbers 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 density array ρ(z),z∈NZ, the array contains the corresponding longitudinal wave velocity value and density value, and the array is arranged in order from small to large according to z. The one-dimensional elastic parameter array after random position processing is consistent with the actual data in the vertical scale, ensuring the scale of the simulation. At the same time, the random increment introduced between the layers simulates the uneven changes in the actual underground lithology and enhances the elastic parameters' ability to characterize the complex underground structure. In the subsequent wave field convolution simulation process, the random size changes of the elastic parameters are equivalent to adding scattering points and reasonable background noise to the simulation process, and the simulation results are more in line with the actual situation.

[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 type of 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 P 0 +2(1-x)xP 1 +x 2 P 2 .

[0148] Among them, x is an independent variable, and the number of its values ​​corresponds to the range of fold deformation. For global transverse fold deformation, the number of x is NX. 0 is the starting point coordinate corresponding to the Bessel basis function, p 1 is the coordinate of the control point corresponding to the Bessel basis function, P 2 is the end point coordinate corresponding to the Bessel basis function.

[0149] After the layer array is transformed by the random Bessel function, the layer array Z with global wrinkles is obtained. g (x), x∈NX. The global fold layer array is introduced, namely the global fold layer array Z g (x) is:

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

[0151] Then, the gradient linear interpolation method is used to calculate the NL global fold layer arrays Z g (x), x∈NX is used for two-dimensional profile interpolation to achieve the continuity of vertical gradient changes in the layer, thereby avoiding the occurrence of step-type anomalies in the layer.

[0152] Specifically, in the process of generating a two-dimensional layer matrix according to multiple global fold layer arrays, the midpoint can be selected as the reference position Z in a single global fold layer array. g (x 0 ), calculate the difference Z between the horizontal layer and the reference position g (x)-Z g (x 0 ). Take the reference position Z of each key layer as g (x 0 ) as the independent variable and the difference as the dependent variable, and the target polynomial is obtained by polynomial fitting:

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

[0154] Among them, t is the independent variable, f x (t) is the dependent variable, a x ,b x ,c x ,d x are the coefficients of the polynomial, which are obtained by dividing the values ​​of the NL independent variables by Z g (x 0 ) into t, dependent variable Z g (x)-Z g (x 0 ) into 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 (x 0 ,z) as the independent variable, interpolation Z g (x,z)-Z g (x 0 ,z) as the dependent variable to fit the 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 the density matrix ρ(x,z),x∈NX,z∈NZ.

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

[0157]

[0158] The value range of the vertical axis subscript z is [1, NZ-1], the value range of the horizontal axis subscript x is [1, NX], and the reflectivity coefficient value of the grid point 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 R 1 (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 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 value 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, in this embodiment, the reflectivity matrix R 1 (x,z),x∈NX,z∈NZ introduces fault characteristics to form a reflectivity model including fault structure, namely the fault reflectivity matrix R 2 (x, z), x∈NX, z∈NZ. In this embodiment, the fault feature can be introduced through the execution process of 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 of the fault starting position [z,x].

[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 indicating 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] Creates a new reflectivity matrix new_reflectivity which 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 new_fault_pos_z of the fault at its current x position.

[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, in this embodiment, after generating the volcano deformation matrix, the volcano deformation matrix can be applied to the reflection coefficient R by using the coordinate mapping method. 2 (x,z),x∈NX,z∈NZ, the reflectivity matrix R containing the volcanic channel structure is formed 3 (x,z),x∈NX,z∈NZ.

[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 profile.

[0195] left_shift and right_shift: Generate random displacement functions for subsequent adjustment of boundary points P 0 and P 2 The vertical axis produces a small disturbance.

[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 P of the left and right boundaries left_bond and right_bond in the current grid 0 and P 2 ;

[0199] P 0 and P 2 The vertical coordinate of the random displacement left_shift and right_shift is increased to get P 0 ′ and P 2 ′.

[0200] Create a random middle position P 1 , as a reference point between the left and right boundaries.

[0201] Calculate the intermediate position P 1 The displacement between the two boundary positions is X 0 and X 2 , according to X 0 and X 2The size of the volcano determines whether the bottom of the volcanic channel has been reached.

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

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

[0204] Determine the displacement X between the left and right boundaries 0 and X 2 Are they all greater than zero:

[0205] If so, calculate two Bezier curves. 1 In P 0 and P 2 The projection of the connecting line is taken as the middle point P k , respectively (P 0 ′,P k ,P 1 ) and (P 1 ,P k ,P 2 ′) 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, in P 0 and P 2 In the corresponding horizontal position, fill the values ​​of the two Bezier curves into the matrix to be assigned.

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

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

[0209] S1084. Coordinate mapping of the fault reflectivity matrix is ​​performed based on the volcanic deformation matrix to obtain the volcanic reflectivity matrix.

[0210] Specifically, after generating the volcanic deformation matrix Q(x,z), the fault reflectivity matrix R can be obtained. 2 (x,z),x∈NX,z∈NZ position coordinate matrix Q 0 (x,z),x∈NX,z∈NZ, then the new coordinate matrix Q with the volcanic channel 1 (x,z) is:

[0211] Q 1 (x,z)=Q 0 (x,z)+Q(x,z).

[0212] Afterwards, this embodiment can be Q 0 (x,z) as independent variables, with Q1 (x,z) is substituted as the dependent variable into the above target polynomial for polynomial fitting. The fitted polynomial will be R 2 Substituting (x,z) into the volcanic channel, we can obtain the reflectivity matrix introduced into the volcanic channel, namely 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 offset seismic profile.

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

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

[0217] The convolution is performed along the longitudinal 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 the seismic data of the relevant technology, volcanic channels have relatively clear identification features. The relatively narrow volcanic channel phase presents a chaotic high-frequency signal, is at a certain angle to the continuous axis of the surrounding strata, and has a relatively obvious upward flow trend. Therefore, the relevant technology can identify the volcanic channels in the offset seismic data through feature recognition methods. However, the volcanic channel phase runs through the volcanic rock development area, and the seismic attribute method is difficult to directly distinguish volcanic rocks of different microphases, and the workload of manual interpretation of volcanic channels is large. The supervised feature recognition method requires a large amount of volcanic channel phase data for feature learning. Therefore, obtaining offset seismic data containing volcanic channels has become an important research issue.

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

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

[0221] A division unit 601 is used to perform grid division in the horizontal direction and the vertical direction on the 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 determination 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 used to construct a zero matrix of corresponding dimension according to the number of horizontal grid rows and the number of vertical grid columns and use it as a 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 to generate the left boundary curve and the right boundary curve corresponding to the volcanic channel;

[0225] A traversal unit 605 is used 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 used to assign values ​​to the elements in the to-be-assigned matrix according to each traversed grid, so as to obtain a matrix after the assignment is completed;

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

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

[0229] It should be noted that the processing processes 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 respectively referred to in 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 values ​​to the elements in the matrix to be assigned, a new matrix to be assigned is obtained, 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] 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;

[0243] 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;

[0244] 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;

[0245] 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.

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

[0247] In the mesh division interface, a first mesh and a second mesh corresponding to the position coordinates of the left boundary point and the position coordinates of the right boundary point are respectively determined;

[0248] 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;

[0249] The zero values ​​at the first element position and the second element position are replaced with the first volcanic deformation value and the second volcanic deformation value respectively to obtain a new matrix to be assigned.

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

[0251] Acquire location data of multiple key layers in the target underground area and elastic parameter data of underground media;

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

[0253] Introducing local fold features and fault features into the reflectivity matrix, the fault reflectivity matrix is ​​obtained;

[0254] Based on the volcanic deformation matrix, coordinate mapping is performed on the fault reflectivity matrix to obtain the volcanic reflectivity matrix;

[0255] The 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 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 sites are expanded to generate a layer site array of each key layer, and the number of elements in each layer site array is equal to the number of horizontal grid rows;

[0259] The elastic parameter value is expanded 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, obtaining 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] Generate a reflectivity matrix based on a 2D elastic parameter matrix.

[0264] The two-dimensional migration seismic profile generation device proposed in this embodiment 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. This embodiment can effectively realize the numerical simulation of the two-dimensional migration seismic profile containing the volcanic channel, and can improve the accuracy of the generated two-dimensional migration seismic profile containing the volcanic channel.

[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 is shown.

[0267] See also Figure 7 , a schematic diagram of the structure 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 a dedicated integrated circuit, a programmable logic device or a combination thereof. The programmable logic device may be a complex programmable logic device, a field programmable logic gate array, a general purpose array logic or any combination thereof.

[0269] The memory 20 stores instructions executable by at least one processor 10, so that at least one processor 10 executes 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, an application required for at least one function. The data storage area may store data created according to 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 devices. In some optional embodiments, the memory 20 may optionally include a memory remotely arranged 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 a volatile memory, such as a random access memory. The memory may also include a non-volatile memory, such as a flash memory, a hard disk or a solid state drive. The memory 20 may also include a combination of the above-mentioned types of memory.

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

[0273] The embodiment of the present invention also provides a computer-readable storage medium. The method according to the embodiment of the present invention can be implemented in hardware, firmware, or can be implemented as a computer code that can be recorded in a storage medium, or can be implemented by downloading through a grid and originally stored in a remote storage medium or a non-temporary machine-readable storage medium and will be stored in a local storage medium, 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 hard disk, etc.; further, the storage medium can also include a combination of the above types of memories. 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 embodiments of the present invention.

Claims

1. A method for generating a two-dimensional migration seismic profile, characterized in that: include: On a two-dimensional vertical section of the target underground area, meshing is performed in the horizontal direction and the vertical direction to obtain a corresponding meshing interface; wherein the two-dimensional vertical section includes a volcanic channel; Determine the number of horizontal grid rows and the number of vertical grid columns in the grid division interface, and construct 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 position distribution of the volcanic channel, determining the vertical starting grid corresponding to the volcanic channel on the grid division interface, and generating the left boundary curve and the right boundary curve corresponding to the volcanic channel; Along the vertical line where the vertical starting grid is located, all the grids between the vertical starting grid and the last grid on the vertical line are traversed in sequence; Assign values ​​to the elements in the matrix to be assigned according to each of the traversed grids, and obtain a matrix after the assignment is completed 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 values ​​according to each of the traversed grids to obtain a matrix after the assignment is completed includes: For the first grid traversed, determine whether the first grid reaches the bottom of the volcanic channel; 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, the elements in the matrix to be assigned are assigned using the volcanic deformation value to obtain a new matrix to be assigned, and the second grid is traversed; 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 The determining whether the first grid reaches the bottom of the volcanic passage comprises: Determine 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; 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 respectively vertically displaced 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 middle 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 reaches the bottom of the volcanic passage.

4. The method according to claim 3, characterized in that The 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; Acquire the horizontal coordinates of the left boundary point and the right boundary point, input the horizontal coordinate of the left boundary point into the first Bessel function to calculate the corresponding first function value, and input the horizontal coordinate of the right boundary point into the second Bessel function to calculate the 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 a whole as the volcanic deformation value corresponding to the first grid.

5. The method according to claim 4, characterized in that The step of using the volcano deformation value to assign values ​​to the elements in the matrix to be assigned to obtain a new matrix to be assigned includes: In the mesh division interface, respectively determining a first mesh and a second mesh corresponding to the position coordinates of the left boundary point and the position coordinates of the right boundary point; 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 with 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 The step of generating a two-dimensional migration seismic profile corresponding to the two-dimensional vertical profile according to the volcanic deformation matrix includes: Acquire location data of multiple key layers in the target underground area and elastic parameter data of underground media; Generate a reflectivity matrix corresponding to the two-dimensional vertical section 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 features and fault features 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 migration seismic profile.

7. The method according to claim 6, characterized in that The position data of each key layer position includes the position coordinates of multiple layer position points, and the elastic parameter data includes the elastic parameter value corresponding to each key layer position; 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 and the number of horizontal grid rows to generate a layer site array for each key layer, 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; Generate a two-dimensional horizon matrix according to the plurality of global fold horizon arrays; Generate 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 division 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, used 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 dimension 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 a vertical starting grid corresponding to the volcanic channel on the grid division interface 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, used for 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; An assignment unit, used for assigning values ​​to the elements in the matrix to be assigned values ​​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 offset 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