共查询到20条相似文献，搜索用时 31 毫秒
1.
在具有垂直对称轴横向各向同性介质中，利用四种参数来确定中间至远偏移距转换波（C-波）动校正。它们是C-波叠加速度Vc2，垂直速度比和有效速度比γ0和γeff以及各向异性参数χeff。我们将这四种参数作为C波叠加速度模型。C-波速度分析的目的就是确定这种叠加速度模型。C-波叠加速度模型Vc2，γ0，γeff，和χeff可以由P-波和C-波反射动校正资料获得。然而错误的传播是C-波反射动校正反演中的严重问题。当前短排列叠加速度由于是从双曲线动校正推算而得，因而其精度不足以为各向异性参数提供有意义的反演值。中间偏移非双曲线动校正不再被人们所勿略，而是可以用一个背景γ加以量化。非双曲线分析通过中间偏移距的γ校正量可以产生Vc2，若数据不含燥音，其误差小于1％。方法稳健，允许γ启始假定值的误差达20％。该方法也适用垂直非均匀各向异性介质。精度的提高使能够用4分量地震资料计算各向异性参数。为此提出了两种工作流程：双扫描和单扫描流程。理论数据和实际数据的应用表明这两种流程得出的结果其精度相似，但是单扫描流程比双扫描更有效。 相似文献
2.
We have developed new basic theories for calculating the conversion point and the travel time of the P-SV converted wave (C-wave)
in anisotropic, inhomogeneous media. This enables the use of conventional procedures such as semblance analysis, Dix-type
model building and Kirchhoff summation, to implement anisotropic processing, and makes anisotropic processing affordable.
Here we present these new developments in two parts: basic theory and application to velocity analysis and parameter estimation.
This part deals with the basic theory, including both conversion-point calculation and moveout analysis.
Existing equations for calculating the PS-wave (C-wave) conversion point in layered media with vertical transverse isotropy
(VTI) are strictly limited to offsets about half the reflector depth (an offset-depth ratio, xlz, of 0.5), and those for calculating the C-wave traveltimes are limited to offsets equal to the reflector depth (x/z=l.0). In contrast, the new equations for calculating the conversion-point extend into offsets about three-times the reflector
depth (x/z=3.0), those for calculating the C-wave traveltimes extend into offsets twice the reflector depth (x/z=2.0). With the improved accuracy, the equations can help in C-wave data processing and parameter estimation in anisotropic,
inhomogeneous media.
This work is funded by the Edinburgh Anisotropy Project (EAP) of the British Geological Survey.
First author:
Xiangyang Li, Mr. Li is currently a professorial research seismologist (Grade 6) and technical director of the Edinburgh Anisotropy Project
in the British Geological Survey. He also holds a honorary professorship in multicomponent seismology at the School of Geosciences,
University of Edinburgh. He received his BSc(1982) in Geophysics from Changchun Geological Institute, China, an MSc (1984)
in applied geophysics from East China Petroleum Institute (now known as the China University of Petroleum), and a PhD (1992)
in seismology from the University of Edinburgh. During 1984–1987, he worked as a lecturer with the East China Petroleum Institute.
Since 1991, he has been employed by the British Geological Survey. His research interests include seismic anisotropy and multicomponent
seismology. 相似文献
3.
Several parameters are needed to describe the converted-wave （C-wave） moveout in processing multi-component seismic data, because of asymmetric raypaths and anisotropy. As the number of parameters increases, the converted wave data processing and analysis becomes more complex. This paper develops a new moveout equation with two parameters for C-waves in vertical transverse isotropy （VTI） media. The two parameters are the C-wave stacking velocity （Vc2） and the squared velocity ratio （7v,i） between the horizontal P-wave velocity and C-wave stacking velocity. The new equation has fewer parameters, but retains the same applicability as previous ones. The applicability of the new equation and the accuracy of the parameter estimation are checked using model and real data. The form of the new equation is the same as that for layered isotropic media. The new equation can simplify the procedure for C-wave processing and parameter estimation in VTI media, and can be applied to real C-wave processing and interpretation. Accurate Vc2 and Yvti can be deduced from C-wave data alone using the double-scanning method, and the velocity ratio model is suitable for event matching between P- and C-wave data. 相似文献
4.
Common‐conversion‐point binning associated with converted‐wave (C‐wave) processing complicates the task of parameter estimation, especially in anisotropic media. To overcome this problem, we derive new expressions for converted‐wave prestack time migration (PSTM) in anisotropic media and illustrate their applications using both 2D and 3D data examples. The converted‐wave kinematic response in inhomogeneous media with vertical transverse isotropy is separated into two parts: the response in horizontally layered vertical transverse isotrophy media and the response from a point‐scatterer. The former controls the stacking process and the latter controls the process of PSTM. The C‐wave traveltime in horizontally layered vertical transverse isotrophy media is determined by four parameters: the C‐wave stacking velocity VC2, the vertical and effective velocity ratios γ0 and γeff, and the C‐wave anisotropic parameter χeff. These four parameters are referred to as the C‐wave stacking velocity model. In contrast, the C‐wave diffraction time from a point‐scatterer is determined by five parameters: γ0, VP2, VS2, ηeff and ζeff, where ηeff and ζeff are, respectively, the P‐ and S‐wave anisotropic parameters, and VP2 and VS2 are the corresponding stacking velocities. VP2, VS2, ηeff and ζeff are referred to as the C‐wave PSTM velocity model. There is a one‐to‐one analytical link between the stacking velocity model and the PSTM velocity model. There is also a simple analytical link between the C‐wave stacking velocities VC2 and the migration velocity VCmig, which is in turn linked to VP2 and VS2. Based on the above, we have developed an interactive processing scheme to build the stacking and PSTM velocity models and to perform 2D and 3D C‐wave anisotropic PSTM. Real data applications show that the PSTM scheme substantially improves the quality of C‐wave imaging compared with the dip‐moveout scheme, and these improvements have been confirmed by drilling. 相似文献
5.
Stacking velocity V C2, vertical velocity ratio γ 0, effective velocity ratio γ eff, and anisotropic parameter χ eff are correlated in the PS-converted-wave (PS-wave) anisotropic prestack Kirchhoff time migration (PKTM) velocity model and are thus difficult to independently determine. We extended the simplified two-parameter (stacking velocity V C2 and anisotropic parameter k eff) moveout equation from stacking velocity analysis to PKTM velocity model updating and formed a new four-parameter (stacking velocity V C2, vertical velocity ratio γ 0, effective velocity ratio γ eff, and anisotropic parameter k eff) PS-wave anisotropic PKTM velocity model updating and process flow based on the simplified two-parameter moveout equation. In the proposed method, first, the PS-wave two-parameter stacking velocity is analyzed to obtain the anisotropic PKTM initial velocity and anisotropic parameters; then, the velocity and anisotropic parameters are corrected by analyzing the residual moveout on common imaging point gathers after prestack time migration. The vertical velocity ratio γ 0 of the prestack time migration velocity model is obtained with an appropriate method utilizing the P- and PS-wave stacked sections after level calibration. The initial effective velocity ratio γ eff is calculated using the Thomsen (1999) equation in combination with the P-wave velocity analysis; ultimately, the final velocity model of the effective velocity ratio γ eff is obtained by percentage scanning migration. This method simplifies the PS-wave parameter estimation in high-quality imaging, reduces the uncertainty of multiparameter estimations, and obtains good imaging results in practice. 相似文献
6.
Song-Yan Song Xue-Song Zhou Chun-Yong Wang Xian-Kang Zhang Jian-Li Song Yi Gong 《地震学报(英文版)》1997,10(1):15-25
On the basis of S wave information from Tai’an-Xinzhou DSS profile and with reference to the results from P-wave interpretation,
the 2-D structures, including S-wave velocity V
s, ratio γ between V
p and V
s; and Poisson’ s ratio σ, are calculated; the structural configuration of the profile is presented and the relevant inferences are drawn from the
above results. Upwarping mantle districts (V
s≈4.30 km/s) and sloping mantle districts (V
s≈4.50 km/s) of the profile with velocity difference about −4% at the top of upper mantle are divided according to the differences
of V
s, γ and σ in different media and structures, also with reference to the information of their neighbouring regions; the existence of
Niujiaqiao-Dongwang high-angle ultra-crustal fault zone is reaffirmed; the properties of low and high velocity blocks (zones)
including the crust-mantle transitionalzone and the boudary indicators of North China rift valley are discussed. A comprehensive
study is conducted on the relation of the interpretation results with earthquakes. It is concluded that the mantle upwarps,
thermal material upwells through the high-angle fault, the primary hypocenter was located at the crust-mantle juncture 30.0∼33.0
km deep, and additional stress excited the M
S=6.8 and M
S=7.2 earthquakes at specific locations around 9.0 km below Niujiaqiao-Dongwang, the earthquakes took place around the high-angle
ultra-crustal fault and centered in the brittle media and rock strata with low γ and low σ values.
This subject is part of the 85-907-02 key project during the “8th Five-Year Plan” from the State Science and Technology Commission. 相似文献
7.
It has been shown in the past that the interval-NMO velocity and the non-ellipticity parameter largely control the P-wave reflection time moveout of VTI media. To invert for these two parameters, one needs either reasonably large offsets, or some structure in the subsurface in combination with relatively mild lateral velocity variation.This paper deals with a simulation of an inversion approach, building on the assumption that accurately measured V
NMO, as defined by small offset asymptotics for a particular reflector, were available. Instead of such measurements we take synthetically computed data. First, an isotropic model is constructed which explains these V
NMO. Subsequently, residual moveout in common image gathers is modelled by ray tracing (replacing real data), along with its sensitivity for changes in the interval-NMO velocity and the non-ellipticity parameter under the constraint that V
NMO is preserved. This enables iterative updating of the non-ellipticity parameter and the interval-NMO velocity in a layer that can be laterally inhomogeneous.This approach is successfully applied for a mildly dipping reflector at the bottom of a layer with laterally varying medium parameters. With the exact V
NMO assumed to be given, lateral inhomogeneity and anisotropy can be distinguished for such a situation. However, for another example with a homogeneous VTI layer overlying a curved reflector with dip up to 30°, there appears to be an ambiguity which can be understood by theoretical analysis. Consistently with existing theory using the NMO-ellipse, the presented approach is successfully applied to the latter example if V
NMO in the strike direction is combined with residual moveout in dip direction. 相似文献
8.
A major complication caused by anisotropy in velocity analysis and imaging is the uncertainty in estimating the vertical velocity and depth scale of the model from surface data. For laterally homogeneous VTI (transversely isotropic with a vertical symmetry axis) media above the target reflector, P‐wave moveout has to be combined with other information (e.g. borehole data or converted waves) to build velocity models for depth imaging. The presence of lateral heterogeneity in the overburden creates the dependence of P‐wave reflection data on all three relevant parameters (the vertical velocity VP0 and the Thomsen coefficients ε and δ) and, therefore, may help to determine the depth scale of the velocity field. Here, we propose a tomographic algorithm designed to invert NMO ellipses (obtained from azimuthally varying stacking velocities) and zero‐offset traveltimes of P‐waves for the parameters of homogeneous VTI layers separated by either plane dipping or curved interfaces. For plane non‐intersecting layer boundaries, the interval parameters cannot be recovered from P‐wave moveout in a unique way. Nonetheless, if the reflectors have sufficiently different azimuths, a priori knowledge of any single interval parameter makes it possible to reconstruct the whole model in depth. For example, the parameter estimation becomes unique if the subsurface layer is known to be isotropic. In the case of 2D inversion on the dip line of co‐orientated reflectors, it is necessary to specify one parameter (e.g. the vertical velocity) per layer. Despite the higher complexity of models with curved interfaces, the increased angle coverage of reflected rays helps to resolve the trade‐offs between the medium parameters. Singular value decomposition (SVD) shows that in the presence of sufficient interface curvature all parameters needed for anisotropic depth processing can be obtained solely from conventional‐spread P‐wave moveout. By performing tests on noise‐contaminated data we demonstrate that the tomographic inversion procedure reconstructs both the interfaces and the VTI parameters with high accuracy. Both SVD analysis and moveout inversion are implemented using an efficient modelling technique based on the theory of NMO‐velocity surfaces generalized for wave propagation through curved interfaces. 相似文献
9.
W.J. VETTER 《Geophysical Prospecting》1987,35(6):700-717
Since the important contributions of Dürbaum and Dix, 30 years ago, velocity profile estimation procedures on horizontally layered and vertically heterogeneous media from seismic probing data have been based largely on hyperbolic moveout models and RMS and stacking velocity concepts. Re-examination of the fundamentals reveals that quantitative velocity heterogeneity and canonical valocity profiles have been implicit factors for moveout modelling and for profile inversion in the use of the Dix procedure. Heterogeneity h is the ratio (and vRMS the geometric or harmonic mean) of the path-average and time-average velocities for a raypath or, in a more restricted sense, for the normal ray belonging to a velocity profile. The canonical profile for a given velocity profile or profile segment is a moveout-equivalent monotonically increasing ramp-like profile. The ramp or constant gradient in depth is the simplest velocity profile approximator which can explicitly accommodate velocity heterogeneity. A ramp model structure is detailed which facilitates moveout simulation and model parameter estimation, and the parametric effects are explored. The horizontal offset range is quantified for which this model can give good moveout approximations. 相似文献
10.
11.
Xiangyang Li Jianxin Yuan British Geological Survey West Mains Road Edinburgh EH LA UK. Formerly at British Geological Survey now at PGS Inc. Richmond Avenue Suite Houston TX USA. 《应用地球物理》2005,2(1)
我们业已研发了计算各向异性、非均质介质中P- SV转换波(C-波)的转换点和旅行时的新理论。据此 可以利用诸如相似性分析、迪克斯模型建模、克契 霍夫求和等常规方法来完成各向异性的处理和各向 异性处理,并使各向异性的处理成为可能。这里将 我们的新发展分作两部分来介绍。第一部分为理 论,第二部分为对速度分析和参数计算的应用。第 一部分理论包括转换点的计算和动校正的分析。 相似文献
12.
Anisotropic P-wave attenuation measured from a multi-azimuth surface seismic reflection survey 总被引：2，自引：0，他引：2
Roger A. Clark Philip M. Benson rew J. Carter Carlos A. Guerrero Moreno 《Geophysical Prospecting》2009,57(5):835-845
A system of aligned vertical fractures produces azimuthal variations in stacking velocity and amplitude variation with offset, characteristics often reported in seismic reflection data for hydrocarbon exploration. Studies of associated attenuation anisotropy have been mostly theoretical, laboratory or vertical seismic profiling based. We used an 11 common‐midpoint‐long portion of each of four marine surface‐seismic reflection profiles, intersecting each other at 45° within circa 100 m of a common location, to measure the azimuthal variation of effective attenuation, Q−1eff and stacking velocity, in a shallow interval, about 100 m thick, in which consistently orientated vertical fracturing was expected due to an underlying salt diapirism. We found qualitative and quantitative consistency between the azimuthal variation in the attenuation and stacking velocity, and published amplitude variation with offset results. The 135° azimuth line showed the least apparent attenuation (1000 Q−1eff= 16 ± 7) and the fastest stacking velocity, hence we infer it to be closest to the fracture trend: the orthogonal 45° line showed the most apparent attenuation (1000Q−1eff= 52 ± 15) and slowest stacking velocity. The variation of Q−1eff with azimuth φ is well fitted by 1000Q−1eff = 34 − 18cos[2(φ+40°)] giving a fracture direction of 140 ± 23° (±1SD, derived from ‘bootstrapping’ fits to all 114 combinations of individual common‐midpoint/azimuth measurements), compared to 134 ± 47° from published amplitude variation with offset data. The effects of short‐window spectral estimation and choices of spectral ratio bandwidth and offset ranges used in attenuation analysis, individually give uncertainties of up to ±13° in fracture direction. This magnitude of azimuthal variation can be produced by credible crack geometries (e.g., dry cracks, radius 6.5 m, aspect ratio 3 × 10−5, crack density 0.2) but we do not claim these to be the actual properties of the interval studied, because of the lack of well control (and its consequences for the choice of theoretical model and host rock physical properties) and the small number of azimuths available here. 相似文献
13.
Reverse time migration of CMP-gathers an effective tool for the determination of interval velocities
One of the most important steps in the conventional processing of reflection seismic data is common midpoint (CMP) stacking. However, this step has considerable deficiencies. For instance the reflection or diffraction time curves used for normal moveout corrections must be hyperbolae. Furthermore, undesirable frequency changes by stretching are produced on account of the dependence of the normal moveout corrections on reflection times. Still other drawbacks of conventional CMP stacking could be listed.One possibility to avoid these disadvantages is to replace conventional CMP stacking by a process of migration to be discussed in this paper. For this purpose the Sherwood-Loewenthal model of the exploding reflector has to be extended to an exploding point model with symmetry to the lineP
EX
M whereP
EX
is the exploding point, alias common reflection point, andM the common midpoint of receiver and source pairs.Kirchhoff summation is that kind of migration which is practically identical with conventional CMP stacking with the exception that Kirchhoff summation provides more than one resulting trace.In this paper reverse time migration (RTM) was adopted as a tool to replace conventional CMP stacking. This method has the merit that it uses the full wave equation and that a direct depth migration is obtained, the velocityv can be any function of the local coordinatesx, y, z. Since the quality of the reverse time migration is highly dependent on the correct choice of interval velocities such interval velocities can be determined stepwise from layer to layer, and there is no need to compute interval velocities from normal moveout velocities by sophisticated mathematics or time consuming modelling. It will be shown that curve velocity interfaces do not impair the correct determination of interval velocities and that more precise velocity values are obtained by avoiding or restricting muting due to non-hyperbolic normal moveout curves.Finally it is discussed how in the case of complicated structures the reverse time migration of CMP gathers can be modified in such a manner that the combination of all reverse time migrated CMP gathers yields a correct depth migrated section. This presupposes, however, a preliminary data processing and interpretation. 相似文献
14.
Velocity model updating in prestack Kirchhoff time migration for PS converted waves: Part I – Theory
A velocity model updating approach is developed based on moveout analysis of the diffraction curve of PS converted waves in prestack Kirchhoff time migration. The diffraction curve can be expressed as a product of two factors: one factor depending on the PS converted‐wave velocity only, and the other factor depending on all parameters. The velocity‐dependent factor represents the hyperbolic behaviour of the moveout and the other is a scale factor that represents the non‐hyperbolic behaviour of the moveout. This non‐hyperbolic behaviour of the moveout can be corrected in prestack Kirchhoff time migration to form an inverse normal‐moveout common‐image‐point gather in which only the hyperbolic moveout is retained. This hyperbolic moveout is the moveout that would be obtained in an isotropic equivalent medium. A hyperbolic velocity is then estimated from this gather by applying hyperbolic moveout analysis. Theoretical analysis shows that for any given initial velocity, the estimated hyperbolic velocity converges by an iterative procedure to the optimal velocity if the velocity ratio is optimal or to a value closer to the optimal velocity if the velocity ratio is not optimal. The velocity ratio (VP/VS) has little effect on the estimation of the velocity. Applying this technique to a synthetic seismic data set confirms the theoretical findings. This work provides a practical method to obtain the velocity model for prestack Kirchhoff time migration. 相似文献
15.
P. HUBRAL 《Geophysical Prospecting》1974,22(4):722-735
The implementation of a stacking filter involves the filtering of each trace with an individual filter and the subsequent summing of all outputs. The actual position of a trace in space as well as certain simultaneous shifts of traces and filter components in time do not influence the process. The resulting output is consequently invariant to various arbitrary coordinate transformations. For a certain useful class of ensembles of non-linear moveout arrival times for signals a particular transformation can be found which transforms a given ensemble into one consisting only of straight lines. It is thus possible to reduce, for instance, the analysis of a stacking filter designed for hyperbola-like moveout curves to the analysis of a velocity filter with linear moveout curves. As the (f—k) transform is a very useful concept to describe a velocity filter, it can consequently be applied to characterize a stacking filter in regard to its performance on input signals with non-linear moveout. 相似文献
16.
M. Malinowski 《Studia Geophysica et Geodaetica》2005,49(4):485-500
The aim of this paper is to show the application of short-period surface waves recorded during deep seismic sounding experiment
for constraining shallow velocity structure of the crust. Phase velocity of fundamental mode Rayleigh waves, observed along
the CELEBRATION 2000 experiment profile CEL09, were obtained by a p-ω method and has been subsequently inverted for one-dimensional
shear velocity models for the top 2 km. Multiple filter technique applied to one shot gather was used to carry out a joint
inversion of phase and group velocity data and to provide γR data to be used for Qβ inversion. Validity of obtained VS and Qβ models was confirmed by the reflectivity method. Noticeably, no clear dispersive wawes were observed in the Tepla-Barrandian
Unit. Quasi-2D model based on the individual 1D VS models is well correlated with the surface geology. Lower VS are observed in the Saxothuringian Zone in comparison to the Moldanubian Zone. In the vicinity of the Central Bohemian and
Moldanubian Plutons, the near-surface VS values are relatively low, but below 1 km depth, they are higher than in surrounding areas. We interpret it as the result
of the weathering and cracks within the granitoid rocks. 相似文献
17.
Common‐midpoint moveout of converted waves is generally asymmetric with respect to zero offset and cannot be described by the traveltime series t2(x2) conventionally used for pure modes. Here, we present concise parametric expressions for both common‐midpoint (CMP) and common‐conversion‐point (CCP) gathers of PS‐waves for arbitrary anisotropic, horizontally layered media above a plane dipping reflector. This analytic representation can be used to model 3D (multi‐azimuth) CMP gathers without time‐consuming two‐point ray tracing and to compute attributes of PS moveout such as the slope of the traveltime surface at zero offset and the coordinates of the moveout minimum. In addition to providing an efficient tool for forward modelling, our formalism helps to carry out joint inversion of P and PS data for transverse isotropy with a vertical symmetry axis (VTI media). If the medium above the reflector is laterally homogeneous, P‐wave reflection moveout cannot constrain the depth scale of the model needed for depth migration. Extending our previous results for a single VTI layer, we show that the interval vertical velocities of the P‐ and S‐waves (VP0 and VS0) and the Thomsen parameters ε and δ can be found from surface data alone by combining P‐wave moveout with the traveltimes of the converted PS(PSV)‐wave. If the data are acquired only on the dip line (i.e. in 2D), stable parameter estimation requires including the moveout of P‐ and PS‐waves from both a horizontal and a dipping interface. At the first stage of the velocity‐analysis procedure, we build an initial anisotropic model by applying a layer‐stripping algorithm to CMP moveout of P‐ and PS‐waves. To overcome the distorting influence of conversion‐point dispersal on CMP gathers, the interval VTI parameters are refined by collecting the PS data into CCP gathers and repeating the inversion. For 3D surveys with a sufficiently wide range of source–receiver azimuths, it is possible to estimate all four relevant parameters (VP0, VS0, ε and δ) using reflections from a single mildly dipping interface. In this case, the P‐wave NMO ellipse determined by 3D (azimuthal) velocity analysis is combined with azimuthally dependent traveltimes of the PS‐wave. On the whole, the joint inversion of P and PS data yields a VTI model suitable for depth migration of P‐waves, as well as processing (e.g. transformation to zero offset) of converted waves. 相似文献
18.
The dispersion of exhaust products of rocket fuel in the direction perpendicular to the motion of a rocket is investigated
in this work. A comparison of the results of numerical calculations with a self-similar approximation of a strong cylindrically
symmetric explosion is fulfilled. It is shown that at sufficiently high rocket velocity V
∞, which exceeds the sum of gas exhaust velocity V
e
from the nozzle and sound speed V
s
(V
∞ > V
e
+V
s
), a gasdynamic hole can arise around the rocket trajectory in the upper atmosphere, inside which the total concentration
of gas becomes less than the equilibrium concentration of gas at a given altitude. The dynamics of the profiles of density
and temperature of the exhaust products inside a rocket plume is calculated. 相似文献
19.
For the problem of matrix compaction and melt segregation a general mush continuity equation is derived, which explicitly
expresses the coupling between the melt percolation and the inelastic matrix deformation and closes the governing equation
set. Besides, a general equation is obtained, which describes the change in the volume of pore space due to all the possible
reasons (inelastic matrix deformation, the phase transitions, and the advection of porosity by the matrix flow). The features
of the isothermal melt segregation inside a partially molten zone are demonstrated using one-dimensional (1D) numerical solutions.
It follows from the solutions that the pattern and the characteristic time of the melt segregation inside a partially molten
zone of thickness L are controlled by the segregation parameter γ
c
= (L/δ
c
)2, where the compaction length δ
c
= k(φ0)η/(φ0μ) depends on the permeability, k, the value of characteristic porosity, φ0, and the viscosities of the matrix, η, and melt, μ. The solutions demonstrate that at any value of γ
c
, layers that are highly enriched in melt compared to the maximum initial porosity are formed in the upper part of the zone.
At the same time, the evolution of the system and the segregation time differ considerably in the limits of γ
c
≪ γ* and γ
c
≫ γ*, where γ* depends on the boundary and initial conditions of the problem, and γ* is about 80 for the problem of melt segregation inside a partially molten zone with the maximum in the initial melt distribution
located in the middle of the zone. At γ
c
≪ γ*, which corresponds to the segregation of low-viscosity ultrabasic melts (kimberlites, carbonatites), all the melt accumulates
to the roof of the zone, and the segregation time does not depend on the matrix permeability and melt viscosity and decreases
with an increase in the thickness of the zone as L
−1. The latter can be the reason for the formation of clusters of the same age and same composition eruptions characteristic
of the kimberlite provinces. In the opposite limiting case, γ
c
≫ γ*, the segregation time does not depend on the matrix viscosity and scales as L with a wave sequence forming in the upper part of the zone, which, probably, elucidates the origin of the rhythmical layering
of the large tholeiitic basalt plutons. 相似文献
20.
转换波四参数速度分析方法在k71地区的应用 总被引：1，自引：0，他引：1
3-D converted-wave data were acquired using digital MEMS （micro-electromechanical system） three component （3C） sensors in the alternating sand and shale sequence in the overburden of the Shengli Ken-71 area. This gives rise to serious non-hyperbolic moveout effects in the converted-wave data due to both the asymmetrical ray path and anisotropic effects. Conventional velocity analysis and moveout correction based on isotropic methods do not flatten reflections events. Here, we use a four-parameter theory to evaluate these effects and process the data. These four parameters include the PS converted wave stacking velocity （Vc2）, the vertical velocity ratio （Y0）, the effective velocity ratio （Yeff）, and the anisotropy parameter （xoff）, The method utilizes the moveout information at different offsets to estimate the different parameters and ensures that the events are properly aligned for stacking, As a result, this four-parameter theory leads to an improvement in imaging quality and correlation between the P-waves and converted-waves. 相似文献