Model construction for improving seismic data density and method for improving seismic data density

Abstract

The invention discloses a model construction method for improving seismic data density and a method for improving seismic data density. The method comprises the following steps: acquiring data of a plurality of seismic sections of a three-dimensional seismic data volume of a target layer section; respectively acquiring corresponding profile matrixes, inserting a plurality of columns into each profile matrix, setting data in the matrix after the columns are inserted, and generating an original data matrix, a to-be-encrypted data matrix, a random data matrix and an identification matrix corresponding to each profile matrix as training samples; inputting each training sample into a pre-built deep learning model for training, calculating a loss function of the deep learning model, adjusting a first parameter of the deep learning model according to a result of the loss function until the loss function reaches a preset convergence condition, and obtaining a trained improved seismic data density model; and the seismic data density improving model outputs a seismic data matrix after density improvement. The resolution of seismic data can be accurately and efficiently improved at low cost.

Description

Technical Field

[0001] The present invention relates to the field of oil and gas exploration, and in particular to a model construction for improving seismic data density and a method for improving seismic data density. Background Art

[0002] With the continuous deepening of oil and gas exploration, the exploration objects are becoming increasingly complex, especially the identification and prediction of thin layers, small faults, cracks, etc., which puts higher requirements on the resolution of seismic data. Improving the resolution of seismic data has always been a difficult problem for geophysicists. At present, there are two ways to improve the resolution of seismic data: increasing the frequency of seismic data and increasing the density of seismic data. Among them, the method of improving the frequency of seismic data is relatively mature, including enhancing effective signals through seismic data processing, expanding the effective frequency band of seismic data, and obtaining high-frequency components of seismic data by frequency division. A number of commercial software have been developed and have achieved good results in actual production applications. However, for the method of improving the density of seismic data, the existing technology mainly focuses on collecting high-density seismic data. Specifically, by increasing the spatial sampling rate to achieve alias-free sampling of the seismic wave field, although the useful signal and noise can be recorded comprehensively and with high fidelity, it is conducive to the subsequent effective signal recovery and noise removal. However, due to the reliance on the acquisition of seismic data, it is often costly and has a long cycle.

[0003] In recent years, with the rapid development of computer technology, artificial intelligence represented by deep learning has been widely used in all aspects of oil and gas exploration and development. Deep learning based on statistics can effectively improve the density and resolution of seismic data, which is of great significance for oil and gas exploration. Summary of the invention

[0004] In view of the above problems, the present invention is proposed to provide a method for constructing a model for improving seismic data density and improving seismic data density, which overcomes the above problems or at least partially solves the above problems.

[0005] In a first aspect, an embodiment of the present invention provides a model construction method for improving seismic data density, comprising:

[0006] Acquire a three-dimensional seismic data volume of a target layer segment, and extract data of multiple seismic sections from the three-dimensional seismic data volume; acquire a corresponding section matrix according to the data of each seismic section, wherein the columns in the section matrix represent the original seismic traces in the corresponding seismic section, and the data in each column represent the real seismic data in the corresponding seismic section;

[0007] Inserting a number of columns into each of the profile matrices, and setting the data in the profile matrix after the insertion of the number of columns, respectively generating an original data matrix, a data matrix to be encrypted, a random data matrix, and an identification matrix corresponding to each profile matrix as training samples, wherein the inserted number of columns characterizes the encrypted seismic traces, and the data on the number of columns characterizes the encrypted seismic data;

[0008] Each of the training samples is input into a pre-built deep learning model for training. During the iterative process of training, the loss function of the deep learning model is calculated. According to the result of the loss function, the first parameter of the deep learning model is adjusted until the loss function reaches a preset convergence condition, thereby obtaining a trained model for improving the density of seismic data; the model for improving the density of seismic data outputs a seismic data matrix with improved density.

[0009] In one embodiment, a plurality of columns are inserted into the profile matrix, and data in the profile matrix after the plurality of columns are inserted is set to generate an original data matrix, a data matrix to be encrypted, a random data matrix, and an identification matrix corresponding to the profile matrix as training samples, including:

[0010] Inserting a preset number of columns between every two adjacent columns in the profile matrix to obtain an extended matrix, wherein the extended matrix includes encrypted columns and original columns, wherein the encrypted columns are inserted columns and the original columns are non-inserted columns;

[0011] Copying the extended matrix to obtain a first copied matrix, and setting the data on the encrypted columns in the first copied matrix to zero, while leaving the data on the original columns unchanged, to obtain a matrix of data to be encrypted;

[0012] Copying the extended matrix to obtain a second copied matrix, and setting the data on the encrypted columns in the second copied matrix to random numbers and the data on the original columns to 0 to obtain a random data matrix;

[0013] Copy the extended matrix to obtain a third copied matrix, and set the data on the encrypted columns in the third copied matrix to 0 and the data on the original columns to 1 to obtain an identification matrix;

[0014] The data on the encrypted column in the matrix after the column insertion is set to be empty, and the data on the original column remains unchanged, to obtain the original data matrix;

[0015] The original data matrix, the data matrix to be encrypted, the random data matrix and the identification matrix are used as training samples corresponding to the profile matrix.

[0016] In one embodiment, the pre-built deep learning model includes a generator and a discriminator;

[0017] Input the training samples into a pre-built deep learning model for training, including:

[0018] Inputting the data matrix to be encrypted, the random data matrix and the identification matrix in the training sample into the generator for encryption, and obtaining the encrypted data matrix corresponding to the training sample;

[0019] The original data matrix, identification matrix and encrypted data matrix of the training samples are input into the discriminator to obtain the probability matrix output by the discriminator. The elements in the probability matrix represent the probability that the data corresponding to the element position in the encrypted data matrix is ​​real seismic data.

[0020] In one embodiment, the generator includes an encoder and a decoder;

[0021] The step of inputting the data matrix to be encrypted, the random data matrix and the identification matrix in the training sample into the generator for encryption to obtain the encrypted data matrix corresponding to the training sample comprises:

[0022] Inputting the to-be-encrypted data matrix, the random data matrix, and the identification matrix in the training sample into the encoder for feature extraction to obtain a first vector;

[0023] Using a semantic coding algorithm, reducing the dimension of the first vector to obtain a second vector;

[0024] The second vector is input into the long short-term memory recurrent neural network module of the decoder to obtain an encrypted data matrix corresponding to the training sample, and the encrypted data matrix has the same dimension as any matrix in the training sample.

[0025] In one embodiment, the encoder includes a convolution module and a bidirectional long short-term memory recurrent neural network module;

[0026] The method of inputting the to-be-encrypted data matrix, the random data matrix, and the identification matrix in the training sample into the encoder for feature extraction to obtain a first vector includes:

