CN105589100B  A kind of microseism hypocentral location and rate pattern Simultaneous Inversion method  Google Patents
A kind of microseism hypocentral location and rate pattern Simultaneous Inversion method Download PDFInfo
 Publication number
 CN105589100B CN105589100B CN201410564587.9A CN201410564587A CN105589100B CN 105589100 B CN105589100 B CN 105589100B CN 201410564587 A CN201410564587 A CN 201410564587A CN 105589100 B CN105589100 B CN 105589100B
 Authority
 CN
 China
 Prior art keywords
 msubsup
 mrow
 msup
 micro
 mfrac
 Prior art date
 Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
 Active
Links
 238000004587 chromatography analysis Methods 0.000 claims abstract description 51
 230000000875 corresponding Effects 0.000 claims description 2
 238000007689 inspection Methods 0.000 claims 1
 238000000034 method Methods 0.000 description 4
 230000004807 localization Effects 0.000 description 3
 238000001514 detection method Methods 0.000 description 2
 238000010586 diagram Methods 0.000 description 2
 238000010276 construction Methods 0.000 description 1
 230000001419 dependent Effects 0.000 description 1
 230000000694 effects Effects 0.000 description 1
 238000005516 engineering process Methods 0.000 description 1
 230000004048 modification Effects 0.000 description 1
 238000006011 modification reaction Methods 0.000 description 1
 238000000926 separation method Methods 0.000 description 1
 230000035939 shock Effects 0.000 description 1
Abstract
The invention provides a kind of microseism hypocentral location and rate pattern Simultaneous Inversion method, belong to seismic source location field in micro seismic monitoring.This method includes：(1) first arrival time of microseismic event, is picked up；(2), using Sonic Logging Data, initial velocity model is built, then obtains initial hypocentral location and the origin time of earthquake of microseismic event；(3) when the theory, based on two spots ray tracing calculating microseismic event is walked, structure joint chromatography equation group；(4) the joint chromatography equation group, is solved, calculates parameter renewal amount δ u, Δ x and Δ τ；(5) judge whether to meet end condition, if it is, being transferred to step (6), if it is not, then being updated using the parameter renewal amount to rate pattern, hypocentral location and the origin time of earthquake, be then back to step (3)；(6) final rate pattern, hypocentral location and the origin time of earthquake, are exported.
Description
Technical field
The invention belongs to seismic source location field in micro seismic monitoring, and in particular to a kind of microseism hypocentral location and rate pattern
Simultaneous Inversion method, information Simultaneous Inversion hypocentral location and rate pattern when being walked using microseismic event.
Background technology
Microseismic is development the 1990s and active new Nontraditional Techniques, and it can be with
Applied to numerous areas such as Reservoir Development, Mine Monitoring, Geological Hazards Monitorings.The position of focus, hair shake in microseismic monitoring
Moment, source level and rate pattern are all unknown, and it is the top priority of microseismic monitoring to determine these parameters.
When microseism positions, because hypocentral location and rate pattern are all unknown quantitys, both are coupled.Microseism is examined
The positioning precision of event is largely dependent on the precision of rate pattern during survey.Conventional localization method is utilized based on well logging
Data and onedimensional rate pattern after perforation corrects are positioned.It is based on when subsurface velocity model has lateral velocity variation
The positioning precision of onedimensional rate pattern will reduce, therefore just need to hypocentral location and rate pattern while carry out inverting.
Thuber (1983,1992) proposes local earthquake's chromatography method, and this method will using Local Subspace inversion technique
Focal shock parameter and speed parameter separation, so as to realize Simultaneous Inversion threedimensional velocity structure and hypocentral location.Domestic and international seismologist
It is widely applied in earthquakes location and velocity inversion.
Zhang (2003) has developed double difference on the basis of double difference positioning mode (Waldhauser and Fllsworth, 2000)
Chromatography method, this has effectively been drawn the advantage of double difference positioning and has taken into account inversion speed.The efficiency of inverse process of zhang methods is better than
Local earthquake's chromatography method.
But Zhang (2003) method ray path when construction chromatographs equation essentially coincides in shallowlayer, is unfavorable for shallow
Interval Velocity Inversion.
The content of the invention
It is an object of the invention to solve problem present in abovementioned prior art, there is provided a kind of microseism hypocentral location and
Rate pattern Simultaneous Inversion method, wave detector double difference chromatography equation is introduced, combine common chromatography equation, focus double difference chromatography side
Journey, three build joint chromatography equation.Compared with zhang double difference CT calculating method, it is beneficial to improve shallowlayer velocity inversion essence
Degree, it can further improve the precision and positioning precision of microseism inversion speed model.
The present invention is achieved by the following technical solutions：
A kind of microseism hypocentral location and rate pattern Simultaneous Inversion method, including：
(1) first arrival time of microseismic event, is picked up；
(2), using Sonic Logging Data, initial velocity model is built, then obtains the initial focus position of microseismic event
Put and the origin time of earthquake；
(3) when the theory, based on two spots ray tracing calculating microseismic event is walked, structure joint chromatography equation group；
(4) the joint chromatography equation group, is solved, calculates parameter renewal amount δ u, Δ x and Δ τ；
(5), judge whether to meet end condition, if it is, step (6) is transferred to, if it is not, then using the parameter more
New amount is updated to rate pattern, hypocentral location and the origin time of earthquake, is then back to step (3)；
(6) final rate pattern, hypocentral location and the origin time of earthquake, are exported.
What the step (3) was realized in：
Establish joint chromatography equation group, including common chromatography equation, focus double difference chromatography equation and wave detector double difference chromatography
Equation, wherein commonly shown in chromatography equation such as formula (2)：
Wherein,Residual error between when being travelled for microseismic event i to receiving point k calculating and when pickup is travelled；δ u are slow
Spend model modification amount；The respectively renewal amount in microseismic event i hypocentral locations x, y, z direction；Δτ^{i}For
Microseismic event i origin time of earthquake renewal amount；Focus partial derivativeFor the slowness vector at focus, it is：
Shown in the focus double difference chromatography equation such as formula (4)：
Wherein,Between when respectively microseismic event i, j to receiving point k calculating are travelled and when pickup is travelled
Residual error；Δτ^{j}For microseismic event j origin time of earthquake renewal amount；It is the same of two events of record for the focus double difference time
The difference of residual error during the travelling of detection：
Shown in the wave detector double difference chromatography equation such as formula (6)：
Wherein,For the wave detector double difference time, the difference of residual error during the travelling for two wave detectors for recording same event：
It is to solve the joint chromatography equation group using damping LSQR method in the step (4).
Rate pattern, hypocentral location and the origin time of earthquake are updated using the parameter renewal amount in the step (5)
It is realized in：
Corresponding renewal amount δ u, Δ x and the Δ τ that parameter in step (2) is obtained plus step (4) respectively.
End condition in the step (5) is as follows：
Residual error when averagely travelling, i.e.,Absolute value arithmetic mean of instantaneous value is less than a certain constant or parameter renewal amount is less than certain
One constant.This four parameter setting threshold values are given respectively, and this four constants are the amount of artificially given very little, can such as be used
0.00001。
Compared with prior art, the beneficial effects of the invention are as follows：
1), inverting stability is stronger.
2), compared with traditional double difference localization method, microseism seismic source location precision is high；
3), compared with double difference chromatography method, inversion result is advantageous to improve shallowlayer velocity inversion precision.
Brief description of the drawings
XY views in Fig. 1 a double differences chromatography positioning result
XZ views in Fig. 1 b double differences chromatography positioning result
YZ views in Fig. 1 c double differences chromatography positioning result
XY views in Fig. 2 a this method positioning results
XZ views in Fig. 2 b this method positioning results
YZ views in Fig. 2 c this method positioning results
Fig. 3 correct velocity models (depth 400m)
Fig. 4 double difference tomographic inversion rate patterns (depth 400m)
Fig. 5 this method inversion speed models (depth 400m)
The step block diagram of Fig. 6 this method.
Embodiment
The present invention is described in further detail below in conjunction with the accompanying drawings：
Joint chromatography equation equation group
1) equation is commonly chromatographed：During microseismic event i to receiving point k actual observation travelling
Wherein, u is slowness field (inverse of speed), τ^{i}For the origin time of earthquake.When microseism positions, focus coordinate (x_{1}, x_{2},
x_{3}), the origin time of earthquake and slowness field be all unknown parameter.
During the arrival that wave detector receives(P ripples or S ripples), using the hypocentral location of experiment, origin time of earthquake and initial
When rate pattern (prior information) carrys out computational theory arrivalBoth travelling when residual error be
Residual error when when calculating travelling and observing travellingCan be approximately linearly with the disturbance of focus and rate pattern
Wherein, focus partial derivativeFor the slowness vector at focus, it is
2) focus double difference chromatography equation：For identical wave detector k, different microseismic event i are received, j is residual during travelling
The difference of difference, can obtain equation below (Zhang, 2003)
Wherein,For focus double difference time (Waldhauser and Fllsworth, 2000), for the same of two events of record
The difference of residual error during the travelling of one detection.
3) wave detector double difference chromatography equation：For same microseismic event i, two different wave detectors, receiving point k's is another
Microseismic event j, also there is similar expression formula；Both, which subtract each other, can obtain equation below
Wherein,For the wave detector double difference time, the difference of residual error during the travelling for two wave detectors for recording same event.
Common chromatography equation, focus double difference chromatography equation and wave detector double difference chromatography equation, i.e. equation (2), (4) and (6),
Simultaneous forms final joint chromatography equation.Furthermore it is possible to give different weight coefficient w for three kinds of chromatography equations_{1}, w_{2}, w_{3},
To control the weight proportion of three.
As shown in fig. 6, the specific implementation step of the inventive method is as follows：
(1) first arrival time of microseismic event, is picked up；
(2) Sonic Logging Data, structure initial velocity model (horizontal layer), are utilized；By artificially given or other
Microseism localization method (such as Gaiger methods, grid data service etc.), obtain initial hypocentral location and the hair of microseismic event
Shake the moment；
(3) when the theory, based on two spots ray tracing calculating microseismic event is walked；Constructed according to formula (2), (4) and (6)
Joint chromatography equation (including common chromatography equation, focus double difference chromatography equation and wave detector double difference chromatography equation)；
(4) joint chromatography equation, is solved using damping LSQR method, parameter renewal amount (δ u, Δ x and Δ τ) is calculated, enters
And renewal speed model, hypocentral location and the origin time of earthquake；
(5), repeat step (3) (4), until meeting end condition, that is, residual error is less than a certain constant or parameter more when travelling
New amount is less than a certain constant；
(6), exported parameters such as final rate pattern, hypocentral location and the origin time of earthquake as output result.
It is specific as follows for examining effectiveness of the invention below with for one embodiment：
Fig. 1 a Fig. 1 c are that double difference chromatographs positioning result, respectively XY views, XZ views and YZ views, and Fig. 2 a Fig. 2 c are this
The positioning result of invention, respectively XY views, XZ views and YZ views；Comparison diagram 1a Fig. 1 c and Fig. 2 a Fig. 2 c can be seen that
The hypocentral location that inverting of the present invention obtains is more accurate.The rate pattern (depth 400m) that Fig. 4 obtains for double difference tomographic inversion, Fig. 5
The rate pattern (depth 400m) obtained for inverting of the present invention, contrast true velocity model are (as shown in Figure 3) as can be seen that this hair
It is bright to be also better than double difference chromatography method in the fine degree of inversion speed model.
The present invention relates to seismic source location method in micro seismic monitoring, is in initial velocity model inaccuracy, improves positioning accurate
Spend the method for Simultaneous Inversion rate pattern.The present invention goes out microseism hypocentral location, hair using joint chromatography equation come Simultaneous Inversion
Shake moment and rate pattern；Joint chromatography equation includes common chromatography equation, focus double difference chromatography equation and wave detector double difference
Chromatograph equation.Difference due to employing travel time residual and travel time residual simultaneously improves seismic source location precision as optimized variable
With the inversion accuracy of rate pattern, especially in the case of rate pattern inaccuracy, the stability of inversion result is also improved.
Abovementioned technical proposal is one embodiment of the present invention, for those skilled in the art, at this
On the basis of disclosure of the invention application process and principle, it is easy to make various types of improvement or deformation, be not limited solely to this
Invent the method described by abovementioned embodiment, therefore previously described mode is simply preferable, and and without limitation
The meaning of property.
Claims (4)
1. a kind of microseism hypocentral location and rate pattern Simultaneous Inversion method, it is characterised in that：Methods described includes：
(1) first arrival time of microseismic event, is picked up；
(2), using Sonic Logging Data, build initial velocity model, then obtain microseismic event initial hypocentral location and
The origin time of earthquake；
(3) when the theory, based on two spots ray tracing calculating microseismic event is walked, structure joint chromatography equation group；
(4) the joint chromatography equation group, is solved, calculates parameter renewal amount δ u, Δ x and Δ τ；
(5) judge whether to meet end condition, if it is, step (6) is transferred to, if it is not, then utilizing the parameter renewal amount
Rate pattern, hypocentral location and the origin time of earthquake are updated, are then back to step (3)；
(6) final rate pattern, hypocentral location and the origin time of earthquake, are exported；
Wherein, the step (3) is realized in：
Joint chromatography equation group, including common chromatography equation, focus double difference chromatography equation and wave detector double difference chromatography equation are established,
Wherein commonly shown in chromatography equation such as formula (2)：
<mrow>
<msubsup>
<mi>r</mi>
<mi>k</mi>
<mi>i</mi>
</msubsup>
<mo>=</mo>
<msup>
<mi>&Delta;&tau;</mi>
<mi>i</mi>
</msup>
<mo>+</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>l</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mn>3</mn>
</munderover>
<mfrac>
<mrow>
<mo>&part;</mo>
<msubsup>
<mi>T</mi>
<mi>k</mi>
<mi>i</mi>
</msubsup>
</mrow>
<mrow>
<mo>&part;</mo>
<msubsup>
<mi>x</mi>
<mi>l</mi>
<mi>i</mi>
</msubsup>
</mrow>
</mfrac>
<msubsup>
<mi>x</mi>
<mi>l</mi>
<mi>i</mi>
</msubsup>
<mo>+</mo>
<msubsup>
<mo>&Integral;</mo>
<mi>i</mi>
<mi>k</mi>
</msubsup>
<mi>&delta;</mi>
<mi>u</mi>
<mi>d</mi>
<mi>s</mi>
<mo></mo>
<mo></mo>
<mo></mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
Wherein,Residual error between when being travelled for microseismic event i to receiving point k calculating and when pickup is travelled；δ u are slowness mould
Type renewal amount；The respectively renewal amount in microseismic event i hypocentral locations x, y, z direction；Δτ^{i}To be micro
Shake event i origin time of earthquake renewal amount；Focus partial derivativeFor the slowness vector at focus, it is:
<mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<msubsup>
<mi>T</mi>
<mi>k</mi>
<mi>i</mi>
</msubsup>
</mrow>
<mrow>
<mo>&part;</mo>
<msubsup>
<mi>x</mi>
<mi>l</mi>
<mi>i</mi>
</msubsup>
</mrow>
</mfrac>
<mo>=</mo>
<mo></mo>
<mfrac>
<mn>1</mn>
<mi>V</mi>
</mfrac>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msubsup>
<mi>dx</mi>
<mi>l</mi>
<mi>i</mi>
</msubsup>
</mrow>
<mrow>
<mi>d</mi>
<mi>s</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mo></mo>
<mo></mo>
<mo></mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
Wherein, V represents speed at focus, and s represents seismic ray length,
Shown in the focus double difference chromatography equation such as formula (4)：
<mrow>
<msubsup>
<mi>dr</mi>
<mi>k</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msubsup>
<mo>=</mo>
<msubsup>
<mi>r</mi>
<mi>k</mi>
<mi>i</mi>
</msubsup>
<mo></mo>
<msubsup>
<mi>r</mi>
<mi>k</mi>
<mi>j</mi>
</msubsup>
<mo>=</mo>
<msup>
<mi>&Delta;&tau;</mi>
<mi>i</mi>
</msup>
<mo>+</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>l</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mn>3</mn>
</munderover>
<mfrac>
<mrow>
<mo>&part;</mo>
<msubsup>
<mi>T</mi>
<mi>k</mi>
<mi>i</mi>
</msubsup>
</mrow>
<mrow>
<mo>&part;</mo>
<msubsup>
<mi>x</mi>
<mi>l</mi>
<mi>i</mi>
</msubsup>
</mrow>
</mfrac>
<msubsup>
<mi>&Delta;x</mi>
<mi>l</mi>
<mi>i</mi>
</msubsup>
<mo>+</mo>
<msubsup>
<mo>&Integral;</mo>
<mi>i</mi>
<mi>k</mi>
</msubsup>
<mi>&delta;</mi>
<mi>u</mi>
<mi>d</mi>
<mi>s</mi>
<mo></mo>
<msup>
<mi>&Delta;&tau;</mi>
<mi>j</mi>
</msup>
<mo></mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>l</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mn>3</mn>
</munderover>
<mfrac>
<mrow>
<mo>&part;</mo>
<msubsup>
<mi>T</mi>
<mi>k</mi>
<mi>j</mi>
</msubsup>
</mrow>
<mrow>
<mo>&part;</mo>
<msubsup>
<mi>x</mi>
<mi>l</mi>
<mi>j</mi>
</msubsup>
</mrow>
</mfrac>
<msubsup>
<mi>&Delta;x</mi>
<mi>l</mi>
<mi>j</mi>
</msubsup>
<mo></mo>
<msubsup>
<mo>&Integral;</mo>
<mi>j</mi>
<mi>k</mi>
</msubsup>
<mi>&delta;</mi>
<mi>u</mi>
<mi>d</mi>
<mi>s</mi>
<mo></mo>
<mo></mo>
<mo></mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
Wherein,It is residual between when respectively microseismic event i, j to receiving point k calculating are travelled and when pickup is travelled
Difference；Δτ^{j}For microseismic event j origin time of earthquake renewal amount；It is the same inspection for recording two events for the focus double difference time
The difference of residual error during the travelling of ripple device：
<mrow>
<msubsup>
<mi>dr</mi>
<mi>k</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msubsup>
<mo>=</mo>
<msubsup>
<mi>r</mi>
<mi>k</mi>
<mi>i</mi>
</msubsup>
<mo></mo>
<msubsup>
<mi>r</mi>
<mi>k</mi>
<mi>j</mi>
</msubsup>
<mo>=</mo>
<msup>
<mrow>
<mo>(</mo>
<msubsup>
<mi>t</mi>
<mi>k</mi>
<mi>i</mi>
</msubsup>
<mo></mo>
<msubsup>
<mi>t</mi>
<mi>k</mi>
<mi>j</mi>
</msubsup>
<mo>)</mo>
</mrow>
<mrow>
<mi>o</mi>
<mi>b</mi>
<mi>s</mi>
</mrow>
</msup>
<mo></mo>
<msup>
<mrow>
<mo>(</mo>
<msubsup>
<mi>t</mi>
<mi>k</mi>
<mi>i</mi>
</msubsup>
<mo></mo>
<msubsup>
<mi>t</mi>
<mi>k</mi>
<mi>j</mi>
</msubsup>
<mo>)</mo>
</mrow>
<mrow>
<mi>c</mi>
<mi>a</mi>
<mi>l</mi>
</mrow>
</msup>
<mo></mo>
<mo></mo>
<mo></mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
</mrow>
Shown in the wave detector double difference chromatography equation such as formula (6)：
<mrow>
<msubsup>
<mi>dr</mi>
<mrow>
<mi>m</mi>
<mi>n</mi>
</mrow>
<mi>i</mi>
</msubsup>
<mo>=</mo>
<msubsup>
<mi>r</mi>
<mi>m</mi>
<mi>i</mi>
</msubsup>
<mo></mo>
<msubsup>
<mi>r</mi>
<mi>n</mi>
<mi>i</mi>
</msubsup>
<mo>=</mo>
<msubsup>
<mo>&Integral;</mo>
<mi>i</mi>
<mi>m</mi>
</msubsup>
<mi>&delta;</mi>
<mi>u</mi>
<mi>d</mi>
<mi>s</mi>
<mo></mo>
<msubsup>
<mo>&Integral;</mo>
<mi>i</mi>
<mi>n</mi>
</msubsup>
<mi>&delta;</mi>
<mi>u</mi>
<mi>d</mi>
<mi>s</mi>
<mo></mo>
<mo></mo>
<mo></mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
Wherein,For the wave detector double difference time, the difference of residual error during the travelling for two wave detectors for recording same event：
<mrow>
<msubsup>
<mi>dr</mi>
<mrow>
<mi>m</mi>
<mi>n</mi>
</mrow>
<mi>i</mi>
</msubsup>
<mo>=</mo>
<msubsup>
<mi>r</mi>
<mi>m</mi>
<mi>i</mi>
</msubsup>
<mo></mo>
<msubsup>
<mi>r</mi>
<mi>n</mi>
<mi>i</mi>
</msubsup>
<mo>=</mo>
<msup>
<mrow>
<mo>(</mo>
<msubsup>
<mi>t</mi>
<mi>m</mi>
<mi>i</mi>
</msubsup>
<mo></mo>
<msubsup>
<mi>t</mi>
<mi>n</mi>
<mi>i</mi>
</msubsup>
<mo>)</mo>
</mrow>
<mrow>
<mi>o</mi>
<mi>b</mi>
<mi>s</mi>
</mrow>
</msup>
<mo></mo>
<msup>
<mrow>
<mo>(</mo>
<msubsup>
<mi>t</mi>
<mi>m</mi>
<mi>i</mi>
</msubsup>
<mo></mo>
<msubsup>
<mi>t</mi>
<mi>n</mi>
<mi>i</mi>
</msubsup>
<mo>)</mo>
</mrow>
<mrow>
<mi>c</mi>
<mi>a</mi>
<mi>l</mi>
</mrow>
</msup>
<mo></mo>
<mo></mo>
<mo></mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
<mo>.</mo>
</mrow>
2. microseism hypocentral location according to claim 1 and rate pattern Simultaneous Inversion method, it is characterised in that：It is described
It is to solve the joint chromatography equation group using damping LSQR method in step (4).
3. microseism hypocentral location according to claim 2 and rate pattern Simultaneous Inversion method, it is characterised in that：It is described
What is be realized in is updated to rate pattern, hypocentral location and the origin time of earthquake using the parameter renewal amount in step (5)：
Corresponding renewal amount δ u, Δ x and the Δ τ that parameter in step (2) is obtained plus step (4) respectively.
4. microseism hypocentral location according to claim 3 and rate pattern Simultaneous Inversion method, it is characterised in that：It is described
End condition in step (5) is as follows：
Residual error when averagely travelling, i.e.,It is normal less than a certain that absolute value arithmetic mean of instantaneous value is less than a certain constant or parameter renewal amount
Number.
Priority Applications (1)
Application Number  Priority Date  Filing Date  Title 

