MIGRATION OF MIXED-MODE VSP WAVEFIELDS1

Two-dimensional VSP surveys are often conducted to provide structural illumination of the subsurface away from the borehole. The illumination is achieved through offsetting the source with respect to the downhole geophone. This inevitably gives rise to mode-conversions in both downgoing and upgoing wavefields. Migration of mixed-mode wavefields is complex because the velocity profile used for wavefield extrapolation is valid only for a particular propagation mode; the other mode always propagates at a different velocity. It is therefore advisable to separate the wave-types (P-wave and SV-wave) prior to migration. This may be achieved through wavemode filtering, a multichannel process which exploits the relation between propagation velocity, slowness of events at the recording array and particle motion. The necessary information about particle motion is available only if the VSP data are acquired with a three-component downhole geophone assembly. The wavemode filter partitions wave-types at the recording array; it provides no information about the various changes of propagation mode experienced by the energy as it travels from source to geophone. For the purpose of migration, the intermediate modes of propagation must be deduced. Much of the energy arriving at the receivers is P-wave which has followed the P-wave velocity profile from the source. It can therefore be imaged by conventional (Kirchhoff) migration. As an example of SV-wave imaging, a common mode-code is P-wave from source to reflector and SV-wave from reflector to geophone. Migration of such data calls for back-propagation of the geophone array wavefield, at SV-wave velocity, to the point in the subsurface where it is time-coincident with the forward propagated downwave, at P-wave velocity.     

FREQUENCY-SELECTIVE DESIGN OF THE KIRCHHOFF MIGRATION OPERATOR1

Integral migration techniques perform a sum over an aperture of input traces to obtain output at a single point. The length of the aperture is limited by a spatial Nyquist criterion, which typically prohibits imaging very steep dips at very high frequencies without generating severe migration artifacts (migration operator aliasing). For time-domain Kirchhoff migration, this can be a fatal shortcoming. The standard way to address this problem is to interpolate traces spatially before migration. This reduces the trace spacing, thereby increasing the frequency content which can be migrated without aliasing at steep dips. An alternative remedy to the operator aliasing problem is to modify the phase response of the Kirchhoff migration operator. This operator is frequency-selective across the migration aperture: it passes all temporal frequencies of the input traces in the innermost portion of the aperture (referring to the shallow dips), and gradually cuts out the higher frequencies as it approaches the outer portion of the aperture. Thus, while all frequencies of the input data contribute to the shallow-dip portion of the migrated image, only the permissible low frequencies of the input data contribute to imaging the steepest dips. Using a simple realization of a frequency-selective Kirchhoff migration operator, this technique is illustrated on a synthetic data set involving greater than vertical dips.     

THE DETERMINATION OF THE SEISMIC QUALITY FACTOR Q FROM VSP DATA: A COMPARISON OF DIFFERENT COMPUTATIONAL METHODS1

Ten methods for the computation of attenuation have been investigated, namely: amplitude decay, analytical signal, wavelet modelling, phase modelling, frequency modelling, rise-time, pulse amplitude, matching technique, spectral modelling and spectral ratio. In particular, we have studied the reliability of each of these methods in estimating correct values of Q using three synthetic VSP seismograms for plane P-waves with different noise contents. The investigations proved that no single method is generally superior. Rather, some methods are more suitable than others in specific situations depending on recording, noise or geology. The analytical signal method has been demonstrated to be superior if true amplitude recordings are available. Otherwise spectral modelling or, in the ‘ noise-free’ case the spectral ratio method, is optimal. Finally, two field VSPs in sediments are investigated. Only in the case of the highest quality VSP can significant information be deduced from the computed attenuation.     

REFLECTION AND POLARIZATION OF TUBE WAVES AS SEEN IN VSP DATA1

P-wave and S-wave data acquired with vertical seismic profiling (VSP) often include tube waves propagating in the borehole, although considerable efforts are generally made to ensure that these waves are not recorded. However, several theoretical studies have indicated that tube waves could provide important information about the rock formation and thus should not be considered as pure noise. In order to study some of these aspects experimentally, tube waves were acquired by VSP in a well in the Paris Basin both before and after casing. A sparker was used as source inside the borehole, which ensured that the data recorded contained high-amplitude tube waves. It is shown that the casing is an obstacle which prevents the study of formation parameters, and thus further tube-wave acquisitions should be carried out in open holes only. The before-casing tube-wave reflection log is compared to a synthetic log computed from the sonic log. The high resolution of the tube waves is of particular interest, revealing layers that are too thin to be detected in body-wave surveys. It has recently been suggested that the projection of the tube-wave polarization in the horizontal plane can be used to determine directions of stress-induced anisotropy in the rock formation. Strong polarization anomalies are observed in the data sets but are attributed to tool problems rather than any rock-formation feature.     

ENHANCED MIGRATION OF SEISMIC DATA1

The so-called ‘enhanced migration’ which uses diffraction tomography as the ‘repair tool’ for correction of amplitudes (reflection coefficients) of migrated sections is discussed. As with any linearized procedure, diffraction tomography requires knowledge of the initial model. It is suggested that the initial model is taken as the migrated image. It will be demonstrated that diffraction tomography applied to the data residuals improves the amplitudes of the migrated images. Migration is redefined as the reconstruction of the wavefront sets of distributions (reflection interfaces), and the inversion process as tomographic correction of migrated images.     

APPLICATION OF ELASTIC MODELLING IN PROCESSING AND INTERPRETATION OF VSP DATA: A CASE HISTORY1

Data from offshore Norway is used to study applications of elastic VSP modelling in detecting shear waves and observing the effects of successive mode conversion in field-recorded VSP data. The shear-wave velocities and densities from log data are used in conjunction with compressional wave velocities determined from surface seismic and log data in the VSP modelling. The time domain non-normal incidence elastic VSP modelling technique of Aminzadeh and Mendel is used as the modelling algorithm. Two surface seismograms are computed first. One is the vertical component and the other is the horizontal component for plane waves that have specified incident angles. A downward continuation method is then applied to generate seismograms at different depth points. The collection of these seismograms constitutes non-normal incidence VSPs. Both vertical and horizontal components of VSP data can be obtained by this procedure. In this paper non-normal incidence VSPs are generated for a 12.5° incident plane wave. The modelling results of layered earth systems of thin layers and thick layers are both compared with field data, and the effect of mode conversions in thin layers is observed. Several events in the field data can be explained by this elastic VSP modelling. Comparison of the model data and field data enabled a probable tube wave or out-of-plane event to be identified, the removal of which significantly improved the final VSP section. This study also shows how the VSP data helped the interpretation of the surface 3D data.     

大地电磁资料偏移成像

在研究Zhdanov的二维电磁场偏移成像的基础上, 针对它在某些方面的不足, 研究并提出了改进算法. 首先, 用一维反演结果确定初始背景电阻率, 然后在延拓计算过程中, 利用迭代算法进行逐层修正、 延拓计算, 从而获得较准确的电阻率; 第二, 研究了有限差分方程的改进算法, 根据大地电磁场变化的特点, 设计了纵向以等比级数变网格的迭代算法. 通过以上两方面技术的改进, 经理论模型和实际资料试算验证, 取得了较好的结果.      

MIGRATION OF VSP DATA BY RAY EQUATION EXTRAPOLATION IN 2-D VARIABLE VELOCITY MEDIA*

The standard Kirchhoff algorithm can be generalized for migration of pre-stack finite-offset data from variable-velocity media. The concentric ellipses over which the data are spread in constant velocity media become significantly distorted (even multi-valued) in the variable velocity case. The specific shapes can be explicitly defined by kinematic extrapolation of the source and recorded wave fields with the ray equation. The use of Kirchhoff migration with a surface source and a subsurface recorder requires that two sets of Kirchhoff loci be superimposed. For each trace, the first set of loci is computed with the source and the actual recorder position as foci; the second set is computed with the source and the virtual recorder position as foci. This dual procedure explicitly incorporates the primary diffracted energy and the free-surface reflections, respectively. Implementation involves the construction of a virtual medium, lying above the free surface, with a velocity distribution that is the mirror image of the actual distribution below the free surface. Ray-equation extrapolation is performed through the real/virtual boundary. The resulting image is produced in a split form, with all the contributions of the primary reflected and diffracted energy lying in the lower ‘real’ half and all the contributions of the energy that was reflected at the free surface lying in the upper ‘virtual’ half. The final image is produced by folding the split image about the free surface and adding the two halves. A practical advantage is that the origin of various contributions (and artifacts) can be more readily identified (for interpretation or removal) in the split images. The ray-equation pre-stack migration algorithm is very general. It is applicable to all source-recorder geometries and variable velocity media and reduces exactly to the standard Kirchhoff algorithm when applied to zero or finite-offset surface survey data. The algorithm is illustrated by application to VSP data. For the VSP geometry, the algorithm does not require any specific trace spacing (in depth) and can be used for data from deviated as well as vertical holes.     