[0027] The data matrix to be encrypted, the random data matrix and the identification matrix in the training sample are input into the convolution module to perform feature extraction and feature fusion to obtain a fused feature vector;

[0028] The fused feature vector is input into the bidirectional long short-term memory recurrent neural network module for feature extraction to obtain a first vector.

[0029] In one embodiment, the original data matrix of the training sample, the identification matrix and the obtained encrypted data matrix are input into the discriminator to obtain the probability matrix output by the discriminator, including:

[0030] The original data matrix of the training sample and the obtained encrypted data matrix are input into the long short-term memory recurrent neural network module of the discriminator for feature extraction to obtain the third vector and the fourth vector respectively;

[0031] A probability matrix is ​​obtained according to the third vector and the fourth vector; the identification matrix is ​​used to indicate whether the data corresponding to each element position in the probability matrix is ​​real seismic data or encrypted seismic data.

[0032] In one embodiment, the loss function of the deep learning model includes a generator loss function and a discriminator loss function;

[0033] In the iterative process of training, the loss function of the deep learning model is calculated, and according to the result of the loss function, the value of the first parameter of the deep learning model is adjusted until the loss function reaches a preset convergence condition, thereby obtaining a trained model for improving seismic data density, including:

[0034] In the training iteration process, the value of the discriminator loss function is calculated until the discriminator loss function reaches a preset convergence condition, and the discriminator is fixed; in the training iteration process after the discriminator is fixed, the value of the generator loss function is calculated until the generator loss function reaches a preset convergence condition, and then the iteration is stopped to obtain a trained seismic data density improvement model; or

[0035] During the training iteration process, the value of the generator loss function is calculated until the generator loss function reaches the preset convergence condition, and the generator is fixed; during the training iteration process after the generator is fixed, the value of the discriminator loss function is calculated until the discriminator loss function reaches the preset convergence condition, then the iteration is stopped to obtain a trained model for improving seismic data density.

[0036] In one embodiment, the method further comprises:

[0037] A test sample is input into the improved seismic data density model completed in each iteration for testing, and the value of the root mean square error of the improved seismic data density model is calculated. According to the value of the root mean square error, the value of the second parameter of the improved seismic data density model is adjusted until the value of the root mean square error is less than or equal to a preset root mean square error threshold. The test sample is determined based on data of other seismic profiles of the three-dimensional seismic data body.

[0038] In a second aspect, an embodiment of the present invention provides a method for improving the density of seismic data, comprising:

[0039] Acquire data of a seismic section to be encrypted of a three-dimensional seismic data volume to be increased in density; acquire a corresponding section matrix to be encrypted according to the data of the seismic section to be encrypted;

[0040] Inserting a number of columns into the section matrix to be encrypted to set the data in the matrix after the inserted columns, generating a data matrix to be encrypted, a random data matrix and an identification matrix corresponding to the section matrix to be encrypted;

[0041] Inputting the to-be-encrypted data matrix, the random data matrix and the identification matrix into a model for increasing the density of seismic data to obtain a seismic data matrix with increased density;

[0042] The model for improving the density of seismic data is obtained through the aforementioned model building method for improving the density of seismic data.

[0043] In one embodiment, the method further comprises: converting the obtained density-enhanced seismic data matrix into a seismic profile.

[0044] In a third aspect, an embodiment of the present invention provides a model building device for improving seismic data density, comprising:

[0045] The first acquisition module is used to acquire a three-dimensional seismic data volume of a target layer segment, and extract data of multiple seismic sections from the three-dimensional seismic data volume; according to the data of each seismic section, a corresponding section matrix is ​​acquired, wherein the columns in the section matrix represent the original seismic traces in the corresponding seismic section, and the data in each column represent the real seismic data in the corresponding seismic section;

[0046] The first generation module is used to insert a number of columns into each of the profile matrices, and set the data in the profile matrix after the insertion of the number of columns, and respectively generate the original data matrix, the data matrix to be encrypted, the random data matrix and the identification matrix corresponding to each profile matrix as training samples, the inserted number of columns characterize the encrypted seismic traces, and the data on the number of columns characterize the encrypted seismic data;

[0047] A training module is used to input each of the training samples into a pre-built deep learning model for training. During the iterative process of training, the loss function of the deep learning model is calculated, and according to the result of the loss function, the first parameter of the deep learning model is adjusted until the loss function reaches a preset convergence condition, thereby obtaining a trained model for improving the density of seismic data; the model for improving the density of seismic data outputs a seismic data matrix with improved density.

[0048] In a fourth aspect, an embodiment of the present invention provides a device for improving seismic data density, comprising:

[0049] The second acquisition module is used to acquire data of the seismic section to be encrypted of the three-dimensional seismic data volume to be increased in density; and acquire the corresponding section matrix to be encrypted according to the data of the seismic section to be encrypted;

[0050] The second generating module is used to insert a number of columns into the section matrix to be encrypted to set the data in the matrix after the inserted columns, and generate a data matrix to be encrypted, a random data matrix and an identification matrix corresponding to the section matrix to be encrypted;

[0051] A density improvement module, used for inputting the to-be-encrypted data matrix, the random data matrix and the identification matrix into a seismic data density improvement model to obtain a seismic data matrix after density improvement;

[0052] The model for improving the density of seismic data is obtained through the aforementioned model building method for improving the density of seismic data.

[0053] In a fifth aspect, an embodiment of the present invention provides a computer storage medium, in which computer executable instructions are stored. When the computer executable instructions are executed by a processor, the aforementioned model construction method for improving seismic data density or the aforementioned method for improving seismic data density is implemented.

[0054] In a sixth aspect, an embodiment of the present invention provides a terminal device, characterized in that it comprises: a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein when the processor executes the program, it implements the aforementioned model construction method for improving the density of seismic data or the aforementioned method for improving the density of seismic data.

[0055] The beneficial effects of the above technical solution provided by the embodiment of the present invention include at least:

