共查询到20条相似文献，搜索用时 31 毫秒
1.
Elastic rock properties can be estimated from prestack seismic data using amplitude variation with offset analysis. P‐wave, S‐wave and density ‘reflectivities’, or contrasts, can be inverted from angle‐band stacks. The ‘reflectivities’ are then inverted to absolute acoustic impedance, shear impedance and density. These rock properties can be used to map reservoir parameters through all stages of field development and production. When P‐wave contrast is small, or gas clouds obscure reservoir zones, multicomponent ocean‐bottom recording of converted‐waves (P to S or Ps) data provides reliable mapping of reservoir boundaries. Angle‐band stacks of multicomponent P‐wave (Pz) and Ps data can also be inverted jointly. In this paper Aki‐Richards equations are used without simplifications to invert angle‐band stacks to ‘reflectivities’. This enables the use of reflection seismic data beyond 30° of incident angles compared to the conventional amplitude variation with offset analysis. It, in turn, provides better shear impedance and density estimates. An important input to amplitude variation with offset analysis is the Vs/Vp ratio. Conventional methods use a constant or a time‐varying Vs/Vp model. Here, a time‐ and space‐varying model is used during the computation of the ‘reflectivities’. The Vs/Vp model is generated using well log data and picked horizons. For multicomponent data applications, the latter model can also be generated from processing Vs/Vp models and available well data. Reservoir rock properties such as λρ, μρ, Poisson's ratio and bulk modulus can be computed from acoustic impedance, shear impedance and density for pore fill and lithology identification. λ and μ are the Lamé constants and ρ is density. These estimations can also be used for a more efficient log property mapping. Vp/Vs ratio or Poisson's ratio, λρ and weighted stacks, such as the one computed from λρ and λ/μ, are good gas/oil and oil/water contact indicators, i.e., pore fill indicators, while μρ mainly indicates lithology. μρ is also affected by pressure changes. Results from a multicomponent data set are used to illustrate mapping of gas, oil and water saturation and lithology in a Tertiary sand/shale setting. Whilst initial log crossplot analysis suggested that pore fill discrimination may be possible, the inversion was not successful in revealing fluid effects. However, rock properties computed from acoustic impedance, shear impedance and density estimates provided good lithology indicators; pore fill identification was less successful. Neural network analysis using computed rock properties provided good indication of sand/shale distribution away from the existing wells and complemented the results depicted from individual rock property inversions. 相似文献
2.
AVO investigations of shallow marine sediments 总被引：2，自引：0，他引：2
Amplitude‐variation‐with‐offset (AVO) analysis is based on the Zoeppritz equations, which enable the computation of reflection and transmission coefficients as a function of offset or angle of incidence. High‐frequency (up to 700 Hz) AVO studies, presented here, have been used to determine the physical properties of sediments in a shallow marine environment (20 m water depth). The properties that can be constrained are P‐ and S‐wave velocities, bulk density and acoustic attenuation. The use of higher frequencies requires special analysis including careful geometry and source and receiver directivity corrections. In the past, marine sediments have been modelled as elastic materials. However, viscoelastic models which include absorption are more realistic. At angles of incidence greater than 40°, AVO functions derived from viscoelastic models differ from those with purely elastic properties in the absence of a critical angle of incidence. The influence of S‐wave velocity on the reflection coefficient is small (especially for low S‐wave velocities encountered at the sea‐floor). Thus, it is difficult to extract the S‐wave parameter from AVO trends. On the other hand, P‐wave velocity and density show a considerably stronger effect. Attenuation (described by the quality factor Q) influences the reflection coefficient but could not be determined uniquely from the AVO functions. In order to measure the reflection coefficient in a seismogram, the amplitudes of the direct wave and the sea‐floor reflection in a common‐midpoint (CMP) gather are determined and corrected for spherical divergence as well as source and streamer directivity. At CMP locations showing the different AVO characteristics of a mud and a boulder clay, the sediment physical properties are determined by using a sequential‐quadratic‐programming (SQP) inversion technique. The inverted sediment physical properties for the mud are: P‐wave velocity α=1450±25 m/s, S‐wave velocity β=90±35 m/s, density ρ=1220±45 kg/m3, quality factor for P‐wave QP=15±200, quality factor for S‐wave QS=10±30. The inverted sediment physical properties for the boulder clay are: α=1620±45 m/s,β=360±200 m/s,ρ=1380±85 kg/m3,QP=790±660,QS=25±10. 相似文献
3.
4.
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. 相似文献
5.
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. 相似文献
6.
土体剪切波速是进行土层地震反应分析的动力学参数，对场地地震动参数确定具有重要意义。基于地质地貌分析，将大同盆地划分为5类典型地质单元。对盆地1429个钻孔剪切波速资料进行分析，探讨VS30与VS20的相关性，研究土体埋深、岩性、地质单元、标贯击数及密实度等地质特征对VS的影响，并基于地质单元、剪切波速比、密实度系数及第四系上部覆盖层厚度相关性分析给出土体VS30预测模型。研究结果表明，基于典型地质特征的VS30预测模型拟合优度R2>0.90，预测精度很高，对于离散性较大、直接拟合估算较差及无剪切波速场地来说，以区分地质单元及土体类型的方式进行VS30分解预测是良好的研究思路。首次在区分地质单元及土体类型的前提下提出剪切波速比及密实度系数，并将其与第四系上部覆盖层厚度综合应用于VS30预测研究。研究结果可为大同盆地城市防震减灾规划、震害预测、区域性地震安全评价提供重要技术支撑。 相似文献
7.
The relationship between Vp and Vs may be used to predict Vs where only Vp is known. Vp/Vs is also used to identify pore fluids from seismic data and amplitude variation with offset analysis. Theoretical, physical, as well as statistical empirical Vp‐Vs relationships have been proposed for reservoir characterization when shear‐wave data are not available. In published work, the focus is primarily on the Vp‐Vs relationship of quartzitic sandstone. In order to broaden the picture we present Vp‐Vs relationships of greensand composed of quartz and glauconite by using data from the Paleocene greensand Nini oil field in the North Sea. A Vp‐Vs relationship derived from modelling is compared with empirical Vp‐Vs regressions from laboratory data as well as from log data. The accuracy of Vs prediction is quantified in terms of root‐mean‐square error. We find that the Vp‐Vs relationship derived from modelling works well for greensand shear‐wave velocity prediction. We model the seismic response of glauconitic greensand by using laboratory data from the Nini field. Our studies here reveal that brine‐saturated glauconitic greensand can have a similar seismic response to that from oil‐saturated quartzitic sandstone and that oil‐saturated strongly cemented greensand can have a similar amplitude variation with offset response to that from brine‐saturated weakly cemented greensand. 相似文献
8.
The use of relaxation mechanisms has recently made it possible to simulate viscoelastic (Q) effects accurately in time-domain numerical computations of seismic responses. As a result, seismograms may now be synthesized for models with arbitrary spatial variations in compressional- and shear-wave quality factors (Q9, and Qs, as well as in density (ρ) and compressional- and shear-wave velocities (Vp, and Vs). Reflections produced by Q contrasts alone may have amplitudes as large as those produced by velocity contrasts. Q effects, including their interaction with Vp, Vs and p, contribute significantly to the seismic response of reservoirs. For band-limited data at typical seismic frequencies, the effects of Q on reflectivity and attenuation are more visible than those on dispersion. Synthetic examples include practical applications to reservoir exploration, evaluation and monitoring. Q effects are clearly visible in both surface and offset vertical seismic profile data. Thus, AVO analyses that neglect Q may produce erroneous conclusions. 相似文献
9.
A practical method is presented for determining three‐dimensional S‐wave velocity (VS) profile from microtremor measurements. Frequency–wave number (f–k) spectral analyses of microtremor array records are combined, for this purpose, with microtremor horizontal‐to‐vertical (H/V) spectral ratio techniques. To demonstrate the effectiveness of the proposed method, microtremor measurements using arrays of sensors were conducted at six sites in the city of Kushiro, Japan. The spectral analyses of the array records yield dispersion characteristics of Rayleigh waves and H/V spectra of surface waves, and joint inversion of these data results in VS profiles down to bedrock at the sites. Conventional microtremor measurements were performed at 230 stations within Kushiro city, resulting in the H/V spectra within the city. Three‐dimensional VS structure is then estimated from inversion of the H/V spectra with the VS values determined from the microtremor array data. This reveals three‐dimensional VS profile of Kushiro city, together with an unknown hidden valley that crosses the central part of the city. The estimated VS profile is consistent with available velocity logs and results of subsequent borings, indicating the effectiveness of the proposed method. Copyright © 2008 John Wiley & Sons, Ltd. 相似文献
10.
The main objective of this work is to establish the applicability of shallow surface‐seismic traveltime tomography in basalt‐covered areas. A densely sampled ～1300‐m long surface seismic profile, acquired as part of the SeiFaBa project in 2003 ( Japsen et al. 2006 ) at Glyvursnes in the Faroe Islands, served as the basis to evaluate the performance of the tomographic method in basalt‐covered areas. The profile is centred at a ～700‐m deep well. VP, VS and density logs, a zero‐offset VSP, downhole‐geophone recordings and geological mapping in the area provided good means of control. The inversion was performed with facilities of the Wide Angle Reflection/Refraction Profiling program package ( Ditmar et al. 1999 ). We tested many inversion sequences while varying the inversion parameters. Modelled traveltimes were verified by full‐waveform modelling. Typically an inversion sequence consists in several iterations that proceed until a satisfactory solution is reached. However, in the present case with high velocity contrasts in the subsurface we obtained the best result with two iterations: first obtaining a smooth starting model with small traveltime residuals by inverting with a high smoothing constraint and then inverting with the lowest possible smoothing constraint to allow the inversion to have the full benefit of the traveltime residuals. The tomogram gives usable velocity information for the near‐surface geology in the area but fails to reproduce the expected velocity distribution of the layered basalt flows. Based on the analysis of the tomogram and geological mapping in the area, a model was defined that correctly models first arrivals from both surface seismic data and downhole‐geophone data. 相似文献
11.
Umberto Spagnolini 《Geophysical Prospecting》1994,42(3):241-259
Minimization of seismic residuals does not guarantee uniqueness of the model, and this implies ambiguities in the inversion. Amplitude vs. offset (AVO) inversion does not lead to a unique solution of single elastic interface parameters unless converted and S-wave or critical angle reflections are available. Given the ambiguity of AVO inversion, this paper discusses the interaction between AVO and velocity estimation. The number of independent parameters necessary to describe an isolated reflection with AVO behaviour and residual velocity error is determined. Statistical analysis allows the establishment of an approximate equivalence of the effects of AVO and slight velocity variations; this equivalence cannot be solved without geological a priori information (kinematic equivalence). The data are then decomposed into compound events (i.e. sequences of N interfaces that follow each other at a fixed time lag). The decomposition is obtained by extrapolating the results of the analysis from narrowband to wideband data. Compound events decomposition demonstrates that AVO inversion is ambiguous, not only in the physical parameter space (P- and S-wave velocities, and density) but also kinematically. As an example of compound event decomposition, a medium is derived. This medium is geologically implausible but is kinematically equivalent. 相似文献
12.
Sensitivity of time-lapse seismic to reservoir stress path 总被引：1，自引：1，他引：1
Colin M. Sayers 《Geophysical Prospecting》2006,54(3):369-380
The change in reservoir pore pressure due to the production of hydrocarbons leads to anisotropic changes in the stress field acting on the reservoir. Reservoir stress path is defined as the ratio of the change in effective horizontal stress to the change in effective vertical stress from the initial reservoir conditions, and strongly influences the depletion‐induced compaction behaviour of the reservoir. Seismic velocities in sandstones vary with stress due to the presence of stress‐sensitive regions within the rock, such as grain boundaries, microcracks, fractures, etc. Since the response of any microcracks and grain boundaries to a change in stress depends on their orientation relative to the principal stress axes, elastic‐wave velocities are sensitive to reservoir stress path. The vertical P‐ and S‐wave velocities, the small‐offset P‐ and SV‐wave normal‐moveout (NMO) velocities, and the P‐wave amplitude‐versus‐offset (AVO) are sensitive to different combinations of vertical and horizontal stress. The relationships between these quantities and the change in stress can be calibrated using a repeat seismic, sonic log, checkshot or vertical seismic profile (VSP) at the location of a well at which the change in reservoir pressure has been measured. Alternatively, the variation of velocity with azimuth and distance from the borehole, obtained by dipole radial profiling, can be used. Having calibrated these relationships, the theory allows the reservoir stress path to be monitored using time‐lapse seismic by combining changes in the vertical P‐wave impedance, changes in the P‐wave NMO and AVO behaviour, and changes in the S‐wave impedance. 相似文献
13.
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. 相似文献
14.
Elastic imaging from ocean bottom cable (OBC) data can be challenging because it requires the prior estimation of both compressional‐wave (P‐wave) and shear‐wave (S‐wave) velocity fields. Seismic interferometry is an attractive technique for processing OBC data because it performs model‐independent redatuming; retrieving ‘pseudo‐sources’ at positions of the receivers. The purpose of this study is to investigate multicomponent applications of interferometry for processing OBC data. This translates into using interferometry to retrieve pseudo‐source data on the sea‐bed not only for multiple suppression but for obtaining P‐, converted P to S‐wave (PS‐wave) and possibly pure mode S‐waves. We discuss scattering‐based, elastic interferometry with synthetic and field OBC datasets. Conventional and scattering‐based interferometry integrands computed from a synthetic are compared to show that the latter yields little anti‐causal response. A four‐component (4C) pseudo‐source response retrieves pure‐mode S‐reflections as well at P‐ and PS‐reflections. Pseudo‐source responses observed in OBC data are related to P‐wave conversions at the seabed rather than to true horizontal or vertical point forces. From a Gulf of Mexico OBC data set, diagonal components from a nine‐component pseudo‐source response demonstrate that the P‐wave to S‐wave velocity ratio (VP/VS) at the sea‐bed is an important factor in the conversion of P to S for obtaining the pure‐mode S‐wave reflections. 相似文献
15.
The azimuthally varying non‐hyperbolic moveout of P‐waves in orthorhombic media can provide valuable information for characterization of fractured reservoirs and seismic processing. Here, we present a technique to invert long‐spread, wide‐azimuth P‐wave data for the orientation of the vertical symmetry planes and five key moveout parameters: the symmetry‐plane NMO velocities, V(1)nmo and V(2)nmo , and the anellipticity parameters, η(1), η(2) and η(3) . The inversion algorithm is based on a coherence operator that computes the semblance for the full range of offsets and azimuths using a generalized version of the Alkhalifah–Tsvankin non‐hyperbolic moveout equation. The moveout equation provides a close approximation to the reflection traveltimes in layered anisotropic media with a uniform orientation of the vertical symmetry planes. Numerical tests on noise‐contaminated data for a single orthorhombic layer show that the best‐constrained parameters are the azimuth ? of one of the symmetry planes and the velocities V(1)nmo and V(2)nmo , while the resolution in η(1) and η(2) is somewhat compromised by the trade‐off between the quadratic and quartic moveout terms. The largest uncertainty is observed in the parameter η(3) , which influences only long‐spread moveout in off‐symmetry directions. For stratified orthorhombic models with depth‐dependent symmetry‐plane azimuths, the moveout equation has to be modified by allowing the orientation of the effective NMO ellipse to differ from the principal azimuthal direction of the effective quartic moveout term. The algorithm was successfully tested on wide‐azimuth P‐wave reflections recorded at the Weyburn Field in Canada. Taking azimuthal anisotropy into account increased the semblance values for most long‐offset reflection events in the overburden, which indicates that fracturing is not limited to the reservoir level. The inverted symmetry‐plane directions are close to the azimuths of the off‐trend fracture sets determined from borehole data and shear‐wave splitting analysis. The effective moveout parameters estimated by our algorithm provide input for P‐wave time imaging and geometrical‐spreading correction in layered orthorhombic media. 相似文献
16.
叠前反演技术已经成为岩性预测和流体识别的关键技术,并将在油气勘探中发挥重要作用,本文提出了模拟退火AVA同步反演方法,该方法以多个角道集地震数据体和先验纵波速度、横波速度、密度模型等为输入,综合利用地震、地质、测井等多种信息为约束,采用模拟退火全局优化方法,直接采用Knott-Zoeppritz方程求解反射系数,同步反演获得纵波阻抗、横波阻抗、纵横波速度、密度、Mu、Lambda等23种岩石弹性参数模型.其主要特点如下：１）采用Zoeppritz方程直接求解,精度高;2) 目标函数选择灵活;３）多角度提取子波;４）可获得深度域的反演结果.实际资料应用证明,该反演方法保留了地震反射振幅随偏移距不同或入射角不同而变化的特征,可有效提高储层预测和流体识别精度.实际应用表明,本方法预测结果可靠,具有较强的实用性. 相似文献
17.
By applying seismic inversion, we can derive rock impedance from seismic data. Since it is an interval property, impedance is valuable for reservoir characterization. Furthermore, the decomposition of the impedance into two fundamental properties, i.e. velocity and density, provides a link to the currently available rock‐physics applications to derive quantitative reservoir properties. However, the decomposition is a challenging task due to the strong influence of noise, especially for seismic data with a maximum offset angle of less than 30°. We present a method of impedance decomposition using three elastic impedance data derived from the seismic inversion of angle stacks, where the far‐stack angle is 23.5°. We discuss the effect of noise on the analysis as being the most significant cause of making the decomposition difficult. As the result, the offset‐consistent component of noise mostly affects the determination of density but not the velocities (P‐ and S‐wave), whereas the effect of the random component of noise occurs equally in the determination of the velocities and density. The effect is controlled by the noise enhancement factor 1/A, which is determined by a combination of stack angles. Based on the results of the analysis, we show an innovative method of decomposition incorporating rock‐physics bounds as constraints for the analysis. The method is applied to an actual data set from an offshore oilfield; we demonstrate the result of analysis for sandbody detection. 相似文献
18.
Haruhisa Nakamichi Satoru Tanaka Hiroyuki Hamaguchi 《Journal of Volcanology and Geothermal Research》2002,116(3-4)
A genetic algorithm inversion of receiver functions derived from a dense seismic network around Iwate volcano, northeastern Japan, provides the fine S wave velocity structure of the crust and uppermost mantle. Since receiver functions are insensitive to an absolute velocity, travel times of P and S waves propagating vertically from earthquakes in the subducting slab beneath the volcano are involved in the inversion. The distribution of velocity perturbations in relation to the hypocenters of the low-frequency (LF) earthquakes helps our understanding of deep magmatism beneath Iwate volcano. A high-velocity region (dVS/VS=10%) exists around the volcano at depths of 2–15 km, with the bottom depth decreasing to 11 km beneath the volcano’s summit. Just beneath the thinning high-velocity region, a low-velocity region (dVS/VS=−10%) exists at depths of 11–20 km. Intermediate-depth LF (ILF) events are distributed vertically in the high-velocity region down to the top of the low-velocity region. This distribution suggests that a magma reservoir situated in the low-velocity region supplies magma to a narrow conduit that is detectable by the hypocenters of LF earthquakes. Another broad low-velocity region (dVS/VS=−5 to −10%) occurs at depths of 17–35 km. Additional clusters of deep LF (DLF) events exist at depths of 32–37 km in the broad low-velocity zone. The DLF and ILF events are the manifestations of magma movement near the Moho discontinuity and in the conduit just beneath the volcano, respectively. 相似文献
19.
We describe a method to invert a walkaway vertical seismic profile (VSP) and predict elastic properties (P‐wave velocity, S‐wave velocity and density) in a layered model looking ahead of the deepest receiver. Starting from Bayes's rule, we define a posterior distribution of layered models that combines prior information (on the overall variability of and correlations among the elastic properties observed in well logs) with information provided by the VSP data. This posterior distribution of layered models is sampled by a Monte‐Carlo method. The sampled layered models agree with prior information and fit the VSP data, and their overall variability defines the uncertainty in the predicted elastic properties. We apply this technique first to a zero‐offset VSP data set, and show that uncertainty in the long‐wavelength P‐wave velocity structure results in a sizable uncertainty in the predicted elastic properties. We then use walkaway VSP data, which contain information on the long‐wavelength P‐wave velocity (in the reflection moveout) and on S‐wave velocity and density contrasts (in the change of reflectivity with offset). The uncertainty of the look‐ahead prediction is considerably decreased compared with the zero‐offset VSP, and the predicted elastic properties are in good agreement with well‐log measurements. 相似文献
20.
Velocities of compressional and shear waves in limestones 总被引：2，自引：1，他引：2
Carbonate rocks are important hydrocarbon reservoir rocks with complex textures and petrophysical properties (porosity and permeability) mainly resulting from various diagenetic processes (compaction, dissolution, precipitation, cementation, etc.). These complexities make prediction of reservoir characteristics (e.g. porosity and permeability) from their seismic properties very difficult. To explore the relationship between the seismic, petrophysical and geological properties, ultrasonic compressional‐ and shear‐wave velocity measurements were made under a simulated in situ condition of pressure (50 MPa hydrostatic effective pressure) at frequencies of approximately 0.85 MHz and 0.7 MHz, respectively, using a pulse‐echo method. The measurements were made both in vacuum‐dry and fully saturated conditions in oolitic limestones of the Great Oolite Formation of southern England. Some of the rocks were fully saturated with oil. The acoustic measurements were supplemented by porosity and permeability measurements, petrological and pore geometry studies of resin‐impregnated polished thin sections, X‐ray diffraction analyses and scanning electron microscope studies to investigate submicroscopic textures and micropores. It is shown that the compressional‐ and shear‐wave velocities (Vp and Vs, respectively) decrease with increasing porosity and that Vp decreases approximately twice as fast as Vs. The systematic differences in pore structures (e.g. the aspect ratio) of the limestones produce large residuals in the velocity versus porosity relationship. It is demonstrated that the velocity versus porosity relationship can be improved by removing the pore‐structure‐dependent variations from the residuals. The introduction of water into the pore space decreases the shear moduli of the rocks by about 2 GPa, suggesting that there exists a fluid/matrix interaction at grain contacts, which reduces the rigidity. The predicted Biot–Gassmann velocity values are greater than the measured velocity values due to the rock–fluid interaction. This is not accounted for in the Biot–Gassmann velocity models and velocity dispersion due to a local flow mechanism. The velocities predicted by the Raymer and time‐average relationships overestimated the measured velocities even more than the Biot model. 相似文献