Correlation analysis of statistical facies fault models

Doklady Earth Sciences - This paper addresses seismic imaging of fault zones and analysis of the seismic data with the use of the fault facies model developed at Uni Research CIPR. Simulated and...     

Atmospheric Corrections Using MODTRAN for TOA and Surface BRDF Characteristics from High Resolution Spectroradiometric/ Angular Measurements from a Helicopter Platform

l. IntroductionThe cIassicaI definition of bi--directional reflectance--distribution function (BRDF) is aderivative, distribution function, relatlng the irradiance incident from one given direction tolts contribution to the radiance renected in another direct1on (N1codemus et al., l977).f r (0,, rp,; 0,, rp, )= dL, (0,, P,; 0,, 9,; E, )/ dE, (0,, 9,) [sr-- 1 ], (l)where 0 (zenith angle) and 9 (azimuth angle) together indicate a direction, the subscr1pt i indi-cates quantities associated wi…     

热液金刚石压腔在地质流体研究中的应用

热液金刚石压腔(HDAC)是专为模拟地壳温压条件下的地质作用而设计的，它尤其适用于观测水或其它流体与地质物质之间的相互作用。HDAC可时-190～1200℃，0～10GPa的热液体系进行实验，并可在实验的温压条件下．用各种先进的光学方法分析样品，更可以把实验的全程录像存档。充满流体的HDAC本身就可当做是一个人工合成的流体包裹体。因此它可以用来研究流体的状态方程和相关系。它又可时其它流体包裹体样品施加外压力，因此在热分析的过程中可免除包裹体的膨胀或爆破的困境。HDAC可应用到颇为宽广的温压范围，它已被广泛地用来观测各种化学体系的临界现象，包括在地质方面特别有用的含水硅酸盐体系。HDAC也可与同步辐射X光源相结合，而取得各种金属或稀土元素水溶液的x射线吸收精细结构(X-ray absorption fine structure；XAFS)光谱，因而时在热液里的金属或稀土元素络合物的组分和结构提供了最基本的资料。然而，X光的强度在透过金刚石时，因绕射和吸收而大大地减弱，因此应用一般的HDAC来获取那些吸收边在10keV以下的元素的XAFS光谱颇为困难。目前已有两种改良式的HDAC解除了这方面的困难，而时在元素周期表上的第一排过渡性金属元素和稀土元素的水溶液，提供清晰的XAFS光谱。这些资料可用来研究金属或稀土元素络合物在地质热液里的特性，及其在元素迁移和成矿作用方面的效应。而这些元素在地质应用方面特别重要。     

3D diffraction imaging of fault and fracture zones via image spectral decomposition of partial images

The paper presents 3D diffraction imaging based on the spectral decomposition of a different combination of selective or partial images. These images are obtained by the pre-stack asymmetric migration procedure, which is weighted data summation. Spectral decomposition is done in the Fourier domain with respect to spatial dip and azimuth angles. Numerical examples with the application of different workflows for the synthetic and real data examples demonstrate detailed reliable reconstruction of the fractured zones and reliable reconstruction of fracture orientation on synthetic and real 3D data examples.     

Finite-difference algorithm with local time-space grid refinement for simulation of waves

This paper presents a new approach to a local time-space grid refinement for a staggered-grid finite-difference simulation of waves. The approach is based on approximation of a wave equation at the interface where two grids are coupled. As no interpolation or projection techniques are used, the finite-difference scheme preserves second order of convergence. We have proved that this approach is low-reflecting, the artificial reflections are about 10^− 4 of an incident wave. We have also shown that if a successive refinement is applied, i.e. temporal and spatial steps are refined at different interfaces, this approach is stable.     

On the interface error analysis for finite difference wave simulation

A common way to construct a finite difference scheme is to satisfy a desired order of approximation, typically as high as possible. For linear wave propagation problem approximation together with stability delivers convergence of the same order as approximation. If a wave propagation proses is considered convergence to a plane wave solution can be derived analytically by means of the dispersion analysis. However, mentioned techniques are applicable only to homogeneous media and provide no knowledge of reflection/transmission coefficients. In this paper we prove that the only way to get second order accuracy of the solution for media with discontinuous parameters is to use a conservative finite difference scheme of the second order, and the only way to do this is to use the arithmetic mean for the density and the harmonic mean for the bulk modulus in the vicinity of the interface.     

True amplitude elastic Gaussian beam imaging of multicomponent walkaway vertical seismic profiling data

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.     

一种新型人工“流体包裹体”：融合二氧化硅毛细管技术

利用融合二氧化硅毛细管技术制作了纯H2O体系、纯CO2体系、H2O-NaCl体系和H2O-CO2体系的人工包裹体样品,并对样品进行了显微测温和激光拉曼光谱测试工作。实验结果显示毛细管样品中的流体成分具有代表性,而且常规的流体包裹体显微测温和显微激光拉曼光谱分析技术完全适于毛细管样品的测试。对样品的显微测温和拉曼光谱研究...     

True amplitude imaging by inverse generalized Radon transform based on Gaussian beam decomposition of the acoustic Green's function

True amplitude migration is one of the most important procedures of seismic data processing. As a rule it is based on the decomposition of the velocity model of the medium into a known macrovelocity component and its sharp local perturbations to be determined. Under this decomposition the wavefield can be considered as the superposition of an incident and reflected/scattered waves. The single scattering approximation introduces the linear integral operator that connects the sharp local perturbations of the macrovelocity model with the multishot/multioffset data formed from reflected/scattered waves. We develop the pseudoinverse of this operator using the Gaussian beam based decomposition of acoustic Green's functions. The computation of this pseudoinverse operator is done pointwise by shooting Gaussian beams from the target area towards the acquisition system. The numerical implementation of the pseudoinverse operator was applied to the synthetic data Sigsbee2A. The results obtained demonstrate the high quality of the true amplitude images computed both in the smooth part of the model and under the salt body.     

Numerical simulation of seismic waves in models with anisotropic formations: coupling Virieux and Lebedev finite-difference schemes

Anisotropy is widespread in the Earth’s interior. However, there is a number of models where anisotropic formations comprise as few as 10–20?% of the volume, and this includes fractured reservoirs, thin-layered packs, etc. while the major part of the medium is isotropic. In this situation, the use of computationally intense anisotropy-oriented approaches throughout the computational domain is prodigal. So this paper presents an original advanced finite-difference algorithm based on the domain decomposition technique with individual scheme used inside subdomains. It means that the standard staggered grid scheme or the Virieux scheme is used in the main part of the model which is isotropic, while the anisotropy-oriented Lebedev scheme is utilized inside domains with anisotropic formations. Finite-difference consistency conditions at the artificial interface where the schemes are coupled are designed to make the artificial reflections as low as possible, namely, for the second-order scheme, the third order of convergence of the reflection coefficients is proved.