CN201410564587.9A CN105589100B (en)  20141021  20141021  A kind of microseism hypocentral location and rate pattern Simultaneous Inversion method 
Applications Claiming Priority (1)
Application Number  Priority Date  Filing Date  Title 

CN201410564587.9A CN105589100B (en)  20141021  20141021  A kind of microseism hypocentral location and rate pattern Simultaneous Inversion method 
Publications (2)
Publication Number  Publication Date 

CN105589100A CN105589100A (en)  20160518 
CN105589100B true CN105589100B (en)  20180309 
Family
ID=55928824
Family Applications (1)
Application Number  Title  Priority Date  Filing Date 

CN201410564587.9A Active CN105589100B (en)  20141021  20141021  A kind of microseism hypocentral location and rate pattern Simultaneous Inversion method 
Country Status (1)
Country  Link 

CN (1)  CN105589100B (en) 
Families Citing this family (13)
Publication number  Priority date  Publication date  Assignee  Title 

CN106772577B (en) *  20160629  20190426  中国石油大学(华东)  Source inversion method based on microseism data and SPSA optimization algorithm 
CN108072892B (en) *  20161109  20200110  中国石油化工股份有限公司  Automatic geological structure constraint chromatography inversion method 
CN106772600B (en) *  20161221  20200828  中国科学技术大学  Doublepair doubledifference seismic positioning method and device 
CN107703540B (en) *  20170626  20180720  河海大学  A kind of microseism positioning and chromatography imaging method 
CN109212594B (en) *  20170701  20200407  中国石油化工股份有限公司  Combined positioning method for longitudinal waves and transverse waves of anisotropic medium 
CN109212593B (en) *  20170701  20200623  中国石油化工股份有限公司  Longitudinal and transverse wave combined positioning method based on multiple perforation double differences 
CN109655892B (en) *  20171011  20200407  中国石油化工股份有限公司  Combined positioning method and system of anisotropic longitudinal and transverse wave travel time ternary polynomial 
CN109655919B (en) *  20171011  20200407  中国石油化工股份有限公司  Anisotropic longitudinal and transverse wave travel time point multiplication positioning method and system 
WO2019071504A1 (en) *  20171012  20190418  南方科技大学  Twopoint ray tracing based seismic travel time tomography inversion method 
CN110716230A (en) *  20180713  20200121  中国石油化工股份有限公司  Wellground combined microseismic positioning method 
CN110764136A (en) *  20180727  20200207  中国石油化工股份有限公司  Combined positioning method for timelapse linear combination and nonlinear combination of anisotropic longitudinal and transverse waves 
CN109521467A (en) *  20181126  20190326  阳泉煤业（集团）股份有限公司  A kind of forward probe method based on projecting coal bed tunnel 
CN111221036B (en) *  20200121  20210330  中南大学  Target area seismic source positioning method and system containing unknown cavity 
Citations (2)
Publication number  Priority date  Publication date  Assignee  Title 

