CN112558147B - Polarization analysis method and system for microseism data in well - Google Patents
Polarization analysis method and system for microseism data in well Download PDFInfo
- Publication number
- CN112558147B CN112558147B CN201910914284.8A CN201910914284A CN112558147B CN 112558147 B CN112558147 B CN 112558147B CN 201910914284 A CN201910914284 A CN 201910914284A CN 112558147 B CN112558147 B CN 112558147B
- Authority
- CN
- China
- Prior art keywords
- data
- well
- scale
- component
- trace
- 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.)
- Active
Links
- 230000010287 polarization Effects 0.000 title claims abstract description 100
- 238000004458 analytical method Methods 0.000 title claims abstract description 31
- 238000000034 method Methods 0.000 claims abstract description 27
- 230000001131 transforming effect Effects 0.000 claims abstract description 7
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 238000001914 filtration Methods 0.000 claims description 3
- 238000012545 processing Methods 0.000 abstract description 6
- 238000012544 monitoring process Methods 0.000 description 8
- 238000011161 development Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000000354 decomposition reaction Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 208000035126 Facies Diseases 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 229910052500 inorganic mineral Inorganic materials 0.000 description 1
- 239000011707 mineral Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012163 sequencing technique Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
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. analysis, for interpretation, for correction
- G01V1/288—Event detection in seismic signals, e.g. microseismics
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
Abstract
The invention provides a polarization analysis method and a polarization analysis system for microseism data in a well, and belongs to the field of microseism data processing. The method comprises the following steps: (1) obtaining three-component in-well microseismic data; (2) Transforming each trace of data in the three-component well microseism data into a multi-scale wavelet domain; (3) Calculating the polarization angle of each trace in the three-component well microseism data in each scale in the wavelet domain of each scale; (4) For each trace of microseismic data in a three-component well, determining the polarization angle of the trace according to the polarization angle of each scale. For low signal-to-noise ratio data, the method has better stability in polarization analysis of the microseism data in the well; compared with a polarization analysis method directly utilizing the original data, the method provided by the invention has higher polarization angle precision.
Description
Technical Field
The invention belongs to the field of microseism data processing, and particularly relates to a polarization analysis method and a polarization analysis system for microseism data in a well.
Background
The microseism monitoring technology is an unconventional technology and is applied to the fields of hydraulic fracturing effect evaluation, mineral product monitoring, geological disaster monitoring and the like. The rapid development of unconventional oil gas, especially shale gas development, plays an important role in optimizing a fracturing scheme and developing well pattern deployment, so that the microseism monitoring technology in oil and gas field development is rapidly developed. Microseism monitoring includes in-well microseism monitoring and ground microseism monitoring; microseism monitoring in wells is widely used due to low cost.
Polarization analysis is a key step in the processing of microseismic data in wells. Polarization analysis extracts the spatial polarization parameters of the wave field by studying the polarization characteristics of the wave. The polarization analysis precision not only affects the positioning precision of the microseism event in the well; in addition, the identification and detection of the seismic facies of microseismic events are also affected.
When the signal-to-noise ratio of the microseism data is high enough, the original data can be directly utilized for polarization analysis; however, the signal-to-noise ratio of the microseismic data is generally low, and a large error exists when the original data is directly utilized for polarization analysis, namely the problem that the polarization angle estimation of three-component data is difficult under the low signal-to-noise ratio exists.
Disclosure of Invention
The invention aims to solve the problems in the prior art and provide a polarization analysis method and a polarization analysis system for microseism data in a well, which improve the accuracy of analyzing the polarization angle of three-component data under low signal-to-noise ratio data.
The invention is realized by the following technical scheme:
a polarization analysis method of microseismic data in a well, comprising the following steps:
(1) Acquiring three-component well microseism data;
(2) Transforming each trace of data in the three-component well microseism data into a multi-scale wavelet domain;
(3) Calculating the polarization angle of each trace in the three-component well microseism data in each scale in the wavelet domain of each scale;
(4) For each trace of microseismic data in a three-component well, determining the polarization angle of the trace according to the polarization angle of each scale.
The operation in step (1) includes:
recording microseism data in the well by using a three-component detector;
obtaining direct longitudinal wave data in the in-well microseism data as three-component in-well microseism data;
the direct longitudinal wave data includes a horizontal component x, a horizontal component y, and a vertical component z.
The operation of step (2) comprises:
and transforming each data in the three-component well microseism data into a multi-scale wavelet domain by adopting a Mallat algorithm to obtain transformed data of each scale.
The step (2) further comprises:
and filtering the transformed data of each scale by adopting a low-pass filter.
The operation of step (3) comprises:
in the wavelet domain of each scale, the following operations are performed on the trace j in the three-component well microseismic data:
(31) Calculating a function
wherein ,
x i (t)、y i (t) respectively representing the transformed data of the horizontal component x and the horizontal component y of the channel j in the microseismic data in the three-component well in the scale i;
x i rot (t)、y i rot (t) represents a horizontal component x and a horizontal component y after the rotation angle θ, and θ is a polarization angle;
(32) Will cause E i θ with the largest value of (θ) is taken as the polarization angle of the trace j at the scale iWhere j represents the number of traces in the three-component well microseismic data, j=1, 2,..n, N is the total number of traces; i denotes the number of scales, i=1, 2.
The operation of step (32) includes:
setting a value range and a step length of theta;
starting from the minimum value of θ, every iterationIncreasing 1 step at a time until the maximum value of theta is reached; the operation of step (31) is performed in each iteration to obtain an E i A value of (θ);
after iteration is completed, E corresponding to different theta is obtained i Value of (θ), all E i Sorting the values of (θ), finding the maximum value thereof;
taking θ corresponding to the maximum value as the polarization angle of trace j in the micro-seismic data in the well in the scale i
The operation of step (4) comprises:
(41) For trace j, its polarization angle at dimension i isSorting the sequence from big to small to obtain a sorted sequence;
(42) Selecting a value at a middle position in the ordered sequence as the polarization angle theta of the track j j 。
The operation of step (42) includes:
if M is an odd number, the value at the middle position in the ordered sequence is taken as the polarization angle theta of the track j j ；
If M is even, taking the average value of two values positioned at the middle position in the ordered sequence as the polarization angle theta of the track j j 。
The invention also provides a polarization analysis system of the microseism data in the well, which comprises the following steps:
the data acquisition unit is used for acquiring microseism data in the three-component well;
the conversion unit is connected with the data acquisition unit and is used for converting each trace of data in the microseism data in the three-component well into a multi-scale wavelet domain;
the calculating unit is connected with the converting unit and is used for calculating the polarization angle of each trace in the three-component well microseism data in each scale in the wavelet domain of each scale;
and the polarization angle determining unit is connected with the calculating unit and is used for determining the polarization angle of each trace of the three-component well microseism data according to the polarization angle of each scale.
The present invention also provides a computer-readable storage medium storing at least one program executable by a computer, which when executed by the computer, causes the computer to perform the steps in a polarization analysis method of microseismic data in a well of the present invention.
Compared with the prior art, the invention has the beneficial effects that:
1) For low signal-to-noise ratio data, the method has better stability in polarization analysis of the microseism data in the well;
2) Compared with a polarization analysis method directly utilizing the original data, the method provided by the invention has higher polarization angle precision.
Drawings
FIG. 1 is a block diagram of the steps of the method of the present invention;
2-1, a horizontal component x of microseismic data in a well simulated in the embodiment;
2-2, a horizontal component y of microseismic data in a well simulated in the embodiment;
polarization angle error at different signal-to-noise ratios in the embodiment of fig. 3;
fig. 4 is a block diagram of the system of the present invention.
Detailed Description
The invention is described in further detail below with reference to the attached drawing figures:
when the microseismic effective signal in the well contains noise, the difference between the noise and the signal can be used for processing the microseismic data because the noise and the signal have different manifestations under different wavelet scales. The original microseism data is transformed into wavelet domains with different scales by utilizing discrete wavelet transformation for processing, so that the interference of noise can be effectively avoided. The method calculates the polarization angle in the multi-scale wavelet domain, increases the useful information amount in the processing process of estimating the polarization angle, and improves the precision of the polarization angle under the condition of low signal to noise ratio.
As shown in fig. 1, the method of the present invention comprises:
the method comprises the steps of obtaining microseism data in a three-component well;
specifically, the three-component in-well microseismic data obtained in the step (1) is direct longitudinal wave (P-wave) data of the microseismic data (the "in-well microseismic data" in the following steps all refer to direct longitudinal wave data, and the direct longitudinal wave data includes a horizontal component x, a horizontal component y and a vertical component z); wherein, the microseismic data in each well is recorded by a three-component detector, which comprises a horizontal component x, a horizontal component y and a vertical component z.
Transforming the data of each trace in the microseism data in the well into a multi-scale wavelet domain;
specifically, three-component data for each trace in the microseismic data in the well is transformed into the multi-scale wavelet domain. Initialization parameters of the multi-scale wavelet transform, such as wavelet base type and maximum wavelet decomposition scale M, etc., may be determined as desired. One skilled in the art can select any suitable wavelet base type (e.g., up to Bei Xixiao waves), maximum decomposition dimension M, etc., as desired.
Transforming the data of each trace in the microseismic data in the well into a multi-scale wavelet domain by adopting the existing Mallat algorithm; recursive computation of microseismic data s by low pass filter H 0 Approximation { s at different levels 1 ,s 2 K }, as follows:
t represents time, s i (l) Is the transformed data of the micro-seismic data in the well in the scale i (namely, the data of each trace in the micro-seismic data in the well is transformed into the data obtained by a multiscale wavelet domain by adopting a Mallat algorithm); h (l) represents an impulse response sequence of the low-pass filter H, and its specific expression form can be determined according to the type of the wavelet base (the specific determination method is wavelet transformation theory, which is not described herein in detail), and l represents an operator length of the filter.
The purpose of the filtering process with the low pass filter H is: 1) Noise influence is reduced, 2) a plurality of scales are calculated simultaneously, and result stability is improved.
Third, in the wavelet domain of each scale, calculating the polarization angle of each trace in the microseismic data in the well at the scaleWhere j represents the number of traces in the microseismic data in the well, j=1, 2,..n, N is the total number of traces in the data; i represents the number of the multi-scale mesoscale, i=1, 2,..m, M is the total number of scales in the multi-scale;
in the wavelet domain of each scale, calculating a horizontal polarization angle for a trace j in the microseismic data in the well specifically comprises the following steps of;
first, a function is calculated
wherein ,
x i (t)、y i (t) represents the transformed data of the horizontal component x and the horizontal component y of any trace j in the microseismic data in the well (namely, the data obtained by the step (2)) in the scale i; x is x i rot (t)、y i rot (t) represents two components processed by the above formula (i.e., rotated by an angle θ), θ being an angle between the polarization direction and the horizontal component x, i.e., a polarization angle, which is a quantity required to be determined for polarization analysis of the horizontal component.
Finally, will cause E i θ with the largest value of (θ) is taken as the polarization angle of each trace in the microseismic data in the well at scale iThe method comprises the following steps:
setting the value range of theta as (-pi, pi) and step length, and calculating the step length according to specific numerical valuesThe calculation accuracy is determined, the values of the horizontal component x and the horizontal component y are determined, different theta are sequentially substituted into the formula, and the E corresponding to different theta is calculated by iteration (in the iteration process, the minimum value of theta and the theta value changed by the step length are utilized, namely, 1 step length is added once per iteration from the minimum value of theta until the maximum value of theta is reached.) i (θ) all E i Ordering (θ), find E i A maximum value of (θ), and θ corresponding to the maximum value is taken as a polarization angle of a trace j in the micro-seismic data in the well in a scale i
Fourth, aiming at each trace of microseismic data in the well, the polarization angles of wavelet domains based on various scales are adoptedThe final polarization angle of the track is determined as follows:
all polarization angles for trace j in the borehole microseismic dataSequencing the array (namely, the array formed by M polarization angles corresponding to M scales) from large to small to obtain a sequenced sequence;
selecting a value at a middle position in the ordered sequence as the polarization angle theta of the track j j Specific: if M is an odd number, the value at the middle position in the ordered sequence is taken as the polarization angle theta of the track j j The method comprises the steps of carrying out a first treatment on the surface of the If M is even, taking the average value of two values positioned at the middle position in the ordered sequence as the polarization angle theta of the track j j 。
The invention also provides a polarization analysis system of the microseismic data in the well, as shown in fig. 4, comprising:
a data acquisition unit 10 for acquiring three-component well microseismic data;
the conversion unit 20 is connected with the data acquisition unit 10 and is used for converting each trace of data in the three-component well microseism data into a multi-scale wavelet domain;
a calculating unit 30, connected to the converting unit 20, for calculating, in the wavelet domain of each scale, the polarization angle of each trace of the microseismic data in the three-component well at the scale;
a polarization angle determining unit 40, connected to the calculating unit 30, for determining, for each trace of the microseismic data in the well, a polarization angle of the trace based on the polarization angles of the respective dimensions.
Examples of polarization analysis of microseismic data in wells using the method of the present invention are as follows:
in this embodiment, the method is validated using simulated microseismic data.
FIGS. 2-1 and 2-2 are simulated microseismic data in a well, FIG. 2-1 being the horizontal component x and FIG. 2-2 being the horizontal component y; the analog data is added with Gaussian white noise, and the signal-to-noise ratio is 10; the true polarization angle of the channel model data was 50 degrees. FIG. 3 is a graph of the polarization angle error of analog data under different signal-to-noise ratios, comparing the polarization angle error of the method of the present invention with the polarization angle error of the method using only the original data; it can be seen from fig. 3 that the polarization angle error calculated by the method of the present invention is smaller. The model test shows that the polarization analysis method has higher precision.
The foregoing technical solution is only one embodiment of the present invention, and various modifications and variations can be easily made by those skilled in the art based on the application methods and principles disclosed in the present invention, not limited to the methods described in the foregoing specific embodiments of the present invention, so that the foregoing description is only preferred and not in a limiting sense.
Claims (7)
1. A polarization analysis method for microseismic data in a well, which is characterized by comprising the following steps of: the method comprises the following steps:
(1) Acquiring three-component well microseism data;
(2) Transforming each trace of data in the three-component well microseism data into a multi-scale wavelet domain;
(3) Calculating the polarization angle of each trace in the three-component well microseism data in each scale in the wavelet domain of each scale;
(4) Determining the polarization angle of each trace according to the polarization angle of each scale aiming at each trace of microseism data in the three-component well;
the operation of step (3) comprises:
in the wavelet domain of each scale, the following operations are performed on the trace j in the three-component well microseismic data:
(31) Calculating a function
wherein ,
x i (t)、y i (t) respectively representing the transformed data of the horizontal component x and the horizontal component y of the channel j in the microseismic data in the three-component well in the scale i; t represents time;
x i rot (t)、y i rot (t) represents a horizontal component x and a horizontal component y after the rotation angle θ, and θ is a polarization angle;
(32) Will cause E i θ with the largest value of (θ) is taken as the polarization angle of the trace j at the scale iWhere j represents the number of traces in the three-component well microseismic data, j=1, 2,..n, N is the total number of traces; i represents the number of scales, i=1, 2,..m, M is the total number of scales;
the operation of step (32) includes:
setting a value range and a step length of theta;
starting from the minimum value of theta, increasing 1 step length every iteration until the maximum value of theta is reached; the operation of step (31) is performed in each iteration to obtain an E i Value of (θ)；
After iteration is completed, E corresponding to different theta is obtained i Value of (θ), all E i Sorting the values of (θ), finding the maximum value thereof;
taking θ corresponding to the maximum value as the polarization angle of trace j in the micro-seismic data in the well in the scale i
The operation of step (4) comprises:
(41) For trace j, its polarization angle at dimension i isSorting the sequence from big to small to obtain a sorted sequence;
(42) Selecting a value at a middle position in the ordered sequence as the polarization angle theta of the track j j 。
2. The method of polarization analysis of microseismic data in a well of claim 1 wherein: the operation in step (1) includes:
recording microseism data in the well by using a three-component detector;
obtaining direct longitudinal wave data in the in-well microseism data as three-component in-well microseism data;
the direct longitudinal wave data includes a horizontal component x, a horizontal component y, and a vertical component z.
3. The method of polarization analysis of microseismic data in a well of claim 1 wherein: the operation of step (2) comprises:
and transforming each data in the three-component well microseism data into a multi-scale wavelet domain by adopting a Mallat algorithm to obtain transformed data of each scale.
4. A method of polarization analysis of microseismic data in a well according to claim 3, wherein: the step (2) further comprises:
and filtering the transformed data of each scale by adopting a low-pass filter.
5. The method of polarization analysis of microseismic data in a well of claim 1 wherein: the operation of step (42) includes:
if M is an odd number, the value at the middle position in the ordered sequence is taken as the polarization angle theta of the track j j ；
If M is even, taking the average value of two values positioned at the middle position in the ordered sequence as the polarization angle theta of the track j j 。
6. A polarization analysis system for microseismic data in a well, characterized by: the system comprises:
the data acquisition unit is used for acquiring microseism data in the three-component well;
the conversion unit is connected with the data acquisition unit and is used for converting each trace of data in the microseism data in the three-component well into a multi-scale wavelet domain;
the calculating unit is connected with the converting unit and is used for calculating the polarization angle of each trace in the three-component well microseism data in each scale in the wavelet domain of each scale;
and the polarization angle determining unit is connected with the calculating unit and is used for determining the polarization angle of each trace according to the polarization angle of each scale aiming at each trace of the microseismic data in the three-component well.
The operations performed by the computing unit include:
in the wavelet domain of each scale, the following operations are performed on the trace j in the three-component well microseismic data:
(31) Calculating a function
wherein ,
x i (t)、y i (t) respectively representing the transformed data of the horizontal component x and the horizontal component y of the channel j in the microseismic data in the three-component well in the scale i; t represents time;
x i rot (t)、y i rot (t) represents a horizontal component x and a horizontal component y after the rotation angle θ, and θ is a polarization angle;
(32) Will cause E i θ with the largest value of (θ) is taken as the polarization angle of the trace j at the scale iWhere j represents the number of traces in the three-component well microseismic data, j=1, 2,..n, N is the total number of traces; i represents the number of scales, i=1, 2,..m, M is the total number of scales;
the operation of step (32) includes:
setting a value range and a step length of theta;
starting from the minimum value of theta, increasing 1 step length every iteration until the maximum value of theta is reached; the operation of step (31) is performed in each iteration to obtain an E i A value of (θ);
after iteration is completed, E corresponding to different theta is obtained i Value of (θ), all E i Sorting the values of (θ), finding the maximum value thereof;
taking θ corresponding to the maximum value as the polarization angle of trace j in the micro-seismic data in the well in the scale i
The polarization angle determining unit performs operations including:
(41) For trace j, its polarization angle at dimension i isOrdering the array from big to smallObtaining a sequenced sequence;
(42) Selecting a value at a middle position in the ordered sequence as the polarization angle theta of the track j j 。
7. A computer-readable storage medium, characterized by: the computer readable storage medium stores at least one program executable by a computer, which when executed by the computer, causes the computer to perform the steps in the polarization analysis method of microseismic data in a well according to any one of claims 1 to 5.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910914284.8A CN112558147B (en) | 2019-09-25 | 2019-09-25 | Polarization analysis method and system for microseism data in well |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910914284.8A CN112558147B (en) | 2019-09-25 | 2019-09-25 | Polarization analysis method and system for microseism data in well |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112558147A CN112558147A (en) | 2021-03-26 |
| CN112558147B true CN112558147B (en) | 2023-10-10 |
Family
ID=75029576
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910914284.8A Active CN112558147B (en) | 2019-09-25 | 2019-09-25 | Polarization analysis method and system for microseism data in well |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112558147B (en) |
Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CA2580355A1 (en) * | 2004-09-17 | 2006-03-23 | Julian Drew | Microseismic event detection and location by continuous map migration |
| WO2010116236A2 (en) * | 2009-04-08 | 2010-10-14 | Schlumberger Technology B.V. | Methods and systems for microseismic mapping |
| WO2014108843A1 (en) * | 2013-01-09 | 2014-07-17 | Westerngeco Seismic Holdings Limited | Noise attenuation of multicomponent microseismic data |
| CN104216005A (en) * | 2013-06-03 | 2014-12-17 | 中国石油化工股份有限公司 | Downhole microseismic three-component rotating method |
| CN104280772A (en) * | 2013-07-12 | 2015-01-14 | 中国石油化工股份有限公司 | Recognition method for microseism phase in well |
| CN106168675A (en) * | 2015-05-18 | 2016-11-30 | 中国石油化工股份有限公司 | Borehole microseismic P/S ripple event recognition method and device |
| CN106547024A (en) * | 2015-09-16 | 2017-03-29 | 中国石油化工股份有限公司 | For the residual static correction amount estimation method and device of microseism perforation data |
| RU2620785C1 (en) * | 2016-06-14 | 2017-05-29 | Общество с ограниченной ответственностью "Макросейс" | Method of determining the location of the microseismic focus |
| CN107367754A (en) * | 2016-05-11 | 2017-11-21 | 中国石油化工股份有限公司 | Microseism first arrival recognition methods and device based on three-component polarization gradient |
| CN108107475A (en) * | 2018-03-05 | 2018-06-01 | 吉林大学 | A kind of borehole microseismic denoising method based on experience wavelet transformation and multi-threshold function |
| CN109655901A (en) * | 2017-10-11 | 2019-04-19 | 中国石油化工股份有限公司 | A kind of adaptive angle of polarization calculation method of frequency domain and system |
-
2019
- 2019-09-25 CN CN201910914284.8A patent/CN112558147B/en active Active
Patent Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CA2580355A1 (en) * | 2004-09-17 | 2006-03-23 | Julian Drew | Microseismic event detection and location by continuous map migration |
| WO2010116236A2 (en) * | 2009-04-08 | 2010-10-14 | Schlumberger Technology B.V. | Methods and systems for microseismic mapping |
| WO2014108843A1 (en) * | 2013-01-09 | 2014-07-17 | Westerngeco Seismic Holdings Limited | Noise attenuation of multicomponent microseismic data |
| CN104216005A (en) * | 2013-06-03 | 2014-12-17 | 中国石油化工股份有限公司 | Downhole microseismic three-component rotating method |
| CN104280772A (en) * | 2013-07-12 | 2015-01-14 | 中国石油化工股份有限公司 | Recognition method for microseism phase in well |
| CN106168675A (en) * | 2015-05-18 | 2016-11-30 | 中国石油化工股份有限公司 | Borehole microseismic P/S ripple event recognition method and device |
| CN106547024A (en) * | 2015-09-16 | 2017-03-29 | 中国石油化工股份有限公司 | For the residual static correction amount estimation method and device of microseism perforation data |
| CN107367754A (en) * | 2016-05-11 | 2017-11-21 | 中国石油化工股份有限公司 | Microseism first arrival recognition methods and device based on three-component polarization gradient |
| RU2620785C1 (en) * | 2016-06-14 | 2017-05-29 | Общество с ограниченной ответственностью "Макросейс" | Method of determining the location of the microseismic focus |
| CN109655901A (en) * | 2017-10-11 | 2019-04-19 | 中国石油化工股份有限公司 | A kind of adaptive angle of polarization calculation method of frequency domain and system |
| CN108107475A (en) * | 2018-03-05 | 2018-06-01 | 吉林大学 | A kind of borehole microseismic denoising method based on experience wavelet transformation and multi-threshold function |
Non-Patent Citations (3)
| Title |
|---|
| 《地面微地震监测资料静校正方法研究》;姜宇东 等;《石油物探》;第52卷(第2期);第132-140页 * |
| 《微地震信号的偏振―位置对比法震相分离技术》;朱卫星 等;《石油地球物理勘探》;第44卷(第4期);第425-429页 * |
| 《微地震资料频域相干―时间域偏振滤波方法》;宋维琪 等;《石油地球物理勘探》;第43卷(第2期);第161-167页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112558147A (en) | 2021-03-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106405651B (en) | Full waveform inversion initial velocity model construction method based on logging matching | |
| WO2016008105A1 (en) | Post-stack wave impedance inversion method based on cauchy distribution | |
| US20170108604A1 (en) | Denoising seismic data | |
| CN109143331B (en) | Seismic wavelet extraction method | |
| CN108020863A (en) | A kind of thin and interbedded reservoir porosity prediction method based on earthquake parity function | |
| CN107179550B (en) | A kind of seismic signal zero phase deconvolution method of data-driven | |
| CN113093272A (en) | Time domain full waveform inversion method based on convolutional coding | |
| US20110166842A1 (en) | Layer stripping method | |
| CN108957554B (en) | Seismic inversion method in geophysical exploration | |
| CN111427096A (en) | Data quality evaluation and filtering processing method for full tensor gravity gradiometer | |
| CN112213782B (en) | Processing method and device for sub-phase seismic data and server | |
| CN113687440A (en) | Milnaciaceae vicky cycle-based ancient water depth quantitative recovery method and storage medium | |
| CN112558147B (en) | Polarization analysis method and system for microseism data in well | |
| CN112649887A (en) | Method and device for quantitatively dividing sequence stratum based on drilling data | |
| CN106291675A (en) | A kind of geological data reconstructing method based on base tracer technique | |
| CN116068644A (en) | Method for improving resolution and noise reduction of seismic data by using generation countermeasure network | |
| CN110749923A (en) | Deconvolution method for improving resolution based on norm equation | |
| CN106842322B (en) | Carbon dioxide flooding monitoring earthquake time difference correction method | |
| Zhang | Ensemble methods of data assimilation in porous media flow for non-Gaussian prior probability density | |
| CN115879020A (en) | Reservoir classification depicting method and device, electronic equipment and storage medium | |
| CN116774293B (en) | Method, system, electronic equipment and medium for automatically picking up same phase shaft | |
| CN113640880B (en) | Sandstone type uranium ore seismic exploration method and device | |
| CN114609668B (en) | High-quality reservoir identification method, device and equipment based on scattering transformation and neural network and storage medium | |
| CN113900140B (en) | Seismic data optimization method and device based on space-time combination | |
| CN113109875B (en) | Inversion method of carbonate rock reservoir under full waveform velocity field constraint |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |