CN112305592A - Microseismic signal filtering method and system by utilizing generalized group sparsity - Google Patents

Microseismic signal filtering method and system by utilizing generalized group sparsity Download PDF

Info

Publication number
CN112305592A
CN112305592A CN202011151372.6A CN202011151372A CN112305592A CN 112305592 A CN112305592 A CN 112305592A CN 202011151372 A CN202011151372 A CN 202011151372A CN 112305592 A CN112305592 A CN 112305592A
Authority
CN
China
Prior art keywords
matrix
value
sparse
iteration
frobenus
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
Application number
CN202011151372.6A
Other languages
Chinese (zh)
Inventor
翟明岳
魏会娜
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Guangdong University of Petrochemical Technology
Original Assignee
Guangdong University of Petrochemical Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Guangdong University of Petrochemical Technology filed Critical Guangdong University of Petrochemical Technology
Priority to CN202011151372.6A priority Critical patent/CN112305592A/en
Publication of CN112305592A publication Critical patent/CN112305592A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/288Event detection in seismic signals, e.g. microseismics
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/30Analysis
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/16Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/30Noise handling
    • G01V2210/32Noise reduction
    • G01V2210/324Filtering

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Remote Sensing (AREA)
  • Mathematical Physics (AREA)
  • Environmental & Geological Engineering (AREA)
  • Mathematical Optimization (AREA)
  • Acoustics & Sound (AREA)
  • Geophysics (AREA)
  • Data Mining & Analysis (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Pure & Applied Mathematics (AREA)
  • Geology (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Analysis (AREA)
  • Computational Mathematics (AREA)
  • Software Systems (AREA)
  • Computing Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Databases & Information Systems (AREA)
  • Algebra (AREA)
  • Business, Economics & Management (AREA)
  • Emergency Management (AREA)
  • Circuit For Audible Band Transducer (AREA)

Abstract

The embodiment of the invention discloses a microseismic signal filtering method and a microseismic signal filtering system by utilizing the sparsity of a generalization group, wherein the method comprises the following steps: step 101, acquiring a signal sequence S acquired according to a time sequence; step 102, initializing parameters of an iterative process; step 103, updating the parameters of the iterative process; step 104, judging the difference value of the two adjacent iteration results and ending the iteration; step 105, obtaining a sparse generalization factor; step 106 finds the signal sequence after noise filtering.

Description

Microseismic signal filtering method and system by utilizing generalized group sparsity
Technical Field
The invention relates to the field of geology, in particular to a microseismic signal filtering method and system.
Background
The hydraulic fracturing microseismic monitoring technology is an important new technology developed in the fields of low-permeability reservoir fracturing, reservoir driving, water-drive leading edges and the like in recent years, and is also an important supporting technology for shale gas development. According to the technology, a multistage three-component detector array is arranged in an adjacent well, a microseismic event generated in a target interval of a fractured well in a hydraulic fracturing process is monitored, and the microseismic event is inverted to obtain parameters such as a seismic source position, so that the geometrical shape and the spatial distribution of crack growth in the hydraulic fracturing process are described, the length, the height, the width and the direction of the crack generated by hydraulic fracturing are provided in real time, and the industrial development of shale gas is realized. The hydraulic fracturing microseismic detection is a hotspot and difficulty of scientific research in the field of current shale gas development. From the social and national demand perspective, the development of the research on the aspect of the microseismic monitoring system is very important, and the microseismic monitoring system has great social and economic values.
An important task in microseismic monitoring systems is the localization of microseismic events. The positioning accuracy is the most important factor affecting the application effect of the microseismic monitoring system, and the accuracy of positioning the microseismic event mainly depends on the related factors such as the accuracy of the fluctuation first-arrival (also called first-arrival) reading. But the problem is that the first arrival pick-up is not as simple as it is imagined. The rock fracture form is very complex under the influence of the mining of ground instruments and geological structures, and then microseismic fluctuation with various forms and energy is generated, the form can be dozens or even hundreds, not only are the dominant frequency, the delay, the energy and the like different, but also the waveform form difference near the first arrival position is huge, and the non-uniformity of the waveform characteristics makes the first arrival picking very difficult. Further studies have also shown that the microseismic source mechanism also affects the first arrival point characteristics: most microseismic fluctuations generated by the shearing action of hard rock have large energy, higher main frequency, short time delay and the position of the maximum peak value closely follows the initial first arrival, and the first arrival point of the waves is clear, the jump-off time delay is short, and the waves are easy to pick up; however, most microseismic fluctuations generated by the stretching action have small energy, low main frequency, long delay time, slow take-off and uniform energy distribution, the amplitude of the waves at the first arrival point is small and is easily submerged by interference signals, the characteristic expressions of the first arrival point are inconsistent, and the first arrival pickup is not easy; the microseismic fluctuation generated by soft rock has concentrated energy distribution, fuzzy initial first arrival points, unobvious boundary lines, is obviously different from hard rock, and is difficult to pick up the first arrival. Meanwhile, according to foreign research, it is found that many algorithms want to certainly consider the first arrival wave as a P wave because the P wave velocity is greater than the S wave velocity, but the fact may be more complicated: the first arrivals may be P-waves, S-waves, and even outliers (outliers). According to the study, 41% of the first arrivals are S-waves, and 10% of the first arrivals are caused by outliers. These all present considerable difficulties for first arrival pick-ups.
In addition to the complexity of first arrival point features, first arrival picking faces another greater challenge: microseismic recordings are mass data. For example, approximately 1 million microseismic events were recorded in a test area of month 1 of 2005. Meanwhile, in order to meet production requirements, the microseismic monitoring system needs to continuously record 24 hours a day. Not only is a significant portion of this data a noise and interference caused by human or mechanical activity, independent of microseisms. The literature further classifies noise into three basic types: high frequency (>200Hz) noise, caused by various job related activities; low frequency noise (<10Hz), typically caused by machine activity far from the recording site, and commercial current (50 Hz). In addition, the microseismic signals themselves are not pure, for example, the professor of sinus name in China considers that the microseismic signals include various signals.
Therefore, how to identify microseismic events and pick up first arrivals from mass data is the basis of microseismic data processing. Compared with the prior art, the production method mostly adopts a manual method, wastes time and labor, has poor precision and reliability, cannot ensure the picking quality, and cannot process mass data. The automatic first arrival pickup is one of the solutions, and the automatic first arrival pickup of the micro-seismic fluctuation is one of the key technologies for processing the micro-seismic monitoring data and is also a technical difficulty for realizing the automatic positioning of the micro-seismic source.
Disclosure of Invention
Due to the complex working environment of water conservancy pressure microseismic detection, microseismic signals are seriously interfered by background noise, the background noise presents obvious non-stationarity and non-Gaussian characteristics, a common low-pass filter cannot achieve an ideal filtering effect in the non-stationarity and non-Gaussian noise environment, the non-stationarity non-Gaussian noise is difficult to filter, and the performance of a microseismic event detection algorithm is seriously influenced.
The invention aims to provide a microseismic signal filtering method and a microseismic signal filtering system utilizing generalized group sparsity. The method has good noise filtering performance and is simple in calculation.
In order to achieve the purpose, the invention provides the following scheme:
a method of filtering microseismic signals using generalized group sparsity, comprising:
step 101, acquiring a signal sequence S acquired according to a time sequence;
step 102 iterates overInitializing program parameters, specifically: the iteration process parameters comprise a group sparse transformation matrix W, a sparse coding vector g and an iteration control parameter k; the initialization value of the group sparse transform matrix W is recorded as W0The initialized value of the sparse coding vector g is recorded as g0(ii) a The calculation formula of the initialization value of the iteration process parameter is as follows:
W0=DCT(STS)
g0=SV
k=0
wherein:
DCT(STs) represents a pair matrix STS performs a two-dimensional discrete cosine transform,
v is a matrix STA left eigenvector matrix of S;
step 103, updating the parameters of the iterative process, specifically:
the value of the iterative control parameter k is increased by 1,
Ak=Wk-1gk-1,
Xk=ηλ[Ak],
Wk=vkuk,
gk=[Wk]TXk,
wherein:
a represents a parameter matrix;
Akthe kth step value of the parameter matrix A is represented;
Figure BDA0002741415940000021
Ak(k,: represents A)kThe k row elements of (1);
||Ak(k,:)||Fis represented by AkThe Frobenus pattern of (k,:);
Figure BDA0002741415940000031
is a Frobenus modulus factor;
m0is the signalThe mean value of the sequence S;
σ0is the mean square error of the signal sequence S;
vkis a matrix [ X ]k]TXkA left feature matrix of (a);
ukis a matrix [ X ]k]TXkA right feature matrix of (a);
Xkis the k-th step value of the sparse vector X;
Wk-1representing the k-1 step value of the group sparse transform matrix W;
Wka kth step value representing the group sparse transform matrix W;
gk-1representing the k-1 step value of the sparse coding vector g;
gka k-th step value representing the sparse coding vector g;
step 104, judging a difference value of two adjacent iteration results and ending the iteration, specifically: if the difference value of the two adjacent iteration results meets | | gk-gk-1||F≥λ||gk||FReturning to the step 103 and the step 104 to continue the iterative updating process; otherwise, the iteration process is ended, and the value of the iteration control parameter K at the moment is given to the iteration step number KoTo obtain the optimal group sparse transformation matrix
Figure BDA0002741415940000034
And an optimal sparse coding matrix
Figure BDA0002741415940000033
Wherein:
||gk-gk-1||Fdenotes gk-gk-1The Frobenus pattern of (1) above,
||gk||Fdenotes gkThe Frobenus moustache of (1);
step 105, obtaining a sparse generalization factor, specifically: the coefficient generalization factor is recorded as theta, and the calculation formula is
Figure BDA0002741415940000032
Wherein:
||Xopt||Frepresenting the optimal sparse coding moment XoptThe Frobenus moustache of (1);
step 106, obtaining a signal sequence after noise filtering, specifically: the signal sequence after noise filtering is recorded as SnewThe formula used is:
Snew=[I+θWopt]-1S
wherein: i denotes an identity matrix.
A microseismic signal filtering system utilizing generalized group sparsity comprising:
the module 201 acquires a signal sequence S acquired in time sequence;
the module 202 initializes parameters of an iterative process, specifically: the iteration process parameters comprise a group sparse transformation matrix W, a sparse coding vector g and an iteration control parameter k; the initialization value of the group sparse transform matrix W is recorded as W0The initialized value of the sparse coding vector g is recorded as g0(ii) a The calculation formula of the initialization value of the iteration process parameter is as follows:
W0=DCT(STS)
g0=SV
k=0
wherein:
DCT(STs) represents a pair matrix STS performs a two-dimensional discrete cosine transform,
v is a matrix STA left eigenvector matrix of S;
module 203 updates the iterative process parameters, specifically:
the value of the iterative control parameter k is increased by 1,
Ak=Wk-1gk-1,
Xk=ηλ[Ak],
Wk=vkuk,
gk=[Wk]TXk,
wherein:
a represents a parameter matrix;
Akthe kth step value of the parameter matrix A is represented;
Figure BDA0002741415940000041
Ak(k,: represents A)kThe k row elements of (1);
||Ak(k,:)||Fis represented by AkThe Frobenus pattern of (k,:);
Figure BDA0002741415940000042
is a Frobenus modulus factor;
m0is the mean of the signal sequence S;
σ0is the mean square error of the signal sequence S;
vkis a matrix [ X ]k]TXkA left feature matrix of (a);
ukis a matrix [ X ]k]TXkA right feature matrix of (a);
Xkis the k-th step value of the sparse vector X;
Wk-1representing the k-1 step value of the group sparse transform matrix W;
Wka kth step value representing the group sparse transform matrix W;
gk-1representing the k-1 step value of the sparse coding vector g;
gka k-th step value representing the sparse coding vector g;
the module 204 determines a difference between two adjacent iteration results and ends the iteration, specifically: if the difference value of the two adjacent iteration results meets | | gk-gk-1||F≥λ||gk||FThen returning to the module 203 and the module 204 to continue the iterative updating process; otherwise, the iteration process is ended, and the iteration at the moment is controlledThe value of the parameter K is given to the number of iteration steps KoTo obtain the optimal group sparse transformation matrix
Figure BDA0002741415940000051
And an optimal sparse coding matrix
Figure BDA0002741415940000052
Wherein:
||gk-gk-1||Fdenotes gk-gk-1The Frobenus pattern of (1) above,
||gk||Fdenotes gkThe Frobenus moustache of (1);
the module 205 finds a sparse generalization factor, specifically: the coefficient generalization factor is recorded as theta, and the calculation formula is
Figure BDA0002741415940000053
Wherein:
||Xopt||Frepresenting the optimal sparse coding moment XoptThe Frobenus moustache of (1);
the module 206 calculates a signal sequence after noise filtering, specifically: the signal sequence after noise filtering is recorded as Snew
The solving formula is as follows:
Snew=[I+θWopt]-1S
wherein: i denotes an identity matrix.
According to the specific embodiment provided by the invention, the invention discloses the following technical effects:
due to the complex working environment of water conservancy pressure microseismic detection, microseismic signals are seriously interfered by background noise, the background noise presents obvious non-stationarity and non-Gaussian characteristics, a common low-pass filter cannot achieve an ideal filtering effect in the non-stationarity and non-Gaussian noise environment, the non-stationarity non-Gaussian noise is difficult to filter, and the performance of a microseismic event detection algorithm is seriously influenced.
The invention aims to provide a microseismic signal filtering method and a microseismic signal filtering system utilizing generalized group sparsity. The method has good noise filtering performance and is simple in calculation.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the embodiments will be briefly described below. It is obvious that the drawings in the following description are only some embodiments of the invention, and that for a person skilled in the art, other drawings can be derived from them without inventive effort.
FIG. 1 is a schematic flow diagram of the process of the present invention;
FIG. 2 is a schematic flow chart of the system of the present invention;
FIG. 3 is a flow chart illustrating an embodiment of the present invention.
Detailed Description
The technical solution in the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings in the embodiments of the present invention. It is to be understood that the described embodiments are merely exemplary of the invention, and not restrictive of the full scope of the invention. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
In order to make the aforementioned objects, features and advantages of the present invention comprehensible, embodiments accompanied with figures are described in further detail below.
FIG. 1 is a schematic flow chart of a microseismic signal filtering method using generalized group sparsity
FIG. 1 is a schematic flow chart of a microseismic signal filtering method using generalized group sparsity according to the present invention. As shown in fig. 1, the microseismic signal filtering method using generalized group sparsity specifically includes the following steps:
step 101, acquiring a signal sequence S acquired according to a time sequence;
step 102, initializing parameters of an iterative process, specifically: the iteration process parameters comprise a group sparse transformation matrix W, a sparse coding vector g and an iteration control parameter k; the initialization value of the group sparse transform matrix W is recorded as W0The initialized value of the sparse coding vector g is recorded as g0(ii) a The calculation formula of the initialization value of the iteration process parameter is as follows:
W0=DCT(STS)
g0=SV
k=0
wherein:
DCT(STs) represents a pair matrix STS performs a two-dimensional discrete cosine transform,
v is a matrix STA left eigenvector matrix of S;
step 103, updating the parameters of the iterative process, specifically:
the value of the iterative control parameter k is increased by 1,
Ak=Wk-1gk-1,
Xk=ηλ[Ak],
Wk=vkuk,
gk=[Wk]TXk,
wherein:
a represents a parameter matrix;
Akthe kth step value of the parameter matrix A is represented;
Figure BDA0002741415940000061
Ak(k,: represents A)kThe k row elements of (1);
||Ak(k,:)||Fis represented by AkThe Frobenus pattern of (k,:);
Figure BDA0002741415940000071
is a Frobenus modulus factor;
m0is the mean of the signal sequence S;
σ0is the mean square error of the signal sequence S;
vkis a matrix [ X ]k]TXkA left feature matrix of (a);
ukis a matrix [ X ]k]TXkA right feature matrix of (a);
Xkis the k-th step value of the sparse vector X;
Wk-1representing the k-1 step value of the group sparse transform matrix W;
Wka kth step value representing the group sparse transform matrix W;
gk-1representing the k-1 step value of the sparse coding vector g;
gka k-th step value representing the sparse coding vector g;
step 104, judging a difference value of two adjacent iteration results and ending the iteration, specifically: if the difference value of the two adjacent iteration results meets | | gk-gk-1||F≥λ||gk||FReturning to the step 103 and the step 104 to continue the iterative updating process; otherwise, the iteration process is ended, and the value of the iteration control parameter K at the moment is given to the iteration step number KoTo obtain the optimal group sparse transformation matrix
Figure BDA0002741415940000073
And an optimal sparse coding matrix
Figure BDA0002741415940000074
Wherein:
||gk-gk-1||Fdenotes gk-gk-1The Frobenus pattern of (1) above,
||gk||Fdenotes gkThe Frobenus moustache of (1);
step 105, obtaining a sparse generalization factor, specifically: the coefficient generalization factor is recorded as theta, and the calculation formula is
Figure BDA0002741415940000072
Wherein:
||Xopt||Frepresenting the optimal sparse coding moment XoptThe Frobenus moustache of (1);
step 106, obtaining a signal sequence after noise filtering, specifically: the signal sequence after noise filtering is recorded as Snew
The solving formula is as follows:
Snew=[I+θWopt]-1S
wherein: i denotes an identity matrix.
FIG. 2 is a schematic diagram of a microseismic signal filtering system utilizing generalized group sparsity
FIG. 2 is a schematic diagram of a microseismic signal filtering system utilizing generalized group sparsity according to the present invention. As shown in fig. 2, the microseismic signal filtering system utilizing generalized group sparsity includes the following structure:
the module 201 acquires a signal sequence S acquired in time sequence;
the module 202 initializes parameters of an iterative process, specifically: the iteration process parameters comprise a group sparse transformation matrix W, a sparse coding vector g and an iteration control parameter k; the initialization value of the group sparse transform matrix W is recorded as W0The initialized value of the sparse coding vector g is recorded as g0(ii) a The calculation formula of the initialization value of the iteration process parameter is as follows:
W0=DCT(STS)
g0=SV
k=0
wherein:
DCT(STs) represents a pair matrix STS performs a two-dimensional discrete cosine transform,
v is a matrix STA left eigenvector matrix of S;
module 203 updates the iterative process parameters, specifically:
the value of the iterative control parameter k is increased by 1,
Ak=Wk-1gk-1,
Xk=ηλ[Ak],
Wk=vkuk,
gk=[Wk]TXk,
wherein:
a represents a parameter matrix;
Akthe kth step value of the parameter matrix A is represented;
Figure BDA0002741415940000081
Ak(k,: represents A)kThe k row elements of (1);
||Ak(k,:)||Fis represented by AkThe Frobenus pattern of (k,:);
Figure BDA0002741415940000082
is a Frobenus modulus factor;
m0is the mean of the signal sequence S;
σ0is the mean square error of the signal sequence S;
vkis a matrix [ X ]k]TXkA left feature matrix of (a);
ukis a matrix [ X ]k]TXkA right feature matrix of (a);
Xkis the k-th step value of the sparse vector X;
Wk-1representing the k-1 step value of the group sparse transform matrix W;
Wka kth step value representing the group sparse transform matrix W;
gk-1representing the k-1 step value of the sparse coding vector g;
gka k-th step value representing the sparse coding vector g;
the module 204 determines a difference between two adjacent iteration results and ends the iteration, specifically: if the difference value of the two adjacent iteration results meets | | gk-gk-1||F≥λ||gk||FThen returning to the module 203 and the module 204 to continue the iterative updating process; otherwise, the iteration process is ended, and the value of the iteration control parameter K at the moment is given to the iteration step number KoTo obtain the optimal group sparse transformation matrix
Figure BDA0002741415940000091
And an optimal sparse coding matrix
Figure BDA0002741415940000092
Wherein:
||gk-gk-1||Fdenotes gk-gk-1The Frobenus pattern of (1) above,
||gk||Fdenotes gkThe Frobenus moustache of (1);
the module 205 finds a sparse generalization factor, specifically: the coefficient generalization factor is recorded as theta, and the calculation formula is
Figure BDA0002741415940000093
Wherein:
||Xopt||Frepresenting the optimal sparse coding moment XoptThe Frobenus moustache of (1);
the module 206 calculates a signal sequence after noise filtering, specifically: the signal sequence after noise filtering is recorded as Snew
The solving formula is as follows:
Snew=[I+θWopt]-1S
wherein: i denotes an identity matrix.
The following provides an embodiment for further illustrating the invention
FIG. 3 is a flow chart illustrating an embodiment of the present invention. As shown in fig. 3, the method specifically includes the following steps:
step 301, acquiring a signal sequence S acquired according to a time sequence;
step 302, initializing parameters of an iterative process, specifically: the iteration process parameters comprise a group sparse transformation matrix W, a sparse coding vector g and an iteration control parameter k; the initialization value of the group sparse transform matrix W is recorded as W0The initialized value of the sparse coding vector g is recorded as g0(ii) a The calculation formula of the initialization value of the iteration process parameter is as follows:
W0=DCT(STS)
g0=SV
k=0
wherein:
DCT(STs) represents a pair matrix STS performs a two-dimensional discrete cosine transform,
v is a matrix STA left eigenvector matrix of S;
step 303 updates the iterative process parameters, specifically:
the value of the iterative control parameter k is increased by 1,
Ak=Wk-1gk-1,
Xk=ηλ[Ak],
Wk=vkuk,
gk=[Wk]TXk,
wherein:
a represents a parameter matrix;
Akthe kth step value of the parameter matrix A is represented;
Figure BDA0002741415940000101
Ak(k,: represents A)kThe k row elements of (1);
||Ak(k,:)||Fis represented by AkThe Frobenus pattern of (k,:);
Figure BDA0002741415940000102
is a Frobenus modulus factor;
m0is the mean of the signal sequence S;
σ0is the mean square error of the signal sequence S;
vkis a matrix [ X ]k]TXkA left feature matrix of (a);
ukis a matrix [ X ]k]TXkA right feature matrix of (a);
Xkis the k-th step value of the sparse vector X;
Wk-1representing the k-1 step value of the group sparse transform matrix W;
Wka kth step value representing the group sparse transform matrix W;
gk-1representing the k-1 step value of the sparse coding vector g;
gka k-th step value representing the sparse coding vector g;
step 304, determining a difference between two adjacent iteration results and ending the iteration, specifically: if the difference value of the two adjacent iteration results meets | | gk-gk-1||F≥λ||gk||FReturning to the step 303 and the step 304 to continue the iterative updating process; otherwise, the iteration process is ended, and the value of the iteration control parameter K at the moment is given to the iteration step number KoTo obtain the optimal group sparse transformation matrix
Figure BDA0002741415940000103
And an optimal sparse coding matrix
Figure BDA0002741415940000104
Wherein:
||gk-gk-1||Fdenotes gk-gk-1The Frobenus pattern of (1) above,
||gk||Fdenotes gkThe Frobenus moustache of (1);
step 305, finding a sparse generalization factor, specifically: the coefficient generalization factor is recorded as theta, and the calculation formula is
Figure BDA0002741415940000111
Wherein:
||Xopt||Frepresenting the optimal sparse coding moment XoptThe Frobenus moustache of (1);
step 306, obtaining the signal sequence after noise filtering, specifically: the signal sequence after noise filtering is recorded as Snew
The solving formula is as follows:
Snew=[I+θWopt]-1S
wherein: i denotes an identity matrix.
The embodiments in the present description are described in a progressive manner, each embodiment focuses on differences from other embodiments, and the same and similar parts among the embodiments are referred to each other. For the system disclosed by the embodiment, the description is simple because the system corresponds to the method disclosed by the embodiment, and the relevant part can be referred to the method part for description.
The principles and embodiments of the present invention have been described herein using specific examples, which are provided only to help understand the method and the core concept of the present invention; meanwhile, for a person skilled in the art, according to the idea of the present invention, the specific embodiments and the application range may be changed. In view of the above, the present disclosure should not be construed as limiting the invention.

