CN114706128A - Microseism data noise suppression method and system, storage medium and seismic information processing equipment - Google Patents
Microseism data noise suppression method and system, storage medium and seismic information processing equipment Download PDFInfo
- Publication number
- CN114706128A CN114706128A CN202210278729.XA CN202210278729A CN114706128A CN 114706128 A CN114706128 A CN 114706128A CN 202210278729 A CN202210278729 A CN 202210278729A CN 114706128 A CN114706128 A CN 114706128A
- Authority
- CN
- China
- Prior art keywords
- model
- noise
- micro
- data
- seismic
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 79
- 230000001629 suppression Effects 0.000 title claims abstract description 61
- 230000010365 information processing Effects 0.000 title claims abstract description 10
- 238000012549 training Methods 0.000 claims abstract description 55
- 238000012360 testing method Methods 0.000 claims abstract description 35
- 238000012795 verification Methods 0.000 claims abstract description 31
- 230000007246 mechanism Effects 0.000 claims abstract description 28
- 230000002123 temporal effect Effects 0.000 claims description 35
- 238000000926 separation method Methods 0.000 claims description 28
- 239000013598 vector Substances 0.000 claims description 28
- 230000008569 process Effects 0.000 claims description 19
- 230000004913 activation Effects 0.000 claims description 15
- 230000006870 function Effects 0.000 claims description 15
- 238000010606 normalization Methods 0.000 claims description 15
- 238000005070 sampling Methods 0.000 claims description 11
- 238000004590 computer program Methods 0.000 claims description 8
- 238000013507 mapping Methods 0.000 claims description 7
- 238000003825 pressing Methods 0.000 claims description 7
- 238000000605 extraction Methods 0.000 claims description 2
- 238000005516 engineering process Methods 0.000 abstract description 4
- 238000012545 processing Methods 0.000 abstract description 3
- 230000008878 coupling Effects 0.000 description 8
- 238000010168 coupling process Methods 0.000 description 8
- 238000005859 coupling reaction Methods 0.000 description 8
- 230000000694 effects Effects 0.000 description 8
- 238000004422 calculation algorithm Methods 0.000 description 6
- 238000011161 development Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 6
- 238000001914 filtration Methods 0.000 description 4
- 238000012800 visualization Methods 0.000 description 4
- 230000008030 elimination Effects 0.000 description 2
- 238000003379 elimination reaction Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 238000000547 structure data Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 238000013527 convolutional neural network Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000002372 labelling Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
- G01V1/364—Seismic filtering
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
A method and a system for suppressing the noise of micro-seismic data, a storage medium and a seismic information processing device belong to the micro-seismic data processing technology and aim to solve the problems that the existing method can not eliminate all noise and can also lose effective signals. The technical points are as follows: collecting microseism data of a target block by using a detector; dividing the acquired micro-seismic data into a training set, a verification set and a test set according to a proportion, inputting the training set data into a time attention mechanism mutual information generation countermeasure network model for learning and outputting a micro-seismic noise suppression pre-training model; inputting verification set data into a trained model, and verifying model parameters by using signal-to-noise ratio to obtain an optimal micro-seismic noise pressure system model hyper-parameter combination; testing the optimal hyper-parameter combination by using the test set data and obtaining a final micro-seismic noise pressure modeling model; and inputting the field micro-seismic data into a final micro-seismic noise pressure control model to obtain de-noised micro-seismic data, and realizing micro-seismic data noise control.
Description
Technical Field
The invention relates to a method and a system for suppressing microseism data noise, a storage medium and seismic information processing equipment, and belongs to the technical field of microseism data processing.
Background
In recent years, the potential of conventional oil and gas exploitation in China is gradually reduced, but the oil and gas demand caused by rapid development of national industry is continuously improved, and the external dependence of oil in China is as high as 73.5% by 2020. The unconventional oil gas in China can be used for detecting about 1075 hundred million tons of reserves, the development potential is huge, and the development of unconventional oil gas exploration and development is a necessary way for ensuring the energy safety in China. The microseism monitoring technology is a key technology for evaluating the hydraulic fracturing effect and acquiring parameters such as the azimuth, the length and the height of a crack, and is an important means for guaranteeing the efficient development of unconventional oil and gas.
In the practical engineering, the ground micro-seismic detector not only collects a large amount of random noise but also contains strong interference caused by complex ground noise sources (such as 50Hz industrial alternating current interference, drilling interference and the like) in the data acquisition process, so that a micro-seismic effective signal is easily submerged, and the identification of a micro-seismic effective event and the accuracy and reliability of seismic source positioning are seriously influenced. Therefore, noise suppression is a core element of microseismic data preprocessing. However, in establishing a noise suppression model, we find that the conventional noise reduction method usually takes subjective expert experience guidance as a main guide when distinguishing effective signals from noise, that is, the noise signal is considered to exist only in the high-frequency part of the whole signal, and then the high-frequency noise is eliminated by using low-pass filtering. However, noise signals contained in the micro-seismic signals have coupling, and the subjective identification of the noise frequency and the corresponding method for elimination cannot eliminate all the noise and also lose effective signals. In addition, the effective event characteristic distribution of the microseisms presents discretization and only exists in local positions of the time series. The microseism signals collected by different detectors also present different characteristics due to geological changes, so that loss of effective events is easily caused in the noise suppression process of microseism data. Therefore, finding new ways to suppress noise in microseismic data without losing significant events has become paramount in unconventional oil and gas exploration development.
Disclosure of Invention
In order to solve the technical problems mentioned in the background art: the traditional denoising method generally takes subjective expert experience guidance as a main guide when distinguishing effective signals and noise, namely, the noise signals are only determined to exist in the high-frequency part of the whole signals, and then the low-pass filtering is utilized to eliminate the high-frequency noise. However, noise signals contained in the micro-seismic signals have coupling, and the noise frequency is subjectively determined and eliminated by adopting a corresponding method, so that not only can all the noise be eliminated, but also effective signals can be lost, and therefore the micro-seismic data noise suppression method and system, the storage medium and the seismic information processing equipment are provided, wherein the micro-seismic data noise suppression method and system are based on the time attention mechanism mutual information generation countermeasure network.
The technical scheme adopted by the invention for solving the technical problems is as follows: a microseism data noise suppression method for generating a countermeasure network based on time attention mechanism mutual information is characterized by comprising the following steps:
step 1: collecting the micro-seismic data of the target block by using a detector to obtain the micro-seismic data X with noise and the non-noise micro-seismic data
Step 2: dividing the microseism data acquired in the step 1 into a training set, a verification set and a test set according to the ratio of 8:1: 1; the training set is represented asWherein xkThe K training trace set samples are K training samples in total; the verification set is represented asWherein xaThe verification trace set sample is the a-th verification trace set sample, and the A-th verification trace set sample is the A-th verification trace set sample; the test set is represented asWherein xbB test trace set samples are the B test trace set samples; the sampling frequency of the training set, the verification set and the test set is f;
and step 3: extraction of the samples collected in step 1Structural semantic feature set ofWherein, cd∈RH×1The d-th micro-seismic structure semantic vector is H;
and 4, step 4: sampling a random noise distribution (such as Gaussian distribution) with a sampling frequency f, and acquiring a random noise vector setWherein z ism∈RN×1The m-th random noise vector is composed of N sampling points;
and 5: subjecting the product obtained in step 3With the product obtained in step 4Dimension combination is carried out to obtain feature vectorsWherein v isu∈R(H+N)×1Is the u-th merged feature vector;
step 6: the method comprises the following steps of utilizing a countermeasure network generated based on time attention mechanism mutual information to learn training set data and output a trained microseism noise suppression pre-training model, wherein the specific path is as follows:
the microseism noise pressing model adopted in the step comprises a generation model, a discrimination model and a separation model;
the generating model is composed of a time attention module and a generating module, wherein the time attention module is specifically defined as follows:
wherein the input to the temporal attention module is v obtained via step 5u,j(ii) a The output of the temporal attention module is the weight αu;σtRepresents the tanh activation function;a jth channel representing a convolution kernel having a size of Z' x 1;is a bias term;
the generation module consists of a Dense layer, a Batch Normalization layer, a ReLU layer, a Deconv layer and a Tanh layer; the discrimination model is composed of a time attention module and a discrimination module, wherein the time attention module is specifically defined as follows:
wherein the input to the temporal attention module is obtained via step 2The output of the temporal attention module is the weight βk;σsRepresents the softmax activation function;a jth channel representing a convolution kernel having a convolution kernel size of F' × 1;is a bias term;
the discrimination module consists of a Dense layer, a Batch Normalization layer, a LeakyReLU layer, a Conv layer and a Softmax layer;
the separation model is composed of a time attention module and a separation module, wherein the time attention module is defined as:
wherein the input to the temporal attention module is G (z, c) obtained via the generative modelh,j(ii) a The output of the temporal attention module is the weight γh;σsigRepresenting a sigmoid activation function;a jth channel representing a convolution kernel having a size of Y' × 1;is a bias term;
the separation module consists of a Dense layer, a Batch Normalization layer, a LeakyReLU layer, a Conv layer and a Sigmoid layer;
subjecting the product obtained in step 5Inputting the data into a generating model to obtain a generated noise-free micro seismic data set; the generated noise-free micro-seismic data set and the data obtained in the step 2Respectively inputting the data into a discrimination model, training the discrimination model and reversely transmitting the gradient to a generation model; inputting the generated noise-free micro-seismic data set into a separation model to extract structural semantic features, calculating mutual information between the generated noise-free micro-seismic data and the extracted structural semantic features, training the separation model, and reversely transmitting the gradient to the generation model; finally, when the maximum iteration times are reached, a pre-training model for suppressing the microseism noise can be obtained;
and 7: verification set data to be obtained via step 2Inputting the data into the noise suppression pre-training model obtained in the step 6, verifying the performance of the trained noise suppression pre-training model by utilizing signal-to-noise ratio, and obtaining the optimal hyper-parameter combination of the micro-seismic noise suppression model;
and 8: utilizing test set data obtained via step 2Testing the optimal hyper-parameter combination of the microseism noise pressing model obtained in the step 7, and obtaining a final microseism data noise pressing model;
and step 9: and (4) inputting the field micro-seismic data into the final micro-seismic noise suppression model obtained in the step (8) to obtain the de-noised micro-seismic data, so as to suppress the noise of the micro-seismic data.
Further, the generative model training process in step 6 is performed according to the following path:
2.1 the generative model is composed of a time attention module and a generation module, wherein the time attention module is defined as:
where the input to the temporal attention module is v obtained via step 5u,j(ii) a The output of the temporal attention module is the weight αu;σtRepresents the tanh activation function;a jth channel representing a convolution kernel having a size of Z' x 1;is the bias term.
2.2 v to be obtained via step 5u,jAnd alpha obtained via equation (1)uIs input into formula (2) to obtain an increaseStrong feature mapping vattSpecifically, the following are defined:
wherein v isattFor enhanced feature mapping of merged feature vectors, vuFor the u-th merged feature vector, αuThe temporal attention weight of the u-th merged feature vector.
2.3 v to be obtained via equation (2)attInputting the data into a generation module to obtain a generated noise-free micro-seismic data setWherein, G (z, c)h∈RS×1 is the noise-free microseism data generated at the h-th; the generating module is composed of a Dense layer, a Batch Normalization layer, a ReLU layer, a Deconv layer and a Tanh layer.
Further, the discriminant model training process in step 6 is performed according to the following path:
3.1 the discriminant model is composed of a time attention module and a discriminant module, wherein the time attention module is defined as:
wherein the input to the temporal attention module is obtained via step 2The output of the temporal attention module is the weight βk;σsRepresents the softmax activation function;a jth channel representing a convolution kernel having a convolution kernel size of F' × 1;is offset byAn item.
3.2 obtaining via step 2With β obtained via equation (3)kInput into formula (4) to obtain an enhanced feature map xattSpecifically, the following are defined:
wherein xattIs composed ofEnhanced feature mapping of (2);for the kth training set sample, βkA temporal attention weight for the kth training set sample;
3.3 x to be obtained via equation (4)attInput into a discrimination module to obtainProbability of belonging to a true sample. The discrimination module is composed of a sense layer, a Batch Normalization layer, a LeakyReLU layer, a Conv layer and a Softmax layer.
Further, the separation model training process in step 6 is performed according to the following path:
4.1 the separation model is composed of a time attention module and a separation module, wherein the time attention module is defined as:
wherein the input to the temporal attention module is G (z, c) obtained via the generative modelh,j(ii) a The output of the temporal attention module is the weight γh;σsigRepresenting a sigmoid activation function;a jth channel representing a convolution kernel having a size of Y' × 1;is the bias term.
4.2 to be obtained via feature 2.3And gamma obtained via the formula (5)hInput into equation (6) to obtain an enhanced feature map qattSpecifically, the following are defined:
wherein q isattEnhanced feature mapping for G (z, c); g (z, c)hFor the h-th generated noiseless microseismic data, gammahTime attention weight for the h-th generated noise-free microseismic data.
4.3 q to be obtained via equation (6)attInputting into a separation module to obtainStructural semantic vector ofWhereinFor the w-th extracted structural semantic feature vector, and calculate G (z, c) and cqMutual information of (2); the separation module is composed of a sense layer, a Batch Normalization layer, a LeakyReLU layer, a Conv layer and a Sigmoid layer.
Further, the model verification process in step 7 is performed according to the following path:
verification set data to be obtained by step 2Inputting the data into the noise suppression pre-training model obtained in the step 6, verifying the performance of the trained noise suppression pre-training model by using the signal-to-noise ratio, and obtaining the optimal hyper-parameter combination of the micro-seismic noise suppression model.
Further, the model testing process in step 8 is performed according to the following path:
test set data to be obtained by step 2And 7, testing the optimal hyper-parameter combination of the micro-seismic noise suppression model obtained in the step 7, and obtaining a final micro-seismic data noise suppression model.
Further, the noise suppression process in step 9 is performed according to the following path: and (4) inputting the field micro-seismic data into the final micro-seismic noise suppression model obtained in the step (8) to obtain the de-noised micro-seismic data, so as to suppress the noise of the micro-seismic data.
A micro-seismic data noise suppression system for generating a countermeasure network based on time attention mechanism mutual information is provided with a program module corresponding to the steps of the technical scheme, and the steps in the micro-seismic data noise suppression method for generating the countermeasure network based on the time attention mechanism mutual information are executed during running.
A computer readable storage medium having stored thereon a computer program configured to, when invoked by a processor, perform the steps of the method for noise suppressing microseismic data based on time attention mechanism mutual information generation countermeasure networks described above.
A seismic information processing apparatus comprising at least one processor, and a memory communicatively coupled to the at least one processor, wherein the memory stores instructions executable by the at least one processor, the instructions being executable by the at least one processor to enable the at least one processor to perform a method of noise suppressing microseismic data based on time attention mechanism mutual information generation countermeasure network as described above.
The invention has the following beneficial technical effects:
the microseism noise suppression method overcomes the defects that the coupling noise cannot be completely removed and the loss of microseism effective events is easily caused by the traditional filtering method. The traditional denoising method generally takes subjective expert experience guidance as a main guide when distinguishing effective signals and noise, namely, the noise signals are only determined to exist in the high-frequency part of the whole signals, and then the low-pass filtering is utilized to eliminate the high-frequency noise. However, noise signals contained in the micro-seismic signals have coupling, and the subjective identification of the noise frequency and the corresponding method for elimination cannot eliminate all the noise and also lose effective signals. In addition, the effective event characteristic distribution of the microseisms presents discretization and only exists in local positions of the time series. The microseism signals collected by different detectors also present different characteristics due to geological changes, so that loss of effective events is easily caused in the noise suppression process of microseism data.
The method comprises the following steps: collecting microseism data of a target block by using a detector; dividing the acquired micro-seismic data into a training set, a verification set and a test set according to a certain proportion, inputting the training set data into time attention mechanism mutual information to generate a confrontation network model for learning, and outputting a micro-seismic noise suppression pre-training model; inputting verification set data into a trained model, verifying model parameters by utilizing signal-to-noise ratio, and obtaining an optimal micro-seismic noise pressure system model hyper-parameter combination; testing the optimal hyper-parameter combination by using the test set data and obtaining a final micro-seismic noise pressure modeling model; and inputting the field micro-seismic data into a final micro-seismic noise pressure control model to obtain de-noised micro-seismic data, and realizing micro-seismic data noise control.
The invention considers that the noise contained in the microseism data has the characteristics of strong coupling and large energy, namely the noise characteristic and the microseism data characteristic are coded in a complex and disordered state in a data space, so that the ideal denoising effect is difficult to achieve by the traditional method. Therefore, the method provided by the invention can be used for learning the noise distribution in the micro-seismic data by introducing the potential noise vector and maximizing the mutual information between the potential noise vector distribution and the generated micro-seismic data so as to achieve the purpose of removing the noise. Finally, the effective event characteristic distribution in the microseism monitoring data is considered to be discretized and only exists in the local position of the time sequence. Therefore, the invention adopts a time attention mechanism to amplify the characteristics of the effective events, avoids the phenomenon that the effective events are lost in the denoising process, and further obtains better denoising effect. The invention provides a generation countermeasure network based on mutual information, which solves the problem that strong coupling noise is difficult to remove in microseism noise suppression, introduces a time attention mechanism to amplify effective events in microseism data, and avoids the phenomenon that the effective events are lost in the denoising process. A set of systematic and scientific microseism data noise suppression method is provided, and the method has operability and practical value.
Therefore, the method can effectively solve the problems and amplify the microseism effective events while realizing microseism noise suppression.
The effectiveness of the method of the invention is reflected in the summary of the invention:
1. the method is built by adopting a deep convolutional neural network, has excellent nonlinear characterization capability, and improves the precision of microseism noise suppression on the network structure level.
2. The method adopts a time attention mechanism to amplify the local characteristics of the micro-seismic data time sequence effective events, can avoid the loss of effective signals in the micro-seismic data in the process of realizing noise suppression, and improves the precision of micro-seismic noise suppression on the characteristic level of the data.
3. The method adopts the calculation mutual information to strengthen the correlation between the effective characteristics of the micro-seismic data and the micro-seismic signals, can effectively suppress the coupling noise contained in the micro-seismic data, and improves the precision of suppressing the micro-seismic noise on the algorithm level.
In addition, the method belongs to an unsupervised learning method based on data driving, dependence on manual labeling is eliminated, and more manpower and material resources are reduced. On-site micro-seismic data can be directly subjected to noise suppression through a trained model, a micro-seismic data processing technology which is easier to operate is provided for on-site constructors, the on-site working efficiency is improved, and the best effect is achieved.
Drawings
FIG. 1 is a block flow diagram of a microseismic noise suppression model (method and system) according to an embodiment of the present invention;
FIG. 2 is a diagram of a temporal attention mechanism mutual information generation countermeasure network core, according to an embodiment of the invention;
FIGS. 3 and 4 are diagrams of examples of noise suppressed visualizations of microseismic data according to embodiments of the invention, FIG. 3 is a diagram of examples of noise suppressed visualizations of microseismic data (noisy microseismic noise data), FIG. 4 is a diagram of examples of noise suppressed visualizations of microseismic data (denoised microseismic noise data); FIGS. 5 and 6 are diagrams of another example of noise-suppressed visualization of microseismic data (noisy microseismic noise data);
in fig. 4 and 6, the four lines after denoising represent the in-phase axes (the in-phase axes are the lines of extreme values of the same vibration phase of each trace on the seismic record), and represent the layering of the stratum; the two groups of the two images have different directions of the same phase axis and different straight potentials, and represent different geological structures. The provided figure shows that the invention can be used for compressing a plurality of different geological structure data, and the practicability of the algorithm of the invention can be better embodied.
Detailed Description
The invention will be further described with reference to the accompanying drawings in which:
in order to make those skilled in the art better understand the technical solutions in the present application, the technical solutions in the embodiments of the present application will be clearly and completely described below with reference to the drawings in the embodiments of the present application.
Fig. 1 is a flow chart of microseismic noise suppression in the technical solution of the present invention, and the whole flow chart is completely implemented by Python language programming (of course, other programming languages may be used to implement the flow chart). According to the figure 1, firstly, on-site micro-seismic noise-free data and noise-containing data are collected, then the data are divided into a training set, a verification set and a test set according to a specific proportion, and the training set data are input into time attention mechanism mutual information to generate a pre-training model for learning in a countermeasure network model and outputting micro-seismic noise suppression; inputting verification set data into a trained model, verifying model parameters by utilizing signal-to-noise ratio, and obtaining an optimal micro-seismic noise pressure system model hyper-parameter combination; testing the optimal hyper-parameter combination by using the test set data, and obtaining a final micro-seismic noise pressure modeling model; and inputting the field micro-seismic data into a final micro-seismic noise pressure making model to obtain de-noised micro-seismic data, thereby realizing micro-seismic data noise suppression.
In the embodiment, the validity of the method is verified by taking field microseism data as an example.
Step 1: and acquiring the micro-seismic data of the target block by using a detector to obtain a micro-seismic data signal A with noise and a micro-seismic data signal B without noise.
Step 2: and (3) dividing the microseism data acquired in the step (1) into a training set, a verification set and a test set according to the ratio of 8:1: 1. The training set is represented asWherein x iskIs the k training set sample; the verification set is represented asWherein x isaIs the a-th verification set sample; the test set is represented asWherein x isbIs the b test set sample. And the lengths of the sampling points of the training set, the verification set and the test set are all S.
And step 3: extracting the structural semantic feature set of the signal B acquired in the step 1Wherein, cd∈RH×1The length of the d micro-seismic structure semantic vector is H.
And 4, step 4: sampling a random noise distribution (e.g., gaussian distribution) and obtaining a set of random noise vectorsWherein z ism∈RN×1And the m-th random noise vector is composed of N sampling points.
And 5: subjecting the product obtained in step 3With the product obtained in step 4Dimension combination is carried out to obtain a combined feature vectorWherein v isu∈R(H+N)×1Is the u-th merged feature vector.
And 6: referring to fig. 2, a countermeasure network is generated by using time attention mechanism mutual information to learn training set data and output a trained microseism noise suppression pre-training model, and the specific path is as follows:
the microseism noise pressing model adopted in the step comprises a generation model, a discrimination model and a separation model.
The generating model is composed of a time attention module and a generating module, wherein the time attention module is specifically defined as:
wherein the input to the temporal attention module is v obtained via step 5u,j(ii) a The output of the temporal attention module is the weight αu;σtRepresents the tanh activation function;a jth channel representing a convolution kernel having a size of Z' x 1;is the bias term.
The generation module is composed of a sense layer, a Batch Normalization layer, a ReLU layer, a Deconv layer and a Tanh layer, and belongs to the field of the prior art.
The discrimination model is composed of a time attention module and a discrimination module, wherein the time attention module is specifically defined as follows:
wherein the input to the temporal attention module is obtained via step 2The output of the temporal attention module is the weight βk;σsRepresents the softmax activation function;a jth channel representing a convolution kernel having a convolution kernel size of F' × 1;is the bias term.
The discrimination module is composed of a Dense layer, a Batch Normalization layer, a LeakyReLU layer, a Conv layer and a Softmax layer.
The separation model is composed of a time attention module and a separation module, wherein the time attention module is defined as follows:
wherein, whenThe input to the inter-attention module is G (z, c) obtained via the generative modelh,j(ii) a The output of the temporal attention module is the weight γh;σsigRepresenting a sigmoid activation function;a jth channel representing a convolution kernel having a size of Y' × 1;is the bias term.
The separation module is composed of a Dense layer, a Batch Normalization layer, a LeakyReLU layer, a Conv layer and a Sigmoid layer.
FIG. 2 illustrates a microseismic noise pressure modeling model: subjecting the product obtained in step 5Inputting the data into a generating model to obtain a generated noise-free micro seismic data set; collecting the generated noise-free micro-seismic data and the data obtained in step 2Respectively inputting the data into a discrimination model, training the discrimination model and reversely transmitting the gradient to a generation model; the generated noise-free micro-seismic data set is input into a separation model to extract structural semantic features, mutual information of the generated noise-free micro-seismic data and the extracted structural semantic features is calculated, the separation model is trained, and the gradients are reversely transmitted to the generation model, so that micro-seismic can be enhanced, the correlation between effective features of the micro-seismic data and micro-seismic signals can be enhanced through the method, and coupling noise in the micro-seismic data can be effectively suppressed.
And 7: verification set data to be obtained via step 2Inputting the noise suppression pre-training model obtained in the step 6, verifying the performance of the trained noise suppression pre-training model by using the signal-to-noise ratio, and obtainingAnd (3) carrying out optimal hyper-parameter combination on the microseism noise suppression model.
And step 8: utilizing test set data obtained via step 2And 7, testing the optimal hyper-parameter combination of the micro-seismic noise suppression model obtained in the step 7, and obtaining a final micro-seismic data noise suppression model.
TABLE 1 comparison of different denoising model effects
And step 9: and (4) inputting the field micro-seismic data into the final micro-seismic noise suppression model obtained in the step (8) to obtain the de-noised micro-seismic data, so as to suppress the noise of the micro-seismic data. The results of the noise suppression experiment for the microseismic data are shown in fig. 3 and 4 and fig. 5 and 6.
The provided figure shows that the method can be used for suppressing a plurality of different geological structure data, and the practicability of the algorithm can be better shown. The method adopts a time attention mechanism to amplify the characteristics of the effective events, avoids the phenomenon that the effective events are lost in the denoising process, further obtains a better denoising effect, and has good universality.
The method provided by the invention verifies the claimed technical effect through simulation experiments and practical application.
The algorithm (method) provided by the invention is the bottom technical kernel of the invention, and various products can be derived based on the algorithm.
The algorithm (method) provided by the invention utilizes a program language to develop a micro-seismic data noise suppression system for generating a confrontation network based on the mutual information of the time attention mechanism, the system is provided with a program module corresponding to the steps of the technical scheme, and the steps in the micro-seismic data noise suppression method for generating the confrontation network based on the mutual information of the time attention mechanism are executed during running.
A computer program of the developed system (software) is stored on a computer readable storage medium, the computer program being configured to implement the steps of the above-described method of noise-suppressing microseismic data for generating a counterpoise network based on time attention mechanism mutual information when invoked by a processor. I.e. the invention may be embodied on a carrier, to form a computer program product.
A seismic information processing apparatus comprising at least one processor and a memory communicatively coupled to the at least one processor, wherein the memory stores instructions executable by the at least one processor, the instructions being executable by the at least one processor to enable the at least one processor to perform a method of noise suppressing microseismic material based on time attention mechanism mutual information generation countermeasure network as described above. The seismic information processing equipment is used as a terminal product applied by the invention.
Various implementations of the systems and techniques described here can be realized in digital electronic circuitry, integrated circuitry, application specific ASICs (application specific integrated circuits), computer hardware, firmware, software, and/or combinations thereof. These various embodiments may include: implemented in one or more computer programs that are executable and/or interpretable on a programmable system including at least one programmable processor, which may be special or general purpose, receiving data and instructions from, and transmitting data and instructions to, a storage system, at least one input device, and at least one output device.
The computer programs (also known as programs, software applications, or code) in the present invention include machine instructions for a programmable processor, and may be implemented using high-level procedural and/or object-oriented programming languages, and/or assembly/machine languages. As used herein, the terms "machine-readable medium" and "computer-readable medium" refer to any computer program product, apparatus, and/or device (e.g., magnetic discs, optical disks, memory, Programmable Logic Devices (PLDs)) used to provide machine instructions and/or data to a programmable processor, including a machine-readable medium that receives machine instructions as a machine-readable signal. The term "machine-readable signal" refers to any signal used to provide machine instructions and/or data to a programmable processor.
It should be understood that various forms of the flows shown above may be used, with steps reordered, added, or deleted. For example, the steps described in the present application may be executed in parallel, sequentially, or in different orders, and are within the scope of the present invention as long as the desired results of the technical solutions disclosed in the present application can be achieved.
Claims (10)
1. A microseism data noise suppression method for generating a countermeasure network based on time attention mechanism mutual information is characterized by comprising the following steps:
step 1: collecting the micro-seismic data of the target block by using a detector to obtain the micro-seismic data X with noise and the non-noise micro-seismic data
And 2, step: dividing the microseism data acquired in the step 1 into a training set, a verification set and a test set according to the ratio of 8:1: 1; the training set is represented asWherein xkFor the kth training trace set sample, K training samples are shared; the verification set is represented asWherein xaThe verification trace set sample is the a-th verification trace set sample, and the A-th verification trace set sample is the A-th verification trace set sample; the test set is represented asWherein x isbB test trace set samples are the B test trace set samples; the sampling frequency of the training set, the verification set and the test set is f;
and step 3: extraction of the samples collected in step 1Structural semantic feature set ofWherein, cd∈RH×1The d-th micro-seismic structure semantic vector is H;
and 4, step 4: sampling random noise distribution with sampling frequency f, and acquiring random noise vector setWherein z ism∈RN×1The m-th random noise vector is composed of N sampling points;
and 5: subjecting the product obtained in step 3And obtained in step 4Dimension combination is carried out to obtain feature vectorsWherein v isu∈R(H+N)×1Merging the feature vectors into the u-th merged feature vector;
step 6: the method comprises the following steps of utilizing a time attention mechanism-based mutual information to generate a confrontation network to learn training set data and output a trained microseism noise suppression pre-training model, wherein the specific path is as follows:
the microseism noise pressing model adopted in the step comprises a generation model, a discrimination model and a separation model;
the generating model is composed of a time attention module and a generating module, wherein the time attention module is specifically defined as:
wherein the input to the temporal attention module is v obtained via step 5u,j(ii) a The output of the temporal attention module is the weight αu;σtRepresents the tanh activation function;a jth channel representing a convolution kernel having a size of Z' x 1;is a bias term;
the generation module consists of a Dense layer, a Batch Normalization layer, a ReLU layer, a Deconv layer and a Tanh layer;
the discrimination model is composed of a time attention module and a discrimination module, wherein the time attention module is specifically defined as follows:
wherein the input to the temporal attention module is obtained via step 2The output of the temporal attention module is the weight βk;σsRepresents the softmax activation function;representing the jth channel of the convolution kernel, wherein the size of the convolution kernel is F' multiplied by 1;is a bias term;
the discrimination module consists of a Dense layer, a Batch Normalization layer, a LeakyReLU layer, a Conv layer and a Softmax layer;
the separation model is composed of a time attention module and a separation module, wherein the time attention module is defined as:
wherein the input to the temporal attention module is G (z, c) obtained via the generative modelh,j(ii) a The output of the temporal attention module is the weight γh;σsigRepresenting a sigmoid activation function;a jth channel representing a convolution kernel having a size of Y' × 1;is a bias term;
the separation module consists of a Dense layer, a Batch Normalization layer, a LeakyReLU layer, a Conv layer and a Sigmoid layer;
subjecting the product obtained in step 5Inputting the data into a generating model to obtain a generated noise-free micro seismic data set; the generated noise-free micro-seismic data set and the data obtained in the step 2Respectively inputting the data into a discrimination model, training the discrimination model and reversely transmitting the gradient to a generation model; inputting the generated noiseless micro-seismic data set into a separation model to extract structural semantic features, calculating mutual information between the generated noiseless micro-seismic data and the extracted structural semantic features, and carrying out separation modelingTraining is carried out, and the gradient is reversely transmitted to a generated model; finally, when the maximum iteration times are reached, a pre-training model for suppressing the microseism noise can be obtained;
and 7: verification set data to be obtained via step 2Inputting the data into the noise suppression pre-training model obtained in the step 6, verifying the performance of the trained noise suppression pre-training model by utilizing signal-to-noise ratio, and obtaining the optimal hyper-parameter combination of the micro-seismic noise suppression model;
and 8: utilizing test set data obtained via step 2Testing the optimal hyper-parameter combination of the microseism noise pressing model obtained in the step 7, and obtaining a final microseism data noise pressing model;
and step 9: and (4) inputting the field micro-seismic data into the final micro-seismic noise suppression model obtained in the step (8) to obtain the de-noised micro-seismic data, so as to suppress the noise of the micro-seismic data.
2. The method of claim 1, wherein the method comprises: the training process of the generative model in the step 6 is performed according to the following path:
2.1 the generative model is composed of a time attention module and a generation module, wherein the time attention module is defined as:
where the input to the temporal attention module is v obtained via step 5u,j(ii) a The output of the temporal attention module is the weight αu;σtRepresents the tanh activation function;a jth channel representing a convolution kernel having a size of Z' x 1;is the bias term.
2.2 v to be obtained via step 5u,jAnd alpha obtained via equation (1)uInputting the data into formula (2) to obtain an enhanced feature map vattSpecifically, the following are defined:
wherein v isattFor enhanced feature mapping of merged feature vectors, vuFor the u-th merged feature vector, αuThe temporal attention weight of the u-th merged feature vector.
2.3 v to be obtained via equation (2)attInputting the data into a generation module to obtain a generated noise-free micro-seismic data setWherein, G (z, c)h∈RS×1Generating noiseless microseismic data for the h-th time; the generating module is composed of a Dense layer, a Batch Normalization layer, a ReLU layer, a Deconv layer and a Tanh layer.
3. The method for suppressing the noise of the microseismic data based on the time attention mechanism mutual information generation countermeasure network as claimed in claim 1 or 2, wherein: the discriminant model training process in the step 6 is performed according to the following path:
3.1 the discriminant model is composed of a time attention module and a discriminant module, wherein the time attention module is defined as:
wherein the input to the temporal attention module is obtained via step 2The output of the temporal attention module is the weight βk;σsRepresents the softmax activation function;a jth channel representing a convolution kernel having a convolution kernel size of F' × 1;is the bias term.
3.2 obtaining via step 2With β obtained via equation (3)kInput into formula (4) to obtain an enhanced feature map xattSpecifically, the following are defined:
wherein x isattIs composed ofEnhanced feature mapping of (2);for the kth training set sample, βkA temporal attention weight for the kth training set sample;
4. The method of claim 3, wherein the method comprises: the training process of the separation model in the step 6 is carried out according to the following path:
4.1 the separation model is composed of a time attention module and a separation module, wherein the time attention module is defined as:
wherein the input to the temporal attention module is G (z, c) obtained via the generative modelh,j(ii) a The output of the temporal attention module is the weight γh;σsigRepresenting a sigmoid activation function;representing the j channel of the convolution kernel, wherein the size of the convolution kernel is Y' × 1;is a bias term.
4.2 to be obtained via feature 2.3And gamma obtained via the formula (5)hInput into equation (6) to obtain an enhanced feature map qattSpecifically, the following are defined:
wherein q isattEnhanced feature mapping for G (z, c); g (z, c)hNoise-free micro for h-th generationSeismic data, gammahTime attention weight for the h-th generated noise-free microseismic data.
4.3 q to be obtained via equation (6)attInputting into a separation module to obtainStructural semantic vector ofWhereinFor the w-th extracted structural semantic feature vector, and calculate G (z, c) and cqMutual information of (2); the separation module is composed of a sense layer, a Batch Normalization layer, a LeakyReLU layer, a Conv layer and a Sigmoid layer.
5. The method of claim 4, wherein the method comprises: the model verification process in step 7 is performed according to the following path:
verification set data to be obtained by step 2Inputting the data into the noise suppression pre-training model obtained in the step 6, verifying the performance of the trained noise suppression pre-training model by using the signal-to-noise ratio, and obtaining the optimal hyper-parameter combination of the micro-seismic noise suppression model.
6. The method of claim 5 for noise suppressing microseismic data based on time attention mechanism mutual information generation countermeasure network, wherein the method comprises the following steps: the model testing process in step 8 is performed according to the following path:
7. The method of claim 6, wherein the method comprises: the noise suppression process in step 9 is performed according to the following path:
and (4) inputting the field micro-seismic data into the final micro-seismic noise suppression model obtained in the step (8) to obtain the de-noised micro-seismic data, so as to suppress the noise of the micro-seismic data.
8. A microseism data noise suppression system for generating a countermeasure network based on time attention mechanism mutual information is characterized in that: the system has program modules corresponding to the steps of any one of claims 1 to 7, and is operable to perform the steps of the method for noise suppressing microseismic data based on time attention mechanism mutual information generation countermeasure network.
9. A computer-readable storage medium characterized by: the computer readable storage medium stores a computer program configured to, when invoked by a processor, implement the steps of the method of noise-suppressing microseismic data based on temporal attention mechanism mutual information generating a counterpoise network of any of claims 1-7.
10. A seismic information processing apparatus characterized by: the seismic information processing apparatus comprises at least one processor and a memory communicatively connected to the at least one processor, wherein the memory stores instructions executable by the at least one processor, the instructions being executable by the at least one processor to enable the at least one processor to perform a method of noise suppressing a microseismic data based on time attention mechanism mutual information generation countermeasure network as claimed in any one of claims 1-7.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210278729.XA CN114706128A (en) | 2022-03-21 | 2022-03-21 | Microseism data noise suppression method and system, storage medium and seismic information processing equipment |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210278729.XA CN114706128A (en) | 2022-03-21 | 2022-03-21 | Microseism data noise suppression method and system, storage medium and seismic information processing equipment |
Publications (1)
Publication Number | Publication Date |
---|---|
CN114706128A true CN114706128A (en) | 2022-07-05 |
Family
ID=82169003
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210278729.XA Pending CN114706128A (en) | 2022-03-21 | 2022-03-21 | Microseism data noise suppression method and system, storage medium and seismic information processing equipment |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114706128A (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116660982A (en) * | 2023-08-02 | 2023-08-29 | 东北石油大学三亚海洋油气研究院 | Full waveform inversion method based on attention convolution neural network |
-
2022
- 2022-03-21 CN CN202210278729.XA patent/CN114706128A/en active Pending
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116660982A (en) * | 2023-08-02 | 2023-08-29 | 东北石油大学三亚海洋油气研究院 | Full waveform inversion method based on attention convolution neural network |
CN116660982B (en) * | 2023-08-02 | 2023-09-29 | 东北石油大学三亚海洋油气研究院 | Full waveform inversion method based on attention convolution neural network |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Wang et al. | Blind source extraction of acoustic emission signals for rail cracks based on ensemble empirical mode decomposition and constrained independent component analysis | |
CN111723329B (en) | Seismic phase feature recognition waveform inversion method based on full convolution neural network | |
CN111950503B (en) | Aviation transient electromagnetic data processing method and device and computing equipment | |
CN109613610B (en) | Automatic picking method for microseism signal arrival time difference | |
CN108710777B (en) | Diversified anomaly detection identification method based on multi-convolution self-coding neural network | |
CN104599292A (en) | Noise-resistant moving target detection algorithm based on low rank matrix | |
CN113687433B (en) | Bi-LSTM-based magnetotelluric signal denoising method and system | |
CN105116442A (en) | Lithologic oil-gas reservoir weak-reflection seismic signal reconstruction method | |
CN110632662A (en) | Algorithm for automatically identifying microseism signals by using DCNN-inclusion network | |
CN111507155A (en) | U-Net + + and UDA combined microseism effective signal first-arrival pickup method and device | |
CN113822201B (en) | Deep learning method for underwater object shape recognition based on flow field velocity component time course | |
CN111948711B (en) | Method and system for extracting low-frequency part of seismic data by deep learning method | |
CN113962244A (en) | Rayleigh wave seismic data noise removal method, storage medium and electronic device | |
CN115291281A (en) | Real-time micro-earthquake magnitude calculation method and device based on deep learning | |
CN114706128A (en) | Microseism data noise suppression method and system, storage medium and seismic information processing equipment | |
RU2412454C2 (en) | Method to process seismic data using discrete wavelet transform | |
CN112464721A (en) | Automatic microseism event identification method and device | |
CN117056680A (en) | Data noise reduction and signal detection method, device and system and storage medium | |
CN116049637A (en) | Gas-liquid two-phase flow pattern recognition method based on flow noise signals | |
Liu et al. | Maximum correntropy criterion-based blind deconvolution and its application for bearing fault detection | |
CN117076858B (en) | Deep learning-based low-frequency geomagnetic strong interference suppression method and system | |
CN115690772A (en) | Gravity fault automatic identification method based on deep learning | |
CN111380680A (en) | Check valve fault feature extraction method based on improved permutation entropy | |
CN115828069A (en) | End-to-end magnetic anomaly signal noise reduction method based on deep learning | |
CN115563480A (en) | Gear fault identification method for screening octave geometric modal decomposition based on kurtosis ratio coefficient |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination |