平面波产生的侧波

一、引言 側波,在地震勘探中常称为折射波,在一般的文献中又常称为首波。“折射波”这名称易和普通的折射波相混,“首波”从其含义說又沒能包括和其本属同一范畴的表面側波。但是在理論研究中,把首波和表面侧波放在一起討論比較方便,故这里就采用“侧波”这名.     

声波测井时裸眼井中弹性波的传播

本文分析了简化模型下声波测井时裸眼井（内充流体一般是泥浆）中弹性波的传播,此波场可看成是轴向传播有不同波矢的分波的叠加。分析表明,对于井外只有单一均匀岩层且无限深的井,仅对如（18）式所示的特征频率成分才存在具有轴向传播最大相速度v_zm=c₂的所谓“首波”;而v_x > c₂及c₂ > v_z>b₂的分波却不存在,此处c₂和b₂分别是井外岩层的纵波和横波速度。这正解释了文献[3]中所报道的裸眼井的“过滤效应”。分析还表明,对于井外具有多层岩层的井,只在一定条件下才存在与上述特性相近的“首波”。本文结果对声速测井显然是重要的。     

地震多次波成像

地震资料含有各种类型多次波,而传统成像方法仅利用地震一次反射波成像,在地震成像前需将多次波去除.然而,多次波携带了丰富的地下结构信息,多次波偏移能够提供除反射波外的额外地下照明.修改传统逆时偏移方法,用包含一次反射波和多次波的原始记录代替震源子波,将SRME方法预测的表面多次波代替一次反射波作为输入数据,可将表面多次波成像.多次波成像的挑战和困难在于大量串扰噪声的产生,针对表面多次波成像中的成像噪声问题,将最小二乘逆时偏移方法与多次波分阶思想结合起来,发展可控阶数的表面多次波反演成像方法,有望初步实现高精度的表面多次波成像.在消除原始记录中的表面多次波后,通过逆散射级数方法预测得到层间多次波,将层间多次波作为逆时偏移方法的输入数据可将其准确归位到地下反射位置.数值实验表明,多次波成像能够有效地为地下提供额外照明,而可控阶表面多次波最小二乘逆时偏移成像方法几乎完全避免成像噪声.     

一阶多次波聚焦变换成像

将多次波转换成反射波并按传统反射波偏移算法成像,是多次波成像的一种方法.聚焦变换能准确的将多次波转换为纵向分辨率更高的新波场记录,其中一阶多次波转换为反射波.本文对聚焦变换提出了两点改进:1)提出局部聚焦变换,以减小存储量和计算量,增强该方法对检波点随炮点移动的采集数据的适应性;2)引入加权矩阵,理论上证明原始记录的炮点比检波点稀疏时,共检波点道集域的局部聚焦变换可以将多次波准确转换成炮点与检波点有相同采样频率的新波场记录.本文在第一个数值实验中对比了对包含反射波与多次波的原始记录做局部聚焦变换和直接对预测的多次波做局部聚焦变换两种方案,验证了第二种方案转换得到的波场记录信噪比更高且避免了第一个方案中切聚焦点这项比较繁杂的工作.第二个数值实验表明:在炮点采样较为稀疏时,该方法能有效的将一阶多次波转换成反射波;转换的反射波能提供更丰富的波场信息,成像结果更均衡、在局部有更高的信噪比,以及较高的纵向分辨率.     

基于波射线路径偏移压制多次波

波射线路径压制多次波的反射波成像是在偏移过程去除多次波同时仅对反射波成像.通过在共炮道集和共检波点道集分别计算炮点射线的入射角和检波点射线的出射角计算射线的路径.从炮点入射的射线与从检波点出射的射线的交点形成的走时,若等于观测走时,可以判断此条射线是反射波;反之,若不相等,则是多次波.数值实验表明此方法可以有效地去掉由于多次波能量产生的假成像点和压制多次波,因此界面可以正确归位,同时去掉由于多次波引起的假成像位置.     

逆子波域消除多次波方法研究

SRMA（与表面相关多次波的衰减）算法包含预测和相减两步.相减算法中,当多次波与反射波同相轴相交时,如何有效减去多次波、保留反射波,是面临的主要问题.通过分析非正交性对滤波器（逆子波）的影响,可以证明：逆子波因非正交性产生的误差呈近似的高斯分布.在此基础上,本文提出了在逆子波域（单道自适应相减滤波的滤波算子的集合）,利用其误差的概率分布特征,对逆子波进行估计,用逆子波的估计对逆子波进行校正来消除多次波的方法.其步骤为：首先用SRMA方法预测出表面多次波,并对每一单炮进行单道自适应相减,得到逆子波,形成逆子波域;其次,在逆子波域采用中值滤波,提取接近真实逆子波的逆子波估计;第三,在逆子波域用逆子波估计对畸变的逆子波进行校正;最后采用校正后的逆子波来衰减多次波.通过简单模型和SMARRT模型的测试,该方法不仅能够有效减去多次波,而且在相交的区域,能够保持反射波同相轴的连续性并恢复其正确的振幅.     

南海深水多次波成像

南海海域水深可达5 km,是我国未来油气勘探最有前景的区域.南海海域地震资料含有各种类型多次波,在地震成像前需将多次波去除.然而与传统的一次反射波偏移比较,多次波携带更多的地下结构信息,多次波偏移可提供更好的地下成像.本文提出了一种修改传统逆时偏移的方法,利用多次波进行成像.该方法用包含一次波和多次波的记录代替震源子波,用预测多次波代替一次反射波作为输入数据.在逆时偏移处理流程中,将地表多次波记录沿时间延拓到各个深度层,通过互相关成像条件将包含一次波和多次波的地震记录沿时间延拓到相同深度层.将多次波偏移结果与一次波偏移结果的差应用目标函数进行匹配,弥补一次波偏移的不足.对Sigsbee2B数据做数值测试的结果表明与传统一次反射波偏移相比,该方法可以得到范围更大、均衡性更好的地下照明度.将该方法应用到南海深水实际资料成像,获得了南海深水地震资料多次波成像剖面.偏移结果表明,多次波成像可以弥补反射波照明的不足,由此说明多次波偏移是未来的一个重要研究方向.     

亚洲中东部地区塑性流动波(“慢波”)与地震活动

根据先前相似模型实验的结果 ,在板块边界驱动下 ,沿岩石圈下层 (含下地壳和岩石圈地幔 )传播的塑性流动波包括“快波”和“慢波” ,二者均属于黏性重力波 ,并分别由主波和辅波叠加而成。“快波”以 10 0 ～ 10 2 km/a量级的速度传播 ,已为塑性流动波控制下地震迁移的研究所证实。文中根据 7级以上强震的条带状分布图像 ,证明在喜马拉雅弧驱动下还存在着波速仅为 10 0 ～ 10 1m/a量级的“慢波” ,其中对 7级以上强震起控制作用的主要是“慢波”的辅波 ,其平均波长为 4 45km ,波速为 0 81～ 2 80m/a ,周期为 0 .16～ 0 .5 5Ma。“慢波”的边界起波时间距今约 1.34～ 4 .6 6Ma ,相当于上新世中期至早更新世中期 ,与喜马拉雅构造运动的主要活跃期 (幕 )之一相吻合。以喜马拉雅弧西段和东段为波源所形成的2个“慢波”系统的波峰带相互重叠 ,为 7级以上强震的发生提供了必要的能量背景     

基于稀疏反演的OBS数据多次波压制方法

自由表面多次波压制是海底地震仪（Ocean Bottom Seismometer,OBS）数据处理和成像中的难点,OBS数据多次波能量强,周期长,严重影响深层一次反射波的处理和成像.不同于常规拖缆观测系统,OBS数据站点一般相隔较远,仅仅利用检波点稀疏的波场信息难以压制OBS数据中的自由表面多次波.本文采用拖缆数据与OBS数据联合,利用稀疏反演估计（Estimation of Primaries and Multiples by Sparse Inversion,EPSI）方法,研究了OBS数据自由表面多次波压制理论,分析了OBS多次波产生的机理,详细推导了拖缆数据与OBS数据联合预测OBS多次波的EPSI方法基本原理.通过利用拖缆数据的信息,实现了OBS检波点稀疏数据多次波的压制问题.EPSI方法通过稀疏反演直接估计一次反射波,避免了SRME（Surface Related Multiple Elimination）方法中自适应相减对有效信号的损害,保真了一次反射有效信号,理论模拟OBS数据验证了方法的有效性.     

