CN106199697A - The elastic wave forward modeling method of simulation microseism - Google Patents
The elastic wave forward modeling method of simulation microseism Download PDFInfo
- Publication number
- CN106199697A CN106199697A CN201610494146.5A CN201610494146A CN106199697A CN 106199697 A CN106199697 A CN 106199697A CN 201610494146 A CN201610494146 A CN 201610494146A CN 106199697 A CN106199697 A CN 106199697A
- Authority
- CN
- China
- Prior art keywords
- microseism
- wave
- simulation
- elastic wave
- forward modeling
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 54
- 238000004088 simulation Methods 0.000 title claims abstract description 31
- 239000011435 rock Substances 0.000 claims abstract description 14
- 239000011159 matrix material Substances 0.000 claims abstract description 6
- 238000004458 analytical method Methods 0.000 claims abstract description 4
- 230000005855 radiation Effects 0.000 claims abstract description 4
- 230000000644 propagated effect Effects 0.000 claims abstract description 3
- 238000005516 engineering process Methods 0.000 claims description 7
- 238000006073 displacement reaction Methods 0.000 claims description 4
- 230000000694 effects Effects 0.000 claims description 4
- 230000008569 process Effects 0.000 claims description 3
- 230000003595 spectral effect Effects 0.000 claims description 3
- 238000004611 spectroscopical analysis Methods 0.000 claims description 3
- 235000007926 Craterellus fallax Nutrition 0.000 abstract description 2
- 240000007175 Datura inoxia Species 0.000 abstract description 2
- 230000004807 localization Effects 0.000 abstract description 2
- 230000035939 shock Effects 0.000 abstract description 2
- 238000004613 tight binding model Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 3
- 230000007547 defect Effects 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000008034 disappearance Effects 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 238000004880 explosion Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 210000003205 muscle Anatomy 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000010189 synthetic 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. for interpretation or for event detection
- G01V1/282—Application of seismic models, synthetic seismograms
-
- 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/30—Analysis
- G01V1/307—Analysis for determining seismic attributes, e.g. amplitude, instantaneous phase or frequency, reflection strength or polarity
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/63—Seismic attributes, e.g. amplitude, polarity, instant phase
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
The present invention provides a kind of elastic wave forward modeling method simulating microseism, including: step 1, according to the real data of geological exploration, analyze the character of rock, determine the pattern of rock rupture, generate the moment tensor matrix expressing microseism source；Step 2, sets up the geology velocity field model of required survey area；Step 3, under wave theory, sets up and is applicable to describe the elastic wave wave equation that microseism wave field is propagated；Step 4, under wave theory, solves microseism elastic wave wave equation, it is achieved the wave field of microseism propagates simulation；Step 5, simulation completes, multi-components record and the energy radiation mode of ripple in length and breadth in output microseism ground, well.The elastic wave forward modeling method of this simulation microseism can be quickly obtained the more accurate analog result of relatively ray theory, and quantity of information more horn of plenty, and provide superior technique means and instrument for follow-up characteristic rule analysis, localization method research and focal shock parameter inverting.
Description
Technical field
The present invention relates to oil field development technical field, especially relate to a kind of elastic wave just side of drilling simulating microseism
Method.
Background technology
Since entering 21 century, unconventional petroleum resources shows huge potentiality to be exploited.At numerous unconventionaloil pool
In exploration and development approach, microseismic plays very important effect.
In order to fully grasp micro-seismic technology, need the direct problem to microseism, focal mechanism carry out deep understanding and
Research, and realizing this target approach the most easily is Seismic wave numerical modeling technology.Conventional Seismic forward technology uses
Seismic wavelet, as focus, is suitable only for hypocenter of the explosion or the controlled source of common seismic exploration, and can not be well to microseism
Focus is described.At present a kind of major way of microseism forward simulation is based on ray theory, give some dielectric model and
Microseismograms is obtained by ray tracing mode under the premise of hypocentral location.But the microseismograms of ray theory Imitating cannot
Reproduce the different response to medium of the different frequency composition, record lacks Mintrop wave information, and may be only available for smooth medium.No
Only such, the microseism source observed in reality has the widest frequency band range and depends on the situation of change of medium, the highest
Frequently the hypothesis of ray has had a strong impact on the scope of application of microseism numerical simulation.Another way, based on wave theory, uses sky
So rock rupture pattern is portrayed by the moment tensor of earthquake, and uses Green function method to carry out elastic wave propagation numerical simulation.
But Green function method needs repeatedly to ask for Green's function during asking for, amount of calculation is relatively big, and asks for Green's function and need it
Carry out approximate expansion, this reduces the precision of simulation.The elastic wave of a kind of new simulation microseism is we have invented for this
Forward modeling method, solves above technical problem.
Summary of the invention
It is an object of the invention to provide a kind of moment tensor concept introducing and portraying focus described in earthquake, it is achieved to pressure
Split the elastic wave forward modeling method of the simulation microseism portrayed of focus.
The purpose of the present invention can be achieved by the following technical measures: the elastic wave forward modeling method of simulation microseism, this mould
The elastic wave forward modeling method intending microseism includes: step 1, according to the real data of geological exploration, analyzes the character of rock, determines
The pattern of rock rupture, generates the moment tensor matrix expressing microseism source；Step 2, sets up the geology speed of required survey area
Field model；Step 3, under wave theory, sets up and is applicable to describe the elastic wave wave equation that microseism wave field is propagated；Step 4,
Under wave theory, solve microseism elastic wave wave equation, it is achieved the wave field of microseism propagates simulation；Step 5, simulation completes,
Multi-components record and the energy radiation mode of ripple in length and breadth in output microseism ground, well.
The purpose of the present invention realizes also by following technical measures:
In step 1, by the analysis to the actual microseism data received in hydraulic fracturing process, it is thus achieved that microseism
Magnitude, thus generate express microseism source moment tensor matrix
Wherein Mij is rock rupture magnitude, and each element in moment tensor represents the effect of a pair power in opposite direction.
In step 2, read the geology velocity field model information of the required survey area pre-build, set up required exploration
The geology velocity field model in region, i.e. completes the modeling of zone velocity field.
In step 3, the microseism elastic wave wave equation of foundation is:
Wherein ui(x, t) and uk(x t) is wave field displacement, cijklFor stiffness coefficient, ρ is Media density, and x is that spatial point is sat
Mark, t is time, under three-dimensional situation, i, j, k, l=x, y, z, xsFor focal point position, S (t) is source wavelet, and M is moment tensor.
In step 3, utilize the moment tensor that step 1 obtains, substitute into microseism wave equation, use equivalent force method the other side
Focus item on the right of journey carries out equivalent dispersion.
In step 4, staggered-mesh high-order finite difference method technology, pseudo-spectrometry, FInite Element, spectral element method or iteration is utilized to ask
Solution solves microseism elastic wave wave equation.
The elastic wave forward modeling method of the simulation microseism in the present invention, overcomes the defect that ray theory high frequency is assumed, keeps away
Exempt from the shortcoming that Green function method is computationally intensive, simulation precision is low, take into full account the microseism shake caused by different pressure break pattern
On the premise of Source Type, use equivalent force method that microseism focus is described, carried out the biography of seismic wave by wave theory
Broadcast simulation, it is possible to be quickly obtained the more accurate analog result of relatively ray theory, and quantity of information more horn of plenty, and be follow-up spy
Levy law-analysing, localization method research and focal shock parameter inverting and superior technique means and instrument are provided.
Accompanying drawing explanation
Fig. 1 is the flow chart of the specific embodiment of the elastic wave forward modeling method of the simulation microseism of the present invention；
Fig. 2 be the present invention a specific embodiment in theoretical record (solid line), analog record of the present invention (chain-dotted line) and profit
The schematic diagram compared with Green Function Method analog record (dotted line) three oscillating curve X-component；
Fig. 3 be the present invention a specific embodiment in theoretical record (solid line), analog record of the present invention (chain-dotted line) and profit
The schematic diagram compared by Green Function Method analog record (dotted line) three oscillating curve Y-component；
Fig. 4 be the present invention a specific embodiment in theoretical record (solid line), analog record of the present invention (chain-dotted line) and profit
The schematic diagram compared with Green Function Method analog record (dotted line) three oscillating curve Z component.
Detailed description of the invention
For making the above and other objects, features and advantages of the present invention to become apparent, cited below particularly go out preferably to implement
Example, and coordinate shown in accompanying drawing, it is described in detail below.
As it is shown in figure 1, the flow chart of the elastic wave forward modeling method of the simulation microseism that Fig. 1 is the present invention.
In step 101, according to the real data of geological exploration, analyze the character of rock, determine the pattern of rock rupture.With
Time, by the analysis to the actual microseism data received in hydraulic fracturing process, it is thus achieved that the magnitude of microseism, thus generate
Express the moment tensor matrix in microseism sourceWherein MijFor rock rupture magnitude.In moment tensor
Each element represents the effect of a pair power in opposite direction.
In step 102, read the geology velocity field model information such as stiffness coefficient of the required survey area pre-build
cijkl, Media density ρ；Set up the geology velocity field model of required survey area, i.e. complete the modeling of zone velocity field.
Under step 103, wave theory, set up and be applicable to describe the elastic wave wave equation that microseism wave field is propagatedWherein ui(x, t) and uk(x t) is wave field displacement, cijkl
For stiffness coefficient, ρ is Media density, and x is spatial point coordinate, and t is the time.Under three-dimensional situation, i, j, k, l=x, y, z, xsFor shake
Source location, S (t) is source wavelet, and M is moment tensor.Utilize the moment tensor that step 101 obtains, substitute into microseism wave equationThe focus item on the right of equivalent force method the other side's journey is used to carry out equivalence
Discrete.
In step 104, under wave theory, utilize in staggered-mesh high-order finite difference method technology solution procedure 103 is micro-
Earthquake wave equation, it is achieved the wave field of microseism propagates simulation；Also can be with pseudo-spectrometry, FInite Element, spectral element method or iterative
Method realizes.
In step 105, simulation completes, multi-components record and the energy radiation mode of ripple in length and breadth in output microseism ground, well
Formula.
Fig. 2-4 is the simulative example utilizing the present invention.Simulation medium uses three-dimensional uniform dielectric, and observation system is earth's surface
Receiving, its source moment tensor is shown in (1) formula
Wherein M0For rock rupture magnitude.Fig. 2-4 illustrates one oscillating curve of three displacement components.Wherein solid line is
Analytic solutions, dotted line and chain-dotted line are the analog result utilizing Green functional based method and the inventive method respectively.It can be seen that it is sharp
The numerical solution obtained by the inventive method is basically identical with analytic solutions, and Green functional based method is owing to existing first approximation, error
Bigger than normal.And owing to Green functional based method needs successive ignition, calculating time to be about three times of the inventive method.Fig. 2 is microseism
The X-component of multi-components record, Fig. 3 is Y-component, and Fig. 4 is Z component.
Present invention introduces the moment tensor concept portraying focus described in earthquake, it is achieved pressure break focus is portrayed.Logical
Cross the definition of moment tensor, it is achieved the microseism fluctuation equation number value analogy method of high-efficiency high-accuracy.This method overcomes ray reason
The defect that opinion high frequency is assumed, it is to avoid the problem of ray simulation lower ripple phenomenon disappearance；Meanwhile, this method uses muscle power mode of equal value
The focus realizing rock rupture pattern realizes, and more conventional kinetics Green's function synthetic method is more flexibly quickly.
Claims (6)
1. simulate the elastic wave forward modeling method of microseism, it is characterised in that the elastic wave forward modeling method of this simulation microseism includes:
Step 1, according to the real data of geological exploration, analyzes the character of rock, determines the pattern of rock rupture, generate expression micro-
The moment tensor matrix of the seismic origin；
Step 2, sets up the geology velocity field model of required survey area；
Step 3, under wave theory, sets up and is applicable to describe the elastic wave wave equation that microseism wave field is propagated；
Step 4, under wave theory, solves microseism elastic wave wave equation, it is achieved the wave field of microseism propagates simulation；
Step 5, simulation completes, multi-components record and the energy radiation mode of ripple in length and breadth in output microseism ground, well.
The elastic wave forward modeling method of simulation microseism the most according to claim 1, it is characterised in that in step 1, pass through
Analysis to the actual microseism data received in hydraulic fracturing process, it is thus achieved that the magnitude of microseism, thus it is micro-to generate expression
The moment tensor matrix of the seismic origin
Wherein Mij is rock rupture magnitude, and each element in moment tensor represents the effect of a pair power in opposite direction.
The elastic wave forward modeling method of simulation microseism the most according to claim 1, it is characterised in that in step 2, reads
The geology velocity field model information of the required survey area pre-build, sets up the geology velocity field model of required survey area,
I.e. complete the modeling of zone velocity field.
The elastic wave forward modeling method of simulation microseism the most according to claim 1, it is characterised in that in step 3, sets up
Microseism elastic wave wave equation be:
Wherein ui(x, t) and uk(x t) is wave field displacement, cijklFor stiffness coefficient, ρ is Media density, and x is spatial point coordinate, t
For time, under three-dimensional situation, i, j, k, l=x, y, z, xsFor focal point position, S (t) is source wavelet, and M is moment tensor.
The elastic wave forward modeling method of simulation microseism the most according to claim 1, it is characterised in that in step 3, utilizes
The moment tensor that step 1 obtains, substitutes into microseism wave equation, uses the focus item on the right of equivalent force method the other side's journey to carry out equivalence
Discrete.
The elastic wave forward modeling method of simulation microseism the most according to claim 1, it is characterised in that in step 4, utilizes
Staggered-mesh high-order finite difference method technology, pseudo-spectrometry, FInite Element, spectral element method or iterative method solve microseism elasticity wave-wave
Dynamic equation.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610494146.5A CN106199697A (en) | 2016-06-29 | 2016-06-29 | The elastic wave forward modeling method of simulation microseism |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610494146.5A CN106199697A (en) | 2016-06-29 | 2016-06-29 | The elastic wave forward modeling method of simulation microseism |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN106199697A true CN106199697A (en) | 2016-12-07 |
Family
ID=57461566
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610494146.5A Pending CN106199697A (en) | 2016-06-29 | 2016-06-29 | The elastic wave forward modeling method of simulation microseism |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106199697A (en) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108983285A (en) * | 2018-07-19 | 2018-12-11 | 中国石油大学(北京) | A kind of a variety of source wavefield analogy methods and device based on moment tensor |
| CN109270575A (en) * | 2018-11-02 | 2019-01-25 | 河南理工大学 | A kind of attenuation of seismic waves model construction method equivalent based on building seismic response |
| CN110598367A (en) * | 2019-10-12 | 2019-12-20 | 中南大学 | Footprint-guided efficient aviation electromagnetic numerical simulation method |
| CN111794727A (en) * | 2020-07-02 | 2020-10-20 | 中国石油大学(北京) | Pump injection frequency selection method and device for pulse circulation hydraulic fracturing |
| CN114660659A (en) * | 2022-03-30 | 2022-06-24 | 中国矿业大学 | Longitudinal and transverse wave decoupling high-precision microseism positioning method |
| CN116953774A (en) * | 2023-07-06 | 2023-10-27 | 四川伟博震源科技有限公司 | Air explosion transverse wave seismic source excitation system and excitation method |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4887244A (en) * | 1988-06-28 | 1989-12-12 | Mobil Oil Corporation | Method for seismic trace interpolation using a forward and backward application of wave equation datuming |
| CN102062875A (en) * | 2010-11-30 | 2011-05-18 | 中国石油集团川庆钻探工程有限公司 | Forward modeling method for fluctuating surface elastic wave equation |
| CN102385066A (en) * | 2010-09-06 | 2012-03-21 | 中国石油天然气股份有限公司 | Pre-stack earthquake quantitative imaging method |
| CN102914799A (en) * | 2012-10-12 | 2013-02-06 | 中国石油天然气股份有限公司 | Forward modeling method and device for nonequivalent wave field |
-
2016
- 2016-06-29 CN CN201610494146.5A patent/CN106199697A/en active Pending
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4887244A (en) * | 1988-06-28 | 1989-12-12 | Mobil Oil Corporation | Method for seismic trace interpolation using a forward and backward application of wave equation datuming |
| CN102385066A (en) * | 2010-09-06 | 2012-03-21 | 中国石油天然气股份有限公司 | Pre-stack earthquake quantitative imaging method |
| CN102062875A (en) * | 2010-11-30 | 2011-05-18 | 中国石油集团川庆钻探工程有限公司 | Forward modeling method for fluctuating surface elastic wave equation |
| CN102914799A (en) * | 2012-10-12 | 2013-02-06 | 中国石油天然气股份有限公司 | Forward modeling method and device for nonequivalent wave field |
Non-Patent Citations (2)
| Title |
|---|
| HUIJIAN LI 等: "Microseismic forward modeling based on different focal mechanisms used by the seismic moment tensor and elastic wave equation", 《JOURNAL OF GEOPHYSICS AND ENGINEERING》 * |
| 张睿璇 等: "波动方程法地震波正演数值模拟研究综述", 《油气地球物理》 * |
Cited By (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108983285A (en) * | 2018-07-19 | 2018-12-11 | 中国石油大学(北京) | A kind of a variety of source wavefield analogy methods and device based on moment tensor |
| CN108983285B (en) * | 2018-07-19 | 2019-12-13 | 中国石油大学(北京) | moment tensor-based multi-seismic source wave field simulation method and device |
| CN109270575A (en) * | 2018-11-02 | 2019-01-25 | 河南理工大学 | A kind of attenuation of seismic waves model construction method equivalent based on building seismic response |
| CN110598367A (en) * | 2019-10-12 | 2019-12-20 | 中南大学 | Footprint-guided efficient aviation electromagnetic numerical simulation method |
| CN111794727A (en) * | 2020-07-02 | 2020-10-20 | 中国石油大学(北京) | Pump injection frequency selection method and device for pulse circulation hydraulic fracturing |
| CN111794727B (en) * | 2020-07-02 | 2021-06-11 | 中国石油大学(北京) | Pump injection frequency selection method and device for pulse circulation hydraulic fracturing |
| CN114660659A (en) * | 2022-03-30 | 2022-06-24 | 中国矿业大学 | Longitudinal and transverse wave decoupling high-precision microseism positioning method |
| CN114660659B (en) * | 2022-03-30 | 2023-03-31 | 中国矿业大学 | Longitudinal and transverse wave decoupling high-precision microseism positioning method |
| CN116953774A (en) * | 2023-07-06 | 2023-10-27 | 四川伟博震源科技有限公司 | Air explosion transverse wave seismic source excitation system and excitation method |
| CN116953774B (en) * | 2023-07-06 | 2024-03-12 | 四川伟博震源科技有限公司 | Air explosion transverse wave seismic source excitation system and excitation method |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106199697A (en) | The elastic wave forward modeling method of simulation microseism | |
| CN108181652B (en) | A kind of subsea node seismic data uplink and downlink wave field numerical method | |
| Vidale et al. | Finite-difference seismograms for SH waves | |
| CN103713315B (en) | A kind of seismic anisotropy parameter full waveform inversion method and device | |
| CN104570082B (en) | Extraction method for full waveform inversion gradient operator based on green function characterization | |
| Manolis et al. | Seismic wave propagation in non-homogeneous elastic media by boundary elements | |
| CN106054248A (en) | Earthquake rock physical inversion method based on large area tight reservoir | |
| CN107894618B (en) | A kind of full waveform inversion gradient preprocess method based on model smoothing algorithm | |
| CN104237937B (en) | Pre-stack seismic inversion method and system thereof | |
| CN102253415A (en) | Method for establishing earthquake response mode based on fracture equivalent medium model | |
| Golubev et al. | Simulation of seismic wave propagation in a multicomponent oil deposit model | |
| CN106324662A (en) | Full waveform inversion method and system aiming at target layer | |
| CN104122578B (en) | A kind of earthquake stimulation analogy method based on rock microstructure and dielectric resilient parameter | |
| CN104732093B (en) | A kind of FCT FDM the Forward Modelings based on disperse viscosity wave equation | |
| Schroder et al. | Elastic waves interacting with buried land mines: a study using the FDTD method | |
| WANG et al. | Multi‐Azimuth Three‐Component Surface Seismic Modeling in Cracked Monoclinic Media | |
| Wang et al. | The seismic responses and its Doppler effects of moving aircraft caused by Acoustic-Seismic Coupling | |
| Huang et al. | Study on dynamics of vibrator baseplate at low and high frequencies | |
| Li et al. | The underground explosion point measurement method based on high-precision location of energy focus | |
| EP2932304A1 (en) | Methods and systems for seismic modeling using multiple seismic source types | |
| CN104977605B (en) | High-resolution multi-scale wave equation inversion method | |
| YAN et al. | High‐Order Finite‐Difference Numerical Modeling of Wave Propagation in Viscoelastic TTI Media Using Rotated Staggered Grid | |
| Track et al. | 3-D cross-well simulations and applications | |
| Petrov et al. | Waveform inversion in Laplace-Fourier domain for estimation of seismic event location and moment tensor | |
| Liu et al. | Numerical modeling of seismic-wave propagation through fractures with nonuniform height and densityin 3d |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20161207 |
|
| WD01 | Invention patent application deemed withdrawn after publication |