Claims (2)

1. A microseismic signal filtering method utilizing generalized group sparsity, comprising:
step 101, acquiring a signal sequence S acquired according to a time sequence;
step 102, initializing parameters of an iterative process, specifically: the iteration process parameters comprise a group sparse transformation matrix W, a sparse coding vector g and an iteration control parameter k; the initialization value of the group sparse transform matrix W is recorded as W0The initialized value of the sparse coding vector g is recorded as g0(ii) a The calculation formula of the initialization value of the iteration process parameter is as follows:
W0=DCT(STS)
g0=SV
k=0
wherein:
DCT(STs) represents a pair matrix STS performs a two-dimensional discrete cosine transform,
v is a matrix STA left eigenvector matrix of S;
step 103, updating the parameters of the iterative process, specifically:
the value of the iterative control parameter k is increased by 1,
Ak=Wk-1gk-1,
Xk=ηλ[Ak],
Wk=vkuk,
gk=[Wk]TXk,
wherein:
a represents a parameter matrix;
Akthe kth step value of the parameter matrix A is represented;
Figure FDA0002741415930000011
Ak(k,: represents A)kThe k row elements of (1);
||Ak(k,:)||Fis represented by AkThe Frobenus pattern of (k,:);
Figure FDA0002741415930000012
is a Frobenus modulus factor;
m0is a stand forThe mean value of the signal sequence S;
σ0is the mean square error of the signal sequence S;
vkis a matrix [ X ]k]TXkA left feature matrix of (a);
ukis a matrix [ X ]k]TXkA right feature matrix of (a);
Xkis the k-th step value of the sparse vector X;
Wk-1representing the k-1 step value of the group sparse transform matrix W;
Wka kth step value representing the group sparse transform matrix W;
gk-1representing the k-1 step value of the sparse coding vector g;
gka k-th step value representing the sparse coding vector g;
step 104, judging a difference value of two adjacent iteration results and ending the iteration, specifically: if the difference value of the two adjacent iteration results meets | | gk-gk-1||F≥λ||gk||FReturning to the step 103 and the step 104 to continue the iterative updating process; otherwise, the iteration process is ended, and the value of the iteration control parameter K at the moment is given to the iteration step number KoTo obtain the optimal group sparse transformation matrix
Figure FDA0002741415930000021
And an optimal sparse coding matrix
Figure FDA0002741415930000022
Wherein:
||gk-gk-1||Fdenotes gk-gk-1The Frobenus pattern of (1) above,
||gk||Fdenotes gkThe Frobenus moustache of (1);
step 105, obtaining a sparse generalization factor, specifically: the coefficient generalization factor is recorded as theta, and the calculation formula is
Figure FDA0002741415930000023
Wherein:
||Xopt||Frepresenting the optimal sparse coding moment XoptThe Frobenus moustache of (1);
step 106, obtaining a signal sequence after noise filtering, specifically: the signal sequence after noise filtering is recorded as SnewThe formula used is:
Snew=[I+θWopt]-1S
wherein: i denotes an identity matrix.
2. A microseismic signal filtering system utilizing generalized group sparsity comprising:
the module 201 acquires a signal sequence S acquired in time sequence;
the module 202 initializes parameters of an iterative process, specifically: the iteration process parameters comprise a group sparse transformation matrix W, a sparse coding vector g and an iteration control parameter k; the initialization value of the group sparse transform matrix W is recorded as W0The initialized value of the sparse coding vector g is recorded as g0(ii) a The calculation formula of the initialization value of the iteration process parameter is as follows:
W0=DCT(STS)
g0=SV
k=0
wherein:
DCT(STs) represents a pair matrix STS performs a two-dimensional discrete cosine transform,
v is a matrix STA left eigenvector matrix of S;
module 203 updates the iterative process parameters, specifically:
the value of the iterative control parameter k is increased by 1,
Ak=Wk-1gk-1,
Xk=ηλ[Ak],
Wk=vkuk,
gk=[Wk]TXk,
wherein:
a represents a parameter matrix;
Akthe kth step value of the parameter matrix A is represented;
Figure FDA0002741415930000031
Ak(k,: represents A)kThe k row elements of (1);
||Ak(k,:)||Fis represented by AkThe Frobenus pattern of (k,:);
Figure FDA0002741415930000032
is a Frobenus modulus factor;
m0is the mean of the signal sequence S;
σ0is the mean square error of the signal sequence S;
vkis a matrix [ X ]k]TXkA left feature matrix of (a);
ukis a matrix [ X ]k]TXkA right feature matrix of (a);
Xkis the k-th step value of the sparse vector X;
Wk-1representing the k-1 step value of the group sparse transform matrix W;
Wka kth step value representing the group sparse transform matrix W;
gk-1representing the k-1 step value of the sparse coding vector g;
gka k-th step value representing the sparse coding vector g;
the module 204 determines a difference between two adjacent iteration results and ends the iteration, specifically: if the difference value of the two adjacent iteration results meets | | gk-gk-1||F≥λ||gk||FThen returning to the module 203 and the module 204 to continue the iterative updating process;otherwise, the iteration process is ended, and the value of the iteration control parameter K at the moment is given to the iteration step number KoTo obtain the optimal group sparse transformation matrix
Figure FDA0002741415930000033
And an optimal sparse coding matrix
Figure FDA0002741415930000034
Wherein:
||gk-gk-1||Fdenotes gk-gk-1The Frobenus pattern of (1) above,
||gk||Fdenotes gkThe Frobenus moustache of (1);
the module 205 finds a sparse generalization factor, specifically: the coefficient generalization factor is recorded as theta, and the calculation formula is
Figure FDA0002741415930000035
Wherein:
||Xopt||Frepresenting the optimal sparse coding moment XoptThe Frobenus moustache of (1);
the module 206 calculates a signal sequence after noise filtering, specifically: the signal sequence after noise filtering is recorded as SnewThe formula used is:
Snew=[I+θWopt]-1S
wherein: i denotes an identity matrix.
CN202011151372.6A 2020-10-25 2020-10-25 Microseismic signal filtering method and system by utilizing generalized group sparsity Pending CN112305592A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011151372.6A CN112305592A (en) 2020-10-25 2020-10-25 Microseismic signal filtering method and system by utilizing generalized group sparsity

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011151372.6A CN112305592A (en) 2020-10-25 2020-10-25 Microseismic signal filtering method and system by utilizing generalized group sparsity