表面多次波最小二乘逆时偏移成像

使用相同的炮记录,多次波偏移能提供比反射波偏移更广的地下照明和更多的地下覆盖但是同时产生很多的串声噪声.相比传统逆时偏移,最小二乘逆时偏移反演的反射波成像结果具有更高的分辨率和更均衡的振幅.我们主要利用最小二乘逆时偏移压制多次波偏移产生的串声噪声.多次波最小二乘逆时偏移通常需要一定的迭代次数以较好地消除串声噪声.若提前将一阶多次波从所有阶数的多次波中过滤出来,使用相同的迭代次数,一阶多次波的最小二乘逆时偏移能够得到具有更高信噪比的成像剖面,而且能够提供与多次波最小二乘逆时偏移相似的有效地下结构成像.     

柴达木盆地低頻地震探測的基岩首波和大角度反射波

1958年6月开始,在柴达木盆地进行低頻折射地震的試驗工作,目的在于探測埋藏很深的基岩,以及基岩以上的主要沉积岩界面。根据地貭資料的推測,盆地中部基岩埋藏的深度可达10公里以上,利用其他的地震勘探方法是不容易得到效果的.所用检波器的固有頻率为13周/秒,最远两个相遇炮点相距約90公里,接收距离达100公里.关于仪器的性能与工作方法,已在前文叙述,对于低頻地震波性貭的研究,仅有簡短的报导.     

柴达木东盆地的基岩首波和反射波

一、引言 目前,研究结晶基底,多采用地震首波的方法,利用反射波只能探测地震地质条件较好的浅层基岩面。由于柴达木盆地区内基岩面埋藏很深,界面性质很复杂,过去曾利用重力、磁力、电测方法和中频地震探测结晶基底,均未能取得满意的结果。也曾对西盆     

RESEARCH ON TIME-DISTANCE INTERPRETATION OF MULTIPLES AND GHOSTS BASED ON MARINE SEISMIC DATA

Currently, the study on the active fault in the land areas is relatively mature, while there is still lack of detection and research on active faults in the sea areas. Marine exploration, which is different from land areas, has a prominent problem due to the existence of strong reflecting interfaces such as water surface and seafloor in the sea, thus the recording is often accompanied by interference of multiples on seafloor reflections. In addition, because of the characteristics of marine seismic exploration, the source exciting in the water and the geophone receiving in the water, ghost wave usually can be recorded simultaneously during the reflected wave propagation. This phenomenon makes it difficult to distinguish the effective waves and the noise, and has always plagued the data and seriously affects the quality of records. In the offshore and other regions of complex structures, such as inclined interfaces, it is difficult to eliminate the interference of multiples accurately by traditional multiples suppression methods, which are based on the horizontal interface assumption. This paper combines the sea area seismic data and its acquisition method, uses simplified model to simulate the multiples based on the time-distance analysis of multiples and their ghost wave in inclined interface. The time-distance characteristics of the multiples and their ghost waves from different interfaces(including the inclined interface)are obtained, and they are consistent with the actual records. The multiples time-distance simulation can help to distinguish the causes of reflected waves, summarize the multiple-wave time-distance characteristics from different interfaces(including inclined interfaces), and analyze the relationship between the characteristics of multiple waves and primary waves. In particular, this simulation has a significant effect on characterizing the internal multiples that are difficult to identify due to inconspicuous periodicity and the multiples of the inclined interface which present the phenomenon that the vertex of the time-distance curve is shifted. On this basis, relying on the time-distance analysis of ghost wave, we analyze the travel time difference characteristics between reflected waves and their accompanying ghost waves. The differences of the travel time characteristics of different orders ghost wave and reflected wave are summarized and the symmetry of the travel time difference between inclined interface and horizontal interface of ghost waves and reflected waves is analyzed. We simulate the distraction of the ghost wave event with the event of the reflected wave and analyze the influence of the ghost wave on the sea area seismic records. These results can improve the practical interpretation of seismic data. At last, the time-distance information is used to synthesize sea area seismic records, which can help us carry out the effective data processing and understand the characteristics of the time-distance and velocity of multiples in different interfaces and the layer artifact caused by multiples. This study combines the time-distance simulation of multiples and their ghost wave with conventional seismic data processing to analyze the pre-stack and post-stack features of multiple waves and their ghost waves in the seismic records of the sea area. The results of this study are conducive to the effective identification of multiples in seismic records in the sea, provide a theoretical basis for multi-wave suppression and prediction, and may facilitate the future study of sea-area seismic activity detection.     

A METHOD TO DETERMINE THE VELOCITIES OF THE SEAFLOOR AND NEAR-SURFACE SEDIMENTS*

Shotpoint gathers from conventional reflection seismic surveys contain both reflected and refracted waves. In this study shot records were processed and analyzed, and the data were modeled with reflected, refracted, and reflected-refracted waves to fit the recorded data. The result is a detailed velocity model. The inverse problem for refracted waves was solved by using the Wiechert-Herglotz inversion. A 500-km-long 26-fold reflection seismic line from the Barents Sea, north of Norway, has been investigated. The data show high velocities, multiple reflections, and various types of noise. To test the method a total of 34 shot gathers were analyzed along this line. The aim of the interpretation was to determine the velocity in the seafloor and the near-surface sediments. It is possible to map the vertical as well as the lateral velocity distribution in detail. Depending on the length of the streamer and the velocity gradient in the sediments, the calculated depth varies between 300 and 500 m below the seafloor. These velocities were also compared to the stacking velocities obtained from the reflection seismic data to see how the velocities determined by different methods were related. The velocity distribution in the sediments is one of the key factors in seismic interpretation. The technique discussed in this paper can contribute to velocity information both in the processing and interpretation of seismic data.     

柴达木东盆地的深层地震反射波和地壳构造

一、引言 地壳和上地幔顶部的构造,对研究地震发生、发展及演变过程的深部背景,以及划分地震活动块体和探讨地震成因有着重要的意义。然而地壳构造的研究,又是与地震学的发展密切相关。1909年,莫霍洛维奇在近震研究中,首先发现了地壳与上地幔分界面的首波,其覆盖层的平均速度 =6.3公里/秒,界面速度V_d=8.0公里/秒（简称M界面）。后于1923年,康拉德（conrad）也是根据天然地震资料鉴别出一个平均速度为=5.4     

REFLECTIONS AND REFRACTIONS FROM CURVED INTERFACES: MODEL STUDY *

Reflections and refractions from curved interfaces were studied on two dimensional scale models. Time of arrival, amplitude and character of reflected, converted, and refracted waves were mainly used for this study. Some reflected refractions, refracted reflections and diffractions were also considered. It was possible to separate PS and SP waves and to study their amplitude and character separately. From the amplitude study of refracted arrivals it was concluded that the refracted ray path penetrates into the high velocity layer rather than propagating along the interface. Although most of the results are interesting from the theoretical point of view, a few applications to exploration problems are suggested.     

SEISMIC MODELLING WITH CHANNEL WAVES IN SEAM STRUCTURES INFLUENCED BY MYLONITE ZONES1

Channel waves generated in coal-seams and their reflections from discontinuities are widely used to indicate the tectonic and stratigraphic features of coal deposits, resulting in greater efficiency and safety in coal-mining. In the mining area of Ibbenbüren (F.R.G.) seam structures sometimes contain so-called mylonite zones, which are crushed coal deposits capable of binding gas. If mining hits a mylonite zone this would probably not only reduce output of the mine, but could even cause gas explosions. To investigate the influence of a mylonite zone on the propagation of channel waves, Rayleigh channel wave measurements for 2D analogue models were performed and synthetic seismograms of Love channel waves were calculated. Analogue modelling of the mylonite zone using Rayleigh seam waves within the seam was carried out using a perforation technique. Calculations were made to obtain an estimate of velocity reduction due to perforation. The results agree well with velocity values measured up to a perforation of 25% in a 2D epoxy resin model. Reflected channel wave energy was found by applying dispersion analysis in the case where the impedance reduction between the mylonite seam structure and the undisturbed seam was approximately 5%. The horizontal width of the mylonite structure was detectable from the time lag between reflected channel wave signals from both in-seam borders of the mylonite zone. Resolution of two discrete borders was possible for a width of 1.5 λ's. The influence of a vertical fault, positioned within the mylonite zone, could only poorly be resolved. Numerical model investigations of Love seam waves were concerned mainly with the variation of the horizontal width of the mylonite zone. Mylonite zones with dimensions of the order of 0.4 λ's allow definite statements about their widths from dispersion and spectral analyses. For zones with smaller widths down to 0.2 λ's, it was found that reflectivity and transmissivity analyses give a qualitative criterion for distinguishing a mylonite structure surrounding a fault from a pure fault.     