[0056] The model building method for improving the density of seismic data provided by the embodiment of the present invention extracts data of multiple seismic profiles from a three-dimensional seismic data volume of a target layer segment, obtains a profile matrix corresponding to each seismic profile according to the data of each seismic profile, each column in the profile matrix represents the original seismic trace in the corresponding seismic profile, and correspondingly, the data in each column in the profile matrix represents the real seismic data in the corresponding seismic profile, inserts several columns into each profile matrix, and sets the data in the profile matrix after inserting several columns, the inserted several columns represent the encrypted seismic trace, the data on the inserted several columns represent the encrypted seismic data, and obtains the original data matrix, the data matrix to be encrypted, the random data matrix and the identification matrix corresponding to each profile matrix, and the corresponding data matrix of each profile matrix is ​​set. The original data matrix, the data matrix to be encrypted, the random data matrix and the identification matrix are used as training samples; each training sample is input into a pre-built deep learning model for training; in the iterative process of training, the loss function of the deep learning model is calculated; according to the result of the loss function, the first parameter of the deep learning model is adjusted until the loss function reaches a preset convergence condition, thereby obtaining a trained model for improving the density of seismic data; based on the deep learning model, the internal relationship and overall characteristics between the seismic data are found, and a nonlinear model of the real seismic data and the encryption target is established, and then the seismic data is reconstructed. The density of the reconstructed seismic data is greater than that of the original seismic data, and the seismic data is efficiently and accurately encrypted to improve the resolution of the seismic data.

[0057] In addition, the encoder part of the generator in the seismic data density model is improved, and the convolution structure is used to fully extract the spatial characteristics of the seismic data to ensure that the characteristics of the horizontal data are fully utilized. At the same time, the temporal characteristics of the seismic data are extracted by combining the bidirectional long short-term memory recurrent neural network structure to ensure that the longitudinal data characteristics and global characteristics are fully utilized, thereby improving the accuracy of the reconstructed seismic data.

[0058] Other features and advantages of the present invention will be described in the following description, and partly become apparent from the description, or understood by practicing the present invention. The purpose and other advantages of the present invention can be realized and obtained by the structures particularly pointed out in the written description, claims, and drawings.

[0059] The technical solution of the present invention is further described in detail below through the accompanying drawings and embodiments. BRIEF DESCRIPTION OF THE DRAWINGS

[0060] The accompanying drawings are used to provide a further understanding of the present invention and constitute a part of the specification. Together with the embodiments of the present invention, they are used to explain the present invention and do not constitute a limitation of the present invention. In the accompanying drawings:

[0061] Figure 1A flow chart of a model building method for improving seismic data density in an embodiment of the present invention;

[0062] Figure 2 is a flow chart of a method for training a model in an embodiment of the present invention;

[0063] Figure 3 A flow chart of a method for generating an encrypted data matrix in an embodiment of the present invention;

[0064] Figure 4 is a flow chart of a method for obtaining a first vector in an embodiment of the present invention;

[0065] Figure 5 A flow chart of a method for improving seismic data density in an embodiment of the present invention;

[0066] Figure 6 A schematic diagram of the structure of a model building device for improving seismic data density in an embodiment of the present invention;

[0067] Figure 7 It is a schematic diagram of the structure of a device for improving seismic data density in an embodiment of the present invention;

[0068] Figure 8 Schematic diagram of the structure of a model for improving seismic data density in an embodiment of the present invention

[0069] Fig. 9 This is a comparison diagram of the relationship between deep learning parameters and RMSE in an embodiment of the present invention;

[0070] Fig.10 This is a schematic diagram of an original seismic section before density enhancement in an embodiment of the present invention;

[0071] Fig.11 This is a schematic diagram of the original seismic profile after density improvement in an embodiment of the present invention. DETAILED DESCRIPTION

[0072] The exemplary embodiments of the present disclosure will be described in more detail below with reference to the accompanying drawings. Although the exemplary embodiments of the present disclosure are shown in the accompanying drawings, it should be understood that the present disclosure can be implemented in various forms and should not be limited by the embodiments set forth herein. On the contrary, these embodiments are provided to enable a more thorough understanding of the present disclosure and to fully convey the scope of the present disclosure to those skilled in the art.

[0073] By using deep learning of artificial intelligence, we can find the internal relationships and overall characteristics between seismic data, establish a nonlinear model of data and targets, and then reconstruct the seismic data. The reconstructed seismic data is denser than the original seismic data, which can effectively encrypt seismic data and improve the resolution of seismic data, which is of great significance for oil and gas exploration.

[0074] In order to solve the problems existing in the prior art, the embodiments of the present invention provide a method for building a model for improving the density of seismic data and improving the density of seismic data.

[0075] In order to facilitate the description of a model construction for improving seismic data density and a method for improving seismic data density provided by an embodiment of the present invention, the process of improving density is also referred to as an encryption process in this embodiment.

[0076] Example

[0077] Embodiment 1 of the present invention provides a model construction method for improving seismic data density, the process of which is as follows: Figure 1 As shown, the following steps are included:

[0078] Step S1: Acquire a three-dimensional seismic data volume of a target layer segment, and extract data of multiple seismic sections from the three-dimensional seismic data volume; acquire a corresponding section matrix according to the data of each seismic section, wherein the columns in the section matrix represent the original seismic traces in the corresponding seismic section, and the data in each column represent the real seismic data in the corresponding seismic section;

[0079] Step S2: inserting several columns into each profile matrix, and setting the data in the profile matrix after inserting several columns, respectively generating the original data matrix, the data matrix to be encrypted, the random data matrix and the identification matrix corresponding to each profile matrix as training samples, the inserted several columns characterize the encrypted seismic traces, and the data on the several columns characterize the encrypted seismic data;

[0080] Step S3: Input each training sample into a pre-built deep learning model for training. During the iterative process of training, calculate the loss function of the deep learning model. According to the result of the loss function, adjust the first parameter of the deep learning model until the loss function reaches a preset convergence condition, and obtain a trained model for improving the density of seismic data; the model for improving the density of seismic data outputs a seismic data matrix with improved density.

[0081] Generally, there are singular values ​​such as maximum and minimum in the data of seismic profiles. The existence of singular values ​​will cause a large span of the data value range, making it difficult for deep learning algorithm models to converge quickly. Therefore, it is necessary to standardize the data of each seismic profile to eliminate singular values ​​such as maximum and minimum, so as to improve the seismic data density model and converge quickly.

[0082] In some optional embodiments, the seismic traces in the seismic profile data are used as the original seismic traces. In the above step S1, the seismic profile is converted into a corresponding matrix according to the original seismic traces in the seismic profile and the sampling interval time in the seismic profile. Referring to Table 1 below, a profile matrix is ​​obtained. In Table 1 below, x 11 represents the amplitude of the first sampling point of the first original seismic trace, xtn represents the amplitude of the tth sampling point of the nth original seismic trace,

[0083] Table 1:

[0084]

[0085] In some optional embodiments, the above step S2 inserts several columns into the profile matrix, and sets the data in the matrix after the columns are inserted, to generate the original data matrix, the data matrix to be encrypted, the random data matrix and the identification matrix corresponding to the profile matrix as training samples, which can be specifically implemented in the following manner:

