CN107167843B - More wave time-domain matching process and device - Google Patents
More wave time-domain matching process and device Download PDFInfo
- Publication number
- CN107167843B CN107167843B CN201710378360.9A CN201710378360A CN107167843B CN 107167843 B CN107167843 B CN 107167843B CN 201710378360 A CN201710378360 A CN 201710378360A CN 107167843 B CN107167843 B CN 107167843B
- Authority
- CN
- China
- Prior art keywords
- longitudinal
- wave
- seed point
- point
- error
- 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
- 238000000034 method Methods 0.000 title claims abstract description 93
- 230000008569 process Effects 0.000 title claims abstract description 34
- 238000009825 accumulation Methods 0.000 claims abstract description 153
- 238000009499 grossing Methods 0.000 claims abstract description 87
- 230000017105 transposition Effects 0.000 claims abstract description 33
- 238000012545 processing Methods 0.000 claims description 32
- 230000001174 ascending effect Effects 0.000 claims description 18
- 238000005070 sampling Methods 0.000 claims description 10
- 230000015572 biosynthetic process Effects 0.000 claims description 7
- 238000005755 formation reaction Methods 0.000 claims description 6
- 239000000203 mixture Substances 0.000 claims description 6
- 235000013399 edible fruits Nutrition 0.000 claims description 5
- 238000012512 characterization method Methods 0.000 claims 1
- 238000010408 sweeping Methods 0.000 claims 1
- 238000010586 diagram Methods 0.000 description 13
- 230000008859 change Effects 0.000 description 7
- 239000010410 layer Substances 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 2
- 230000001419 dependent effect Effects 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 238000003384 imaging method Methods 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 238000003672 processing method Methods 0.000 description 2
- 238000012952 Resampling Methods 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 230000036039 immunity Effects 0.000 description 1
- 239000011229 interlayer Substances 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000002922 simulated annealing Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 230000006641 stabilisation Effects 0.000 description 1
- 238000011105 stabilization Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000003786 synthesis reaction Methods 0.000 description 1
- 230000002194 synthesizing effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/303—Analysis for determining velocity profiles or travel times
-
- 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/303—Analysis for determining velocity profiles or travel times
- G01V1/305—Travel times
-
- 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/62—Physical property of subsurface
- G01V2210/622—Velocity, density or impedance
- G01V2210/6222—Velocity; travel time
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 embodiment of the present invention provides a kind of more wave time-domain matching process and device, wherein this method comprises: each dimension in migrated section data determines default value seed point upwards;For the migrated section data of longitudinal wave, the smoothing error accumulation results of longitudinal wave and shear wave at longitudinal seed point are obtained, the smoothing error accumulation results of all longitudinal direction seed points form the first smooth seismic data;By the first 90 degree of smooth earthquake data transposition, the smoothing error accumulated result of longitudinal wave and shear wave at lateral seed point is obtained, the smoothing error accumulated result of all transverse direction seed points forms the second smooth seismic data;By the second 90 degree of smooth earthquake data transposition, for other dimensions other than vertical and horizontal to each seed point, obtain the smoothing error accumulation results of longitudinal wave and shear wave at each seed point;The seismic data that transposition obtains is returned and is chased after, the longitudinal wave of longitudinal wave time domain and the time difference of shear wave are obtained;According to the speed ratio field of the longitudinal wave of time difference field computation longitudinal wave time domain and shear wave.
Description
Technical field
The present invention relates to seismic exploration technique field, in particular to a kind of more wave time-domain matching process and device.
Background technique
Due to containing more abundant subsurface information in multi-wave seismic data, in inverting reservoir properties and petroleuon-gas prediction etc.
Aspect shows good application prospect, and multi-wave seismic survey technology has been used widely since the nineties in last century.Phase
For the reservoir information that single longitudinal wave exploration can not accurately identify, the seismic data of more waves can save and excavate longitudinal wave and turn
Change wave reflection information, structure imaging, lithology estimation, anisotropic analysis, fluid detection and storage description in terms of
It has a wide range of applications.Wherein, the longitudinal and shear wave matching process research in multi-wave seismic data is current puzzlement multi-wave seismic data
One critical issue of processing and subsequent inverting.Because of the difference of longitudinal wave and shear wave propagation speed, conversion wave propagation time is wanted
Longitudinal wave propagation time more corresponding than its is long.And if wanting to make full use of multi-wave seismic data that must just match two kinds of waves.
Therefore, good longitudinal and shear wave high-precision matching process is that more waves is made full use of to carry out joint inversion, the important foundation of explanation and guarantor
Card, has very important significance.
In existing conventional longitudinal and shear wave time match method, since underground inhomogeneous medium can produce incident longitudinal wave
Raw different response, generates reflected P-wave and reflection wave.So the P-wave And S matching of early stage is dependent on phase in wave profile in length and breadth
As imaging characteristics, the axis of reflection on same stratum is calibrated by handling or explaining personnel in longitudinal and shear wave seismic profile, and
Compressing and converting wave profile is to corresponding longitudinal wave time domain.This method is suitable for the obvious area of reflectance signature.And layer position
Calibration process can pass through two kinds of approach and complete: the first is by handling and explaining that similar reflection is found out in the observation of personnel
Axis;Second method is to make longitudinal and shear wave synthetic seismogram according to well-log information to demarcate longitudinal and shear wave seismic profile.But
It is this matching process dependent on artificial geology experiences, matching result varies with each individual, and has biggish random and error.
Assuming that longitudinal and shear wave has similar reflective indexes and identical wavelet frequency and wavelength, Gaiser scans following formula
(1) longitudinal and shear wave similarity factor shown in carries out time match, available Mean aspect wave velocity ratio, and in turn to more wave data
Interlayer speed ratio is obtained,
Wherein PP (t1) and PS (t2) it is longitudinal wave and shear wave data.Fomel and Backus carries out PP wave using least square method
With the automatic lineups matching of PS wave time-domain, reduce the artificial interference for picking up layer position.This method needs one preferable first
Beginning model avoids convergence result from falling into local minimum.Maximum comparability objective function based on longitudinal and shear wave, J.X.Yuan etc. are utilized
Simulated annealing matches longitudinal and shear wave.The problem of this method can be to avoid local minimum, and initial model is missed
Difference has a well adapting to property.But based on this method of longitudinal and shear wave normalized crosscorrelation for Mean aspect wave velocity ratio
There is certain requirement.Gaiser points out, this optimal method based on correlation method require its average speed ratio have one more
Smooth variation tendency, variation adaptability biggish to P-S wave velocity ratio are poor.So this optimization based on cross-correlation
Method still has a series of problem, is typically necessary a prior information, and to the speed ratio change rate of data information, noise
Than etc. have higher requirement.
Summary of the invention
The embodiment of the invention provides a kind of more wave time-domain matching process, to solve the longitudinal and shear wave time in the prior art
The method of completing the square technical problem poor to variation adaptability biggish in P-S wave velocity ratio.This method comprises: obtaining longitudinal wave and turning
The migrated section data of wave are changed, the migrated section data are multi-dimensional earthquake data;In each of the migrated section data
Dimension determines default value seed point upwards, evens up longitudinal wave lineups, calculates the superposition road envelope in longitudinal wave earthquake road, in longitudinal direction
On time dimension, longitudinal seed point is that energy value reaches preset value and characterizes the reflection point of subsurface formations information;On lateral ground
On the Spatial Dimension for shaking road, pre-set space interval is kept between adjacent transverse seed point, so that all transverse direction seed point compositions
Line completely runs through in the seismic data;For the migrated section data of longitudinal wave, for each sampled point, scanning obtains the sampling
The squared magnitude of the longitudinal wave and shear wave between sampled point in the scanning range of point shear wave corresponding with itself is poor, as the sampling
The scanning error of point, the scanning range are that longitudinal wave earthquake records time span;For seed point longitudinally in each, according to longitudinal wave
Shake records time ascending direction, misses to the scanning of the sampled point between the longitudinal direction seed point and a upper longitudinal seed point
Difference is added up, and obtains the positive error accumulation results of the longitudinal direction seed point, and descending according to the longitudinal wave earthquake record time
Direction, the scanning error of the sampled point between the longitudinal direction seed point and next longitudinal seed point is added up, is somebody's turn to do
The reversed error accumulation of longitudinal seed point is as a result, by the positive error accumulation results and reversed error accumulation at the longitudinal direction seed point
Results added, then the corresponding scanning error of sampled point itself at the longitudinal direction seed point is subtracted, obtain longitudinal wave at the longitudinal direction seed point
The first smooth earthquake number is formed with the smoothing error accumulation results of the smoothing error accumulation results of shear wave, all longitudinal direction seed points
According to;It is ascending according to seismic channel on seismic survey lines for each lateral seed point by the first 90 degree of smooth earthquake data transposition
Direction, the scanning error of the sampled point between the transverse direction seed point and a upper lateral seed point is added up, is somebody's turn to do
The positive error accumulation results of lateral seed point, and the direction descending according to the seismic channel, to the transverse direction seed point under
The scanning error of sampled point between one lateral seed point is added up, and the reversed error accumulation knot of the transverse direction seed point is obtained
Fruit by the positive error accumulation results and reversed error accumulation results added at the transverse direction seed point, then subtracts the transverse direction seed
The corresponding scanning error of sampled point itself at point, obtains the smoothing error accumulation results of longitudinal wave and shear wave at the transverse direction seed point,
The smoothing error accumulation results of all transverse direction seed points form the second smooth seismic data;By the second smooth earthquake data transposition 90
Degree, for other dimensions other than vertical and horizontal to each seed point, obtain longitudinal wave and shear wave at each seed point
Smoothing error accumulation results, until obtain all dimensions of the migrated section data to each seed point at longitudinal wave and shear wave
Smoothing error accumulation results, and finally obtained smooth earthquake data transposition is obtained into the arrangement shape with the migrated section data
The consistent seismic data of formula;The seismic data obtained to transposition is returned and is chased after, obtain longitudinal wave and shear wave at each longitudinal seed point when
Between error, obtain the longitudinal wave of longitudinal wave time domain and time difference of shear wave;According to the time difference of the longitudinal wave of longitudinal wave time domain and shear wave
Calculate the longitudinal wave of longitudinal wave time domain and the speed ratio field of shear wave.
In one embodiment, for seed point longitudinally in each, time ascending direction is recorded according to longitudinal wave earthquake,
The scanning error of sampled point between the longitudinal direction seed point and a upper longitudinal seed point is added up, the longitudinal direction seed is obtained
The positive error accumulation results of point, comprising: obtain the longitudinal direction seed point and a upper longitudinal seed according to the method for linear interpolation
Point between the different time differences straight path, sampled point corresponding to different straight paths, according to longitudinal wave earthquake record the time by
It is small to arrive big direction, it is scanned the accumulative of error from first sampled point to a last sampled point, obtains the longitudinal direction seed
The positive error accumulation results of point.
In one embodiment, for seed point longitudinally in each, time descending direction is recorded according to longitudinal wave earthquake,
The scanning error of sampled point between the longitudinal direction seed point and next longitudinal seed point is added up, the longitudinal direction seed is obtained
The reversed deviation accumulation result of point, comprising: obtain the longitudinal direction seed point and next longitudinal seed according to the method for linear interpolation
Point between the different time differences straight path, sampled point corresponding to different straight paths, according to longitudinal wave earthquake record the time by
Small direction is arrived greatly, and the accumulative of error is scanned from first sampled point to a last sampled point, obtains the longitudinal direction seed
The reversed error accumulation result of point.
In one embodiment, the seismic data that transposition obtains is returned and is chased after, obtain longitudinal wave and cross at each longitudinal seed point
The time error of wave obtains the longitudinal wave of longitudinal wave time domain and the time difference of shear wave, comprising: chases after to obtain by following formula time each
The time error of longitudinal wave and shear wave at longitudinal seed point:
I=N-1, N-2 ... 1
Wherein, argmin refers to the path found in all searching routes and make d [i-1, l] the smallest u [i-1];u
[i-1] is the time error of smoothing error accumulation results smallest point corresponding longitudinal wave and shear wave at longitudinal seed point i-1, be by
It gradually returns according to time descending direction and to chase after;D [i-1, l] is at longitudinal seed point i-1 according to the scanning model
Enclose the smoothing error accumulation results being accumulated by scanning error;N is the number of longitudinal seed point;U (0:N-1) is from most
The latter longitudinal direction seed point to first longitudinal seed point reversely seek as a result, representing from the last one longitudinal kind
Son o'clock longitudinal wave of longitudinal wave time domain and time error of shear wave at first longitudinal seed point;It is indulged according at longitudinal seed point
The time error of wave and shear wave carries out linear interpolation, obtains the time error of longitudinal wave and shear wave at each sampled point, obtains longitudinal wave
The longitudinal wave of time-domain and the time difference of shear wave.
In one embodiment, by following formula according to the longitudinal wave of multidimensional migrated section data and the time difference meter of shear wave
Calculate the speed ratio field of multidimensional migrated section data:
Wherein, misfit (Tpp) be longitudinal wave time domain longitudinal wave and time difference of shear wave, γi(tpp) it is longitudinal wave time domain
The speed ratio field of longitudinal wave and shear wave.
The embodiment of the invention also provides a kind of more wave time-domain coalignments, to solve the longitudinal and shear wave time in the prior art
The matching process technical problem poor to variation adaptability biggish in P-S wave velocity ratio.The device includes: data acquisition mould
Block, for obtaining the migrated section data of longitudinal wave and converted wave, the migrated section data are multi-dimensional earthquake data;Seed point is true
Cover half block determines default value seed point for each dimension in the migrated section data upwards, evens up the same phase of longitudinal wave
Axis calculates the superposition road envelope in longitudinal wave earthquake road, and on longitudinal time dimension, longitudinal seed point is that energy value reaches preset value
And characterize the reflection point of subsurface formations information;On the Spatial Dimension of lateral seismic channel, kept between adjacent transverse seed point
Pre-set space interval, so that the line of all transverse direction seed point compositions completely runs through in the seismic data;First error is smoothly located
Module is managed, for being directed to the migrated section data of longitudinal wave, for each sampled point, it is corresponding with itself that scanning obtains the sampled point
The squared magnitude of the longitudinal wave and shear wave between sampled point in the scanning range of shear wave is poor, as the scanning error of the sampled point,
The scanning range is that longitudinal wave earthquake records time span;For seed point longitudinally in each, according to longitudinal wave earthquake record the time by
It is small to arrive big direction, the scanning error of the sampled point between the longitudinal direction seed point and a upper longitudinal seed point is added up,
The positive error accumulation results of the longitudinal direction seed point, and the direction descending according to the longitudinal wave earthquake record time are obtained, to this
The scanning error of sampled point between longitudinal seed point and next longitudinal seed point is added up, and the longitudinal direction seed point is obtained
Reversed error accumulation as a result, by the positive error accumulation results and reversed error accumulation results added at the longitudinal direction seed point, then
The corresponding scanning error of sampled point itself at the longitudinal direction seed point is subtracted, the smooth of longitudinal wave and shear wave at the longitudinal direction seed point is obtained
Error accumulation is as a result, the smoothing error accumulation results of all longitudinal direction seed points form the first smooth seismic data;Second error is flat
Sliding processing module, for by the first 90 degree of smooth earthquake data transposition, for each lateral seed point, according to seismic survey lines Shangdi
The ascending direction in road is shaken, the scanning error of the sampled point between the transverse direction seed point and a upper lateral seed point is carried out
It is accumulative, the positive error accumulation results of the transverse direction seed point, and the direction descending according to the seismic channel are obtained, to the transverse direction
The scanning error of sampled point between seed point and next lateral seed point is added up, and the reversed of the transverse direction seed point is obtained
Error accumulation is as a result, by the positive error accumulation results and reversed error accumulation results added at the transverse direction seed point, then subtracts
The corresponding scanning error of sampled point itself at the transverse direction seed point, obtains the smoothing error of longitudinal wave and shear wave at the transverse direction seed point
The smoothing error accumulation results of accumulation results, all transverse direction seed points form the second smooth seismic data;Data processing module is used
In by the second 90 degree of smooth earthquake data transposition, for other dimensions other than vertical and horizontal to each seed point, obtain
The smoothing error accumulation results of longitudinal wave and shear wave at each seed point, until obtain all dimensions of the migrated section data to
Each seed point at longitudinal wave and shear wave smoothing error accumulation results, and finally obtained smooth earthquake data transposition is obtained
With the consistent seismic data of spread pattern of the migrated section data;Time difference data processing module, for what is obtained to transposition
Seismic data is returned and is chased after, and obtains the time error of longitudinal wave and shear wave at each longitudinal seed point, obtain longitudinal wave time domain longitudinal wave and
The time difference of shear wave;Speed ratio data processing module, for vertical according to the longitudinal wave of longitudinal wave time domain and the time difference field computation of shear wave
The longitudinal wave of wave time-domain and the speed ratio field of shear wave.
In one embodiment, the first error smoothing module, comprising: positive error smoothing processing unit is used
The straight line rail of the different time differences between the longitudinal direction seed point and upper one longitudinal seed point is obtained in the method according to linear interpolation
Mark, sampled point corresponding to different straight paths are sampled according to the direction that the longitudinal wave earthquake record time is ascending from first
Point is scanned the accumulative of error to a last sampled point, obtains the positive error accumulation results of the longitudinal direction seed point.
In one embodiment, the first error smoothing module, further includes: reversed error smoothing processing unit,
The straight line of the different time differences between the longitudinal direction seed point and next longitudinal seed point is obtained for the method according to linear interpolation
Track, sampled point corresponding to different straight paths are adopted according to the direction that the longitudinal wave earthquake record time is descending from first
Sampling point is scanned the accumulative of error to a last sampled point, obtains the reversed error accumulation result of the longitudinal direction seed point.
In one embodiment, the time difference data processing module, comprising: time difference computing unit at seed point, for leading to
It crosses following formula and returns the time error for chasing after to obtain longitudinal wave and shear wave at each longitudinal seed point:
I=N-1, N-2 ... 1
Wherein, argmin refers to the path found in all searching routes and make d [i-1, l] the smallest u [i-1];u
[i-1] is the time error of smoothing error accumulation results smallest point corresponding longitudinal wave and shear wave at longitudinal seed point i-1, be by
It gradually returns according to time descending direction and to chase after;D [i-1, l] is at longitudinal seed point i-1 according to the scanning model
Enclose the smoothing error accumulation results being accumulated by scanning error;N is the number of longitudinal seed point;U (0:N-1) is from most
The latter longitudinal direction seed point to first longitudinal seed point reversely seek as a result, representing from the last one longitudinal kind
Son o'clock longitudinal wave of longitudinal wave time domain and time error of shear wave at first longitudinal seed point;The time difference calculates single at sampled point
Member, for obtaining indulging at each sampled point according to the time error of the longitudinal wave and shear wave progress linear interpolation at longitudinal seed point
The time error of wave and shear wave obtains the longitudinal wave of longitudinal wave time domain and the time difference of shear wave.
In one embodiment, the speed passes through following formula according to multidimensional migrated section data than data processing module
Longitudinal wave and shear wave time difference field computation multidimensional migrated section data speed ratio field:
Wherein, misfit (Tpp) be longitudinal wave time domain longitudinal wave and time difference of shear wave, γi(tpp) it is longitudinal wave time domain
The speed ratio field of longitudinal wave and shear wave.
In embodiments of the present invention, for multidimensional migrated section data, each ties up upward seed point, is obtained by scanning
The positive error accumulation results of longitudinal wave and shear wave at each seed point are obtained, meanwhile, pass through scanning and obtains longitudinal wave at each seed point
Reversed error accumulation with shear wave as a result, by each seed point positive error accumulation results and reversed error accumulation result tire out
Add, then subtract the corresponding scanning error of sampled point at the seed point itself, obtains the smooth of longitudinal wave and shear wave at each seed point
Error accumulation is as a result, since the selection of seed point is to maintain pre-set space interval, and due to longitudinal wave and shear wave at each seed point
Smoothing error accumulation results be by track carry out square-error accumulation obtain, the stabilization of the inverting time difference can be increased
Property;Meanwhile the speed ratio of the longitudinal wave and shear wave according to the longitudinal wave of longitudinal wave time domain and the time difference field computation longitudinal wave time domain of shear wave
, available change rate it is lesser in length and breadth with the speed ratio of shear wave.
Detailed description of the invention
The drawings described herein are used to provide a further understanding of the present invention, constitutes part of this application, not
Constitute limitation of the invention.In the accompanying drawings:
Fig. 1 is a kind of flow chart of more wave time-domain matching process provided in an embodiment of the present invention;
Fig. 2 is a kind of one-dimensional longitudinal wave (PP wave) provided in an embodiment of the present invention and shear wave (PS wave) model data;
Fig. 3 is the schematic diagram of one-dimensional model longitudinal wave time domain longitudinal and shear wave time difference provided in an embodiment of the present invention a kind of;
Fig. 4 is a kind of schematic diagram of one-dimensional model longitudinal wave time domain P-S wave velocity ratio provided in an embodiment of the present invention;
Fig. 5 is a kind of schematic illustration of smooth dynamic time adjusting method provided in an embodiment of the present invention;
Fig. 6 is a kind of schematic diagram for synthesizing longitudinal wave and converted wave poststack section provided in an embodiment of the present invention;
Fig. 7 be a kind of longitudinal wave time domain provided in an embodiment of the present invention the longitudinal and shear wave time difference and speed than field schematic diagram;
Fig. 8 is a kind of schematic diagram of longitudinal and shear wave synthetic seismic data containing noise provided in an embodiment of the present invention;
Fig. 9 is a kind of dynamic smoothing time difference longitudinal and shear wave time match side using in the application provided in an embodiment of the present invention
Schematic diagram of the longitudinal and shear wave time difference that method obtains with speed than field;
Figure 10 is a kind of schematic diagram of practical longitudinal wave poststack data provided in an embodiment of the present invention;
Figure 11 is a kind of schematic diagram of practical converted wave poststack data provided in an embodiment of the present invention;
Figure 12 is the schematic diagram of seed point distribution situation in a kind of longitudinal wave data provided in an embodiment of the present invention;
Figure 13 is a kind of dynamic smoothing time difference longitudinal and shear wave time match side using the application provided in an embodiment of the present invention
Schematic diagram of the longitudinal and shear wave time difference that method obtains with speed than field;
Figure 14 is the schematic diagram of the conversion wave seismic data after a kind of time match provided in an embodiment of the present invention.
Figure 15 is a kind of structural block diagram of more wave time-domain coalignments provided in an embodiment of the present invention.
Specific embodiment
To make the objectives, technical solutions, and advantages of the present invention clearer, right below with reference to embodiment and attached drawing
The present invention is described in further details.Here, exemplary embodiment and its explanation of the invention is used to explain the present invention, but simultaneously
It is not as a limitation of the invention.
In embodiments of the present invention, a kind of more wave time-domain matching process are provided, as shown in Figure 1, this method comprises:
Step 101: obtaining the migrated section data of longitudinal wave and converted wave, the migrated section data are multi-dimensional earthquake number
According to;
Step 102: determining default value seed point upwards in each dimension of the migrated section data, even up longitudinal wave
Lineups calculate the superposition road envelope in longitudinal wave earthquake road, and on longitudinal time dimension, longitudinal seed point is that energy value reaches pre-
If value and the reflection point for characterizing subsurface formations information;On the Spatial Dimension of lateral seismic channel, between adjacent transverse seed point
Pre-set space interval is kept, so that the line of all transverse direction seed point compositions completely runs through in the seismic data;
Step 103: for the migrated section data of longitudinal wave, for each sampled point, scanning obtains the sampled point and itself
The squared magnitude of the longitudinal wave and shear wave between sampled point in the scanning range of corresponding shear wave is poor, the scanning as the sampled point
Error, the scanning range are that longitudinal wave earthquake records time span;For seed point longitudinally in each, when being recorded according to longitudinal wave earthquake
Between ascending direction, the scanning error of the sampled point between the longitudinal direction seed point and a upper longitudinal seed point is carried out tired
Meter obtains the positive error accumulation results of the longitudinal direction seed point, and the direction descending according to the longitudinal wave earthquake record time, right
The scanning error of sampled point between the longitudinal direction seed point and next longitudinal seed point is added up, and the longitudinal direction seed point is obtained
Reversed error accumulation as a result, by the positive error accumulation results and reversed error accumulation results added at the longitudinal direction seed point,
The corresponding scanning error of sampled point itself at the longitudinal direction seed point is subtracted again, obtains the flat of longitudinal wave and shear wave at the longitudinal direction seed point
Sliding error accumulation is as a result, the smoothing error accumulation results of all longitudinal direction seed points form the first smooth seismic data;
Step 104: by the first 90 degree of smooth earthquake data transposition, for each lateral seed point, according on seismic survey lines
The ascending direction of seismic channel, to the scanning error of the sampled point between the transverse direction seed point and a upper lateral seed point into
Row is accumulative, the positive error accumulation results of the transverse direction seed point, and the direction descending according to the seismic channel is obtained, to the cross
Added up to the scanning error of the sampled point between seed point and next lateral seed point, obtains the anti-of the transverse direction seed point
To error accumulation as a result, by the positive error accumulation results and reversed error accumulation results added at the transverse direction seed point, then subtract
The corresponding scanning error of sampled point itself at the transverse direction seed point is gone, the smooth mistake of longitudinal wave and shear wave at the transverse direction seed point is obtained
The smoothing error accumulation results of poor accumulation results, all transverse direction seed points form the second smooth seismic data;
Step 105: by the second 90 degree of smooth earthquake data transposition, for other dimensions other than vertical and horizontal to
Each seed point obtains the smoothing error accumulation results of longitudinal wave and shear wave at each seed point, until obtaining the migrated section
All dimensions of data to each seed point at longitudinal wave and shear wave smoothing error accumulation results, and by it is finally obtained smoothly
Shake data transposition obtains the consistent seismic data of spread pattern with the migrated section data;
Step 106: the seismic data that transposition obtains being returned and is chased after, the time of longitudinal wave and shear wave at each longitudinal seed point is obtained
Error obtains the longitudinal wave of longitudinal wave time domain and the time difference of shear wave;
Step 107: according to the longitudinal wave and shear wave of the longitudinal wave of longitudinal wave time domain and the time difference field computation longitudinal wave time domain of shear wave
Speed ratio field.
Process as shown in Figure 1 is it is found that in embodiments of the present invention, obtain longitudinal wave and cross at each sampled point by scanning
The scanning error of wave is tieed up at each for multidimensional migrated section data and determines default value seed point upwards, pass through seed
The track that point linear interpolation determines is missed according to the ascending corresponding scanning of directional cumulation sampled point of longitudinal wave earthquake record time
Difference obtains the positive error accumulation results of seed point, meanwhile, it is adopted according to longitudinal wave earthquake record time descending directional cumulation
The corresponding scanning error of sampling point obtains the reversed error accumulation of seed point as a result, the positive error at each seed point is accumulated
As a result the scanning error for adding up with reversed error accumulation result and subtracting seed point itself, obtains longitudinal wave and cross at each seed point
The smoothing error accumulation results of wave, this process are referred to as two-way smoothing process, since the selection of seed point is to maintain default sky
Between be spaced, and due at each seed point the smoothing error accumulation results of longitudinal wave and shear wave be by carrying out square-error to track
Accumulation obtain, the stability of the inverting time difference can be increased;Meanwhile according to the time difference of the longitudinal wave of longitudinal wave time domain and shear wave
Calculate longitudinal wave time domain longitudinal wave and shear wave speed ratio field, available change rate it is smaller also more stable in length and breadth with shear wave
Speed ratio.
When it is implemented, smooth time difference in order to obtain, the dynamic time alignment of present inventor in the prior art
It is improved on the basis of method (DTW), obtains smooth dynamic time adjusting method (Smoothing Dynamic Time
Warping, SDTW).
Specifically, dynamic time alignment method principle is as follows:
For having the seismic signal f (i) and g (i) of error on two waveforms there are similitude and on the time, we can make
The time error u (0:n-1) of seismic signal f (i) and g (i) are calculated with dynamic time alignment method.
Wherein, u (0:N-1) represents the optimal solution of a time match time difference.N is the length of seismic signal f (i) He g (i)
Degree.When two consecutive points correspond to two non-conterminous points in g (i) in u [j]-u [j-1]=1, f (i).And work as u [j]-u [j-
1] two consecutive points correspond to a point in g (i) in=- 1, f (i).Because two similar seismic waves will not change very
Acutely, so the limitation of formula (2) has actual geological Significance.The process carried out below is all the limit in formula (2)
The lower progress of system.
The first step of DTW is accumulation.We define scanning error are as follows:
Wherein, l [i] indicates that the amount that the time error of two seismic signals is scanned at ith sample point, scanning range are
L, e [i, l] are the ith sample points of longitudinal wave, and corresponding to the i-th+l sampled points of converted wave, the error that the time difference is l scans knot
Fruit;
D [i, l] is the amount summed from 0 to the Sao Miao error accumulation of ith sample point.
And when one-dimensional DTW method is applied to the data of multidimensional, if being referred to as formula (4) from first sampled point
It is to the rollforward process of the N-1 sampled point, then corresponding to be missed from the last one sampled point to the scanning of first sampled point
The reversed accumulation of difference are as follows:
The second step of DTW is back to chase after,
I=N-1, N-2 ... 1 (6)
Wherein, " argmin " represents the road found in all searching routes and make d [i-1, l] the smallest u [i-1]
Diameter.This time has the heels of journey and accumulation on the contrary, carrying out from N-1 to 0.The time that the u (0:N-1) acquired represents two seismic signals misses
The optimal solution of difference.
The smooth dynamic time adjusting method principle that present inventor proposes is as follows:
We design a model, and the consistent random reflection series sequence of Depth Domain is transformed into longitudinal wave and converted wave time
Domain, and synthesis longitudinal wave as shown in Figure 2 (PP wave) and converted wave (PS wave) data are obtained with Ricker wavelet convolution,
The longitudinal and shear wave time difference of time-domain and corresponding speed ratio are as shown in Figure 3 and Figure 4.It can be seen that the time difference of longitudinal and shear wave is in longitudinal wave
Between domain be a smooth increase process.The thought " with straight Dai Qu " is used, we can be fitted this smooth with segmented linear
Time difference.Using formula (3) in DTW algorithm, it is as shown in Figure 5 to obtain its scanning error.If we expand in DTW method
The time difference range of the accumulation of each step is extended to dx sampled point from a sampled point, according to new interval dx sampled point conduct
The length of segmented linear, two endpoint is as new " sampled point " in seismic channel, we can be called seed point here.
It is consistent with DTW method, we can assume that time difference amplitude of variation be 1, from the time difference of each seed point be k from, to upper one kind
The time difference k to k-dx of son point, the different tracks defined as seed point linear interpolation (as shown in solid black lines in Fig. 5) are to k to k-
The scanning error of sampled point is accumulated within the scope of dx, and the available time difference is the seed point error accumulation result at k.Due to kind
The time difference at sub- point is the integral multiple sampling interval, and sampled point of the linear track at non-seed point is all located at greatly among grid,
But it was verified that the linear interpolation of simple adjacent mesh or directly rounding calculate and also can satisfy most of need
It asks.Finally, seed point is such as regarded as to the sparse sampled point of resampling, then using formula (6) acquire corresponding seed point when
Between error amount, and linear interpolation identical with accumulation track is finally used to obtain the time error value of each sampled point.By
In the flatness of time error, this straight path can be used as a kind of rule method to increase the stability of the inverting time difference.
This method is known as smooth dynamic time adjusting method (Smoothing Dynamic time warping), abbreviation SDTW method.
For the seismic data of multidimensional, it would be desirable to first be done to one-dimensional longitudinal wave earthquake road and corresponding converted wave seismic channel
One two-way smoothing processing, for seed point longitudinally in each, according to the direction that the longitudinal wave earthquake record time is ascending, to this
The scanning error of sampled point between longitudinal seed point and a upper longitudinal seed point is added up, and the longitudinal direction seed point is obtained
Positive error accumulation results, comprising: according to the method for linear interpolation obtain the longitudinal direction seed point and a upper longitudinal seed point it
Between the different time differences straight path, sampled point corresponding to different straight paths, according to longitudinal wave earthquake record the time by it is small to
Big direction is scanned the accumulative of error from first sampled point to a last sampled point, obtains the longitudinal direction seed point
Positive error accumulation results, i.e. realization rollforward process;
For seed point longitudinally in each, according to the direction that the longitudinal wave earthquake record time is descending, to the longitudinal direction seed point
The scanning error of sampled point between next longitudinal seed point is added up, and the reversed error for obtaining the longitudinal direction seed point is tired
Count result, comprising: when obtaining the difference between the longitudinal direction seed point and next longitudinal seed point according to the method for linear interpolation
The straight path of difference, sampled point corresponding to different straight paths record time descending direction according to longitudinal wave earthquake, from
First sampled point is scanned the accumulative of error to a last sampled point, obtains the reversed error accumulation of the longitudinal direction seed point
As a result, realizing reversely accumulative process.
For example, we are positive to carrying out from first sampled point to the scanning error of the N-1 sampled point using formula (4)
It is accumulative, obtain the positive error accumulation results of seed pointThen we seismic data can be regarded as a propagation time from
Arrive small inverse process greatly, carry out the reversed smooth accumulation of a reversed, identical seed point, i.e., using formula (5) to from
The scanning error of the last one sampled point to first sampled point is reversely added up, and the reversed error accumulation knot of seed point is obtained
FruitAnd by following formula (7) by the positive error accumulation results and reversed error accumulation results added at seed point,
The corresponding scanning error e [i, l] of sampled point at the seed point itself is subtracted again, obtains the smooth of longitudinal wave and shear wave at the seed point
Deviation accumulation is as a result, complete a two-way smoothing processing:
The corresponding scanning error e [i, l] of seed point itself is cut in formula (7) from accumulation results, is corresponded to each
Summation process of the scanning error of a sampled point, avoids the accumulation twice at each sampled point.It should be noted that
It is, although formula (7) has carried out a two-way smothing filtering, for one-dimensional data, if we are straight by sharpening result
It taps into go to return as shown in formula (7) and chases after processing, we will be apparent that have any improvement to result.This is because formula
(7) carried out only one it is one-dimensional smooth, for global optimizing method this for SDTW, substantially with an one direction
Accumulation results be not different.And in multidimensional data, this two-way smoothing method need different dimensions alternately,
This method is for the good noise immunity of multidimensional data and the accuracy of solution.
Specifically, carrying out rollforward since the application is in carrying out two-way smoothing process and reversely adding up
Cheng Zhong is first to obtain between the longitudinal direction seed point and next (or upper one) longitudinal seed point according to the method for linear interpolation
The straight path of the different time differences, and then the scanning error of the corresponding sampled point of different straight paths is added up, it is not to two
Sampled point between a seed point is added up one by one, i.e. rollforward in the application and reversed accumulative process is added to the time difference
Track restriction, it is corresponding to the different straight paths between the longitudinal seed point of longitudinal seed point and next (or upper one)
The scanning error of sampled point is added up, and the direct variation tendency of sampled point can be constrained and be consistent.
When it is implemented, obtain migrated section data all dimensions to each seed point at longitudinal wave and shear wave smooth mistake
After poor accumulation results, finally obtained smooth earthquake data transposition is obtained consistent with the spread pattern of the migrated section data
Seismic data, and the seismic data obtained using above-mentioned formula (6) to transposition is returned and is chased after, and longitudinal wave at each longitudinal seed point is obtained
Each sampling is obtained further according to the time difference of the longitudinal wave and shear wave progress linear interpolation at longitudinal seed point with the time difference of shear wave
The time difference of longitudinal wave and shear wave at point obtains the longitudinal wave of longitudinal wave time domain and the time difference of shear wave.At this point, in formula (6)
Argmin refers to the path found in all searching routes and make d [i-1, l] the smallest u [i-1];U [i-1] is at longitudinal kind
The time error of smoothing error accumulation results smallest point corresponding longitudinal wave and shear wave, is descending according to the time at sub- point i-1
Direction gradually return and to chase after;D [i-1, l] be at longitudinal seed point i-1 according to the scanning range to scanning error into
The smoothing error accumulation results that row is accumulated by;N is the number of longitudinal seed point;U (0:N-1) is from the last one longitudinal seed point
To first longitudinal seed point reversely seek as a result, representing from the last one longitudinal seed point to first longitudinal direction
The time error of the longitudinal wave of longitudinal wave time domain and shear wave at seed point.
When it is implemented, this method has no essential distinction for two and three dimensions seismic data, thus herein we with two
It is illustrated for dimension data.For the P wave data with the road M and N number of sampled point 2-D data, we first carry out it
(the flattening processing method is referred to method in the prior art, for example, Lomark etc., proposes for 2006 for flattening processing
Method;Fomel, the method proposed in 2010), and be overlapped, seek superposition road envelope, and press the condition of selected seed point
Longitudinal seed point location and its number K in selective stacking roadZ.The inverse process of sharp flattening obtains longitudinal seed point and is entirely cuing open
Distribution in the seismic data of face.Then, for seed point longitudinally in each, to each pair of seismic channel number in longitudinal and shear wave seismic data
According to an one-dimensional two-way smoothing processing (only carrying out two-way accumulation, reversed return does not chase after) is carried out, seed point longitudinally in each is obtained
Locate the smoothing error accumulation results of longitudinal wave and shear wave, the smoothing error accumulation results of all longitudinal direction seed points form three-dimensional size
For M*KZ* L sharpening result data.Sharpening result data are subjected to 90 degree of transposition processing, regard a new K asZRoad and M sampling
Point, size KZ* the seismic data error scanning result of M*L, for this new seismic data error scanning result, reference
The two-way smoothing process of initial data longitudinal direction carries out corresponding lateral two-way smooth place for each lateral seed point
Reason process obtains the smoothing error accumulated result of longitudinal wave and shear wave at each lateral seed point.It is noted herein that due to
Stacked seismic data construction feature is in lateral spread more, and poor in vertical continuity from above, generally can be with to reduce memory requirements
The lateral seed point distribution for (waiting seismic channels interval) at equal intervals is chosen, without in addition calculating its envelope size.We can be by M
A sampled point is equidistant to be divided into different line segments, chooses corresponding lateral seed point, number KX.It will finally be calculated
Sharpening result re-start one -90 degree transposition processing, obtain a size of data be KX*KZ* L three-dimensional sharpening result.
Finally, being K to final sizeX*KZ* L sharpening result data, we are to KXIt is each where a transverse direction seed point
The size of seismic channel is the smoothing error accumulated result of Kz*L, seeks its corresponding longitudinal seed point using the one-dimensional method that chases after of returning
Time difference value, obtain the time difference field data at longitudinal seed point of two dimensional cross-section, size of data KX*KZ.It was accumulated using meeting
The bilinear interpolation of journey, the time difference of each sampled point of our available initial data, size of data M*N.It can be with root
Corresponding P-S wave velocity ratio field distribution is obtained according to the time difference field data at longitudinal seed point of two dimensional cross-section.When it is implemented,
By following formula according to the speed of the longitudinal wave of multidimensional migrated section data and the time difference field computation multidimensional migrated section data of shear wave
Degree is than field:
Wherein, misfit (Tpp) be longitudinal wave time domain longitudinal wave and time difference of shear wave, γi(tpp) it is longitudinal wave time domain
The speed ratio field of longitudinal wave and shear wave.For underground medium, longitudinal P-S wave velocity ratio change rate is smaller, longitudinal and shear wave
Time difference misfit (Tpp) (longitudinal wave time domain) only corresponding P-S wave velocity ratio γi(tpp) related, thus its time difference has
Slowly varying trend.
Illustrate above-mentioned more wave time match methods below in conjunction with specific example.
Example 1. we establish depth domain model and corresponding longitudinal wave and converted wave be transformed into according to longitudinal wave and shear wave velocity
Time-domain, obtains corresponding longitudinal wave and converted waves data is as shown in Figure 6.The longitudinal and shear wave time difference of corresponding longitudinal wave time domain and speed
It is more as shown in Figure 7 than being distributed.For the stability for verifying the more wave time match methods of the application, random noise is added to initial data,
Make signal-to-noise ratio 2:1, as shown in Figure 8.Using more wave time match methods of the application, we obtain corresponding containing noise data
The longitudinal and shear wave time difference and inversion speed are than as shown in Figure 9.From the figure, it can be seen that this method has the data containing stronger noise
There is good stability, the time difference asked meets model data, and the available good speed for meeting underground layer distribution
Degree is than field result.
Example 2. we test this method using actual P wave data (Figure 10) and converted waves data after treatment
(Figure 11).Flattening processing is carried out to P wave data, obtains its seed point white wire in the distribution situation such as Figure 12 in section
It is shown.Using more wave time match methods of the application, its available longitudinal and shear wave time difference and speed are more as shown in figure 13 than field.Root
According to the longitudinal and shear wave time difference, time adjustment is carried out to converted wave, obtains result such as Figure 14, therefrom it can be seen that, longitudinal wave and converted wave can
To obtain good automatic time matching, Human disturbance is avoided.
Based on the same inventive concept, a kind of more wave time match devices are additionally provided in the embodiment of the present invention, it is such as following
Described in embodiment.Since the principle that more wave time match devices solve the problems, such as is similar to more wave time match methods, more waves
The implementation of time match device may refer to the implementation of more wave time match methods, and overlaps will not be repeated.It is following to be used
, the combination of the software and/or hardware of predetermined function may be implemented in term " unit " or " module ".Although following embodiment institute
The device of description preferably realized with software, but the combined realization of hardware or software and hardware be also may and quilt
Conception.
Figure 15 is a kind of structural block diagram of more wave time-domain coalignments of the embodiment of the present invention, and as shown in figure 15, this is more
Wave time match device includes:
Data acquisition module 1501, for obtaining the migrated section data of longitudinal wave and converted wave, the migrated section data
For multi-dimensional earthquake data;
Seed point determining module 1502 determines default value for each dimension in the migrated section data upwards
Seed point evens up longitudinal wave lineups, calculates the superposition road envelope in longitudinal wave earthquake road, on longitudinal time dimension, longitudinal seed
Point is that energy value reaches preset value and characterizes the reflection point of subsurface formations information;On the Spatial Dimension of lateral seismic channel, phase
Pre-set space interval is kept between adjacent transverse direction seed point, so that the line of all transverse direction seed point compositions is completely applied in earthquake number
In;
First error smoothing module 1503 is swept for the migrated section data for longitudinal wave for each sampled point
Retouch the squared magnitude of the longitudinal wave and shear wave between the sampled point in the scanning range for obtaining sampled point shear wave corresponding with itself
Difference, as the scanning error of the sampled point, the scanning range is that longitudinal wave earthquake records time span;For seed longitudinally in each
Point adds up the scanning error of the sampled point between the longitudinal direction seed point and a upper longitudinal seed point, obtains the longitudinal direction
The positive error accumulation results of seed point, and the direction descending according to the longitudinal wave earthquake record time, to the longitudinal direction seed point
The scanning error of sampled point between next longitudinal seed point is added up, and the reversed error for obtaining the longitudinal direction seed point is tired
Product is as a result, by the positive error accumulation results and reversed error accumulation results added at the longitudinal direction seed point, then subtracts the longitudinal direction
The corresponding scanning error of sampled point itself at seed point obtains the smoothing error of longitudinal wave and shear wave accumulation knot at the longitudinal direction seed point
The smoothing error accumulation results of fruit, all longitudinal direction seed points form the first smooth seismic data;
Second error smoothing module 1504, for by the first 90 degree of smooth earthquake data transposition, for each transverse direction
Seed point, according to the direction that seismic channel on seismic survey lines is ascending, to the transverse direction seed point and a upper lateral seed point it
Between the scanning error of sampled point added up, obtain the positive error accumulation results of the transverse direction seed point, and according to the earthquake
The descending direction in road carries out the scanning error of the sampled point between the transverse direction seed point and next lateral seed point tired
Meter obtains the reversed error accumulation of the transverse direction seed point as a result, by the positive error accumulation results at the transverse direction seed point and instead
To error accumulation results added, then the corresponding scanning error of sampled point itself at the transverse direction seed point is subtracted, obtains this lateral kind
The smoothing error accumulation results of longitudinal wave and shear wave at sub- point, it is flat that the smoothing error accumulation results of all transverse direction seed points form second
Sliding seismic data;
Data processing module 1505, for by the second 90 degree of smooth earthquake data transposition, in addition to vertical and horizontal it
Other outer dimensions to each seed point, the smoothing error accumulation results of longitudinal wave and shear wave at each seed point are obtained, until obtain
The migrated section data all dimensions to each seed point at longitudinal wave and shear wave smoothing error accumulation results, and will most
The smooth earthquake data transposition obtained eventually obtains the consistent seismic data of spread pattern with the migrated section data;
Time difference data processing module 1506, the seismic data for obtaining to transposition, which returns, to be chased after, and obtains each longitudinal seed point
The time error for locating longitudinal wave and shear wave obtains the longitudinal wave of longitudinal wave time domain and the time difference of shear wave;
Speed is than data processing module 1507, for according to the longitudinal wave of longitudinal wave time domain and the time difference field computation longitudinal wave of shear wave
The longitudinal wave of time-domain and the speed ratio field of shear wave.
In one embodiment, the first error smoothing module, comprising: positive error smoothing processing unit is used
The straight line rail of the different time differences between the longitudinal direction seed point and upper one longitudinal seed point is obtained in the method according to linear interpolation
Mark, sampled point corresponding to different straight paths are sampled according to the direction that the longitudinal wave earthquake record time is ascending from first
Point is scanned the accumulative of error to a last sampled point, obtains the positive error accumulation results of the longitudinal direction seed point.
In one embodiment, the first error smoothing module, further includes: reversed error smoothing processing unit,
The straight line of the different time differences between the longitudinal direction seed point and next longitudinal seed point is obtained for the method according to linear interpolation
Track, sampled point corresponding to different straight paths are adopted according to the direction that the longitudinal wave earthquake record time is descending from first
Sampling point is scanned the accumulative of error to a last sampled point, obtains the reversed error accumulation result of the longitudinal direction seed point.
In one embodiment, the time difference data processing module, comprising: time difference computing unit at seed point, for leading to
It crosses following formula and returns the time error for chasing after to obtain longitudinal wave and shear wave at each longitudinal seed point:
I=N-1, N-2 ... 1
Wherein, argmin refers to the path found in all searching routes and make d [i-1, l] the smallest u [i-1];u
[i-1] is the time error of smoothing error accumulation results smallest point corresponding longitudinal wave and shear wave at longitudinal seed point i-1;d[i-
It 1, l] is the smoothing error accumulation knot being accumulated by according to the scanning range to scanning error at longitudinal seed point i-1
Fruit;N is the number of longitudinal seed point;U (0:N-1) is carried out instead from the last one longitudinal seed point to first longitudinal seed point
To seeking as a result, representing the longitudinal wave from the last one longitudinal seed point longitudinal wave time domain to from first longitudinal seed point
With the time error of shear wave;Time difference computing unit at sampled point, for the time according to longitudinal wave and shear wave at longitudinal seed point
Error carry out linear interpolation, obtain the time error of longitudinal wave and shear wave at each sampled point, obtain longitudinal wave time domain longitudinal wave and
The time difference of shear wave.
In one embodiment, the speed passes through following formula according to multidimensional migrated section data than data processing module
Longitudinal wave and shear wave time difference field computation multidimensional migrated section data speed ratio field:
Wherein, misfit (Tpp) be longitudinal wave time domain longitudinal wave and time difference of shear wave, γi(tpp) it is longitudinal wave time domain
The speed ratio field of longitudinal wave and shear wave.
In embodiments of the present invention, the scanning error that longitudinal wave and shear wave at each sampled point are obtained by scanning, for more
Migrated section data are tieed up, is tieed up at each and determines default value seed point upwards, the rail determined by seed point linear interpolation
Mark is obtaining seed point just according to the corresponding scanning error of directional cumulation sampled point that the longitudinal wave earthquake record time is ascending
To error accumulation as a result, simultaneously, being missed according to the descending corresponding scanning of directional cumulation sampled point of longitudinal wave earthquake record time
Difference, obtain the reversed error accumulation of seed point as a result, by each seed point positive error accumulation results and reversed error tire out
The scanning error that product result add up and subtracts seed point itself obtains the smoothing error of longitudinal wave and shear wave accumulation at each seed point
As a result, this process is referred to as two-way smoothing process, since the selection of seed point is to maintain pre-set space interval, and due to each
The smoothing error accumulation results of longitudinal wave and shear wave are obtained by the accumulation for carrying out square-error to track at seed point, can be with
Increase the stability of the inverting time difference;Meanwhile according to the time difference field computation longitudinal wave time domain of the longitudinal wave of longitudinal wave time domain and shear wave
The speed ratio field of longitudinal wave and shear wave, available change rate it is smaller also more stable in length and breadth with the speed ratio of shear wave.
Obviously, those skilled in the art should be understood that each module of the above-mentioned embodiment of the present invention or each step can be with
It is realized with general computing device, they can be concentrated on a single computing device, or be distributed in multiple computing devices
On composed network, optionally, they can be realized with the program code that computing device can perform, it is thus possible to by it
Store and be performed by computing device in the storage device, and in some cases, can be held with the sequence for being different from herein
The shown or described step of row, perhaps they are fabricated to each integrated circuit modules or will be multiple in them
Module or step are fabricated to single integrated circuit module to realize.In this way, the embodiment of the present invention be not limited to it is any specific hard
Part and software combine.
The foregoing is only a preferred embodiment of the present invention, is not intended to restrict the invention, for the skill of this field
For art personnel, the embodiment of the present invention can have various modifications and variations.All within the spirits and principles of the present invention, made
Any modification, equivalent substitution, improvement and etc. should all be included in the protection scope of the present invention.
Claims (10)
1. a kind of more wave time-domain matching process characterized by comprising
The migrated section data of longitudinal wave and converted wave are obtained, the migrated section data are multi-dimensional earthquake data;
Default value seed point is determined upwards in each dimension of the migrated section data, is evened up longitudinal wave lineups, is calculated
The superposition road envelope in longitudinal wave earthquake road, on longitudinal time dimension, longitudinal seed point is that energy value reaches preset value and characterization
The reflection point of subsurface formations information;Default sky is kept on the Spatial Dimension of lateral seismic channel, between adjacent transverse seed point
Between be spaced so that it is all transverse direction seed point composition line completely run through in the seismic data;
For the migrated section data of longitudinal wave, for each sampled point, scanning obtains sampled point shear wave corresponding with itself
The squared magnitude of the longitudinal wave and shear wave between sampled point in scanning range is poor, described to sweep as the scanning error of the sampled point
Retouching range is that longitudinal wave earthquake records time span;It is ascending according to the longitudinal wave earthquake record time for seed point longitudinally in each
Direction, the scanning error of the sampled point between the longitudinal direction seed point and a upper longitudinal seed point is added up, is somebody's turn to do
The positive error accumulation results of longitudinal seed point, and the direction descending according to the longitudinal wave earthquake record time, to this longitudinal kind
The scanning error of sampled point between son point and next longitudinal seed point is added up, and the reversed mistake of the longitudinal direction seed point is obtained
Poor accumulation results by the positive error accumulation results and reversed error accumulation results added at the longitudinal direction seed point, then subtract this
It is tired to obtain the smoothing error of longitudinal wave and shear wave at the longitudinal direction seed point for the corresponding scanning error of sampled point itself at longitudinal seed point
Product is as a result, the smoothing error accumulation results of all longitudinal direction seed points form the first smooth seismic data;
By the first 90 degree of smooth earthquake data transposition, for each lateral seed point, according to seismic channel on seismic survey lines by it is small to
Big direction adds up the scanning error of the sampled point between the transverse direction seed point and a upper lateral seed point, obtains
The positive error accumulation results of the transverse direction seed point, and the direction descending according to the seismic channel, to the transverse direction seed point with
The scanning error of sampled point between next transverse direction seed point is added up, and the reversed error accumulation of the transverse direction seed point is obtained
As a result, by positive error accumulation results and reversed error accumulation results added at the transverse direction seed point, then subtract this lateral kind
The corresponding scanning error of sampled point itself at sub- point obtains the smoothing error of longitudinal wave and shear wave accumulation knot at the transverse direction seed point
The smoothing error accumulation results of fruit, all transverse direction seed points form the second smooth seismic data;
By the second 90 degree of smooth earthquake data transposition, for other dimensions other than vertical and horizontal to each seed point,
The smoothing error accumulation results of longitudinal wave and shear wave at each seed point are obtained, until obtaining all dimensions of the migrated section data
To each seed point at longitudinal wave and shear wave smoothing error accumulation results, and finally obtained smooth earthquake data transposition is obtained
To the consistent seismic data of spread pattern with the migrated section data;
The seismic data that transposition obtains is returned and is chased after, the time error of longitudinal wave and shear wave at each longitudinal seed point is obtained, is indulged
The longitudinal wave of wave time-domain and the time difference of shear wave;
According to the speed ratio field of the longitudinal wave and shear wave of the longitudinal wave of longitudinal wave time domain and the time difference field computation longitudinal wave time domain of shear wave.
2. more wave time-domain matching process as described in claim 1, which is characterized in that it is directed to seed point longitudinally in each, according to
Longitudinal wave earthquake records time ascending direction, to the sampled point between the longitudinal direction seed point and a upper longitudinal seed point
Scanning error is added up, and the positive error accumulation results of the longitudinal direction seed point are obtained, comprising:
The straight line of the different time differences between the longitudinal direction seed point and upper one longitudinal seed point is obtained according to the method for linear interpolation
Track, sampled point corresponding to different straight paths are adopted according to the direction that the longitudinal wave earthquake record time is ascending from first
Sampling point is scanned the accumulative of error to a last sampled point, obtains the positive error accumulation results of the longitudinal direction seed point.
3. more wave time-domain matching process as claimed in claim 2, which is characterized in that it is directed to seed point longitudinally in each, according to
Longitudinal wave earthquake records time descending direction, to the sampled point between the longitudinal direction seed point and next longitudinal seed point
Scanning error is added up, and the reversed deviation accumulation result of the longitudinal direction seed point is obtained, comprising:
The straight line of the different time differences between the longitudinal direction seed point and next longitudinal seed point is obtained according to the method for linear interpolation
Track, sampled point corresponding to different straight paths are adopted according to the direction that the longitudinal wave earthquake record time is descending from first
Sampling point is scanned the accumulative of error to a last sampled point, obtains the reversed error accumulation result of the longitudinal direction seed point.
4. more wave time-domain matching process as described in claim 1, which is characterized in that returned to the seismic data that transposition obtains
It chases after, obtains the time error of longitudinal wave and shear wave at each longitudinal seed point, obtain the longitudinal wave of longitudinal wave time domain and the time difference of shear wave
, comprising:
The time error for chasing after to obtain longitudinal wave and shear wave at each longitudinal seed point is returned by following formula:
Wherein, arg min refers to the path found in all searching routes and make d [i-1, l] the smallest u [i-1];u[i-1]
The time error of smoothing error accumulation results smallest point corresponding longitudinal wave and shear wave at longitudinal seed point i-1, be according to when
Between descending direction gradually return and chase after;D [i-1, l] is at longitudinal seed point i-1 according to the scanning range pair
The smoothing error accumulation results that scanning error is accumulated by;N is the number of longitudinal seed point;U (0:N-1) is from last
A longitudinal direction seed point to first longitudinal seed point reversely seek as a result, representing from the last one longitudinal seed point
The longitudinal wave of longitudinal wave time domain and the time error of shear wave at first longitudinal seed point;
According at longitudinal seed point the time error of longitudinal wave and shear wave carry out linear interpolation, obtain at each sampled point longitudinal wave and
The time error of shear wave obtains the longitudinal wave of longitudinal wave time domain and the time difference of shear wave.
5. more wave time-domain matching process according to any one of claims 1 to 4, which is characterized in that pass through following formula
According to the speed ratio field of the time difference field computation multidimensional migrated section data of the longitudinal wave of multidimensional migrated section data and shear wave:
Wherein, misfit (Tpp) be longitudinal wave time domain longitudinal wave and time difference of shear wave, γi(tpp) be longitudinal wave time domain longitudinal wave
With the speed ratio field of shear wave.
6. a kind of more wave time-domain coalignments characterized by comprising
Data acquisition module, for obtaining the migrated section data of longitudinal wave and converted wave, the migrated section data are for multidimensional
Shake data;
Seed point determining module determines default value seed point for each dimension in the migrated section data upwards,
Longitudinal wave lineups are evened up, the superposition road envelope in longitudinal wave earthquake road is calculated, on longitudinal time dimension, longitudinal seed point is energy
Value reaches preset value and characterizes the reflection point of subsurface formations information;On the Spatial Dimension of lateral seismic channel, adjacent transverse kind
Pre-set space interval is kept between son point, so that the line of all transverse direction seed point compositions completely runs through in the seismic data;
First error smoothing module, for being directed to the migrated section data of longitudinal wave, for each sampled point, scanning is somebody's turn to do
The squared magnitude of the longitudinal wave and shear wave between sampled point in the scanning range of sampled point shear wave corresponding with itself is poor, as this
The scanning error of sampled point, the scanning range are that longitudinal wave earthquake records time span;For seed point longitudinally in each, according to vertical
Wave earthquake record time ascending direction, to sweeping for the sampled point between the longitudinal direction seed point and a upper longitudinal seed point
It retouches error to be added up, obtains the positive error accumulation results of the longitudinal direction seed point, and according to the longitudinal wave earthquake record time by big
To small direction, the scanning error of the sampled point between the longitudinal direction seed point and next longitudinal seed point is added up, is obtained
To the longitudinal direction seed point reversed error accumulation as a result, by the positive error accumulation results and reversed error at the longitudinal direction seed point
Accumulation results are added, then subtract the corresponding scanning error of sampled point itself at the longitudinal direction seed point, are obtained at the longitudinal direction seed point
The smoothing error accumulation results of the smoothing error accumulation results of longitudinal wave and shear wave, all longitudinal direction seed points form the first smooth earthquake
Data;
Second error smoothing module, for pressing the first 90 degree of smooth earthquake data transposition for each lateral seed point
According to the direction that seismic channel on seismic survey lines is ascending, to the sampled point between the transverse direction seed point and a upper lateral seed point
Scanning error added up, obtain the positive error accumulation results of the transverse direction seed point, and descending according to the seismic channel
Direction, the scanning error of the sampled point between the transverse direction seed point and next lateral seed point is added up, is somebody's turn to do
The reversed error accumulation of lateral seed point is as a result, by the positive error accumulation results and reversed error accumulation at the transverse direction seed point
Results added, then the corresponding scanning error of sampled point itself at the transverse direction seed point is subtracted, obtain longitudinal wave at the transverse direction seed point
The second smooth earthquake number is formed with the smoothing error accumulation results of the smoothing error accumulation results of shear wave, all transverse direction seed points
According to;
Data processing module, for by the second 90 degree of smooth earthquake data transposition, for other other than vertical and horizontal
Tie up to each seed point, obtain the smoothing error accumulation results of longitudinal wave and shear wave at each seed point, until obtain it is described partially
Move all dimensions of cross-sectional data to each seed point at longitudinal wave and shear wave smoothing error accumulation results, and will be finally obtained
Smooth earthquake data transposition obtains the consistent seismic data of spread pattern with the migrated section data;
Time difference data processing module, seismic data for obtaining to transposition are returned and are chased after, obtain at each longitudinal seed point longitudinal wave and
The time error of shear wave obtains the longitudinal wave of longitudinal wave time domain and the time difference of shear wave;
Speed ratio data processing module, for according to the time difference field computation longitudinal wave time domain of the longitudinal wave and shear wave of longitudinal wave time domain
The speed ratio field of longitudinal wave and shear wave.
7. more wave time-domain coalignments as claimed in claim 6, which is characterized in that the first error smoothing processing mould
Block, comprising:
Positive error smoothing processing unit obtains the longitudinal direction seed point and upper one longitudinal kind for the method according to linear interpolation
The straight path of the different time differences between son point, sampled point corresponding to different straight paths record the time according to longitudinal wave earthquake
Ascending direction is scanned the accumulative of error from first sampled point to a last sampled point, obtains this longitudinal kind
The positive error accumulation results of son point.
8. more wave time-domain coalignments as claimed in claim 7, which is characterized in that the first error smoothing processing mould
Block, further includes:
Reversed error smoothing processing unit, obtains the longitudinal direction seed point and next longitudinal kind for the method according to linear interpolation
The straight path of the different time differences between son point, sampled point corresponding to different straight paths record the time according to longitudinal wave earthquake
Descending direction is scanned the accumulative of error from first sampled point to a last sampled point, obtains this longitudinal kind
The reversed error accumulation result of son point.
9. more wave time-domain coalignments as claimed in claim 6, which is characterized in that the time difference data processing module, packet
It includes:
Time difference computing unit at seed point chases after to obtain longitudinal wave and shear wave at each longitudinal seed point for returning by following formula
Time error:
Wherein, arg min refers to the path found in all searching routes and make d [i-1, l] the smallest u [i-1];u[i-1]
The time error of smoothing error accumulation results smallest point corresponding longitudinal wave and shear wave at longitudinal seed point i-1, be according to when
Between descending direction gradually return and chase after;D [i-1, l] is at longitudinal seed point i-1 according to the scanning range pair
The smoothing error accumulation results that scanning error is accumulated by;N is the number of longitudinal seed point;U (0:N-1) is from last
A longitudinal direction seed point to first longitudinal seed point reversely seek as a result, representing from the last one longitudinal seed point
The longitudinal wave of longitudinal wave time domain and the time error of shear wave at first longitudinal seed point;
Time difference computing unit at sampled point, for basis, the time error of longitudinal wave and shear wave is linearly inserted at longitudinal seed point
Value, obtains the time error of longitudinal wave and shear wave at each sampled point, obtains the longitudinal wave of longitudinal wave time domain and the time difference of shear wave.
10. more wave time-domain coalignments as described in any one of claim 6 to 9, which is characterized in that the speed ratio number
It is cutd open by following formula according to the time difference field computation multidimensional offset of the longitudinal wave and shear wave of multidimensional migrated section data according to processing module
The speed ratio field of face data:
Wherein, misfit (Tpp) be longitudinal wave time domain longitudinal wave and time difference of shear wave, γi(tpp) be longitudinal wave time domain longitudinal wave
With the speed ratio field of shear wave.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710378360.9A CN107167843B (en) | 2017-05-25 | 2017-05-25 | More wave time-domain matching process and device |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710378360.9A CN107167843B (en) | 2017-05-25 | 2017-05-25 | More wave time-domain matching process and device |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107167843A CN107167843A (en) | 2017-09-15 |
CN107167843B true CN107167843B (en) | 2019-05-07 |
Family
ID=59821661
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710378360.9A Active CN107167843B (en) | 2017-05-25 | 2017-05-25 | More wave time-domain matching process and device |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107167843B (en) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108426689B (en) * | 2017-12-02 | 2020-05-12 | 青岛理工大学 | Earthquake motion selection method based on earthquake resistance specification |
CN108363099A (en) * | 2018-02-11 | 2018-08-03 | 中国石油化工股份有限公司 | The longitudinal and shear wave correlation technique analyzed based on simulated annealing and multiple dimensioned ridge ripple |
CN110297269B (en) * | 2018-03-23 | 2020-12-04 | 中国石油化工股份有限公司 | Seismic data bidirectional prediction interpolation method based on improved BP neural network |
CN108761533B (en) * | 2018-05-17 | 2019-10-11 | 中国石油天然气集团有限公司 | A kind of method, apparatus and system of determining P-S wave velocity ratio |
CN110618448B (en) * | 2018-06-19 | 2021-02-09 | 中国石油化工股份有限公司 | Multi-wave dynamic matching method and system |
CN111596347B (en) * | 2019-02-21 | 2023-08-22 | 中国石油天然气集团有限公司 | Method and device for rapidly obtaining surface layer longitudinal and transverse wave speed ratio |
US11933930B2 (en) | 2021-02-25 | 2024-03-19 | Saudi Arabian Oil Company | Automatic image registration of multicomponent seismic data |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4686657A (en) * | 1984-09-24 | 1987-08-11 | Mobil Oil Corporation | Matching P wave and shear wave sections in seismic exploration |
CN103091709A (en) * | 2012-12-25 | 2013-05-08 | 中国石油天然气集团公司 | Method and device for acquiring longitudinal wave and converted wave seismic data time matching relation |
CN103472482A (en) * | 2013-09-03 | 2013-12-25 | 中国石油天然气集团公司 | Multi-wave seismic data time domain matching method and system based on gene ordering system |
CN104570079A (en) * | 2013-10-29 | 2015-04-29 | 中国石油化工股份有限公司 | Time matching method of longitudinal wave and converted shear wave seismic data |
CN105182420A (en) * | 2015-10-13 | 2015-12-23 | 中国石油天然气集团公司 | Dynamic-matching dynamic-correction method |
CN106646611A (en) * | 2016-12-08 | 2017-05-10 | 中国石油天然气集团公司 | Longitudinal wave and converted wave seismic data time matching method and device |
-
2017
- 2017-05-25 CN CN201710378360.9A patent/CN107167843B/en active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4686657A (en) * | 1984-09-24 | 1987-08-11 | Mobil Oil Corporation | Matching P wave and shear wave sections in seismic exploration |
CN103091709A (en) * | 2012-12-25 | 2013-05-08 | 中国石油天然气集团公司 | Method and device for acquiring longitudinal wave and converted wave seismic data time matching relation |
CN103472482A (en) * | 2013-09-03 | 2013-12-25 | 中国石油天然气集团公司 | Multi-wave seismic data time domain matching method and system based on gene ordering system |
CN104570079A (en) * | 2013-10-29 | 2015-04-29 | 中国石油化工股份有限公司 | Time matching method of longitudinal wave and converted shear wave seismic data |
CN105182420A (en) * | 2015-10-13 | 2015-12-23 | 中国石油天然气集团公司 | Dynamic-matching dynamic-correction method |
CN106646611A (en) * | 2016-12-08 | 2017-05-10 | 中国石油天然气集团公司 | Longitudinal wave and converted wave seismic data time matching method and device |
Non-Patent Citations (1)
Title |
---|
"P-SV与P-P波的时间匹配";陈双全 等;《石油地球物理勘探》;20100831;第45卷(第4期);第497-501页 |
Also Published As
Publication number | Publication date |
---|---|
CN107167843A (en) | 2017-09-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107167843B (en) | More wave time-domain matching process and device | |
Kosloff et al. | Velocity and interface depth determination by tomography of depth migrated gathers | |
US5930730A (en) | Method and apparatus for seismic signal processing and exploration | |
US7379854B2 (en) | Method of conditioning a random field to have directionally varying anisotropic continuity | |
CN102109616B (en) | Wave impedance inversion method for deposition cycle restriction | |
CN108802812A (en) | Well-seismic fusion stratum lithology inversion method | |
CN102053263B (en) | Method for inspecting surface structure | |
CN102053270A (en) | Sedimentary formation unit-based seismic facies analysis method | |
CN103033846A (en) | Seismic inversion system and seismic inversion method controlled by geologic facies | |
Shragge et al. | Low-frequency ambient distributed acoustic sensing (DAS): case study from Perth, Australia | |
CN107065013B (en) | A kind of interval velocity under earthquake scale determines method and device | |
CN109188520A (en) | Thin reservoir thickness prediction method and device | |
CN104570116A (en) | Geological marker bed-based time difference analyzing and correcting method | |
CA2387760A1 (en) | Transfer function method of seismic signal processing and exploration | |
CN104216013B (en) | The method of C3 coherent body based on wide-azimuth data | |
CN111487692A (en) | Method for predicting seismic response characteristics and reservoir thickness of salt shale oil rhythm layer | |
CN111624649B (en) | Method and device for establishing speed model of transverse speed change layer by using zero offset VSP | |
EA030770B1 (en) | System and method for seismic adaptive optics | |
CN106842297A (en) | Borehole restraint unstable state method for correcting phase | |
CN107272056B (en) | A method of initial model is established based on the stratum Duo Jing transformation factor | |
Bakulin et al. | Seismic imaging of vertical array data acquired using smart DAS uphole acquisition system | |
CN109143335B (en) | Method, system, medium and apparatus for producing synthetic seismic record | |
US6965830B1 (en) | System for estimating thickness of thin subsurface strata | |
CN107678065B (en) | The guarantor for improving seismic resolution constructs well control space the Method of Deconvolution and device | |
CN108363739A (en) | A kind of seismic data low-and high-frequency expanding method based on sparse acquisition |
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 |