COMPUTING FIELD STATICS WITH THE HELP OF SEISMIC TOMOGRAPHY*

Field static corrections in general need be applied to all onshore seismic reflection data to eliminate the disturbing effects a weathering layer or near-surface low velocity zone has on the continuity of deep seismic reflections. The traveltimes of waves refracted at the bottom of the low velocity zone (or intermediate refracting interfaces) can often be observed as first breaks on shot records and used to develop a laterally inhomogeneous velocity model for this layer, from which the field static corrections can then be obtained. A simple method is described for computing accurate field statics from first breaks. It is based on a linearization principal for traveltimes and leads to the algorithms that are widely and successfully applied within the framework of seismic tomography. We refine an initial model for the low velocity layer (estimated by a standard traveltime inversion technique) by minimizing the errors between the observed first arrivals on field records and those computed by ray theory through an initial model of the low velocity layer. Thus, one can include more lateral velocity variations within the low velocity layers, which are important to obtain good field static corrections. Traditional first break traveltime inversion methods cannot, in general, provide such refined velocity values. The technique is successfully applied to seismic data from the Amazon Basin. It is based on a simple model for the low velocity layer that consists of an undulating earth surface and one planar horizontal refractor overlain by a laterally changing velocity field.     

MULTIPLE REFLECTIONS AND HEAD WAVES IN THE GULF OF SUEZ*

Inferior reflection quality in the Gulf of Suez at the target depth interval is attributable in part to surficial multiple reflections. An excellent example of the latter is observed on a typical seismic line in the northern portion of the Gulf. An increase in prominence of the multiple reflections appears associated with decreasing depth to a dipping highvelocity layer. Inversion of a second-order polynomial time-distance function, fitted to the observed refraction onset time-distance values, gives the velocity-depth function for sediments between the water bottom and a high-velocity layer. Velocities thus determined increase non-linearly with depth from a value near water velocity at the water bottom. Depths to the high-velocity layer are obtained from the associated head-wave linear time-distance function and by ray tracing in the overlying sediments. As the high-velocity layer approaches the water bottom from sub-water depths exceeding 0.6 km to a depth of 56 m, intensity of the multiple reflections increases to the extent of completely dominating individual records to a time of at least 3 s. The estimated plane-wave normal-incident reflection coefficient at the top of the high-velocity layer increases with decreasing depth to this layer, approaching 0.5 at the shallowest depth. This strong reflection coefficient further substantiates the existence of multiple reflections between the high-velocity layer and water layer. However, existence of water-layer multiples cannot be ruled out. The estimated water-bottom reflection coefficient is approximately 0.3, a substantial value. Multiple reflections of considerably less intensity are apparent where the high-velocity layer is deepest, and it is likely that such are waterlayer multiple reflections. Unfortunately, water-layer multiple reflections and multiple reflections between the water surface and high-velocity layer cannot be separated by their coincidence with time-distance (normal moveout) curves, the configuration of each visibly matching the curves equally well.     

我国西北地区地壳中的高速夹层

在我国西北地区的柴达木盆地东部和甘肃地区,在距离炮点40互100公里处,能够接收到不少能量较强的地壳深界面反射波。另外还发现一种与一般反射波性质不同的波,其视速度特大,视速度随距离的变化不大,而且有较明显的终点;其吋距曲线与一般深界面反射波的时距曲线相交。根据它的特征可以判断地壳中存在具有速度梯度的高速夹层.求得的夹层参数为: 甘肃地区柴达木盆地东部覆盖层厚度 18.8公里 30.5公里覆盖层平均速度 5.5公里/秒 5.3公里/秒夹层厚度 6.0公里 3.2公里夹层速度 7.5-8.5公里/秒 7.5-8.0公里/秒夹层的上下界面均为强反射面,可以产生多次反射波。分別利用相邻两个反射波可以求得各层参数,并能避免射线折射的影响。甘肃地区和柴达木盆地东部的地壳厚度分別为51和52公里。地壳中有高速夹层的存在,可以更好地说明P~*速度分散的原因,而且也能够解释Lg波的传播机制。