共查询到20条相似文献，搜索用时 15 毫秒
1.
Most amplitude versus offset (AVO) analysis and inversion techniques are based on the Zoeppritz equations for plane‐wave reflection coefficients or their approximations. Real seismic surveys use localized sources that produce spherical waves, rather than plane waves. In the far‐field, the AVO response for a spherical wave reflected from a plane interface can be well approximated by a plane‐wave response. However this approximation breaks down in the vicinity of the critical angle. Conventional AVO analysis ignores this problem and always utilizes the plane‐wave response. This approach is sufficiently accurate as long as the angles of incidence are much smaller than the critical angle. Such moderate angles are more than sufficient for the standard estimation of the AVO intercept and gradient. However, when independent estimation of the formation density is required, it may be important to use large incidence angles close to the critical angle, where spherical wave effects become important. For the amplitude of a spherical wave reflected from a plane fluid‐fluid interface, an analytical approximation is known, which provides a correction to the plane‐wave reflection coefficients for all angles. For the amplitude of a spherical wave reflected from a solid/solid interface, we propose a formula that combines this analytical approximation with the linearized plane‐wave AVO equation. The proposed approximation shows reasonable agreement with numerical simulations for a range of frequencies. Using this solution, we constructed a two‐layer three‐parameter least‐squares inversion algorithm. Application of this algorithm to synthetic data for a single plane interface shows an improvement compared to the use of plane‐wave reflection coefficients. 相似文献
2.
The effect of surface phenomena occurring at the interfaces between immiscible fluids and a solid on the seismic attributes of partially saturated rocks has not yet been fully studied. Meanwhile, over the past two decades considerable progress has been made in the physics of wetting to understand effects such as contact line friction, contact line pinning, contact angle hysteresis, and equilibrium contact angle. In this paper, we developed a new rock physics model considering the aforementioned effects on seismic properties of the rock with a partially saturated plane-strain crack. We demonstrated that for small wave-induced stress perturbations, the contact line of the interface meniscus will remain pinned, while the meniscus will bulge and change its shape through the change of the contact angles. When the stress perturbation is larger than a critical value, the contact line will move with advancing or receding contact angle depending on the direction of contact line motion. A critical stress perturbation predicted by our model can be in the range of ∼102−104 Pa, that is typical for linear seismic waves. Our model predicts strong seismic attenuation in the case when the contact line is moving. When the contact line is pinned, the attenuation is negligibly small. Seismic attenuation is associated with the hysteresis of loading and unloading bulk moduli, predicted by our model. The hysteresis is large when the contact line is moving and negligibly small when the contact line is pinned. Furthermore, we demonstrate that the bulk modulus of the rock with a partially saturated crack depends also on the surface tension and on the contact angle hysteresis. These parameters are typically neglected during calculation of the effecting fluid moduli by applying different averaging techniques. We demonstrate that contact line friction may be a dominant seismic attenuation mechanism in the low frequency limit (<∼10 Hz) when capillary forces dominate over viscous forces during wave-induced two-phase fluid flow. 相似文献
3.
Contact angle is a principal control of the flow of multiple fluid phases through porous media; however its measurement on other than flat surfaces remains a challenge. A new method is presented for the measurement of the contact angle between immiscible fluids at the pore scale at reservoir conditions (10 MPa and 50 °C) inside a quarry limestone through the use of X-ray microtomography. It is applied to a super-critical CO2–brine–carbonate system by resampling the micro-CT data onto planes orthogonal to the contact lines, allowing for vectors to be traced along the grain surface and the CO2–brine interface. A distribution of contact angles ranging from 35° to 55° is observed, indicating that the CO2–brine–carbonate system is weakly water-wet. This range of contact angles can be understood as the result of contact angle hysteresis and surface heterogeneity on a range of length scales. Surface heterogeneity is examined by comparison of micro-CT results with optical thin sections and SEM images. 相似文献
4.
In this paper, we derive an approximation of the SS-wave reflection coefficient and the expression of S-wave ray elastic impedance (SREI) in terms of the ray parameter. The SREI can be expressed by the S-wave incidence angle or P-wave reflection angle, referred to as SREIS and SREIP, respectively. Our study using elastic models derived from real log measurements shows that SREIP has better capability for lithology and fluid discrimination than SREIS and conventional S-wave elastic impedance (SEI). We evaluate the SREIP feasibility using 25 groups of samples from Castagna and Smith (1994). Each sample group is constructed by using shale, brine-sand, and gas-sand. Theoretical evaluation also indicates that SREIP at large incident angles is more sensitive to fluid than conventional fluid indicators. Real seismic data application also shows that SREIP at large angles calculated using P-wave and S-wave impedance can efficiently characterize tight gas-sand. 相似文献
5.
The analytical study of the reflection and transmission coefficient of the seafloor interface is essential for the characterization of the ocean bottom in marine seismic exploration. Based on the boundary conditions of the seafloor interface, the analytical expression of the reflection and transmission coefficient at the submarine interface is derived in this study by using the steady-state wave solution of the elastic wave in a homogeneous, isotropic medium. With this analytical expression, the characteristics of the reflection and transmission coefficient at the submarine interface are analysed and discussed using critical angles. The results show that the change in the reflection and transmission coefficient with the incidence angle presents a “segmented” characteristic, in which the critical angle is the dividing point. The amplitude value and phase angle of the coefficient at the submarine interface change dramatically at the critical angle, which is related to the P- and S-wave velocities in the seabed layer. Compared with the stiff seabed, the soft seabed has a larger P-wave critical angle and an absence of the converted S-wave critical angle, owing to the low P- and S-wave velocities in the solid seabed layer. By analysing and discussing the special changes that occur in the coefficient values at the critical angle, the reflection and transmission characteristics of the different incident angles are obtained. Synthetic models of both stiff and soft seafloors are provided in this study to verify the analytical results. Finally, we compared our synthetic results with real data from the Gulf of Mexico, which enabled the validation of our conclusions. 相似文献
6.
A quasi-static scheme based on pore space spatial statistics is presented to simulate pore-scale two-phase capillary-dominant displacement processes. The algorithm is coupled with computational fluid dynamics in order to evaluate saturation functions. Wettability heterogeneity in partial and fractional/mixed-wet media is implemented using a contact angle map. The simulation process is pixel-wised and performed directly on binary images. Bypassing and snap-off are tackled as non-wetting phase trapping mechanisms. Post-processing results include residual saturations, effective permeability and capillary pressure curves for drainage and imbibition scenarios. The primary advantages of the proposed workflow are eliminating pore space skeletisation/ discretization, superior time efficiency and minimal numerical drawbacks when compared to other direct or network-based simulation techniques. 相似文献
7.
Manuel Carbajal-Romero Esteban Flores-Mendez Norberto Flores-Guzmán Jaime Núñez-Farfán Enrique Olivera-Villaseñor Francisco José Sánchez-Sesma 《Geofísica Internacional》2013,52(1):21-30
The present work shows the propagation of Scholte interface waves at the boundary of a fluid in contact with an elastic solid, for a broad range of solid materials. It has been demonstrated that by an analysis of diffracted waves in a fluid it is possible to infer the mechanical properties of the elastic solid medium, specifically, its propagation velocities. For this purpose, the diffracted wave field of pressures and displacements, due to an initial wave of pressure in the fluid, are expressed using boundary integral representations, which satisfy the equation of motion. The source in the fluid is represented by a Hankel’s function of second kind and zero order. The solution to this wave propagation problem is obtained by means of the Indirect Boundary Element Method, which is equivalent to the well-known Somigliana representation theorem. The validation of the results is carried out by using the Discrete Wave Number Method and the Spectral Element Method. Firstly, we show spectra of pressures that illustrate the behavior of the fluid for each solid material considered, then, we apply the Fast Fourier Transform to show results in time domain. Snapshots to exemplify the emergence of Scholte’s waves are also included. 相似文献
8.
A novel semi-analytical model for computation of capillary entry pressures and associated fluid configurations in arbitrary, potentially non-convex, 2D pore space geometries at uniform wettability is developed. The model computes all possible centre positions of circular arcs, and physically sound criteria are implemented to determine the set of these arcs that correspond to geometrically allowed interfaces. Interfaces and pore boundary segments are connected to form closed boundaries of identified geometrical regions. These regions are classified as either oil regions, located in the wider parts of the pore space, or as water regions located in pore space constrictions. All possible region combinations are identified and evaluated for each radius value in an iterative procedure to determine the favourable entry radius and corresponding configuration based on minimisation of free energy. The model has been validated by comparison with known analytical solutions in idealised pore geometries. In cases where different analytical solutions are geometrically possible, the model generates several oil and water regions, and the valid solution is determined by the region combination that corresponds to the most favourable entry pressure, consistent with the analytical solution. Entry pressure radii and configurations are computed in strongly non-convex pore spaces extracted from an image of Bentheimer sandstone, which demonstrates that the model captures successfully well-known characteristics of capillary behaviour at different wetting conditions. The computations also demonstrate the importance of selecting the fluid configuration of minimum change in free energy. In some cases, a merged region formed by a combination of oil and water regions corresponds to the favourable entry configuration of oil, whereas in other cases, an individual oil region may correspond to the favourable oil entry configuration. It is also demonstrated that oil entry configurations may constitute merged regions for weakly water-wet conditions and individual oil regions for strongly water-wet conditions in the same pore space. The computations show that the ratio of pore area to perimeter multiplied by the cosine of the contact angle under-predicts the entry pressure radii in Bentheimer sandstone pore spaces. An alternative formula is proposed for prediction of entry radii for nonzero contact angles based on the entry radii obtained for zero contact angles. 相似文献
9.
Seismic velocity analysis in the scattering-angle/azimuth domain 总被引：2，自引：0，他引：2
Migration velocity analysis is carried out by analysing the residual moveout and amplitude variations in common image point gathers (CIGs) parametrized by scattering angle and azimuth. The misfit criterion in the analysis is of the differential-semblance type. By using angles to parametrize the imaging we are able to handle and exploit data with multiple arrivals, although artefacts may occur in the CIGs and need to be suppressed. The CIGs are generated by angle migration, an approach based on the generalized Radon transform (GRT) inversion, and they provide multiple images of reflectors in the subsurface for a range of scattering angles and azimuths. Within the differential semblance applied to these CIGs, we compensate for amplitude versus angle (AVA) effects. Thus, using a correct background velocity model, the CIGs should have no residual moveout nor amplitude variation with angles, and the differential semblance should vanish. If the velocity model is incorrect, however, the events in the CIGs will appear at different depths for different angles and the amplitude along the events will be non-uniform. A standard, gradient-based optimization scheme is employed to develop a velocity updating procedure. The model update is formed by backprojecting the differential semblance misfits through ray perturbation kernels, within a GRT inverse. The GRT inverse acts on the data, subject to a shift in accordance with ray perturbation theory. The performance of our algorithm is demonstrated with two synthetic data examples using isotropic elastic models. The first one allows velocity variation with depth only. In the second one, we reconstruct a low-velocity lens in the model that gives rise to multipathing. The velocity model parametrization is based upon the eigentensor decomposition of the stiffness tensor and makes use of B-splines. 相似文献
10.
沿裂隙发育的天然气水合物是印度深水盆地细粒沉积物中水合物的重要产出方式,水合物以结核状或脉状充填在高角度裂隙中.天然气水合物主要沿着构造主应力方向生成,由于裂隙的存在,含水合物的沉积物层呈现各向异性.利用孔隙介质中水合物呈均匀分布的速度模型计算的NGHP01-10D井水合物饱和度高达40%,而压力取芯表明水合物饱和度占孔隙空间的20%左右.为了研究水合物饱和度差异,基于层状介质的各向异性模型计算了裂隙充填型水合物的饱和度.在垂直井孔中,由于波入射角与裂隙倾角有关,考虑裂隙倾角变化,利用纵波和横波速度同时反演水合物饱和度和裂隙倾角.利用层状介质模型计算的水合物占孔隙空间的15%～25%,裂隙的倾角在60°～90°,多为高角度裂隙.在NGHP01-10D井中,纵横波速度联合计算的饱和度与压力取芯结果吻合更好. 相似文献
11.
Full waveform inversion algorithms are widely used in the construction of subsurface velocity models. In the following study, we propose a Laplace–Fourier-domain waveform inversion algorithm that uses both Laplace-domain and Fourier-domain wavefields to achieve the reconstruction of subsurface velocity models. Although research on the Laplace–Fourier-domain waveform inversion has been published recently that study is limited to fluid media. Because the geophysical targets of marine seismic exploration are usually located within solid media, waveform inversion that is approximated to acoustic media is limited to the treatment of properly identified submarine geophysical features. In this study, we propose a full waveform inversion algorithm for isotropic fluid–solid media with irregular submarine topography comparable to a real marine environment. From the fluid–solid system, we obtained P and S wave velocity models from the pressure data alone. We also suggested strategies for choosing complex frequency bands constructed of frequencies and Laplace coefficients to improve the resolution of the restored velocity structures. For verification, we applied our Laplace–Fourier-domain waveform inversion for fluid–solid media to synthetic data that were reconstructed for fluid–solid media. Through this inversion test, we successfully restored reasonable velocity structures. Furthermore, we successfully extended our algorithm to a field data set. 相似文献
12.
13.
Ute Mann Daniel J. Frost David C. Rubie 《Physics of the Earth and Planetary Interiors》2008,167(1-2):1-7
The wetting characteristics of liquid Fe–Si alloys in a matrix of the respective predominating stable silicate mantle mineral (forsterite or silicate perovskite) at pressures of 2–5 and 25 GPa and temperatures of 1600–2000 °C were studied by determining the liquid metal–solid silicate contact angles. The median angle values from texturally equilibrated samples were found to be independent of pressure, temperature, silicate mineralogy and the Si content in the metal fraction and range between 130° and 140° which is far above the critical wetting boundary of 60°. This shows that within the studied range of conditions dissolved Si does not lower the surface energies between Fe-rich liquids and silicate mantle grains. As a consequence, under reducing conditions the presence of Si in the metal phase of planetary bodies would not have enhanced percolative flow as an effective metal–silicate separation process. 相似文献
14.
We present a semi-analytical, combinatorial approach to compute three-phase capillary entry pressures for gas invasion into pore throats with constant cross-sections of arbitrary shapes that are occupied by oil and/or water. For a specific set of three-phase capillary pressures, geometrically allowed gas/oil, oil/water and gas/water arc menisci are determined by moving two circles in opposite directions along the pore/solid boundary for each fluid pair such that the contact angle is defined at the front circular arcs. Intersections of the two circles determine the geometrically allowed arc menisci for each fluid pair. The resulting interfaces are combined systematically to allow for all geometrically possible three-phase configuration changes. The three-phase extension of the Mayer and Stowe – Princen method is adopted to calculate capillary entry pressures for all determined configuration candidates, from which the most favorable gas invasion configuration is determined. The model is validated by comparing computed three-phase capillary entry pressures and corresponding fluid configurations with analytical solutions in idealized triangular star-shaped pores. It is demonstrated that the model accounts for all scenarios that have been analyzed previously in these shapes. Finally, three-phase capillary entry pressures and associated fluid configurations are computed in throat cross-sections extracted from segmented SEM images of Bentheim sandstone. The computed gas/oil capillary entry pressures account for the expected dependence of oil/water capillary pressure in spreading and non-spreading fluid systems at the considered wetting conditions. Because these geometries are irregular and include constrictions, we introduce three-phase displacements that have not been identified previously in pore-network models that are based on idealized pore shapes. However, in the limited number of pore geometries considered in this work, we find that the favorable displacements are not generically different from those already encountered in network models previously, except that the size and shape of oil layers that are surrounded by gas and water are described more realistically. The significance of the results for describing oil connectivity in porous media accurately can only be evaluated by including throats with more complex cross-sections in three-phase pore-network models. 相似文献
15.
本文针对中国石油集团测井有限公司(CPL)三维阵列感应测井仪器(Three Dimensional Induction Tool,TDIT)的真实结构与参数,利用2.5维数值模式匹配技术研究建立水平层状非均质TI地层中仪器偏心情况下三维感应测井响应的高效数值模拟算法.首先,利用2.5维数值模式匹配技术给出水平层状非均质TI模型中偏心磁流源并矢Green函数的半解析解.在此基础上,结合三维阵列感应测井仪器参数以及测量过程中出现的仪器偏心与旋转等实际情况,通过叠加原理推导出该仪器响应的有效计算公式.最后,通过数值模拟结果分析考察仪器偏心、仪器自旋角等参数对测井响应的影响以及变化特征,为多分量感应资料处理与解释提供理论基础. 相似文献
16.
O. HOLBERG 《Geophysical Prospecting》1988,36(2):99-114
Numerical wavefield extrapolation represents the backbone of any algorithm for depth migration pre- or post-stack. For such depth imaging techniques to yield reliable and interpretable results, the underlying wavefield extrapolation algorithm must propagate the waves through inhomogeneous media with a minimum of numerically induced distortion, over a range of frequencies and angles of propagation. A review of finite-difference (FD) approximations to the acoustic one-way wave equation in the space-frequency domain is presented. A straightforward generalization of the conventional FD formulation leads to an algorithm where the wavefield is continued downwards with space-variant symmetric convolutional operators. The operators can be precomputed and made accessible in tables such that the ratio between the temporal frequency and the local velocity is used to determine the correct operator at each grid point during the downward continuation. Convolutional operators are designed to fit the desired dispersion relation over a range of frequencies and angles of propagation such that the resulting numerical distortion is minimized. The optimization is constrained to ensure that evanescent energy and waves propagating at angles higher than the maximum design angle are attenuated in each extrapolation step. The resulting operators may be viewed as optimally truncated and bandlimited spatial versions of the familiar phase shift operator. They are unconditionally stable and can be applied explicitly. This results in a simple wave propagation algorithm, eminently suited for implementation on pipelined computers and on large parallel computing systems. 相似文献
17.
Joint facies and rock properties Bayesian amplitude‐versus‐offset inversion using Markov random fields 
下载免费PDF全文
下载免费PDF全文 Seismic reflection pre‐stack angle gathers can be simultaneously inverted within a joint facies and elastic inversion framework using a hierarchical Bayesian model of elastic properties and categorical classes of rock and fluid properties. The Bayesian prior implicitly supplies low frequency information via a set of multivariate compaction trends for each rock and fluid type, combined with a Markov random field model of lithotypes, which carries abundance and continuity preferences. For the likelihood, we use a simultaneous, multi‐angle, convolutional model, which quantifies the data misfit probability using wavelets and noise levels inferred from well ties. Under Gaussian likelihood and facies‐conditional prior models, the posterior has simple analytic form, and the maximum a‐posteriori inversion problem boils down to a joint categorical/continuous non‐convex optimisation problem. To solve this, a set of alternative, increasingly comprehensive optimisation strategies is described: (i) an expectation–maximisation algorithm using belief propagation, (ii) a globalisation of method (i) using homotopy, and (iii) a discrete space approach using simulated annealing. We find that good‐quality inversion results depend on both sensible, elastically separable facies definitions, modest resolution ambitions, reasonably firm abundance and continuity parameters in the Markov random field, and suitable choice of algorithm. We suggest usually two to three, perhaps four, unknown facies per sample, and usage of the more expensive methods (homotopy or annealing) when the rock types are not strongly distinguished in acoustic impedance. Demonstrations of the technique on pre‐stack depth‐migrated field data from the Exmouth basin show promising agreements with lithological well data, including prediction accuracy improvements of 24% in and twofold in density, in comparison to a standard simultaneous inversion. Much clearer and extensive recovery of the thin Pyxis gas field was evident using stronger coupling in the Markov random field model and use of the homotopy or annealing algorithms. 相似文献
18.
Experimental characterization of the impact of temperature and humidity on the breakdown of soil water repellency in sandy soils and composts 下载免费PDF全文
下载免费PDF全文 Soil water repellency is a widespread phenomenon with the capacity to alter hydrological and geomorphological processes. Water repellency decays with time, and the consequences are only of concern during the timescale at which the water repellency persists. This study aimed to characterize the influence of temperature and humidity on the breakdown of water repellency. Apparent contact angle measurements were carried out on samples consisting of sand treated with stearic acid as well as naturally repellent dune sands and composts. Temperature and humidity were controlled using a cooled incubator and a purpose designed enclosed box in which humidity could be raised or lowered. Results showed the contact angle of the stearic‐acid‐treated sands decayed with time and that there was a significant increase with stearic acid concentration. For all samples, the decay in apparent contact angle could be described with a continuous breakdown model. The stearic‐acid‐treated sands showed a significant increase in contact angle with relative humidity at a temperature of 10 and 20 °C. These differences diminished with increasing temperature. Similar results were seen for the dune sands and composts. Despite the influence of temperature and humidity on contact angles, there was no significant change in the rate at which the contact angle decayed in any sample. Absolute humidity was found to provide a more relevant indicator than relative humidity when assessing the influence of humidity on repellency over a range of temperatures. The contact angle initially increased with absolute humidity before plateauing owing to the confounding effect of temperature. Copyright © 2014 John Wiley & Sons, Ltd. 相似文献
19.
地下岩石由岩石骨架和孔隙流体组成,通常流体含黏性.地震波在地下介质中传播时受岩石骨架和黏性流体的影响会呈现出复杂的变化.本文将流、固体位移和应力连续作为边界条件,推导出含黏性流体孔隙介质分界面上反透射系数方程;通过建立上层为饱油、下层为饱盐水的砂岩孔隙介质模型,开展反透射系数特征研究,分别分析不同频率、不同黏滞系数条件下,含黏性流体孔隙介质分界面上反透射系数随入射角的变化.研究表明,孔隙介质分界面上和等效介质分界面上的反透射系数分别随入射角的变化趋势基本一致,说明方程推导和数值计算的正确性;快纵波反透射系数受频率、流体黏性的影响较小,而快横波反透射系数在一定入射角范围内受频率、流体黏性的影响比较大;由于黏性孔隙流体的作用,慢纵波和慢横波的反透射系数受入射角、频率及流体黏性的影响都很大. 相似文献
20.
We develop the true‐amplitude prestack migration of multicomponent data based on the use of elastic Gaussian beams for walkaway vertical seismic profile (VSP) acquisition systems. It consists in a weighted summation of multishot data with specific weights, computed by tracing elastic Gaussian beams from each imaging point of the target area towards the sources and receivers. Each pair of beams may be connected with either a pair of P‐rays (PP‐image) or the P‐ray towards sources and the S‐ray to receivers (PS‐image) and is uniquely determined by dip (the angle of the bisector between the rays and the vertical direction) and opening (the angle between the rays) angles. Shooting from the bottom towards the acquisition system helps to avoid well‐known troubles, in particular multipathing for the imaging conditions in complex velocity models. The ability to fix the dip angle and implement summation over opening angles leads to the so‐called selective images that contain mostly interfaces with desired slopes. On the other hand, a set of images computed for a range of opening angles by summation over all available dip angles is used as input of an AVO‐like inversion procedure for the recovery of elastic parameters. The feasibility of this imaging procedure is verified by synthetic data for 2D realistic elastic models. 相似文献