Publications (1)

Publication Number Publication Date
CN112305592A true CN112305592A (en) 2021-02-02

Family

ID=74330447

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011151372.6A Pending CN112305592A (en) 2020-10-25 2020-10-25 Microseismic signal filtering method and system by utilizing generalized group sparsity

Country Status (1)

Country Link
CN (1) CN112305592A (en)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108827605A (en) * 2018-03-20 2018-11-16 南京航空航天大学 A kind of mechanical breakdown characteristic automatic extraction method based on improvement sparseness filtering
CN109783739A (en) * 2019-01-23 2019-05-21 北京工业大学 A kind of collaborative filtering recommending method based on the sparse noise reduction self-encoding encoder enhancing of stacking
US20190346522A1 (en) * 2018-05-10 2019-11-14 Siemens Healthcare Gmbh Method of reconstructing magnetic resonance image data
US20200217979A1 (en) * 2019-01-08 2020-07-09 King Fahd University Of Petroleum And Minerals Observation-driven method based on iir wiener filter for microseismic data denoising
CN111679321A (en) * 2020-05-29 2020-09-18 广东石油化工学院 Microseismic signal reconstruction method and system by using generalized gradient

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108827605A (en) * 2018-03-20 2018-11-16 南京航空航天大学 A kind of mechanical breakdown characteristic automatic extraction method based on improvement sparseness filtering
US20190346522A1 (en) * 2018-05-10 2019-11-14 Siemens Healthcare Gmbh Method of reconstructing magnetic resonance image data
US20200217979A1 (en) * 2019-01-08 2020-07-09 King Fahd University Of Petroleum And Minerals Observation-driven method based on iir wiener filter for microseismic data denoising
CN109783739A (en) * 2019-01-23 2019-05-21 北京工业大学 A kind of collaborative filtering recommending method based on the sparse noise reduction self-encoding encoder enhancing of stacking
CN111679321A (en) * 2020-05-29 2020-09-18 广东石油化工学院 Microseismic signal reconstruction method and system by using generalized gradient

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
付丽华;李宏伟;刘智慧;赵浩岚;: "基于多核模型的地震信号高效稀疏分解", 石油地球物理勘探, no. 03 *

