CN110824562A - Microseism signal correction value calculation method and system - Google Patents
Microseism signal correction value calculation method and system Download PDFInfo
- Publication number
- CN110824562A CN110824562A CN201810896900.7A CN201810896900A CN110824562A CN 110824562 A CN110824562 A CN 110824562A CN 201810896900 A CN201810896900 A CN 201810896900A CN 110824562 A CN110824562 A CN 110824562A
- Authority
- CN
- China
- Prior art keywords
- correction amount
- signal
- perforation
- correction value
- energy profile
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000012937 correction Methods 0.000 title claims abstract description 133
- 238000004364 calculation method Methods 0.000 title claims abstract description 31
- 238000000034 method Methods 0.000 claims abstract description 20
- 230000008569 process Effects 0.000 claims description 9
- 238000004590 computer program Methods 0.000 claims description 3
- 238000012544 monitoring process Methods 0.000 abstract description 7
- 238000011161 development Methods 0.000 abstract description 6
- 230000000694 effects Effects 0.000 abstract description 3
- 238000004519 manufacturing process Methods 0.000 abstract description 2
- 230000008901 benefit Effects 0.000 description 3
- 230000015572 biosynthetic process Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 239000003245 coal Substances 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000001902 propagating effect Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 1
- 230000008676 import Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
- G01V1/362—Effecting static or dynamic corrections; Stacking
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
A method and system for calculating the correction value of microseism signal are disclosed. The method can comprise the following steps: calculating the dynamic correction value of each fracturing section; performing signal identification based on the dynamic correction value, and selecting a strong energy profile; obtaining a residual correction value based on the selected strong energy profile; a total correction amount for each fracture zone is calculated based on the calculated motional correction amount and the remaining correction amount. According to the method, the microseism signal correction value calculation method is researched and developed according to the requirement of quick and automatic pickup of the microseism signals, the effective fracture signals are ensured to be quickly identified, the guiding effect of microseism monitoring on actual production is achieved, and powerful support is provided for the exploration and development of unconventional oil and gas.
Description
Technical Field
The invention relates to the field of exploration and development of unconventional oil and gas such as shale gas and coal bed gas, in particular to a method and a system for calculating a microseismic signal correction value.
Background
The microseismic monitoring technology is a commonly used monitoring technology in the development of compact reservoir oil and gas fields, and the change of fractures in the fracturing and oil and gas development processes is analyzed by processing detection signals, so that the reservoir transformation effect is dynamically evaluated, and technical support is provided for unconventional oil and gas development.
Microseismic monitoring requires precise, real-time reversal of the location of the microseismic source. The seismic source position inversion is to determine the occurrence position of an event by using seismic wave information collected by the detectors, namely to determine the coordinates of fracture points caused by fracturing, and further to determine the extending angle and length of the fracture formed by fracturing. The real-time requirement of field fracturing on fracturing monitoring processing results is higher and higher, in order to quickly provide positioning results, an important step in the processing flow is to quickly and automatically pick up useful fracture signals in a large amount of collected data so as to provide monitoring results for fracturing personnel in real time, timely master the extension condition of underground fracturing fractures and provide a basis for optimizing a fracturing process.
Because the detectors for collecting signals are distributed at different positions, the same signal generated by underground fracturing can reach all the detectors at different time. In order to quickly and automatically pick up a useful burst signal, signals received by different detectors are first corrected to the same time (signal burst time), wherein the signal profile is corrected by using the correction amount. Therefore, it is necessary to develop a method and a system for calculating the correction amount of the micro-seismic signal.
The information disclosed in this background section is only for enhancement of understanding of the general background of the invention and should not be taken as an acknowledgement or any form of suggestion that this information forms the prior art already known to a person skilled in the art.
Disclosure of Invention
The invention provides a microseism signal correction value calculation method and a microseism signal correction value calculation system, which can calculate the correction value of a microseism fracture signal through a certain technical means so as to quickly and automatically pick up an effective microseism signal.
According to an aspect of the present invention, a microseismic signal correction amount calculation method is provided. The method may include:
1) calculating the dynamic correction value of each fracturing section;
2) performing signal identification based on the dynamic correction value, and selecting a strong energy profile;
3) obtaining a residual correction value based on the strong energy profile selected in the step 2);
4) calculating a total correction amount for each fracture section based on the kinetic correction amount calculated in step 1) and the remaining correction amount obtained in step 3).
Preferably, in step 1), the time for the perforation signal to propagate to the detector is calculated according to the perforation coordinates and the perforation-to-detector speed, and the dynamic correction value of each fracturing section is obtained.
Preferably, the strong energy profile in step 2) is an energy profile with a perforation signal at each detector.
Preferably, in step 3), the remaining correction amount is obtained by picking up the jitter of the perforation waveform and applying further straightening on the perforation signal.
Preferably, in step 4), the total correction value of each fracture section is obtained by summing the motion correction value of each fracture section and the residual correction value.
According to another aspect of the present invention, there is provided a microseismic signal correction amount calculation system having a computer program stored thereon, wherein the program when executed by a processor implements the steps of:
step 1: calculating the dynamic correction value of each fracturing section;
step 2: performing signal identification based on the dynamic correction value, and selecting a strong energy profile;
and step 3: obtaining a residual correction value based on the strong energy profile selected in the step 2;
and 4, step 4: calculating a total correction amount for each fracture section based on the motional correction amount calculated in step 1 and the remaining correction amount obtained in step 3.
Preferably, in step 1, the time for the perforation signal to propagate to the detector is calculated according to the perforation coordinates and the perforation-to-detector speed, and the dynamic correction value of each fracture section is obtained.
Preferably, in step 2, the energy profile with signals in each detector is selected as the strong energy profile in the signal identification process.
Preferably, in step 3, the remaining correction amount is obtained by picking up the jitter of the perforation waveform and applying a further straightening on the perforation signal.
Preferably, in step 4, the motion correction amount of each fracturing section obtained by the motion correction amount calculation module and the residual correction amount obtained by the residual correction amount calculation module are summed to obtain a total correction amount of each fracturing section.
The present invention has other features and advantages which will be apparent from or are set forth in detail in the accompanying drawings and the following detailed description, which are incorporated herein, and which together serve to explain certain principles of the invention.
Drawings
The above and other objects, features and advantages of the present invention will become more apparent by describing in more detail exemplary embodiments thereof with reference to the attached drawings, in which like reference numerals generally represent like parts.
FIG. 1 is a flow chart illustrating the steps of a method of calculating a microseismic signal correction in accordance with the present invention;
FIG. 2 is a schematic diagram of a window of a microseismic signal correction amount calculation system according to the present invention;
fig. 3 shows a schematic diagram of a selected intense energy profile.
Detailed Description
The invention will be described in more detail below with reference to the accompanying drawings. While the preferred embodiments of the present invention are shown in the drawings, it should be understood that the present invention may be embodied in various forms and should not be limited to the embodiments set forth herein. Rather, these embodiments are provided so that this disclosure will be thorough and complete, and will fully convey the scope of the invention to those skilled in the art.
FIG. 1 is a flow chart showing the steps of a microseismic signal correction amount calculation method of the present invention.
In this embodiment, the microseismic signal correction amount calculation method according to the present invention may include:
step 101, calculating a dynamic correction value of each fracturing section;
in one example, in step 101, the time for the perforation signal to propagate to the geophone is calculated according to the perforation coordinates and the perforation-to-geophone velocity, and the dynamic correction value of each fracture section is obtained.
Specifically, since the receiving points (detectors) are different from the breaking point in distance, the correction amount is also different; microseismic fracturing is typically a multi-stage fracture, each stage being at a different distance from the same detector. Therefore, the correction amount for different fracture stages is different for the same detector, and the correction amount is changed and is called as a dynamic correction amount (NMO). The amount of dynamic correction needs to be recalculated for each segment of signal identification at each segment of microseismic fracture.
102, carrying out signal identification based on the dynamic correction value, and selecting a strong energy profile;
in one example, the strong energy profile is an energy profile with a perforation signal at each receiver.
Specifically, under normal conditions, an accurate correction can be obtained if the perforation signal is relatively strong. However, the perforation signals are generally lost greatly when propagating through the formation to the detectors, and not all the detectors can clearly receive the perforation signals, so that the Residual correction amount Residual cannot be calculated, and the picked-up fracture signals are not positioned accurately. Therefore, the dynamic correction value obtained in step 101 is applied to a signal identification process, that is, a process of combining different signal processing methods to detect an effective signal, so that a section where each detector has a signal is selected as a strong energy section, and a residual correction value is picked up on the strong energy section, so that the signals identified later can be accurately positioned.
103, obtaining a residual correction value based on the strong energy profile selected in the step 102;
in an exemplary embodiment, the remaining correction amount is obtained in step 103 by picking up the jitter of the perforation waveform and applying further straightening on the perforation signal.
Specifically, the dynamic correction value is applied to the perforation signal, and the waveforms of the perforation signals received by the detectors are basically in a straight line. However, because the velocity model cannot be built accurately, the corrected perforation signal waveform will have small jitter and errors will occur in locating such perforation signals. The remaining correction can be picked up by picking up the jitter of the perforation waveform and applying a further straightening on the perforation signal. Wherein, the signal straightening is to draw the signal with the waveform of arc or shaking along the horizontal direction to be a straight line in the horizontal direction by using the correction amount. Meanwhile, because each fracturing stage is not far away, the residual correction amount calculated for the first time can be commonly used in adjacent stages.
And 104, calculating a total correction amount of each fracturing section based on the dynamic correction amount calculated in the step 101 and the residual correction amount obtained in the step 103.
In an exemplary embodiment, the amount of the per fracture zone dynamics correction and the remaining correction are summed to obtain a total correction per fracture zone in step 104.
And further applying the total correction value in a microseism signal detection process so as to automatically detect effective signals.
Application example
To facilitate understanding of the solution of the embodiments of the present invention and the effects thereof, a specific application example is given below. It will be understood by those skilled in the art that this example is merely for the purpose of facilitating an understanding of the present invention and that any specific details thereof are not intended to limit the invention in any way.
The invention also proposes a microseismic signal correction amount calculation system having a computer program stored thereon, wherein the program when executed by a processor implements the steps of:
step 1: calculating the dynamic correction value of each fracturing section;
step 2: performing signal identification based on the dynamic correction value, and selecting a strong energy profile;
and step 3: obtaining a residual correction value based on the strong energy profile selected in the step 2;
and 4, step 4: calculating a total correction amount for each fracture section based on the motional correction amount calculated in step 1 and the remaining correction amount obtained in step 3.
In step 1, the time of the perforation signal propagating to the detector is calculated according to the perforation coordinate and the perforation-to-detector speed, and the dynamic correction value of each fracturing section is obtained.
In step 2, the energy profile of each detector with signals is selected as a strong energy profile in the signal identification process.
In step 3, the remaining correction is obtained by picking up the jitter of the perforation waveform and applying further straightening on the perforation signal.
In step 4, the dynamic correction value of each fracturing section obtained by the dynamic correction value calculation module and the residual correction value obtained by the residual correction value calculation module are summed to obtain the total correction value of each fracturing section.
As shown in fig. 2, the actual window interface of the microseismic signal correction amount calculation system according to the present invention shows all interfaces of the microseismic signal correction amount calculation system according to the present invention, through which the dynamic correction amount, the residual correction amount, and the total correction amount of each fracture section solved step by step in the microseismic signal correction amount calculation system can be realized, and the function of saving and applying the calculation results can be realized. In the figure, ShootX represents a perforation X coordinate, ShootY represents a perforation Y coordinate, ShootZ represents a perforation Z coordinate, ShootT represents a perforation time, and VelFile represents a formation velocity file path; LoadVel denotes an import speed file, the formation speed file can be imported directly for calculation, Total denotes a Total correction amount, NMO denotes a dynamic correction amount, Residual denotes a remaining correction amount, Ele denotes a correction amount Residual, WaveLen denotes a wavelength, CorHaf denotes a coordinate interval, TraNum denotes a number of tracks, where WaveLen wavelength and CorHaf coordinate interval and TraNum number are default parameters.
The data values of the amount of fluctuation calculated by the system are shown in table 1.
TABLE 1 amount of dynamics correction
| Road number | 1 | 2 | 3 | …… | N |
| NMO (millisecond) | 550.322 | 555.809 | 553.719 | …… | 551.968 |
As shown in fig. 3, for the strong energy profile selected through the signal identification process, it can be seen from the data in the figure that when the length of the profile is 2120ms, a strong energy valid signal that can be detected by all the detectors appears, that is, the profile is selected as the strong energy profile. The ordinate of fig. 2 represents the cross-sectional length in units: milliseconds (ms), the lower abscissa indicates the seismic trace number and the upper abscissa indicates the line number.
In the energy profile of fig. 3, the remaining correction values calculated as shown in table 2 were selected.
TABLE 2 remaining correction amount
| Road number | 1 | 2 | 3 | …… | N |
| Residual(ms) | 12.571 | -13.571 | 10.571 | …… | -9.571 |
Those of ordinary skill in the art will understand that: all or a portion of the steps of implementing the above-described method embodiments may be performed by hardware associated with program instructions. The program may be stored in a computer-readable storage medium. When executed, the program performs steps comprising the method embodiments described above; and the aforementioned storage medium includes: various media that can store program codes, such as ROM, RAM, magnetic or optical disks.
In conclusion, the invention can calculate the correction value of the microseism fracture signal so as to quickly and automatically pick up the effective microseism signal, can ensure the real-time guiding significance of microseism monitoring on the actual production by utilizing the method technology, and provides powerful support for the exploration and development service of unconventional oil gas and coal bed gas.
It will be appreciated by persons skilled in the art that the above description of embodiments of the invention is intended only to illustrate the benefits of embodiments of the invention and is not intended to limit embodiments of the invention to any examples given.
Having described embodiments of the present invention, the foregoing description is intended to be exemplary, not exhaustive, and not limited to the embodiments disclosed. Many modifications and variations will be apparent to those of ordinary skill in the art without departing from the scope and spirit of the described embodiments.
Claims (10)
1. A method of calculating a microseismic signal correction amount, comprising:
1) calculating the dynamic correction value of each fracturing section;
2) performing signal identification based on the dynamic correction value, and selecting a strong energy profile;
3) obtaining a residual correction value based on the strong energy profile selected in the step 2);
4) calculating a total correction amount for each fracture section based on the kinetic correction amount calculated in step 1) and the remaining correction amount obtained in step 3).
2. The microseismic signal correction amount calculation method according to claim 1 wherein in step 1), the time for the perforation signal to propagate to the geophone is calculated from the perforation coordinates and the perforation-to-geophone velocity, and the dynamic correction amount for each fracture section is obtained.
3. The microseismic signal correction amount calculation method of claim 1 wherein the strong energy profile in step 2) is an energy profile with a perforation signal at each receiver.
4. The microseismic signal correction amount calculation method of claim 1 wherein in step 3), the remaining correction amount is obtained by picking up the jitter of the perforation waveform and applying further straightening on the perforation signal.
5. The microseismic signal correction amount calculation method of claim 4 wherein in step 4) the motional correction amount for each fracture zone and the remaining correction amount are summed to obtain a total correction amount for each fracture zone.
6. A microseismic signal correction amount calculation system having a computer program stored thereon, wherein the program when executed by a processor implements the steps of:
step 1: calculating the dynamic correction value of each fracturing section;
step 2: performing signal identification based on the dynamic correction value, and selecting a strong energy profile;
and step 3: obtaining a residual correction value based on the strong energy profile selected in the step 2;
and 4, step 4: calculating a total correction amount for each fracture section based on the motional correction amount calculated in step 1 and the remaining correction amount obtained in step 3.
7. The microseismic signal correction amount calculation system of claim 6 wherein the time for the perforation signal to propagate to the geophone is calculated from the perforation coordinates and the perforation-to-geophone velocity in step 1 to obtain a dynamic correction amount for each fracture section.
8. The microseismic signal correction amount calculation system of claim 6 wherein in step 2, the strong energy profile is selected by selecting the energy profile with signals for each detector in the signal identification process.
9. The microseismic signal correction amount calculation system of claim 6 wherein in step 3, a residual correction amount is obtained by picking up the jitter of the perforation waveform and applying further straightening on the perforation signal.
10. The microseismic signal correction amount calculation system of claim 6 wherein in step 4, the motion correction amount for each fracture zone obtained by the motion correction amount calculation module and the remaining correction amount obtained by the remaining correction amount calculation module are summed to obtain a total correction amount for each fracture zone.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810896900.7A CN110824562A (en) | 2018-08-08 | 2018-08-08 | Microseism signal correction value calculation method and system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810896900.7A CN110824562A (en) | 2018-08-08 | 2018-08-08 | Microseism signal correction value calculation method and system |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN110824562A true CN110824562A (en) | 2020-02-21 |
Family
ID=69540705
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810896900.7A Pending CN110824562A (en) | 2018-08-08 | 2018-08-08 | Microseism signal correction value calculation method and system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110824562A (en) |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP0430688A2 (en) * | 1989-11-30 | 1991-06-05 | Conoco Inc. | Common offset depth migration with residual moveout correction |
| EP2592439A2 (en) * | 2011-11-10 | 2013-05-15 | PGS Geophysical AS | Method and system for separating seismic sources in marine simultaneous shooting acquisition |
| CN104216008A (en) * | 2013-06-05 | 2014-12-17 | 中国石油天然气集团公司 | Downhole fracturing microseismic event identification method |
| CN106154321A (en) * | 2016-07-20 | 2016-11-23 | 中国石油天然气集团公司 | The detection method of perforation signal and device |
| CN106547024A (en) * | 2015-09-16 | 2017-03-29 | 中国石油化工股份有限公司 | For the residual static correction amount estimation method and device of microseism perforation data |
| WO2017123240A1 (en) * | 2016-01-15 | 2017-07-20 | Landmark Graphics Corporation | Semblance-based anisotropy parameter estimation using isotropic depth-migrated common image gathers |
-
2018
- 2018-08-08 CN CN201810896900.7A patent/CN110824562A/en active Pending
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP0430688A2 (en) * | 1989-11-30 | 1991-06-05 | Conoco Inc. | Common offset depth migration with residual moveout correction |
| EP2592439A2 (en) * | 2011-11-10 | 2013-05-15 | PGS Geophysical AS | Method and system for separating seismic sources in marine simultaneous shooting acquisition |
| CN104216008A (en) * | 2013-06-05 | 2014-12-17 | 中国石油天然气集团公司 | Downhole fracturing microseismic event identification method |
| CN106547024A (en) * | 2015-09-16 | 2017-03-29 | 中国石油化工股份有限公司 | For the residual static correction amount estimation method and device of microseism perforation data |
| WO2017123240A1 (en) * | 2016-01-15 | 2017-07-20 | Landmark Graphics Corporation | Semblance-based anisotropy parameter estimation using isotropic depth-migrated common image gathers |
| CN106154321A (en) * | 2016-07-20 | 2016-11-23 | 中国石油天然气集团公司 | The detection method of perforation signal and device |
Non-Patent Citations (3)
| Title |
|---|
| 刁瑞，等: "三维地震与地面微地震联合校正方法", 《地球物理学报》 * |
| 姜宇东，等: "地面微地震监测资料静校正方法研究", 《石油物探》 * |
| 程磊磊，等: "基于强事件约束的微地震剩余静校正量估算方法", 《石油物探》 * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106353792B (en) | Method suitable for positioning micro-seismic source of hydraulic fracturing | |
| CN100340871C (en) | Methods and apparatus for processing acoustic waveforms received in a borehole | |
| CN104216008B (en) | Downhole fracturing microseismic event identification method | |
| CN110529087B (en) | Method and device for evaluating hydraulic fracturing effect of stratum | |
| US10345463B2 (en) | Methods and systems for using known source events in seismic data processing | |
| AU2014345427B2 (en) | Method and device for processing seismic signals | |
| CN110687592B (en) | Microseism seismic phase identification first arrival picking method, device and storage medium | |
| Zhu et al. | Accurate identification of microseismic P-and S-phase arrivals using the multi-step AIC algorithm | |
| CN105974407A (en) | Ground penetrating radar underground horizon detection method | |
| CN105204065A (en) | Method and device for picking up preliminary wave | |
| CN110398775B (en) | Tunnel water burst disaster micro-seismic event signal fluctuation first arrival pickup method and system | |
| CN106646592A (en) | Single shot record real-time stacking method | |
| CN106353794B (en) | Microseism velocity model correction method based on relative first arrival matching error | |
| CN110824562A (en) | Microseism signal correction value calculation method and system | |
| CN106291699B (en) | The inclined detection method of big gun and the correct method for blowing out position of positioning in seismic prospecting | |
| CN116466396A (en) | Nonlinear optimal positioning method for microseism focus | |
| CN107664771B (en) | A kind of microseism Full wave shape localization method based on likeness coefficient | |
| CN114415231A (en) | Microseismic positioning method based on EDT surface probability distribution function of station | |
| CN105319592B (en) | The adaptive drawing method of ground Active noise in microseism data | |
| CN105093276B (en) | Knowledge method is sentenced in the digitlization of microseism useful signal | |
| CN116338774A (en) | Underground microseism inversion method and system based on distributed optical fiber sensor | |
| CN109507721B (en) | Method and system for evaluating quality of microseism event in well | |
| CN110967762B (en) | Microseism stratum velocity calibration method and system | |
| CN104749627A (en) | Microseismic signal highlighting method based on similarity | |
| CN111312272A (en) | Products, methods and systems for reducing noise signals in near-wellbore acoustic data sets |
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 | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20200221 |
|
| RJ01 | Rejection of invention patent application after publication |