CN111158045A - Reservoir transformation microseism event scattered point clustering analysis method and system - Google Patents
Reservoir transformation microseism event scattered point clustering analysis method and system Download PDFInfo
- Publication number
- CN111158045A CN111158045A CN202010012257.4A CN202010012257A CN111158045A CN 111158045 A CN111158045 A CN 111158045A CN 202010012257 A CN202010012257 A CN 202010012257A CN 111158045 A CN111158045 A CN 111158045A
- Authority
- CN
- China
- Prior art keywords
- event
- microseism
- oil reservoir
- error ellipse
- 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.)
- Granted
Links
- 238000004458 analytical method Methods 0.000 title claims abstract description 21
- 230000009466 transformation Effects 0.000 title claims abstract description 19
- 238000000034 method Methods 0.000 claims abstract description 21
- 239000011159 matrix material Substances 0.000 claims description 18
- 238000004364 calculation method Methods 0.000 claims description 16
- 238000007621 cluster analysis Methods 0.000 claims description 14
- 238000004422 calculation algorithm Methods 0.000 claims description 10
- 239000006185 dispersion Substances 0.000 claims description 6
- 230000006872 improvement Effects 0.000 claims description 2
- 238000005457 optimization Methods 0.000 abstract description 7
- 230000000694 effects Effects 0.000 abstract description 6
- 230000008569 process Effects 0.000 abstract description 2
- 239000003921 oil Substances 0.000 description 35
- 238000012544 monitoring process Methods 0.000 description 8
- 239000007789 gas Substances 0.000 description 5
- 230000008859 change Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 239000007924 injection Substances 0.000 description 2
- 238000002347 injection Methods 0.000 description 2
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 238000011084 recovery Methods 0.000 description 2
- 230000007480 spreading Effects 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 239000000295 fuel oil Substances 0.000 description 1
- 239000003345 natural gas Substances 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 239000003209 petroleum derivative Substances 0.000 description 1
- 238000010223 real-time analysis Methods 0.000 description 1
- 238000002407 reforming Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 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/288—Event detection in seismic signals, e.g. microseismics
Landscapes
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Environmental & Geological Engineering (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Emergency Management (AREA)
- Business, Economics & Management (AREA)
- Acoustics & Sound (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention provides a scattered point clustering analysis method and a scattered point clustering analysis system for micro-seismic events in oil reservoir reconstruction, wherein the method comprises the following steps: step 1, acquiring micro-seismic event scatter points of oil reservoir reconstruction in real time; step 2, calculating error ellipses of all microseism events; step 3, performing weighting operation according to the moment magnitude of the microseism event; step 4, calculating the geometric center of the microseism event within the error ellipse range; and 5, minimizing the objective function through continuous iteration, and accurately describing and finely depicting the artificial fracture for oil reservoir reconstruction. The method and the system for analyzing the scattered point clustering of the micro-seismic events for oil reservoir transformation can effectively evaluate the oil reservoir transformation effect and guide the optimization and perfection of the transformation process, and have an important effect on the efficient exploration and development of unconventional oil and gas resources.
Description
Technical Field
The invention relates to the technical field of micro-seismic monitoring for petroleum and natural gas reservoir transformation, in particular to a scattered point clustering analysis method and a scattered point clustering analysis system for micro-seismic events of reservoir transformation.
Background
The oil reservoir transformation operations such as hydraulic fracturing, water injection, gas injection, heavy oil thermal recovery and the like can induce micro-seismic events, and the micro-seismic activities have important influence on the exploration and development of unconventional oil and gas fields. In order to effectively evaluate and optimize oil and gas yield increasing operation measures, the oil reservoir transformation effect needs to be rapidly and timely mastered, so that the process technology and the transformation parameters are optimized, and micro-seismic monitoring is one of the most effective technologies for evaluating the oil reservoir transformation effect. At present, due to the problems of longitudinal wave (P wave) or transverse wave (S wave) first arrival time error, uncertainty of a velocity model, source vector angle error, unsatisfactory space-time distribution of a detector and the like of a microseism event, the reliability of a microseism monitoring result is low, the error is large, the problem of divergence of microseism event points exists, the framework structure of an artificial fracture for oil deposit reconstruction is difficult to effectively extract, and the popularization and the application of a microseism monitoring technology in unconventional oil and gas exploration and development are restricted. Therefore, a novel method and a novel system for analyzing scattered point clustering of microseism events in oil reservoir transformation are invented, and the technical problems are solved.
Disclosure of Invention
The invention aims to provide a scattered point clustering analysis method and a scattered point clustering analysis system for micro-seismic events of oil reservoir reconstruction, which can quickly acquire an artificial fracture skeleton of oil reservoir reconstruction through a scattered point clustering analysis method and clearly display dynamic distribution and real-time change conditions of the artificial fracture skeleton of oil reservoir reconstruction.
The object of the invention can be achieved by the following technical measures: the scattered point cluster analysis method for the oil reservoir reconstruction micro-seismic events comprises the following steps: step 1, acquiring micro-seismic event scatter points of oil reservoir reconstruction in real time; step 2, calculating error ellipses of all microseism events; step 3, performing weighting operation according to the moment magnitude of the microseism event; step 4, calculating the geometric center of the microseism event within the error ellipse range; and 5, minimizing the objective function through continuous iteration, and accurately describing and finely depicting the artificial fracture for oil reservoir reconstruction.
The object of the invention can also be achieved by the following technical measures:
in step 1, the microseism event is obtained by positioning a plurality of detectors arranged in a well or on the ground to record waveforms, and comprises the occurrence time, the spatial position and the seismic source parameter information of the microseism event.
In step 2, error ellipses of all the micro-seismic event scatter points are calculated and used as constraint conditions of scatter point cluster analysis, namely as a movable geometric range of the micro-seismic event.
In step 2, calculating an error ellipse after acquiring the microseism scatter points, wherein the error ellipse calculation formula of the microseism event scatter points is as follows:
wherein: equation of time difference r with sigma value as residuei＝ti-t0-Ti(x0,y0,z0) Variance of partial differential matrix, t0Is the initial time of occurrence of the microseismic event, (x)0,y0,z0) Is the spatial position of the occurrence of the microseism event, i is the serial number of the detector, i is 1,2, …, N is the total number of the detectors, tiIs the observed microseismic event arrival time, T, of the ith geophoneiIs the theoretical analog arrival time of the ith detector; the B value is represented by the residual time difference equation ri＝ti-t0-Ti(x0,y0,z0) Calculating a partial differential matrix;
is a χ with a confidence of α and a degree of freedom of β2Distribution, α is a confidence parameter, β is a degree of freedom parameter, i.e., contains independent variablesNumber of (1), χ2Is a chi-square distribution.
In step 2, the value B in the error ellipse calculation formula is represented by the residual time difference equation ri＝ti-t0-Ti(x0,y0,z0) Is defined by the partial differential matrix a, the operational formula is as follows:
B＝ATA
wherein:
is formed by the equation of residual time difference ri＝ti-t0-Ti(x0,y0,z0) Is calculated from the partial differential matrix of (a),
is the partial differential of the theoretical simulated time function T; a. theTIs the transpose of matrix a.
In step 2, the sigma value in the error ellipse calculation formula is represented by the residual time difference equation ri＝ti-t0-Ti(x0,y0,z0) Is defined by the variance of the partial differential matrix a, the formula is as follows:
wherein: r isiIs residual time, i ═ 1,2, …, N; n is the total number of detectors.
In step 3, weighting operation is carried out according to the moment magnitude of the microseism event, the weight is obtained by calculating the moment magnitude of the microseism event, and the specific calculation formula is as follows:
Mw＝a·log10(M0)+b
wherein: a is a constant factor, taking value 2/3; b is a constant factor, and takes the value of-6; m0The seismic moment is determined by the amplitude low-frequency component of the seismic wave, and reflects the size of the fracture energy at the seismic source point, and the larger the fault area is, the larger the energy of the excited long-period seismic wave is, and the longer the period is.
In step 4, an Euclidean distance algorithm is adopted for analysis, whether the microseism event scatter point can move to the geometric center in the error ellipse range is judged, and the Euclidean distance is calculated through the following formula:
wherein: (xx)j,yyj,zzj) Is the spatial three-dimensional coordinate of the microseismic event scatter point in the error ellipse; (xx)0,yy0,zz0) Is the geometric center three-dimensional coordinate of the error ellipse; m is the number of microseismic event scatter points.
The object of the invention can also be achieved by the following technical measures:
the microseism incident scatter cluster analysis system is reformed transform to oil reservoir, and this microseism incident scatter cluster analysis system is reformed transform to oil reservoir includes: a microseismic event acquisition module: acquiring a microseism event, wherein the microseism event is obtained by positioning waveforms recorded by a plurality of detectors distributed in a well or on the ground, and the dispersion point of the microseism event is positioned in the oil reservoir transformation range and comprises the occurrence time, the spatial position and the seismic source point parameter information of the microseism event; an error ellipse calculation module: calculating an error ellipse of a scattered point of the microseism event, and taking the error ellipse as a constraint condition of scattered point clustering analysis, namely as a movable geometric range of the microseism event; moment-magnitude weighting module: calculating the weighting weight of the dispersion point of the microseism event, wherein the value of the weighting weight is obtained by calculating the moment magnitude of the microseism event; an error ellipse judgment module: calculating the geometric center of the microseism event in the error ellipse range, analyzing by adopting an Euclidean distance algorithm, and judging whether the microseism event can move to the geometric center in the error ellipse range; a microseismic event cluster iteration module: and minimizing the objective function through continuous iteration, and accurately describing and finely depicting the artificial fracture for oil reservoir reconstruction.
According to the scattered point cluster analysis method and system for the oil reservoir reconstruction micro-seismic events, scattered point cluster analysis is carried out on a micro-seismic event point set, and through the application of the method and system, the framework structure of the oil reservoir reconstruction artificial fracture can be rapidly obtained, and the dynamic spreading and real-time change conditions of the oil reservoir reconstruction artificial fracture framework can be clearly displayed, so that the reliability of the oil reservoir reconstruction micro-seismic monitoring result is greatly improved, and the function of the micro-seismic monitoring result in improving the oil gas recovery ratio is exerted. The scattered point clustering analysis method and the scattered point clustering analysis system for the micro-seismic events in oil reservoir reconstruction solve the problem of divergence of the micro-seismic events caused by P-wave or S-wave first arrival time-arrival errors, uncertainty of a velocity model, source vector angle errors and unsatisfactory space-time distribution of a detector. This patent can draw oil reservoir transformation artificial crack skeleton fast, and the clear artificial crack skeleton of showing oil reservoir transformation is opened up and is distributed and the situation of change to improve the reliability of oil reservoir transformation microseism monitoring result by a wide margin, the effect of performance microseism monitoring result.
Drawings
FIG. 1 is a flow chart of an embodiment of a method for reservoir reformation microseismic event scatter cluster analysis of the present invention;
FIG. 2 is a block diagram of an embodiment of a reservoir reformation microseismic event scatter cluster analysis system of the present invention.
Detailed Description
In order to make the aforementioned and other objects, features and advantages of the present invention comprehensible, preferred embodiments accompanied with figures are described in detail below.
As shown in fig. 1, fig. 1 is a flow chart of a reservoir reformation microseismic event scatter cluster analysis method of the present invention.
In step 101, the dispersion point of the micro-seismic event of oil reservoir reconstruction is obtained in real time, wherein the micro-seismic event is obtained by positioning the recording waveforms of a plurality of detectors arranged in a well or on the ground and comprises the occurrence time, the spatial position and the seismic source parameter information of the micro-seismic event.
Reservoir reformation over a period of time may produce tens or hundreds of microseismic events that are distributed over the reformed formation. The micro-seismic event occurrence time can reflect the opening sequence of the artificial cracks, the micro-seismic event space position can reflect the expansion condition and the spreading characteristic of the artificial cracks, and the micro-seismic event seismic source parameters can reflect the opening and sliding characteristics of the artificial cracks.
In step 102, error ellipses of all micro-seismic event scatter points are calculated and used as constraint conditions of scatter point cluster analysis, namely as a geometric range in which the micro-seismic events can move.
Calculating an error ellipse after acquiring the microseism scatter points, wherein the error ellipse calculation formula of the microseism event scatter points is as follows:
wherein: equation of time difference r with sigma value as residuei＝ti-t0-Ti(x0,y0,z0) Variance of partial differential matrix, t0Is the initial time of occurrence of the microseismic event, (x)0,y0,z0) Is the spatial position of the occurrence of the microseism event, i is the serial number of the detector, i is 1,2, …, N is the total number of the detectors, tiIs the observed microseismic event arrival time, T, of the ith geophoneiIs the theoretical analog arrival time of the ith detector; the B value is represented by the residual time difference equation ri＝ti-t0-Ti(x0,y0,z0) Calculating a partial differential matrix;
is a χ with a confidence of α and a degree of freedom of β2Distribution, α is a confidence parameter, β is a degree of freedom parameter, i.e., the number of independent variables, χ2Is a chi-square distribution.
The value B in the error ellipse calculation formula is represented by the residual time difference equation ri＝ti-t0-Ti(x0,y0,z0) Is defined by the partial differential matrix a, the operational formula is as follows:
B＝ATA
wherein:
is formed by the equation of residual time difference ri＝ti-t0-Ti(x0,y0,z0) Is calculated from the partial differential matrix of (a),
is the partial differential of the theoretical simulated time function T; a. theTIs the transpose of matrix a.
Residual time difference equation r of sigma value in error ellipse calculation formulai＝ti-t0-Ti(x0,y0,z0) Is defined by the variance of the partial differential matrix a, the formula is as follows:
wherein: r isiIs residual time, i ═ 1,2, …, N; n is the total number of detectors.
At step 103, a weighting operation is performed based on the magnitude of the microseismic event. The weight is obtained by calculating the moment magnitude of the microseismic event, and the specific calculation formula is as follows:
Mw＝a·log10(M0)+b
wherein: a is a constant factor, taking value 2/3; b is a constant factor, and takes the value of-6; m0The seismic moment is determined by the amplitude low-frequency component of the seismic wave, and reflects the size of the fracture energy at the seismic source point, and the larger the fault area is, the larger the energy of the excited long-period seismic wave is, and the longer the period is.
In step 104, the geometric center of the microseismic event within the error ellipse is calculated, and the euclidean distance algorithm is used for analysis to determine whether the microseismic event scatter points can move towards the geometric center within the error ellipse.
The euclidean distance is calculated by the following formula:
wherein: (xx)j,yyj,zzj) Is a micro in the error ellipseSpatial three-dimensional coordinates of seismic event scatter; (xx)0,yy0,zz0) Is the geometric center three-dimensional coordinate of the error ellipse; m is the number of microseismic event scatter points.
In step 105, through continuous optimization iteration of a nonlinear optimization algorithm, the optimization method enables the objective function value to be minimum, and finally the skeleton of the artificial fracture is generated through calculation, so that accurate description and fine drawing of the skeleton of the artificial fracture for oil reservoir reconstruction are achieved.
The nonlinear optimization algorithm is the combination of the least square optimization algorithm and the particle swarm optimization algorithm, has the advantages of fast convergence of a target function, high operation efficiency and small iteration error, and can meet the requirement of real-time dynamic analysis of artificial fractures in an oil reservoir reconstruction field.
A reservoir reforming micro-seismic event scatter cluster analysis system, as shown in fig. 2, comprising the following modules:
(1) a microseismic event acquisition module: the microseism event is obtained by positioning the waveform recorded by a plurality of detectors arranged in a well or on the ground, and the dispersion point of the microseism event is positioned in the oil reservoir transformation range and comprises information such as the occurrence time, the spatial position, the seismic source point parameter and the like of the microseism event;
(2) an error ellipse calculation module: the method comprises the steps of calculating an error ellipse of a microseism event scatter point, and using the error ellipse as a constraint condition of scatter point clustering analysis, namely as a movable geometric range of the microseism event;
(3) moment-magnitude weighting module: the weighting weight is used for calculating the scatter point of the microseism event, and the value of the weighting weight is calculated by the moment magnitude of the microseism event;
(4) an error ellipse judgment module: calculating the geometric center of the microseism event in the error ellipse range, analyzing by adopting an Euclidean distance algorithm, and judging whether the microseism event can move to the geometric center in the error ellipse range;
(5) a microseismic event cluster iteration module: the objective function is minimized through continuous iteration, so that the generated microseism event clusters are gathered and compacted as much as possible, the purpose of obtaining an artificial fracture framework is achieved, and accurate description and fine drawing of the oil reservoir reconstruction artificial fracture are achieved.
In the embodiment of the invention, the scattered-point clustering analysis method and the scattered-point clustering analysis system for the micro-seismic events of oil reservoir reconstruction are provided, the framework shape of the artificial fracture of oil reservoir reconstruction can be rapidly and clearly displayed, and the real-time analysis capability of the dynamic extension of the artificial fracture of oil reservoir reconstruction is greatly improved.
The above description is only a preferred embodiment of the present invention, and is not intended to limit the present invention, and various modifications and changes may be made to the embodiment of the present invention by those skilled in the art. Any modification, equivalent replacement, or improvement made within the spirit and principle of the present invention should be included in the protection scope of the present invention.
Claims (9)
1. The method for analyzing the scattered point clustering of the oil reservoir reconstruction micro-seismic events is characterized by comprising the following steps of:
step 1, acquiring micro-seismic event scatter points of oil reservoir reconstruction in real time;
step 2, calculating error ellipses of all microseism events;
step 3, performing weighting operation according to the moment magnitude of the microseism event;
step 4, calculating the geometric center of the microseism event within the error ellipse range;
and 5, minimizing the objective function through continuous iteration, and accurately describing and finely depicting the artificial fracture for oil reservoir reconstruction.
2. The method of claim 1, wherein in step 1, the microseismic events are obtained by positioning the recording waveforms of a plurality of detectors deployed in the well or on the ground, including the occurrence time, spatial position and source parameter information of the microseismic events.
3. The method of claim 1, wherein in step 2, the error ellipses of all the microseismic event scatters are calculated and used as constraints of the scatterer cluster analysis, namely, as the movable geometric range of the microseismic event.
4. The method for clustering microseism event scatter analysis of oil reservoir improvement according to claim 3, wherein in step 2, after the microseism scatter is obtained, an error ellipse is calculated, and an error ellipse calculation formula of the microseism event scatter is as follows:
wherein: equation of time difference r with sigma value as residuei＝ti-t0-Ti(x0,y0,z0) Variance of partial differential matrix, t0Is the initial time of occurrence of the microseismic event, (x)0,y0,z0) Is the spatial position of the occurrence of the microseism event, i is the serial number of the detector, i is 1,2, …, N is the total number of the detectors, tiIs the observed microseismic event arrival time, T, of the ith geophoneiIs the theoretical analog arrival time of the ith detector; the B value is represented by the residual time difference equation ri＝ti-t0-Ti(x0,y0,z0) Calculating a partial differential matrix;
is a χ with a confidence of α and a degree of freedom of β2Distribution, α is a confidence parameter, β is a degree of freedom parameter, i.e., the number of independent variables, χ2Is a chi-square distribution.
5. The method of claim 4, wherein in step 2, the B value in the error ellipse calculation formula is determined by the residual time difference equation ri＝ti-t0-Ti(x0,y0,z0) Is defined by the partial differential matrix a, the operational formula is as follows:
B＝ATA
wherein:
is formed by the equation of residual time difference ri＝ti-t0-Ti(x0,y0,z0) Is calculated from the partial differential matrix of (a),
is the partial differential of the theoretical simulated time function T; a. theTIs the transpose of matrix a.
6. The method of claim 4, wherein in step 2, the sigma value in the error ellipse calculation formula is represented by the residual equation ri＝ti-t0-Ti(x0,y0,z0) Is defined by the variance of the partial differential matrix a, the formula is as follows:
wherein: r isiIs residual time, i ═ 1,2, …, N; n is the total number of detectors.
7. The method for performing scattered clustering analysis on the microseismic events for oil reservoir reconstruction as claimed in claim 1, wherein in step 3, the weighting operation is performed according to the moment magnitude of the microseismic events, the weighting is obtained by calculating the moment magnitude of the microseismic events, and the specific calculation formula is as follows:
Mw＝a·log10(M0)+b
wherein: a is a constant factor, taking value 2/3; b is a constant factor, and takes the value of-6; m0The seismic moment is determined by the amplitude low-frequency component of the seismic wave, and reflects the size of the fracture energy at the seismic source point, and the larger the fault area is, the larger the energy of the excited long-period seismic wave is, and the longer the period is.
8. The method of claim 1, wherein in step 4, the method uses Euclidean distance algorithm to determine whether the microseismic event scatter moves towards the geometric center within the error ellipse, and the Euclidean distance is calculated by the following formula:
wherein: (xx)j,yyj,zzj) Is the spatial three-dimensional coordinate of the microseismic event scatter point in the error ellipse; (xx)0,yy0,zz0) Is the geometric center three-dimensional coordinate of the error ellipse; m is the number of microseismic event scatter points.
9. Little earthquake event scatter cluster analysis system is reformed transform to oil reservoir, its characterized in that, this little earthquake event scatter cluster analysis system of oil reservoir transformation includes:
a microseismic event acquisition module: acquiring a microseism event, wherein the microseism event is obtained by positioning waveforms recorded by a plurality of detectors distributed in a well or on the ground, and the dispersion point of the microseism event is positioned in the oil reservoir transformation range and comprises the occurrence time, the spatial position and the seismic source point parameter information of the microseism event;
an error ellipse calculation module: calculating an error ellipse of a scattered point of the microseism event, and taking the error ellipse as a constraint condition of scattered point clustering analysis, namely as a movable geometric range of the microseism event;
moment-magnitude weighting module: calculating the weighting weight of the dispersion point of the microseism event, wherein the value of the weighting weight is obtained by calculating the moment magnitude of the microseism event;
an error ellipse judgment module: calculating the geometric center of the microseism event in the error ellipse range, analyzing by adopting an Euclidean distance algorithm, and judging whether the microseism event can move to the geometric center in the error ellipse range;
a microseismic event cluster iteration module: and minimizing the objective function through continuous iteration, and accurately describing and finely depicting the artificial fracture for oil reservoir reconstruction.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010012257.4A CN111158045B (en) | 2020-01-06 | 2020-01-06 | Reservoir transformation microseism event scattered point clustering analysis method and system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010012257.4A CN111158045B (en) | 2020-01-06 | 2020-01-06 | Reservoir transformation microseism event scattered point clustering analysis method and system |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111158045A true CN111158045A (en) | 2020-05-15 |
| CN111158045B CN111158045B (en) | 2022-02-22 |
Family
ID=70561740
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010012257.4A Active CN111158045B (en) | 2020-01-06 | 2020-01-06 | Reservoir transformation microseism event scattered point clustering analysis method and system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111158045B (en) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112379423A (en) * | 2020-11-02 | 2021-02-19 | 中国石油天然气集团有限公司 | Method and device for determining moment-magnitude of microseism event |
| CN113253343A (en) * | 2021-05-12 | 2021-08-13 | 中油奥博(成都)科技有限公司 | Method for identifying fault activity of underground gas storage based on microseism monitoring technology |
| CN113687411A (en) * | 2021-09-03 | 2021-11-23 | 西南石油大学 | Earth stress azimuth prediction method based on microseism |
| CN113919500A (en) * | 2021-12-14 | 2022-01-11 | 矿冶科技集团有限公司 | Microseism event space-time fractal dimension calculation method and device and electronic equipment |
| CN114063177A (en) * | 2021-11-16 | 2022-02-18 | 核工业北京地质研究院 | Magnetotelluric data denoising method and system |
| CN114779330A (en) * | 2022-04-26 | 2022-07-22 | 中国矿业大学 | Mining working face main fracture azimuth analysis and prediction method based on microseismic monitoring |
Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101825723A (en) * | 2004-03-17 | 2010-09-08 | 维斯特恩格科地震控股有限公司 | Offshore earthquake measurement method and system |
| US20140334261A1 (en) * | 2011-08-29 | 2014-11-13 | Jonathan S. Abel | Method and system for microseismic event location error analysis and display |
| CN104182651A (en) * | 2014-09-12 | 2014-12-03 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Automatic microseism event azimuth angle quality control method used for three-component detector reception |
| CN104618137A (en) * | 2014-12-25 | 2015-05-13 | 中国海洋石油总公司 | Method and device for realizing towing exploration control network |
| CN105510880A (en) * | 2014-09-23 | 2016-04-20 | 中国石油化工股份有限公司 | Microseism focus positioning method based on double-difference method |
| US20160109590A1 (en) * | 2014-10-17 | 2016-04-21 | Pgs Geophysical As | Method and system of calculating hypocenters of microseismic events |
| CN107479093A (en) * | 2017-09-18 | 2017-12-15 | 中南大学 | A kind of micro-seismic event denoising and clustering method based on potential function |
| US9945970B1 (en) * | 2011-08-29 | 2018-04-17 | Seismic Innovations | Method and apparatus for modeling microseismic event location estimate accuracy |
| CN109188515A (en) * | 2018-10-31 | 2019-01-11 | 中国石油化工股份有限公司 | Micro-seismic monitoring crack focal point position calculating method and system |
-
2020
- 2020-01-06 CN CN202010012257.4A patent/CN111158045B/en active Active
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101825723A (en) * | 2004-03-17 | 2010-09-08 | 维斯特恩格科地震控股有限公司 | Offshore earthquake measurement method and system |
| US20140334261A1 (en) * | 2011-08-29 | 2014-11-13 | Jonathan S. Abel | Method and system for microseismic event location error analysis and display |
| US9945970B1 (en) * | 2011-08-29 | 2018-04-17 | Seismic Innovations | Method and apparatus for modeling microseismic event location estimate accuracy |
| CN104182651A (en) * | 2014-09-12 | 2014-12-03 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Automatic microseism event azimuth angle quality control method used for three-component detector reception |
| CN105510880A (en) * | 2014-09-23 | 2016-04-20 | 中国石油化工股份有限公司 | Microseism focus positioning method based on double-difference method |
| US20160109590A1 (en) * | 2014-10-17 | 2016-04-21 | Pgs Geophysical As | Method and system of calculating hypocenters of microseismic events |
| CN104618137A (en) * | 2014-12-25 | 2015-05-13 | 中国海洋石油总公司 | Method and device for realizing towing exploration control network |
| CN107479093A (en) * | 2017-09-18 | 2017-12-15 | 中南大学 | A kind of micro-seismic event denoising and clustering method based on potential function |
| CN109188515A (en) * | 2018-10-31 | 2019-01-11 | 中国石油化工股份有限公司 | Micro-seismic monitoring crack focal point position calculating method and system |
Non-Patent Citations (7)
| Title |
|---|
| DIAORUI ET AL.: "Application of spectrum simulation high frequency compensation in time frequency domain", 《2018 CPS/SEG INTERNATIONAL GEOPHYSICAL CONFERENCE》 * |
| M. HALLO ET AL.: "Microseismic Surface Monitoring Network Design - Sensitivity and Accuracy", 《74TH EAGE CONFERENCE & EXHIBITION INCORPORATING SPE EUROPEC 2012》 * |
| SUDIPTA SARKAR ET AL.: "Uncertainty in microseismic event source locations associated with a fracture stimulation in a tight sand", 《ISTANBUL INTERNATIONAL GEOPHYSICAL CONFERENCE AND OIL & GAS EXHIBITION》 * |
| 刁瑞等: "三维地震与地面微地震联合校正方法", 《地球物理学报》 * |
| 巩思园等: "煤矿微震监测台网优化布设研究", 《地球物理学报》 * |
| 延斯.哈夫斯科夫等著，刘俊民等译: "《地震学中的例行数据处理——实例数据、练习和软件》", 31 March 2017 * |
| 李应平等: "微震分析水压致裂的破裂过程", 《地震学报》 * |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112379423A (en) * | 2020-11-02 | 2021-02-19 | 中国石油天然气集团有限公司 | Method and device for determining moment-magnitude of microseism event |
| CN113253343A (en) * | 2021-05-12 | 2021-08-13 | 中油奥博(成都)科技有限公司 | Method for identifying fault activity of underground gas storage based on microseism monitoring technology |
| CN113687411A (en) * | 2021-09-03 | 2021-11-23 | 西南石油大学 | Earth stress azimuth prediction method based on microseism |
| CN114063177A (en) * | 2021-11-16 | 2022-02-18 | 核工业北京地质研究院 | Magnetotelluric data denoising method and system |
| CN114063177B (en) * | 2021-11-16 | 2023-09-29 | 核工业北京地质研究院 | Method and system for denoising magnetotelluric data |
| CN113919500A (en) * | 2021-12-14 | 2022-01-11 | 矿冶科技集团有限公司 | Microseism event space-time fractal dimension calculation method and device and electronic equipment |
| CN114779330A (en) * | 2022-04-26 | 2022-07-22 | 中国矿业大学 | Mining working face main fracture azimuth analysis and prediction method based on microseismic monitoring |
| CN114779330B (en) * | 2022-04-26 | 2022-12-27 | 中国矿业大学 | Mining working face main fracture azimuth analysis and prediction method based on microseismic monitoring |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111158045B (en) | 2022-02-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111158045B (en) | Reservoir transformation microseism event scattered point clustering analysis method and system | |
| CA2743611C (en) | Methods and systems for monitoring and modeling hydraulic fracturing of a reservoir field | |
| Wassermann | Volcano seismology | |
| Chouet | Volcano seismology | |
| US8121792B2 (en) | Integration of geomechanics and seismic analysis for passive seismic feasibility analysis | |
| Wang et al. | Current developments on micro-seismic data processing | |
| CN111175815B (en) | Method and system for solving micro-seismic monitoring crack seismic source mechanism in oil reservoir transformation | |
| WO2012139082A1 (en) | Event selection in the image domain | |
| Watson et al. | Volcano infrasound: Progress and future directions | |
| Yin et al. | Source time function clustering reveals patterns in earthquake dynamics | |
| Li et al. | Real-time characterization of finite rupture and its implication for earthquake early warning: application of FinDer to existing and planned stations in Southwest China | |
| CN111175816B (en) | Method and device for constructing micro-seismic fracture network in real time in oil reservoir reconstruction | |
| Jena et al. | A model for visual assessment of fault plane solutions and active tectonics analysis using the global centroid moment tensor catalog | |
| Battaglia et al. | Precise relocation of earthquakes following the 15 June 1991 eruption of Mount Pinatubo (Philippines) | |
| Atkinson et al. | On the mN, M relation for eastern North American earthquakes | |
| Jones et al. | Accurate location of synthetic acoustic emissions and location sensitivity to relocation methods, velocity perturbations, and seismic anisotropy | |
| Dahm | Rupture dimensions and rupture processes of fluid-induced microcracks in salt rock | |
| Zhang et al. | Velocity structure building and ground motion simulation of the 2014 Ludian M s 6.5 Earthquake | |
| Johnson | Stress modulation of earthquakes: A study of long and short period stress perturbations and the crustal response | |
| Wells | Examining Shallow Crustal Structure in the Intermountain West Using Portable Seismic Arrays | |
| CN110146921B (en) | Microseismic event detection method and system based on Dirac distribution probability | |
| Yang | A Seismologic Study of the Northern Mississippi Embayment | |
| Arrowsmith et al. | Infrasound signal characteristics from small earthquakes | |
| Yang | Seismic imaging of the lithosphere beneath the southern Illinois Basin and its tectonic implications | |
| Bunger et al. | Hydraulic Fracturing Mechanisms Leading to Self‐Organization Within Dyke Swarms |
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 |