Similar Documents

Publication Publication Date Title
CN110146918B (en) Grouping-based microseismic event detection method and system
CN110146920B (en) Microseismic event detection method and system based on relative change of amplitude
CN110794456A (en) Microseismic signal reconstruction method and system by using Gaussian model
CN111596362A (en) Microseismic signal filtering method and system by utilizing Lagrange factor
CN111856563A (en) Microseismic signal reconstruction method and system by using conversion sparsity
CN111679321A (en) Microseismic signal reconstruction method and system by using generalized gradient
CN111596367A (en) Microseismic signal reconstruction method and system based on subspace learning optimization
CN110703319B (en) Microseismic event detection method and system based on Khichin-Einstein theorem
CN110703321B (en) Microseismic event detection method and system using dictionary theory
CN110361779B (en) Microseismic event detection method and system based on chi-square distribution
CN112305592A (en) Microseismic signal filtering method and system by utilizing generalized group sparsity
CN112257560B (en) Microseismic event detection method and system by utilizing cumulative similarity
CN112394403B (en) Microseismic event detection method and system by using edge detection
CN110161560B (en) Method and device for detecting microseismic event
CN110703324A (en) Microseismic event detection method and system represented by random dictionary
CN110146919B (en) Microseismic event detection method and system based on orthogonal projection
CN111596361A (en) Microseismic signal filtering method and system using local limit point
CN110333530B (en) Microseismic event detection method and system
CN112257565B (en) Microseismic event detection method and system using maximum Hull distance
CN111614344A (en) Microseismic signal filtering method and system by using near-end gradient
CN110146921B (en) Microseismic event detection method and system based on Dirac distribution probability
CN110685665A (en) Microseismic event detection method and system based on boundary detection
CN111679323A (en) Microseismic signal filtering method and system based on error vector optimization
CN112285769A (en) Microseismic signal filtering method and system by means of data high-dimensional embedding
CN112180438A (en) Microseismic signal reconstruction method and system by using selection matrix

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
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20210202

WD01 Invention patent application deemed withdrawn after publication