ATTENUATION OF RANDOM NOISE BY 2-D AND 3-D CDP STACKING AND KIRCHHOFF MIGRATION*

Three-dimensional seismic surveys, in general, do not need the same high degree of CDP coverage as 2-D surveys to achieve a certain signal-to-noise ratio after migration. This can be shown theoretically for Kirchhoff migration and laterally uncorrelated noise. More precisely, there exists a formal relationship between the multiplicity of CDP coverage of a 3-D survey and that of a 2-D survey with the same signal-to-uncorrelated-noise ratio. Frequency and aperture are parameters in the corresponding expression. Heuristically the relationship can be obtained by applying the concept of the Fresnel zone. Though the mathematics in this paper refer to laterally uncorrelated noise, the underlying concepts can probably also be used for weakly correlated noise, e.g., for multiple reflections and for the low-frequency remnants of surface waves.     

有限差分格式与物理量守恒性

为了成功地开展大气数值模拟,如何构造相应的离散格式是一项十分重要的工作.本文以大气球面浅水方程为例,从不同形式的微分方程组出发,分别构造了能量守恒格式、准辛格式和多守恒格式,并作相应的数值检验和比较,从中可以清楚地看到三类格式的优劣.这三类格式构造方法不同,所保持的物理守恒性不完全相同,数值计算时所需机时也不同,这些比较结果为实际工作的不同使用需求提供了选择依据.     

对波达波夫和Pride震电波方程组的对比分析

用Biot介质参数说明了波达波夫震电波方程组中弹性动力学 参数的含义，解释了第一类和第二类震电效应的意义，在忽略第一类震电效应条件下将该方 程组与Pride方程组进行比较，说明了二者在描述第二类震电效应方面的异同点. 同时指出 ：波达波夫方程组忽略了流体与固体的耦合质量；方程中的黏性耗散项丢掉了一个孔隙度因 子，依据该方程组计算出的弹性波和转换电场的幅度将偏大；边界条件之一存在错误，会影 响对波在界面上的反射透射规律的描述.     

阿克苏跨断层形变观测资料分析及其与地震关系的研究

结合阿克苏台周边5级以上中强地震和100 km范围内的4级地震,对阿克苏台跨断层形变资料进行分析.分析结果表明,断层的垂直、水平位移形变异常变化与中强地震活动有明显的对应关系.在监控区350 km范围内的中强地震前,垂直形变异常形态表现为次波效应,而水平形变在100 km范围内的4级地震前则出现明显的伸长或缩短的异常特征.     

A STUDY OF DIFFERENT METHODS OF WAVEFIELD SEPARATION FOR APPLICATION TO VSP DATA1

During the last couple of years there has been much research in the area of wavefield separation of borehole seismic data, and several articles have been published on various separation techniques. Methods involving the application of two-dimensional Fourier transformation, the Radon transformation, multi-level median filters or optimal filters, are all suggested as possible approaches to the wavefield separation problem. This paper compares some of these methods commonly used in the industry. The theories of the chosen methods are described to see how they are related. Using the different methods on synthetic and real data, we show how this theoretical relation is reflected in the relatively similar results obtained. We also show how the different filters treat coherent and random noise.     

绝对和相对输入能量谱对比及延性系数的影响研究

利用266条强震记录,在研究绝对输入能量谱和相对输入能量谱衰减规律的基础上,对由衰减关系所确立的两种输入能量谱进行对比分析,讨论了延性系数对这两种输入能量谱的影响. 研究发现,在弹性情况下,两种输入能量谱在周期0.5~1.0 s范围内相差不大; 在非弹性情况下 两种输入能量谱在周期0.5 s处相差不大. 周期较小时,绝对输入能量谱要大于相对输入能量谱;周期较大时,绝对输入能量谱小于相对输入能量谱. 延性系数对这两种输入能量谱影响均比较大, 对绝对输入能量谱而言,周期小于0.3 s时,随着延性系数的增大, 能量谱升高;周期大于0.3 s时,随着延性系数的增大能量谱降低. 不同延性系数的绝对输入能量谱在分界点0.3 s左右相等. 相对输入能量谱受延性系数影响与绝对输入能量谱相似, 但分界点在0.5 s左右. 与绝对输入能量谱相比, 相对输入能量谱在短周期段受延性系数的影响较大, 特别是当场地较软时更为明显.      