CN101561512A (en) *  20080418  20091021  中国石油化工股份有限公司  Multiscale crosshole SIRT tomography method 
CN102053269A (en) *  20091027  20110511  中国石油化工股份有限公司  Analysis method of speed in seismic data 

2014
 20141021 CN CN201410564587.9A patent/CN105589100B/en active Active
Patent Citations (2)
Publication number  Priority date  Publication date  Assignee  Title 

CN101561512A (en) *  20080418  20091021  中国石油化工股份有限公司  Multiscale crosshole SIRT tomography method 
CN102053269A (en) *  20091027  20110511  中国石油化工股份有限公司  Analysis method of speed in seismic data 
NonPatent Citations (2)
Title 

DoubleDifference Tomography: The Method and Its Application to the Hayward Fault，California;Haijiang Zhang;《Bulletin of the Seismological Society of America》;20031031;第93卷(第5期);全文 * 
龙门山断裂带精细速度结构的双差层析成像研究;邓文泽;《中国优秀硕士学位论文全文数据库 基础科学辑》;20140515;第二章第二节、第三章第三节 * 
Also Published As
Publication number  Publication date 

CN105589100A (en)  20160518 
Similar Documents
Publication  Publication Date  Title 

CN102239429B (en)  Multiple anisotropic parameter inversion for a tti earth model  
US7663972B2 (en)  Method for three dimensional seismic travel time tomography in transversely isotropic media  
RU2331089C2 (en)  Methods of identification of bed and well parametres by means of tomography of fresnel volume  
CN106353792B (en)  Method suitable for positioning microseismic source of hydraulic fracturing  
CN105388518B (en)  A kind of centroid frequency and earthquake inversion of quality factor method in the united well of Frequency spectrum ratio  
US7315783B2 (en)  Traveltime calculation in three dimensional transversely isotropic (3D TI) media by the fast marching method  
US7558153B2 (en)  Radial profiling of slowness: methods and apparatus for nearwellbore alteration estimation  
CN104391323B (en)  A kind of method utilizing lower wave number composition in reflected wave information inversion speed field  
CN102967882B (en)  The modeling method of the interval velocity model on stratum  
CN102937721B (en)  Limited frequency tomography method for utilizing preliminary wave travel time  
CN102841373B (en)  Microseism positioning method based on azimuth angle constraint  
CN107884822B (en)  A method of improving mining microseism seismic source location precision  
WO2017133361A1 (en)  Coal collapse column identification method and apparatus  
CN101086535A (en)  A process and program for characterising evolution of an oil reservoir over time  
CN106094032B (en)  A kind of method for building formation velocity model  
US9720130B2 (en)  Guided bayesian experimental design  
CN102841379B (en)  Method for analyzing prestack time migration and speed based on common scatter point channel set  
CN105510880A (en)  Microseism focus positioning method based on doubledifference method  
CN106405651B (en)  Full waveform inversion initial velocity model construction method based on logging matching  
CN104678435A (en)  Method for extracting Rayleigh surface wave frequency dispersion curve  
CN104536043B (en)  A kind of Depth Domain bulk velocity Model Fusion method and device  
CN104730579B (en)  A kind of joint static correcting method of ripple in length and breadth based on calculation of near surface shear velocity inverting  
WO2017167191A1 (en)  Method and device for processing seismic data  
CN103250072A (en)  Fracture characterization by interferometric drillbit imaging, time reversal imaging of fractures using drill bit seismics, and monitoring of fracture generation via time reversed acoustics and electroseismics  
CN102508303A (en)  Advanced detection method for focusing chromatography induced polarization of underground engineering 
Legal Events
Date  Code  Title  Description 

C06  Publication  
PB01  Publication  
C10  Entry into substantive examination  
SE01  Entry into force of request for substantive examination  
GR01  Patent grant  
GR01  Patent grant 