共查询到20条相似文献，搜索用时 78 毫秒
1.
The automatic detection of geological features such as faults and channels is a challenging problem in today's seismic exploration industry. Edge detection filters are generally applied to locate features. It is desirable to reduce noise in the data before edge detection. The application of smoothing or low‐pass filters results in noise suppression, but this causes edge blurring as well. Edge‐preserving smoothing is a technique that results in simultaneous edge preservation and noise suppression. Until now, edge‐preserving smoothing has been carried out on rectangular sampled seismic data. In this paper, an attempt has been made to detect edges by applying edge‐preserving smoothing as a pre‐processing step in the hexagonally sampled seismic‐data spatial domain. A hexagonal approach is an efficient method of sampling and has greater symmetry than a rectangular approach. Here, spiral architecture has been employed to handle the hexagonally sampled seismic data. A comparison of edge‐preserving smoothing on both rectangular and hexagonally sampled seismic data is carried out. The data used were provided by Saudi Aramco. It is shown that hexagonal processing results in well‐defined edges with fewer computations. 相似文献
2.
3.
三维复杂山地条件下的各种地震波型的走时计算技术,可以直接用于复杂山地区域地震波运动学特性的分析、地震数据采集观测系统的设计以及直接基于三维复杂地表的地震数据处理技术的研发.为了在三维复杂地表条件下准确、灵活且稳定地计算各种地震波型的走时,提出一种多级次群推进迎风混合法.该算法利用不等距迎风差分法简洁稳定地处理三维复杂地表及附近的局部走时计算问题,利用计算精度不错的迎风双线性插值法处理绝大部分均匀正方体网格中的局部走时计算问题,利用群推进法模拟三维复杂地表条件下地震波前的扩展问题,利用多级次算法处理各种类型的地震波的走时计算问题.算法分析和计算实例表明:新方法具有很好的计算精度与效率,且能灵活稳定地处理三维复杂地表复杂介质条件下的多波型走时计算问题. 相似文献
4.
Numerical simulation of 2D wave propagation on unstructured grids using explicit differential operators 总被引：2，自引：0，他引：2
We present a numerical method of simulating seismic wave propagation on unstructured 2D grids. The algorithm is based on the velocity–stress formulation of the elastic wave equation and therefore uses a staggered grid approach. Unlike finite-element or spectral-element methods, which can also handle flexible unstructured grids, we use explicit differential operators for the calculation of spatial derivatives in each time step. As shown in previous work, three types of these operators are used, and their particular performance is analysed and compared with standard explicit finite-difference operators on regular quadratic and hexagonal grids. Our investigations are especially focused on the influence of grid irregularity, sampling rate (i.e. gridpoints per wavelength) and numerical anisotropy on the accuracy of numerical seismograms. The results obtained from the various methods are therefore compared with analytical solutions. The algorithm is then applied to a number of models that are difficult to handle using (quasi-)regular grid methods. Such alternative techniques may be useful in modelling the full wavefield of bodies with complex geometries (e.g. cylindrical bore-hole samples, 2D earth models) and, because of their local character, they are well suited for parallelization. 相似文献
5.
回顾了城市活断层探测现状,指出了采用二维地震探测在活断层精确定位中存在中浅层小构造控制程度较差,在地层倾角较大地区断层归位不够准确,受建筑物和其它障碍物限制致使测线布设与断层走向斜交造成断层定位不准,以及由于测网密度限制造成的断层交汇处、断层分段点、端点位置控制程度较差等主要问题．从理论上分析了三维地震勘探的优势,认为三维地震勘探可以较好地解决由于二维地震勘探本身的技术缺陷和地表环境、地质条件影响所造成的断层定位不准的问题．通过实例剖析,从野外采集难点, 仪器选择, 观测系统定义, 优化施工设计, 野外质量监控, 特殊数据处理手段, 以及取得的良好地质效果等方面进行详细阐述,论证了利用三维地震探测在城市进行活断层精确定位中的可行性及其意义． 相似文献
6.
加权近似解析离散化(WNAD) 方法是近年发展的一种在粗网格步长条件下能有效压制数值频散的数值模拟技术. 在地震勘探的实际应用中, 不是所有情况都适合使用空间大网格步长. 为适应波场模拟的实际需要, 本文给出了求解波动方程的非一致网格上的WNAD算法. 这种方法在低速区、介质复杂区域使用细网格, 在其他区域采用粗网格计算. 在网格过渡区域, 根据近似解析离散化方法的特点, 采用了新的插值公式, 使用较少的网格点得到较高的插值精度. 数值算例表明, 非一致网格上的WNAD方法能够有效压制数值频散, 显著减少计算内存需求量和计算时间, 进一步提高了地震波场的数值模拟效率. 相似文献
7.
The processing of a small-scale, very high resolution (VHR) shallow marine 3D data volume is described. The data were acquired over a small clay diapir, on the river Schelde, in 1990. Using an array of 12 dual-channel microstreamers towed from a catamaran, a network of 1 m × 1 m bins could be produced over an area of 50 m × 180 m (< 100 m penetration). Positioning was performed with an auto-tracking laser ranging system, assuring an absolute accuracy of a few decimetres.
Preliminary processing steps included tidal correction and multiple removal. An important step concerned the application of 3D prestack Kirchhoff depth migration. Indeed this processing allows easy handling of the exact positions of both source and receivers as the latter were not set out on a conventional regular grid due to navigation difficulties. Because of the restricted data volume and the more-or-less stratified medium, a 1D velocity model could be used. This allowed a considerable simplification of the migration algorithm, based on summation. Traveltimes were calculated only once, using a 2D time grid with 0.1 m intervals.
This migration method proved very efficient, greatly improving the seismic image, and involved only limited CPU time on a small computer (Sparc 10 workstation). It clearly demonstrates that advanced seismic processing can form a valuable and economically feasible tool for VHR shallow subsurface 3D seismics, as long as the velocity field is not too complex. This method should therefore no longer be restricted to large computers and hydrocarbon exploration, but should also become a routine for VHR 3D shallow seismic work. 相似文献
Preliminary processing steps included tidal correction and multiple removal. An important step concerned the application of 3D prestack Kirchhoff depth migration. Indeed this processing allows easy handling of the exact positions of both source and receivers as the latter were not set out on a conventional regular grid due to navigation difficulties. Because of the restricted data volume and the more-or-less stratified medium, a 1D velocity model could be used. This allowed a considerable simplification of the migration algorithm, based on summation. Traveltimes were calculated only once, using a 2D time grid with 0.1 m intervals.
This migration method proved very efficient, greatly improving the seismic image, and involved only limited CPU time on a small computer (Sparc 10 workstation). It clearly demonstrates that advanced seismic processing can form a valuable and economically feasible tool for VHR shallow subsurface 3D seismics, as long as the velocity field is not too complex. This method should therefore no longer be restricted to large computers and hydrocarbon exploration, but should also become a routine for VHR 3D shallow seismic work. 相似文献
8.
Reconstruction of seismic data is routinely used to improve the quality and resolution of seismic data from incomplete acquired seismic recordings. Curvelet‐based Recovery by Sparsity‐promoting Inversion, adapted from the recently‐developed theory of compressive sensing, is one such kind of reconstruction, especially good for recovery of undersampled seismic data. Like traditional Fourier‐based methods, it performs best when used in conjunction with randomized subsampling, which converts aliases from the usual regular periodic subsampling into easy‐to‐eliminate noise. By virtue of its ability to control gap size, along with the random and irregular nature of its sampling pattern, jittered (sub)sampling is one proven method that has been used successfully for the determination of geophone positions along a seismic line. In this paper, we extend jittered sampling to two‐dimensional acquisition design, a more difficult problem, with both underlying Cartesian and hexagonal grids. We also study what we term separable and non‐separable two‐dimensional jittered samplings. We find hexagonal jittered sampling performs better than Cartesian jittered sampling, while fully non‐separable jittered sampling performs better than separable jittered sampling. Two other 2D randomized sampling methods, Poisson Disk sampling and Farthest Point sampling, both known to possess blue‐noise spectra, are also shown to perform well. 相似文献
9.
Fei Li Tao Xu Minghui Zhang Zhenbo Wu Chenglong Wu Zhongjie Zhang Jiwen Teng 《地震科学（英文版）》2014,27(2):127-136
Seismic traveltime tomographic inversion has played an important role in detecting the internal structure of the solid earth. We use a set of blocks to approximate geologically complex media that cannot be well described by layered models or cells. The geological body is described as an aggregate of arbitrarily shaped blocks, which are separated by triangulated interfaces. We can describe the media as homogenous or heterogeneous in each block. We define the velocities at the given rectangle grid points for each block, and the heterogeneous velocities in each block can be calculated by a linear interpolation algorithm. The parameters of the velocity grid positions are independent of the model parameterization, which is advantageous in the joint inversion of the velocities and the node depths of an interface. We implement a segmentally iterative ray tracer to calculate traveltimes in the 3D heterogeneous block models. The damped least squares method is employed in seismic traveltime inversion, which includes the partial derivatives of traveltime with respect to the depths of nodes in the triangulated interfaces and velocities defined in rectangular grids. The numerical tests indicate that the node depths of a triangulated interface and homogeneous velocity distributions can be well inverted in a stratified model. 相似文献
10.
三维声波方程相比二维声波方程能够更好的模拟三维空间的地震波传播,模拟标量近似下的弹性波在三维复杂介质的传播过程.基于非规则网格的正演模拟方法的格子法可以处理很好的刻画起伏地表、速度间断面等复杂构造,但是这类方法需要大量的几何描述来描述网格.本文提出了三维六面体双重网格的格子法来模拟声波方程,一方面该方法继承了格子法能够灵活处理自由表面和速度间断面的特性.另一方面,该方法通过双重网格的实现极大的减少了几何描述文件的大小,可以最大的实现GPU加速,实现粗粒度并行,在节省了几何描述空间的同时达到了很高的加速比. 相似文献
11.
在许多地震反演和偏移成像方法中,都要涉及到射线路径和旅行时的计算.本文将波前面三角形网格剖分和三维波前重建法射线追踪技术结合使用,实现了射线路径和旅行时的准确快速计算.三维波前重建法射线追踪过程中可以保证稳定合理的射线密度,克服了常规射线追踪方法存在阴影区的问题.波前面三角形网格剖分在描述和拆分波前面时更加准确有效,而且不需太多的网格数目,从而提高了射线追踪的精度和效率.该方法在三维复杂构造成像方面有独特的优势,目前在实际的Kirchhoff 偏移中的已经有相关应用. 相似文献
12.
3-D acoustic wave equation forward modeling with topography 总被引：1，自引：1，他引：1
In order to model the seismic wave field with surface topography, we present a method of transforming curved grids into rectangular grids in two different coordinate systems. Then the 3D wave equation in the transformed coordinate system is derived. The wave field is modeled using the finite-difference method in the transformed coordinate system. The model calculation shows that this method is able to model the seismic wave field with fluctuating surface topography and achieve good results. Finally, the energy curves of the direct and reflected waves are analyzed to show that surface topography has a great influence on the seismic wave＇s dynamic properties. 相似文献
13.
表面多次波是海洋地震勘探中的主要问题.目前,二维数据驱动的表面多次波压制技术(SRME)已经比较成熟,并且已经成为工业界压制海洋表面多次波的主流方法.但是由于二维SRME算法没有考虑横测线方向上多次波的贡献,导致在处理实际三维海洋资料时存在比较大的误差.将二维SRME算法扩展到三维空间后可以得到三维SRME算法,但是由于目前实际采集的三维海洋资料的观测系统存在拖缆漂移,而且横测线方向采样过于稀疏,直接应用三维SRME算法无法准确预测表面多次波.本文提出的通过数据规则化配合稀疏反演的三维表面多次波压制方法能够解决这种实际资料和三维SRME算法之间的矛盾.本文通过研究数据规则化与反规则化技术,使得数据分布满足三维SRME的要求;通过研究稀疏反演技术,有效解决了横测线方向采样稀疏对于多次波预测的影响,三维实际海洋资料的应用结果验证了方法的有效性和可行性. 相似文献
14.
A crucial point in the processing of 3D seismic data is the migration step, both because of its 3D nature and the computational cost involved. The efficiency and accuracy of 3D migration are determined by the wavefield extrapolation technique employed. Wavefield extrapolation based on second-order differential operators of variable-length is very efficient and accurate at the same time. Compared to migration based on the McClellan transform and operator splitting, the use of variable-length second-order differential operators offers significant advantages. The 3D migration operator has an almost perfect circular symmetry. No positioning errors in the 45° azimuth between the in-line and cross-line directions are evident. The method is, in practice, only limited by spatial aliasing and does not require expensive interpolation of data to reduce numerical artifacts. This reduces the computational cost of 3D one-pass depth migration by a large factor. 相似文献
15.
Mark Grasmueck Miquel Coll Quintà Kenri Pomar Gregor P. Eberli 《Geophysical Prospecting》2013,61(5):907-918
Vertical fractures with openings of less than one centimetre and irregular karst cause abundant diffractions in Ground‐Penetrating Radar (GPR) records. GPR data acquired with half‐wavelength trace spacing are uninterpretable as they are dominated by spatially undersampled scattered energy. To evaluate the potential of high‐density 3D GPR diffraction imaging a 200 MHz survey with less than a quarter wavelength grid spacing (0.05 m × 0.1 m) was acquired at a fractured and karstified limestone quarry near the village of Cassis in Southern France. After 3D migration processing, diffraction apices line up in sub‐vertical fracture planes and cluster in locations of karstic dissolution features. The majority of karst is developed at intersections of two or more fractures and is limited in depth by a stratigraphic boundary. Such high‐resolution 3D GPR imaging offers an unprecedented internal view of a complex fractured carbonate reservoir model analogue. As seismic and GPR wave kinematics are similar, improvements in the imaging of steep fractures and irregular voids at the resolution limit can also be expected from high‐density seismic diffraction imaging. 相似文献
16.
3D地震数据不规则采样缺失重建是地震勘探数据处理流程中的重要问题.本文提出了一种基于具有保幅特性的非均匀高阶抛物Radon变换（NHOPRT）地震数据重建方法.在最小二乘反演方程中引入Delaunay三角网格剖分来计算空间不规则加权系数,从而获得最接近完整规则数据的高阶抛物Radon变换域系数.在用SVD求解反演方程过程中,利用高阶抛物Radon变换算子在频率域为指数函数,具有线性可分解特性,将二维空间的高阶抛物Radon变换算子分解为两个独立的一维空间变换算子,减小了变换算子的矩阵大小,从而很大程度地提高了计算效率.理论模型和实际地震数据重建测试证明了本文方法的有效性以及实用性. 相似文献
17.
3‐D shallow‐water seismic survey evaluation and design using the focal‐beam method: a case study offshore Abu Dhabi 
下载免费PDF全文
下载免费PDF全文 For 3‐D shallow‐water seismic surveys offshore Abu Dhabi, imaging the target reflectors requires high resolution. Characterization and monitoring of hydrocarbon reservoirs by seismic amplitude‐versus‐offset techniques demands high pre‐stack amplitude fidelity. In this region, however, it still was not clear how the survey parameters should be chosen to satisfy the required data quality. To answer this question, we applied the focal‐beam method to survey evaluation and design. This subsurface‐ and target‐oriented approach enables quantitative analysis of attributes such as the best achievable resolution and pre‐stack amplitude fidelity at a fixed grid point in the subsurface for a given acquisition geometry at the surface. This method offers an efficient way to optimize the acquisition geometry for maximum resolution and minimum amplitude‐versus‐offset imprint. We applied it to several acquisition geometries in order to understand the effects of survey parameters such as the four spatial sampling intervals and apertures of the template geometry. The results led to a good understanding of the relationship between the survey parameters and the resulting data quality and identification of the survey parameters for reflection imaging and amplitude‐versus‐offset applications. 相似文献
18.
陕北富县地区以前的地震勘探沿沟布设弯曲测 线,平面展布呈不规则树枝状,难以形成闭合回 路,很难对二维地震资料在平面上进行储层反 演。2001年在陕北富县探区进行野外采集方法和 地震处理技术攻关,通过采用极其灵活的可变线 元观测系统实现直测线过塬,使测线成网连片; 资料处理中采用基于模型反演的折射波静校正、 地表一致性振幅补偿和反褶积、叠前噪音衰减等 技术,为后续的储层预测、低幅度构造解释和综 合地质研究提供了高保真的地震资料;由于测线 的连片成网,通过伪测井的方法对无井测线进行 控制,最终实现储层反演的连片解释。 相似文献
19.
三维起伏地表条件下的地震波走时计算技术是研究三维起伏地表地区很多地震数据处理技术的基础性工具.为了获得适应于任意三维起伏地表且计算精度高的走时算法,提出三维不等距迎风差分法.该方法采用不等距网格剖分三维起伏地表模型,通过在迎风差分格式中引入不等距差分格式、Huygens原理及Fermat原理来建立地表附近的局部走时计算公式,并通过在窄带技术中设定新的网格节点类型来获得三维起伏地表条件下算法的整体实现步骤.精度及算例分析表明:三维不等距迎风差分法具有很高的计算精度且能够适应于任意三维起伏地表模型. 相似文献
20.
New alternative formulations of exact boundary conditions for arbitrary three-dimensional (3D) free-surface topographies on seismic media have been derived. They are shown to be equivalent to previously published formulations, thereby verifying the validity of each set of formulations. The top of a curved grid represents the free-surface topography while the interior of the grid represents the physical medium. We assume the velocity–stress version of the viscoelastic wave equations to be valid in this grid before transforming the equations to a rectangular grid. In order to perform the numerical discretization we apply the latter version of the equations for seismic wave propagation simulation in the medium. The numerical discretization of the free-surface topography boundary conditions by second-order finite differences (FDs) is shown, as well as the spatially unconditional stability of the resulting system of equations. The FD order is increased by two for each point away from the free surface up to eight, which is the order used in the interior. We use staggered grids in both space and time and the second-order leap-frog and Crank– Nicholson methods for wavefield time propagation. An application using parameters typical of teleseismic earthquakes and explosions is presented using a 200 × 100 km2 area of real topography from southwestern Norway over a homogeneous medium. A dipping plane wave simulates a teleseismic P-wave incident on the surface topography. Results show conversion from P- to Rg- (short period fundamental mode Rayleigh) waves in the steepest and/or roughest topography, as well as attenuated waves in valleys and fjords. The codes are parallelized for simulation on fast supercomputers and PC-clusters to model high frequencies and/or large areas. 相似文献