银川台DSQ浮子水管仪数字与其模拟记录观测资料的对比分析

通过对银川台形变DSQ浮子水管仪数字化记录观测资料和模拟记录观测资料的数据各项参数进行对比分析,认为数字化记录资料的观测精度、信息量、异常拾取量均优于模拟观测;数字化记录观测资料和模拟记录观测资料的各项对比分析结果都有很好吻合,但数字化记录资料数据利用率及数据处理方法有待进一步提高。数字化记录资料人机结合的处理方式和结果是提高资料精度的关键。     

COMPARISON BETWEEN RESULTS OF MODEL TESTS AND OBSERVED DATA FOR THE FLUCTUATING BACKWATER REGION OF THE GEZHOUBA RESERVOIR

The fluctuationg backwater region of a reservoir has a twofold property of a river as well as a reservoir and undergoes scouring and deposition alternately. Phys- ical model simulating the fluvial processes in a fluctuating backwater region has to ensure the similarities of both scouring and deposition. This causes some difficulties in model design and the selection of model sediment. This paper describes the model study of the fluctuating backwater region of the Gezhouba Project on the Yangtze River. Plastic particles with specific gravity of 1 .056 are used. The physical. model covers a reach of mountainous river 130 km in length. The experimental results show that the total amount of deposit, its distribution along the reach, the size distribution of the deposit, and lengths of the backwater region for various discharges all agree very well with the measured data in reservoir after 5 years of operation. This fact indicates that the model design, including the selection of model sediment, is quite acceptable.     

格尔木台数字与模拟地磁资料的一致性分析

对格尔木台正在运行的地磁相对记录数字资料及模拟资料进行了对比分析,总结出数字记录资料要明显优于模拟资料,数字仪器取代模拟仪器是必然趋势。     

山东地区流动磁测与地震关系的探讨

利用山东地区1991—2007年的野外地磁流动监测资料,通过分区计算地磁场总强度均值,并绘制均值时序图,分析了地磁场异常,为地震预测提供参考依据。     

A COMPARISON OF STATISTICAL AND DETERMINISTIC WIENER DECONVOLUTION OF DEEP-TOW SEISMIC DATA*

Wiener ‘spiking’ deconvolution of seismic traces in the absence of a known source wavelet relies upon the use of digital filters, which are optimum in a least-squares error sense only if the wavelet to be deconvolved is minimum phase. In the marine environment in particular this condition is frequently violated, since bubble pulse oscillations result in source signatures which deviate significantly from minimum phase. The degree to which the deconvolution is impaired by such violation is generally difficult to assess, since without a measured source signature there is no optimally deconvolved trace with which the spiked trace may be compared. A recently developed near-bottom seismic profiler used in conjunction with a surface air gun source produces traces which contain the far-field source signature as the first arrival. Knowledge of this characteristic wavelet permits the design of two-sided Wiener spiking and shaping filters which can be used to accurately deconvolve the remainder of the trace. In this paper the performance of such optimum-lag filters is compared with that of the zero-lag (one-sided) operators which can be evaluated from the reflected arrival sequence alone by assuming a minimum phase source wavelet. Results indicate that the use of zero-lag operators on traces containing non-minimum phase wavelets introduces significant quantities of noise energy into the seismic record. Signal to noise ratios may however be preserved or even increased during deconvolution by the use of optimum-lag spiking or shaping filters. A debubbling technique involving matched filtering of the trace with the source wavelet followed by optimum-lag Wiener deconvolution did not give a higher quality result than can be obtained simply by the application of a suitably chosen Wiener shaping filter. However, cross correlation of an optimum-lag spike filtered trace with the known ‘actual output’ of the filter when presented with the source signature is found to enhance signal-to-noise ratio whilst maintaining improved resolution.     

地震层析成象技术在测井和地面反射资料集中的应用研究

本文以射线的折射和反射定律为基础,采用阻尼最小二乘法（在CT技术中,常称为SIRT算法）和单井单侧多偏移距的反射资料以及地面反射资料,研究反射波层析成象技术对二维弯曲界面的成层结构进行重构的情况。数字模拟试验表明,本文提出的逐层成象方法对透镜体和尖灭地层的成象颇为有效,反演结果对误差扰动是稳定的。 本文所述的成象方法及处理技术可推广到只用地面反射资料进行地层重构的情形。