[0086] (1) inserting a preset number of columns between every two adjacent columns in the profile matrix to obtain an extended matrix, wherein the extended matrix includes encrypted columns and original columns, wherein the encrypted columns are inserted columns and the original columns are non-inserted columns;

[0087] Specifically, a column is inserted between every two adjacent columns in the profile matrix to obtain an extended matrix, as shown in Table 2 below. Table 2 has a total of n columns of data, representing n seismic traces, of which the first and third columns correspond to the first column (x 11 The column where it is located) and the second column (x 12 The column where the first and second columns are located is called the original column, representing the original seismic traces. The second and fourth columns in Table 2, that is, the inserted columns, are called the encrypted columns, representing the encrypted seismic traces.

[0088] Table 2:

[0089]

[0090] (2) Copy the extended matrix to obtain a first copy matrix, and set the data on the encrypted columns in the first copy matrix to zero, while leaving the data on the original columns unchanged, to obtain a matrix of data to be encrypted;

[0091] Specifically, the extended matrix shown in Table 2 is copied to obtain a first copied matrix, and the data on the encrypted column in the first copied matrix is ​​set to 0, and the data on the original column remains unchanged to obtain a matrix of data to be encrypted, which is recorded as Refer to Table 3 below:

[0092] Table 3:

[0093]

[0094] (3) Copy the extended matrix to obtain a second copy matrix, and set the data on the encrypted columns in the second copy matrix to random numbers and the data on the original columns to 0 to obtain a random data matrix;

[0095] Specifically, the extended matrix shown in Table 2 is copied to obtain a second copy matrix, the data on the encrypted column in the second copy matrix is ​​supplemented with a random number z, and the data on the original column is set to 0 to obtain a random data matrix, denoted as Z, as shown in Table 4 below:

[0096] Table 4:

[0097]

[0098]

[0099] (4) Copy the extended matrix to obtain a third copy matrix, and set the data on the encrypted columns in the third copy matrix to 0 and the data on the original columns to 1 to obtain an identification matrix;

[0100] Specifically, the extended matrix shown in Table 2 is copied to obtain a third copy matrix, and the data on the encrypted column in the third copy matrix is ​​set to 0, and the data on the original column is set to 1, to obtain an identification matrix, denoted as M. Referring to Table 5 below, the identification matrix is ​​used to record the position of missing data, 1 indicates that the data at the position exists, and 0 indicates that the data at the position is missing:

[0101] Table 5:

[0102]

[0103] (5) Set the data on the encrypted column in the matrix after the insertion column to empty, and keep the data on the original column unchanged, to obtain the original data matrix, denoted as X. Specifically, refer to Table 6,

[0104] Table 6:

[0105]

[0106] Where n is the number of seismic traces, t is the sampling time, then x tn is the sampling point data of the nth channel when the sampling time is t, where N represents the total number of traces in a profile, and T represents the total number of sampling points in each seismic trace.

[0107] (6) The original data matrix, the data matrix to be encrypted, the random data matrix and the identification matrix are used as training samples corresponding to the profile matrix.

[0108] In some optional embodiments, the pre-built deep learning model includes a generator and a discriminator. In the above step S3, the training sample is input into the pre-built deep learning model for training. Specifically, Figure 2 As shown, this can be achieved by:

[0109] Step S21: input the data matrix to be encrypted, the random data matrix and the identification matrix in the training sample into the generator for encryption, and obtain the encrypted data matrix corresponding to the training sample;

[0110] Step S22: Input the original data matrix of the training sample, the identification matrix and the obtained encrypted data matrix into the discriminator to obtain the probability matrix output by the discriminator. The elements in the probability matrix represent the probability that the data corresponding to the element position in the encrypted data matrix is ​​the real earthquake data.

[0111] In the above step S21, the data matrix to be encrypted The random data matrix Z and the identification matrix M are input into the generator, and the output result is the earthquake estimation matrix, that is, the encrypted data matrix, denoted as X′, as shown in Table 7 below:

[0112] Table 7:

[0113]

[0114] In some optional embodiments, the generator includes an encoder and a decoder, and the above step S21, reference Figure 3 As shown, this can be achieved by:

[0115] Step S31: inputting the to-be-encrypted data matrix, the random data matrix and the identification matrix in the training sample into the encoder for feature extraction to obtain a first vector;

[0116] Step S32: using a semantic coding algorithm to reduce the dimension of the first vector to obtain a second vector;

[0117] Step S33: Input the second vector into the long short-term memory recurrent neural network module of the decoder to obtain an encrypted data matrix corresponding to the training sample, and the encrypted data matrix has the same dimension as any matrix in the training sample.

[0118] The encoder and decoder are connected by a low-dimensional semantic coding vector. Specifically, the output data of the generator (first vector) is mapped into a low-dimensional semantic coding vector (second vector) through the semantic encoder, and the task of the decoder is to reconstruct the low-dimensional semantic coding vector into complete seismic data.

[0119] In some optional embodiments, the encoder includes a convolution module and a bidirectional long short-term memory recurrent neural network module. The above step S31, specifically, refers to Figure 4 As shown, this can be achieved by:

[0120] Step S41: input the to-be-encrypted data matrix, the random data matrix, and the identification matrix in the training sample into the convolution module, perform feature extraction and feature fusion, and obtain a fused feature vector;

[0121] The data matrix to be encrypted, the random data matrix and the identification matrix are input into the convolution module. The embodiment of the present invention adopts a multi-scale convolution structure. In the multi-scale convolution, the embodiment of the present invention adopts a three-channel mode: the first channel has only one convolution layer, the convolution kernel size is 1×1, the padding is 0, and the step length is 1; the second channel contains two convolution layers, the first convolution kernel size is 1×1, the padding is 1, and the step length is 1; the second convolution kernel size is 3×3, the padding is 1, the step length is 1, and the activation function is the ReLU function; the third channel contains three convolution layers, the first convolution kernel size is 1×1, the padding is 1, and the step length is 1, the second convolution kernel size is 3×3, the padding is 1, the step length is 1, and the third convolution kernel size is 5×5, and two 3×3 convolutions are used instead. The features extracted by the three channels are fused with the input layer after a 1×1 convolution. The spatial characteristics of the seismic data are fully extracted.

[0122] Specifically, feature extraction is achieved through the following formula:

[0123]

[0124] In formula (1), conv represents the convolution operation.

[0125] The features extracted by convolution at different scales are fused using the following formula:

[0126]

[0127] Step S42: Input the fused feature vector into the bidirectional long short-term memory recurrent neural network module for feature extraction to obtain a first vector.

[0128] Specifically, the fused features are input into BiLSTM, and the first vector is obtained using the following formula:

[0129]

[0130] In some optional embodiments, the decoder is composed of 64 LSTMs, which are used to restore the data dimension. Specifically, it is implemented using the following formula:

[0131]

[0132] In formula (4), C is the encoder operation, and X′ is the encrypted seismic data and identification matrix output by the encoder G.

[0133] In some optional embodiments, the above step S22 can be implemented in the following way:

[0134] (1) Inputting the original data matrix of the training sample and the obtained encrypted data matrix into the long short-term memory recurrent neural network module of the discriminator for feature extraction to obtain the third vector and the fourth vector respectively;

[0135] (2) A probability matrix is ​​obtained based on the third vector and the fourth vector; the identification matrix is ​​used to indicate whether the data corresponding to each element position in the probability matrix is ​​real earthquake data or encrypted earthquake data.

[0136] The discriminator is composed of a long short-term memory network (LSTM) module and a fully connected layer (FC) module. The LSTM module extracts the feature vectors of the encrypted seismic data X' and the original seismic data X, and the identification matrix M locates the position. The fully connected layer FC reduces the feature vector from high dimension to one dimension, and then calculates the probability of each encrypted data element being true or false through the activation function (Sigomid function). Finally, the obtained probability data is located by the identification matrix to form a probability matrix

[0137] In some optional embodiments, the loss function of the deep learning model includes a generator loss function and a discriminator loss function, and the above step S3 can be implemented in the following manner:

[0138] In the training iteration process, the value of the discriminator loss function is calculated until the discriminator loss function reaches a preset convergence condition, and the discriminator is fixed; in the training iteration process after the discriminator is fixed, the value of the generator loss function is calculated until the generator loss function reaches a preset convergence condition, and then the iteration is stopped to obtain a trained model for improving seismic data density; or

[0139] During the training iteration process, the value of the generator loss function is calculated until the generator loss function reaches the preset convergence condition, and the generator is fixed; during the training iteration process after the generator is fixed, the value of the discriminator loss function is calculated until the discriminator loss function reaches the preset convergence condition, then the iteration is stopped to obtain a trained model for improving seismic data density.

[0140] The error between the encrypted data matrix generated by the generator and the original data matrix is ​​taken as the discriminant error to construct the discriminator loss function; the error between the data matrix to be encrypted and the encrypted data matrix is ​​taken as the reconstruction error to construct the generator loss function. The generator loss function and the discriminator loss function are used to form the loss function of the seismic data density model through the following formula:

[0141] MGB LOSS =-(1-λ)D(X+Z)+λ‖M⊙(G(X+Z)-X)‖2 (5)

[0142] In formula (5), λ represents the number of the generator loss function, D represents the discriminator operation, G represents the generator operation, M represents the identification matrix, X represents the original data matrix, Z represents the random matrix, and ‖…‖ represents the l2 norm.

[0143] The generator generates an encrypted data matrix for training the discriminator, updates the parameters of the discriminator until the discriminator loss function reaches the preset convergence point, fixes the discriminator, updates the parameters of the generator through the output of the discriminator, and iterates repeatedly until the loss function represented by formula (5) reaches the convergence condition, stops iteration, and obtains a trained model for improving seismic data density.

[0144] In some optional embodiments, a test sample is input into the improved seismic data density model completed in each iteration for testing, and the root mean square error value of the improved seismic data density model is calculated. According to the value of the root mean square error, the value of the second parameter of the improved seismic data density model is adjusted until the value of the root mean square error is less than or equal to a preset root mean square error threshold. The test sample is determined based on data of other seismic profiles of the three-dimensional seismic data body.

[0145] In order to make the seismic data density improvement model converge faster, the seismic data density improvement model can be trained once per iteration, and the seismic data density improvement model can be tested with a test sample. Specifically, the test sample is generated by generating a training sample. The generation method has been described in the foregoing, and the embodiment of the present invention will not be repeated here. The test sample is input into the seismic data density improvement model, and the root mean square error of the model is calculated. The root mean square error is used as an evaluation index, as shown in the following formula:

[0146]

[0147] In formula (6), y t The real seismic data of the sampling point at time t in the encrypted seismic data; y t ′ represents the seismic data generated by the sampling point at time t in the encrypted seismic data,

[0148] If the loss function of the seismic data density model has converged and the value of the root mean square error is greater than the preset root mean square error threshold, the value of the second parameter in the seismic data density model is adjusted, for example, the second parameter is the learning rate, the number of hidden layer nodes, the parameter λ, the number of iterations and other parameters, until the loss function of the seismic data density model reaches the preset convergence condition and the value of the root mean square error is less than the preset root mean square error threshold.

[0149] Based on the aforementioned model building method for improving seismic data density, the embodiment of the present invention further provides a method for improving seismic data density, the process of which is as follows: Figure 5 As shown, the following steps are included:

[0150] Step S61: obtaining data of the seismic section to be encrypted of the three-dimensional seismic data volume to be increased in density; obtaining a corresponding section matrix to be encrypted according to the data of the seismic section to be encrypted;

[0151] Step S62: inserting a number of columns into the section matrix to be encrypted to set the data in the matrix after the inserted columns, and generating a data matrix to be encrypted, a random data matrix and an identification matrix corresponding to the section matrix to be encrypted;

[0152] Step S63: inputting the data matrix to be encrypted, the random data matrix and the identification matrix into the seismic data density improvement model to obtain the seismic data matrix after density improvement;

[0153] The model for improving the density of seismic data is obtained through the aforementioned model building method for improving the density of seismic data.

[0154] In some optional embodiments, after obtaining the seismic data matrix with increased density in the above step S13, the obtained seismic data matrix with increased density can also be converted into a seismic profile.

[0155] Based on the same inventive concept, the embodiment of the present invention also provides a model building device for improving the density of seismic data, the structure of which is as follows: Figure 6 As shown, including:

[0156] The first acquisition module 71 is used to acquire the three-dimensional seismic data volume of the target layer segment, and extract the data of multiple seismic sections from the three-dimensional seismic data volume; according to the data of each seismic section, the corresponding section matrix is ​​acquired, and the columns in the section matrix represent the original seismic traces, and the data in each column represent the real seismic data;

[0157] The first generating module 72 is used to insert a number of columns into each profile matrix, and set the data in the matrix after the columns are inserted, and respectively generate the original data matrix, the data matrix to be encrypted, the random data matrix and the identification matrix corresponding to each profile matrix as training samples, the inserted number of columns characterize the encrypted seismic traces, and the data on the number of columns characterize the encrypted seismic data;

[0158] The training module 73 is used to input each training sample into a pre-built deep learning model for training. During the iterative process of training, the loss function of the deep learning model is calculated, and the first parameter of the deep learning model is adjusted according to the result of the loss function until the loss function reaches a preset convergence condition, thereby obtaining a trained model for improving the density of seismic data; the model for improving the density of seismic data outputs a seismic data matrix with improved density.

[0159] Regarding the device for improving the density of seismic data in the above-mentioned embodiment, the specific manner in which each module performs operations has been described in detail in the embodiment of the method, and will not be elaborated here.

[0160] Based on the same inventive concept, the embodiment of the present invention further provides a device for improving the density of seismic data. The structure of the device is as follows: Figure 7 As shown, including:

[0161] The second acquisition module 81 is used to acquire data of the seismic section to be encrypted of the three-dimensional seismic data volume to be increased in density; and acquire the corresponding section matrix to be encrypted according to the data of the seismic section to be encrypted;

[0162] The second generating module 82 is used to insert a number of columns into the section matrix to be encrypted to set the data in the matrix after the inserted columns, and generate a data matrix to be encrypted, a random data matrix and an identification matrix corresponding to the section matrix to be encrypted;

[0163] The density increasing module 83 is used to input the data matrix to be encrypted, the random data matrix and the identification matrix into the seismic data density increasing model to obtain the seismic data matrix after density increasing;

[0164] The model for improving the density of seismic data is obtained through the aforementioned model building method for improving the density of seismic data.

[0165] Regarding the device for improving the density of seismic data in the above-mentioned embodiment, the specific manner in which each module performs operations has been described in detail in the embodiment of the method, and will not be elaborated here.

[0166] The above-mentioned model construction method and method for improving seismic data density are described by a specific embodiment. Taking the Cambrian platform margin belt in Gucheng area of ​​Tarim Basin in northwest China as an example, the specific steps are as follows:

[0167] 1. Data preprocessing

[0168] Preprocess the 3D seismic data volume of the target layer in the study area; extract 1000 sections of seismic data with equal spacing (1×1) between the main survey line and the connecting survey line, of which 800 are used as training data and 200 are used as test data. Each seismic section contains 100 seismic traces, each containing 500 data, and the sampling interval is 4ms. On this basis, each seismic section in the 800 training data is converted into the original data matrix, the data matrix to be encrypted, the random matrix and the identification matrix;

[0169] 2. Build the model generator

[0170] The generator of the model consists of three parts: encoder, decoder and low-dimensional semantic encoding. The encoder part consists of multi-scale convolution and BiLSTM. The encoder and decoder are connected by a low-dimensional semantic encoding vector. The generator input is mapped to a low-dimensional semantic encoding vector through the encoder, and the decoder's task is to reconstruct the vector into complete seismic data. The decoder consists of 64 LSTMs, which are used to restore the data dimension and generate the input of the discriminator;

[0171] 3. Build the discriminator of the model

[0172] The discriminator of the model is mainly composed of LSTM modules and FC modules. Among them, the LSTM module consists of 64 LSTM units, which is responsible for processing the generated seismic data and obtaining the historical memory vector from the encrypted data matrix and the original data matrix. The fully connected layer maps the historical memory vector into a probability output, that is, the Sigmoid function is used to calculate the probability that the input data is true. The training goal of the discriminator is to identify the real samples as true as much as possible and the generated samples as false. The generator and the discriminator form a model to improve the density of seismic data. Figure 8 As shown;

[0173] 4. Input the training samples into the generator, and the generator generates encrypted seismic data for training the discriminator and updating the parameters of the discriminator. Fix the discriminator and update the parameters of the generator based on the results of the discriminator. Repeat the iteration until the MC-GAN-BiLSTM model reaches the optimal value.

[0174] 5. Use the test samples to test the model for improving seismic data density. By continuously adjusting the learning rate, number of hidden layer nodes, parameter λ, number of iterations and other parameters in deep learning, and monitoring the changes of the evaluation index RMSE in real time, until the value of RMSE is less than or equal to the preset root mean square error threshold, the trained model for improving seismic data density is obtained. Fig. 9 As shown, Fig. 9 (a) is the learning rate curve; (b) is the hidden layer node number curve; (c) is the parameter λ curve; (d) is the number of iterations. When the learning rate is 0.001 and the number of hidden layer nodes is 16, the RMSE is the lowest and the model performance is the best. Therefore, the number of hidden layer nodes and the learning rate of the BiLSTM layer are set to 16 and 0.001 respectively, and then the λ parameter of the model is adjusted. When the parameter λ is close to 1 or 0, the encryption effect is poor. This means that the two parts of the generator loss function do not play their due role. When the parameter λ is set to 0.3, the RMSE is the lowest and the completion effect is the best. When the number of iterations is 1000, the RMSE is basically stable. Therefore, the maximum number of iterations is set to 1000. The learning rate, the number of hidden layer nodes, the parameter λ and the number of iterations are set to 0.001, 16, 0.3 and 1000 respectively.

[0175] 6. Use the trained model to improve the density of seismic data to encrypt the seismic data to be encrypted, such as Fig.10 As shown in the figure, it is the original seismic section before encryption. Fig.11 As shown in the figure, the seismic profile after density improvement, the original seismic data trace spacing is 25 meters, and the trace spacing after doubling the seismic data density is equivalent to 12.5 meters. The improved seismic profile has a clearer structure of the progradational mound and shoal body in the Cambrian platform margin belt in the target area, with a higher resolution, and can better reflect the distribution of multi-stage platform margin mound and shoal.

[0176] Due to the development of carbonate platform margin mound-shoal deposits in the Cambrian in the Gucheng area, multiple wells have revealed the development of multi-stage mound-shoal bodies in the platform margin. Through well-seismic calibration, the multi-stage mound-shoal bodies appear as the reflection characteristics of the superimposed foreset mounds on the seismic profile. However, due to the large burial depth of the target area and the low main frequency of the seismic data, the low-resolution seismic data is not obvious for the foreset reflection characteristics of the mound-shoal bodies. How to identify the multi-stage mound-shoal bodies in the platform margin belt through seismic data is a difficulty in the target area. The identification of mound-shoal bodies by increasing the frequency of seismic data through conventional explanatory processing such as frequency extension or frequency division of seismic data often creates many illusions and has poor results. The present invention is used to increase the density of seismic data and thus increase the resolution of seismic data to identify the Cambrian mound-shoal bodies in the Gucheng area to achieve a better effect.

[0177] Based on the same inventive concept, an embodiment of the present invention also provides a computer storage medium, in which computer executable instructions are stored. When the computer executable instructions are executed by a processor, the aforementioned model construction method for improving seismic data density or the aforementioned method for improving seismic data density is implemented.

[0178] Based on the same inventive concept, an embodiment of the present invention also provides a terminal device, including: a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein when the processor executes the program, the aforementioned model construction method for improving the density of seismic data or the aforementioned method for improving the density of seismic data is implemented.

[0179] Unless otherwise specifically stated, terms such as processing, computing, calculating, determining, displaying, etc. may refer to the actions and / or processes of one or more processing or computing systems, or similar devices, which operate and convert data represented as physical (e.g., electronic) quantities within registers or memories of a processing system into other data similarly represented as physical quantities within memories, registers, or other such information storage, transmission, or display devices of the processing system. Information and signals may be represented using any of a variety of different techniques and methods. For example, data, instructions, commands, information, signals, bits, symbols, and chips mentioned throughout the above description may be represented by voltages, currents, electromagnetic waves, magnetic fields or particles, light fields or particles, or any combination thereof.

[0180] It should be understood that the specific order or hierarchy of steps in the disclosed process is an example of an exemplary method. Based on design preferences, it should be understood that the specific order or hierarchy of steps in the process can be rearranged without departing from the scope of protection of the present disclosure. The attached method claims present the elements of the various steps in an exemplary order and are not intended to be limited to the specific order or hierarchy described.

[0181] In the above detailed description, various features are grouped together in a single embodiment to simplify the disclosure. This method of disclosure should not be interpreted as reflecting an intention that the embodiments of the claimed subject matter require more features than are clearly stated in each claim. On the contrary, as reflected in the appended claims, the invention is in a state of having less than all the features of the disclosed individual embodiments. Therefore, the appended claims are hereby expressly incorporated into the detailed description, with each claim standing on its own as a separate preferred embodiment of the invention.

[0182] Those skilled in the art will also appreciate that the various illustrative logic blocks, modules, circuits, and algorithmic steps described in conjunction with the embodiments herein can all be implemented as electronic hardware, computer software, or a combination thereof. In order to clearly illustrate the interchangeability between hardware and software, various illustrative components, blocks, modules, circuits, and steps are generally described above around their functions. Whether such functions are implemented as hardware or software depends on specific applications and the design constraints imposed on the entire system. A skilled person can implement the described functions in an alternative manner for each specific application, but such implementation decisions should not be interpreted as departing from the scope of protection of the present disclosure.

[0183] The steps of the method or algorithm described in conjunction with the embodiments herein may be directly embodied as hardware, a software module executed by a processor, or a combination thereof. The software module may be located in a RAM memory, a flash memory, a ROM memory, an EPROM memory, an EEPROM memory, a register, a hard disk, a mobile disk, a CD-ROM, or any other form of storage medium well known in the art. An exemplary storage medium is connected to the processor so that the processor can read information from the storage medium and can write information to the storage medium. Of course, the storage medium may also be an integral part of the processor. The processor and the storage medium may be located in an ASIC. The ASIC may be located in a user terminal. Of course, the processor and the storage medium may also be present in a user terminal as discrete components.

[0184] For software implementation, the techniques described in this application can be implemented with modules (e.g., procedures, functions, etc.) that perform the functions described in this application. These software codes can be stored in a memory unit and executed by a processor. The memory unit can be implemented within the processor or outside the processor. In the latter case, it is coupled to the processor in a communication manner via various means, which are well known in the art.

[0185] The above description includes examples of one or more embodiments. Of course, it is impossible to describe all possible combinations of components or methods for the purpose of describing the above embodiments, but it should be recognized by those skilled in the art that the various embodiments may be further combined and arranged. Therefore, the embodiments described herein are intended to cover all such changes, modifications and variations that fall within the scope of protection of the appended claims. In addition, with respect to the term "comprising" used in the specification or claims, the word is covered in a manner similar to the term "including", just as "including," is explained as a transitional word in the claims. In addition, any term "or" used in the specification of the claims is intended to mean "non-exclusive or".

Claims

1. A model construction method for improving seismic data density, characterized in that: include: Acquire a three-dimensional seismic data volume of a target layer segment, and extract data of a plurality of seismic sections from the three-dimensional seismic data volume; According to the data of each seismic profile, a corresponding profile matrix is ​​obtained, wherein the columns in the profile matrix represent the original seismic traces in the corresponding seismic profile, and the data in each column represent the real seismic data in the corresponding seismic profile; Inserting a number of columns into each of the profile matrices, and setting the data in the profile matrix after the insertion of the number of columns, respectively generating an original data matrix, a data matrix to be encrypted, a random data matrix, and an identification matrix corresponding to each profile matrix as training samples, wherein the inserted number of columns characterizes the encrypted seismic traces, and the data on the number of columns characterizes the encrypted seismic data; Each of the training samples is input into a pre-built deep learning model for training. During the iterative process of training, the loss function of the deep learning model is calculated. According to the result of the loss function, the first parameter of the deep learning model is adjusted until the loss function reaches a preset convergence condition, thereby obtaining a trained model for improving the density of seismic data; the model for improving the density of seismic data outputs a matrix of seismic data with improved density.

2. The method according to claim 1, characterized in that Inserting a number of columns into the profile matrix, and setting the data in the profile matrix after inserting the number of columns, generating an original data matrix, a data matrix to be encrypted, a random data matrix, and an identification matrix corresponding to the profile matrix as training samples, including: Inserting a preset number of columns between every two adjacent columns in the profile matrix to obtain an extended matrix, wherein the extended matrix includes encrypted columns and original columns, wherein the encrypted columns are inserted columns and the original columns are non-inserted columns; Copying the extended matrix to obtain a first copied matrix, and setting the data on the encrypted columns in the first copied matrix to zero, while leaving the data on the original columns unchanged, to obtain a matrix of data to be encrypted; Copying the extended matrix to obtain a second copied matrix, and setting the data on the encrypted columns in the second copied matrix to random numbers and the data on the original columns to 0 to obtain a random data matrix; Copy the extended matrix to obtain a third copied matrix, and set the data on the encrypted columns in the third copied matrix to 0 and the data on the original columns to 1 to obtain an identification matrix; The data on the encrypted column in the matrix after the column insertion is set to be empty, and the data on the original column remains unchanged, to obtain the original data matrix; The original data matrix, the data matrix to be encrypted, the random data matrix and the identification matrix are used as training samples corresponding to the profile matrix.

3. The method according to claim 2, characterized in that The pre-built deep learning model includes a generator and a discriminator; Input the training samples into a pre-built deep learning model for training, including: Inputting the data matrix to be encrypted, the random data matrix and the identification matrix in the training sample into the generator for encryption, and obtaining the encrypted data matrix corresponding to the training sample; The original data matrix, identification matrix and encrypted data matrix of the training samples are input into the discriminator to obtain the probability matrix output by the discriminator. The elements in the probability matrix represent the probability that the data corresponding to the element position in the encrypted data matrix is ​​real seismic data.

4. The method according to claim 3, characterized in that The generator includes an encoder and a decoder; The step of inputting the data matrix to be encrypted, the random data matrix and the identification matrix in the training sample into the generator for encryption to obtain the encrypted data matrix corresponding to the training sample comprises: Inputting the to-be-encrypted data matrix, the random data matrix, and the identification matrix in the training sample into the encoder for feature extraction to obtain a first vector; Using a semantic coding algorithm, reducing the dimension of the first vector to obtain a second vector; The second vector is input into the long short-term memory recurrent neural network module of the decoder to obtain an encrypted data matrix corresponding to the training sample, and the encrypted data matrix has the same dimension as any matrix in the training sample.

5. The method according to claim 4, characterized in that The encoder includes a convolution module and a bidirectional long short-term memory recurrent neural network module; The method of inputting the to-be-encrypted data matrix, the random data matrix, and the identification matrix in the training sample into the encoder for feature extraction to obtain a first vector includes: The data matrix to be encrypted, the random data matrix and the identification matrix in the training sample are input into the convolution module to perform feature extraction and feature fusion to obtain a fused feature vector; The fused feature vector is input into the bidirectional long short-term memory recurrent neural network module for feature extraction to obtain a first vector.

6. The method according to claim 3, characterized in that The original data matrix, the identification matrix and the obtained encrypted data matrix of the training sample are input into the discriminator to obtain the probability matrix output by the discriminator, including: The original data matrix of the training sample and the obtained encrypted data matrix are input into the long short-term memory recurrent neural network module of the discriminator for feature extraction to obtain the third vector and the fourth vector respectively; A probability matrix is ​​obtained according to the third vector and the fourth vector; the identification matrix is ​​used to indicate whether the data corresponding to each element position in the probability matrix is ​​real seismic data or encrypted seismic data.

7. The method according to claim 3, characterized in that The loss function of the deep learning model includes a generator loss function and a discriminator loss function; In the iterative process of training, the loss function of the deep learning model is calculated, and according to the result of the loss function, the value of the first parameter of the deep learning model is adjusted until the loss function reaches a preset convergence condition, thereby obtaining a trained model for improving seismic data density, including: In the training iteration process, the value of the discriminator loss function is calculated until the discriminator loss function reaches a preset convergence condition, and the discriminator is fixed; in the training iteration process after the discriminator is fixed, the value of the generator loss function is calculated until the generator loss function reaches a preset convergence condition, and then the iteration is stopped to obtain a trained seismic data density improvement model; or During the training iteration process, the value of the generator loss function is calculated until the generator loss function reaches the preset convergence condition, and the generator is fixed; during the training iteration process after the generator is fixed, the value of the discriminator loss function is calculated until the discriminator loss function reaches the preset convergence condition, then the iteration is stopped to obtain a trained model for improving seismic data density.

8. The method according to claim 1, characterized in that The method further comprises: A test sample is input into the improved seismic data density model completed in each iteration for testing, and the value of the root mean square error of the improved seismic data density model is calculated. According to the value of the root mean square error, the value of the second parameter of the improved seismic data density model is adjusted until the value of the root mean square error is less than or equal to a preset root mean square error threshold. The test sample is determined based on data of other seismic profiles of the three-dimensional seismic data body.

9. A method for increasing seismic data density, characterized in that: include: Acquire data of a seismic section to be encrypted of a three-dimensional seismic data volume to be increased in density; acquire a corresponding section matrix to be encrypted according to the data of the seismic section to be encrypted; Inserting a number of columns into the section matrix to be encrypted to set the data in the matrix after the inserted columns, generating a data matrix to be encrypted, a random data matrix and an identification matrix corresponding to the section matrix to be encrypted; Inputting the to-be-encrypted data matrix, the random data matrix and the identification matrix into a model for increasing the density of seismic data to obtain a seismic data matrix with increased density; The model for improving the density of seismic data is obtained by the model construction method for improving the density of seismic data as described in any one of claims 1-8.

10. The method according to claim 9, characterized in that The method further comprises: The obtained seismic data matrix with increased density is converted into a seismic profile.

11. A model building device for improving seismic data density, characterized in that: include: A first acquisition module is used to acquire a three-dimensional seismic data volume of a target layer segment, and extract data of multiple seismic sections from the three-dimensional seismic data volume; According to the data of each seismic profile, a corresponding profile matrix is ​​obtained, wherein the columns in the profile matrix represent the original seismic traces in the corresponding seismic profile, and the data in each column represent the real seismic data in the corresponding seismic profile; The first generation module is used to insert a number of columns into each of the profile matrices, and set the data in the profile matrix after the insertion of the number of columns, and respectively generate the original data matrix, the data matrix to be encrypted, the random data matrix and the identification matrix corresponding to each profile matrix as training samples, the inserted number of columns characterize the encrypted seismic traces, and the data on the number of columns characterize the encrypted seismic data; A training module is used to input each of the training samples into a pre-built deep learning model for training. During the iterative process of training, the loss function of the deep learning model is calculated, and according to the result of the loss function, the first parameter of the deep learning model is adjusted until the loss function reaches a preset convergence condition, thereby obtaining a trained model for improving the density of seismic data; the model for improving the density of seismic data outputs a seismic data matrix with improved density.

12. A device for increasing the density of seismic data, characterized in that: include: The second acquisition module is used to acquire data of the seismic section to be encrypted of the three-dimensional seismic data volume to be increased in density; According to the data of the seismic section to be encrypted, obtaining the corresponding section matrix to be encrypted; The second generating module is used to insert a number of columns into the section matrix to be encrypted to set the data in the matrix after the inserted columns, and generate a data matrix to be encrypted, a random data matrix and an identification matrix corresponding to the section matrix to be encrypted; A density improvement module, used for inputting the to-be-encrypted data matrix, the random data matrix and the identification matrix into a seismic data density improvement model to obtain a seismic data matrix after density improvement; The model for improving the density of seismic data is obtained by the model construction method for improving the density of seismic data as described in any one of claims 1-8.

13. A computer storage medium, characterized in that: The computer storage medium stores computer executable instructions, which, when executed by a processor, implement the model building method for improving seismic data density as described in any one of claims 1 to 8 or the method for improving seismic data density as described in claim 9.

14. A terminal device, characterized in that: include: A memory, a processor, and a computer program stored in the memory and executable on the processor, wherein when the processor executes the program, the model building method for improving the density of seismic data as described in any one of claims 1 to 8 or the method for improving the density of seismic data as described in claim 9 is implemented.

Images