共查询到20条相似文献，搜索用时 31 毫秒
1.
在城市活断层探测中 ,浅层结构常常表现为强烈的非均匀性 ,界面横向强烈起伏 ,层内速度变化较大 ,传统的基于平界面均匀层模型的折射资料处理方法不能适用。研究开发能应用于复杂介质结构中折射资料处理的方法就显得十分必要。文中基于惠更斯原理 ,用波前扩张法对波场作正演计算 ,根据哈格多恩折射波前成像原理 ,在lecomte算法和Hole有限差分计算程序的基础上 ,开发出 1种复杂介质结构中折射资料的处理方法与软件 ,并用此方法处理了福州城市活断层折射探测试验中在义序完成的 2条折射剖面资料。结果表明 :探测区浅层为 3层结构 ,分别为盖层、强风化层和基岩。基岩顶界面的埋深约为 5 8～ 5 2m ,盖层P波速度变化较大 相似文献
2.
速度是地震偏移成像准确与否的关键所在.全波形反演综合利用地震波场运动学和动力学信息,能够得到相比传统速度建模方法更高频的成分.全波形反演的理论比较成熟,但实际应用成功的例子相对较少,特别是对于陆上地震资料.塔里木盆地地震地质条件复杂,为了实现缝洞型储层的准确成像,本文开展了针对目标靶区的全波形反演精细速度建场研究.采用一种时间域分层多尺度全波形反演流程:首先通过层析成像建立初始速度模型;其次利用折射波反演浅层速度模型;最后利用反射波反演中深层速度模型.偏移成像结果表明基于全波形反演的速度建模技术能有效改善火成岩下伏构造的成像精度,显示了全波形反演在常规陆上采集资料的应用潜力. 相似文献
3.
Dennis Woodward 《Journal of Applied Geophysics》1994,31(1-4)
The US Geological Survey, in cooperation with the National Drilling Company of Abu Dhabi, is conducting a 4-year study of the fresh and slightly saline groundwater resources of the eastern Abu Dhabi Emirate. Most of this water occurs in a shallow aquifer, generally less than 150 m deep, in the Al Ain area. A critical part of the Al Ain area coincides with a former petroleum concession area where about 2780 km of vibroseis data were collected along 94 seismic lines during 1981–1983. Field methods, acquistion parameters, and section processing were originally designed to enhance reflections expected at depths ranging from 5000 to 6000 m, and subsurface features directly associated with the shallow aquifer system were deleted from the original seismic sections. The original field tapes from the vibroseis survey were reprocessed in an attempt to extract shallow subsurface information (depths less than 550 m) for investigating the shallow aquifer.A unique sequence of reproccessing parameters was established after reviewing the results from many experimental tests. Many enhancements to the resolution of shallow seismic reflections resulted from: (1) application of a 20-Hz, low-cut filter; (2) recomputation of static corrections to a datum nearer the land surface; (3) intensive velocity analyses; and (4) near-trace muting analyses. The number, resolution, and lateral continuity of shallow reflections were greatly enhanced on the reprocessed sections, as was the delineation of shallow, major faults. Reflections on a synthetic seismogram, created from a borehole drilled to a depth of 786 m on seismic line IQS-11, matcheddprecisely with shallow reflections on the reprocessed section. The 33 reprocessed sections were instrumental in preparing a map showing the major structural features that affect the shallow aquifer system. Analysis of the map provides a better understanding of the effect of these shallow features on the regional occurrence, movement, and quality of groundwater in the concession area. Results from this study demonstrate that original seismic field tapes collected for deep petroleum exploration can be reprocessed to explore for groundwater. 相似文献
4.
利用广州地区详尽的深、浅部地震勘探资料,大深度钻孔和测速等资料,运用自主开发的建模程序和多功能商业软件的二次开发功能高速处理海量数据进行建模,采用开发的二元三次B—Spline函数建模方法解决大城市圈建模数据空间分布不均匀的问题并进行模型优化,建立了广州地区符合地质学和地形学理论三维不均匀地下速度结构模型,为基于活断层探测和地震动危险性评价基础上的强地震动预测（模拟）奠定了基础。 相似文献
5.
B Damotte 《Journal of Applied Geophysics》1993,29(3-4)
The French Ecors program was launched in 1983 by a cooperation agreement between universities and petroleum companies. Crustal surveys have tried to find explanations for the formation of geological features, such as rifts, mountains ranges or subsidence in sedimentary basins. Several seismic surveys were carried out, some across areas with complex geological structures. The seismic techniques and equipment used were those developed by petroleum geophysicists, adapted to the depth aimed at (30–50 km) and to various physical constraints encountered in the field. In France, Ecors has recorded 850 km of deep seismic lines onshore across plains and mountains, on various kinds of geological formations. Different variations of the seismic method (reflection, refraction, long-offset seismic) were used, often simultaneously. Multiple coverage profiling constitutes the essential part of this data acquisition. Vibrators and dynamite shots were employed with a spread generally 15 km long, but sometimes 100 km long.Some typical seismic examples show that obtaining crustal reflections essentialy depends on two factors: (1) the type and structure of shallow formations, and (2) the sources used. Thus, when seismic energy is strongly absorbed across the first kilometers in shallow formations, or when these formations are highly structured, standard multiple-coverage profiling is not able to provide results beyond a few seconds. In this case, it is recommended to simultaneously carry out long-offset seismic in low multiple coverage.Other more methodological examples show: how the impact on the crust of a surface fault may be evaluated according to the seismic method implemented (vibroseis 96-fold coverage or single dynamite shot); that vibrators make it possible to implement wide-angle seismic surveying with an offset 80 km long; how to implement the seismic reflection method on complex formations in high mountains.All data were processed using industrial seismic software, which was not always appropriate for records at least 20 s long. Therefore, a specific procedure adapted to deep seismic surveys was developed for several processing steps. The long duration of the vibroseis sweeps often makes it impossible to perform correlation and stack in the recording truck in the field. Such field records were first preprocessed, in order to be later correlated and stacked in the processing center. Because of the long duration of the recordings and the great length of the spread, several types of final sections were replayed, such as: (1) detailed surface sections (0–5 s), (2) entire sections (0–20 s) after data compression, (3) near-trace sections and far-trace sections, which often yield complementary information.Standard methods of reflection migration gave unsatisfactory results. Velocities in depth are inaccurate, the many diffractions do not all come from the vertical plane of the line, and the migration software is poorly adapted to deep crustal reflections. Therefore, migration is often performed graphically from arrivals picked in the time section. Some line-drawings of various onshore lines, especially those across the Alps and the Pyrenees, enable to judge the results obtained by Ecors. 相似文献
6.
Derecke Palmer 《Geophysical Prospecting》2010,58(4):561-575
The tau‐p inversion algorithm is widely employed to generate starting models with many computer programs that implement refraction tomography. However, this algorithm can frequently fail to detect even major lateral variations in seismic velocities, such as a 50 m wide shear zone, which is the subject of this study. By contrast, the shear zone is successfully defined with the inversion algorithms of the generalized reciprocal method. The shear zone is confirmed with a 2D analysis of the head wave amplitudes, a spectral analysis of the refraction convolution section and with numerous closely spaced orthogonal seismic profiles recorded for a later 3D refraction investigation. Further improvements in resolution, which facilitate the recognition of additional zones with moderate reductions in seismic velocity, are achieved with a novel application of the Hilbert transform to the refractor velocity analysis algorithm. However, the improved resolution also requires the use of a lower average vertical seismic velocity, which accommodates a velocity reversal in the weathering. The lower seismic velocity is derived with the generalized reciprocal method, whereas most refraction tomography programs assume vertical velocity gradients as the default. Although all of the tomograms are consistent with the traveltime data, the resolution of each tomogram is comparable only with that of the starting model. Therefore, it is essential to employ inversion algorithms that can generate detailed starting models, where detailed lateral resolution is the objective. Non‐uniqueness can often be readily resolved with head wave amplitudes, attribute processing of the refraction convolution section and additional seismic traverses, prior to the acquisition of any borehole data. It is concluded that, unless specific measures are taken to address non‐uniqueness, the production of a single refraction tomogram that fits the traveltime data to sufficient accuracy does not necessarily demonstrate that the result is either correct, or even the most probable. 相似文献
7.
《Journal of Applied Geophysics》2001,46(1):31-44
During two Antarctic summers (1996–1997 and 1997–1998), five seismic refraction and two reflection profiles were acquired on the Johnsons Glacier (Livingston Island, Antarctica) in order to obtain information about the structure of the ice, characteristics of the ice-bed contact and basement topography. An innovative technique has been used for the acquisition of reflection data to optimise the field survey schedule. Different shallow seismic sources were used during each field season: Seismic Impulse Source System (SISSY) for the first field survey and low-energy explosives (pyrotechnic noisemakers) during the second one. A comparison between these two shallow seismic sources has been performed, showing that the use of the explosives is a better seismic source in this ice environment. This is one of the first studies where this type of source has been used. The analysis of seismic data corresponding to one of the reflection profiles (L3) allows us to delineate sectors with different glacier structure (accumulation and ablation zones) without using glaciological data. Moreover, vertical discontinuities were detected by the presence of back-scattered energy and the abrupt change in frequency content of first arrivals shown in shot records. After the raw data analysis, standard processing led us to a clear seismic image of the underlying bed topography, which can be correlated with the ice flow velocity anomalies. The information obtained from seismic data on the internal structure of the glacier, location of fracture zones and the topography of the ice-bed interface constrains the glacial dynamics of Johnsons Glacier. 相似文献
8.
Fiona M. Campbell A. Kaiser H. Horstmeyer A.G. Green F. Ghisetti A.R. Gorman M. Finnemore D.C. Nobes 《Journal of Applied Geophysics》2010,70(4):332-342
In an attempt to understand the structure of active faults as they emerge from bedrock into shallow semi-consolidated and unconsolidated sediments, we have recorded a comprehensive high-resolution seismic reflection/refraction data set across the Ostler Fault zone on the central South Island of New Zealand. This fault zone, which absorbs 1–2 mm/yr of compression associated with oblique convergence of the Pacific and Australian tectonic plates, consists of a series of surface-rupturing N–S trending, west-dipping reverse faults that offset a thick sequence of Quaternary glacial outwash and late Neogene fluvio-lacustrine sediments of the Mackenzie Basin. Our study focuses on a region of the basin where two non-overlapping fault segments are separated by a transfer zone. Deformation in this area is accommodated by offsets on multiple small faults and by folding in their hanging walls. The seismic data with source and receiver spacing of 6 and 3 m and nominal CMP fold of 60 was acquired along twelve 1.2 km long lines orthogonal to fault strike and an additional 1.6 km long tie-line parallel to fault strike. The combination of active deformation and shallow glacial outwash sediments results in particularly complicated seismic data, such that application of relatively standard processing schemes yields only poor quality images. We have designed a pre- and post-stack reflection/refraction processing scheme that focuses on minimising random and source-generated noise, determining appropriate static corrections and resolving contrasting reflection dips. Application of this processing scheme to the Ostler Fault data provides critical information on fault geometry and offset and on sedimentary structures from the surface to ～ 800 m depth. Our preliminary interpretation of one of the lines includes complex deformation structures with folding and multiple subsidiary fault splays on either side of a ～ 50° west-dipping primary fault plane. 相似文献
9.
针对近地表物质非均质极强、各向异性明显及地形复杂等特点,系统阐述和讨论了近地表折射和反射法的国内外研究与应用进展,认为:1综合利用纵、横波的优势,开展多波多分量联合勘探对提高浅层地震勘探的精度和分辨率具有重要作用;2现在的浅层地震勘探主要是对地震剖面进行解释,容易忽略一些隐含的地质异常现象,属性提取技术是充分提取地震信息,进行全面综合解释的有效手段;3开展多层折射介质的观测系统和解释方法研究,尤其是折射层析成像研究,是提高多层折射介质成像精度的途径;4开展黏弹性、双相和各向异性介质的地震反射与折射波法研究是提高近地表地震勘探成像和物性参数提取精度的新思路. 相似文献
10.
E. BRÜCKL 《Geophysical Prospecting》1987,35(9):973-992
We consider multiply covered traveltimes of first or later arrivals which are gathered along a refraction seismic profile. The two-dimensional distribution of these traveltimes above a coordinate frame generated by the shotpoint axis and the geophone axis or by the common midpoint axis and the offset axis is named a traveltime field. The application of the principle of reciprocity to the traveltime field implies that for each traveltime value with a negative offset there is a corresponding equal value with positive offset. In appendix A procedures are demonstrated which minimize the observational errors of traveltimes inherent in particular traveltime branches or complete common shotpoint sections. The application of the principle of parallelism to an area of the traveltime field associated with a particular refractor can be formulated as a partial differential equation corresponding to the type of the vibrating string. The solution of this equation signifies that the two-dimensional distribution of these traveltimes may be generated by the sum of two one-dimensional functions which depend on the shotpoint coordinate and the geophone coordinate. Physically, these two functions may be interpreted as the mean traveltime branches of the reverse and the normal shot. In appendix B procedures are described which compute these two functions from real traveltime observations by a least-squares fit. The application of these regressed traveltime field data to known time-to-depth conversion methods is straightforward and more accurate and flexible than the use of individual traveltime branches. The wavefront method, the plus-minus method, the generalized reciprocal method and a ray tracing method are considered in detail. A field example demonstrates the adjustment of regressed traveltime fields to observed traveltime data. A time-to-depth conversion is also demonstrated applying a ray tracing method. 相似文献
11.
海州—韩山断裂是连云港地区一条重要断裂, 属于海泗断裂带的西边界断裂。 断裂隐伏于较浅的覆盖层之下且控制了基岩岩性分界。 浅层地震反射法作为断层探测的首选方法, 对海州—韩山断裂进行探测时仅能识别出基岩顶面反射波, 难以实现对断裂的准确判别。 而折射层析成像法适用于速度横向差异大的区域, 可获取断层两盘岩性速度差异信息, 进而判定断层位置, 弥补反射法的不足。 本文在跨海州—韩山断裂同一位置上联合应用浅层地震反射和折射层析成像两种探测方法, 进行了钻孔联合剖面探测和合成地震记录验证。 研究表明, 浅层地震反射和折射层析成像两种方法联合探测海州—韩山断裂, 较单一方法可对断层实现更精准的定位, 获取更丰富的断层信息, 为类似地质条件下的断层探测提供了思路。 相似文献
12.
A combined reflection/refraction (wide-angle) seismic survey was conducted on the continental shelf north-west of Britain, using a conventional streamer with an airgun source, and static ocean-bottom seismometers (OBS) to record wide-angle energy. The shallow structure down to a basaltic layer was reasonably well imaged on the stacked reflection section. The basalts, however, proved to be opaque to the conventional reflection method and prevented the imaging of deeper horizons, where an important velocity inversion was anticipated. This paper reports on the processing, modelling and interpretation of the densely sampled wide-angle OBS data that were coincident with the reflection profile. Eleven OBS instruments were deployed along a 75 km line and recorded signal from a powerful 149 litre (9100 in.3) airgun array fired every 50 m. Data processing was performed using a standard industrial reflection seismic software package prior to first-arrival picking. Processing steps included geometry definition, trace summation and display of the data using various scaling algorithms. An initial model was constructed from 1D velocity-time profiles digitized every 4 km along the stacked section. First arrival traveltime modelling rapidly converged to a detailed model of the structure of the top 5 km of the crust. Modelling revealed the existence of a buried low-velocity Mesozoic sedimentary basin, of a prominent basement horst and of a normal fault penetrating to the basement. 相似文献
13.
为提高过障碍变观设计的效率,利用C++编程语言开发出一个浅层地震过障碍变观设计软件,以提高过障碍变观的方便性和效率。当勘探中遇到地表障碍时可以利用此软件灵活变观,使地震测线跨越障碍物(江河、高速公路等),以保证反射同相轴能连续追踪对比。经实际应用,该软件效果较好,能够辅助数据采集人员在采集现场快捷地进行地震勘探过障碍变观设计,从而保证地震资料的质量。 相似文献
14.
DERECKE PALMER 《Geophysical Prospecting》1991,39(8):1031-1060
The depth to the surface of a refractor and the seismic velocity within the refractor are very often intimately related. In the shallow environment, increased thicknesses of weathering occur in areas of jointing, shearing or lithological variations, and these zones of deeper weathering can have lower subweathering refractor velocities. This association is important in geotechnical investigations and in the measurement of weathering thicknesses and sub-weathering velocities for statics corrections for reflection seismic surveys. Algorithms, which employ forward and reverse traveltime data and which explicitly accommodate the offset distance through the process known as refraction migration, are necessary if detailed structure on a refractor and rapid lateral variations of the seismic velocity within it are to be resolved. These requirements are satisfied with wavefront construction techniques, Hales’ method and the generalized reciprocal method (GRM). However, these methods employ refraction migration in fundamentally different manners. Most methods compute an offset distance with an often imprecise knowledge of the seismic velocities of the overlying layers. In contrast, the GRM uses a range of offset distances from less than to greater than the optimum value, with the optimum value being selected with a minimum-variance criterion. The approach of the GRM is essential where there are undetected layers and where there are rapid variations in the depth to a refractor and the seismic velocity within it. In the latter situations the offset distance necessary to define the seismic velocities can differ considerably from the value required to define depths. The efficacy of the GRM in resolving structure and seismic velocity is demonstrated with three model studies and two field examples. 相似文献
15.
《Journal of Applied Geophysics》2010,72(4):149-156
Characterization of shallow structures was performed by using different approaches analysing both P- and S-wave seismic data with different resolution. The refraction tomography provided P and S velocity models of the first 80 m, while the reflection seismic processing gives a reasonable stacking velocity field until 300 m depth for both P- and S-wave data. So, we estimated the Vp/Vs ratio and an empirical relationship between the two velocities. We characterised the shallow layers using tomographic velocity models and the deeper layers using seismic images with different resolution. The seismic images were obtained by conventional CMP reflection seismic processing and by a novel multi-refractor imaging technique. 相似文献
16.
Bengt Sjgren 《Geophysical Prospecting》2000,48(5):815-834
An analysis of the generalized reciprocal method (GRM), developed by Palmer for the interpretation of seismic refraction investigations, has been carried out. The aim of the present study is to evaluate the usefulness of the method for geotechnical investigations in connection with engineering projects. Practical application of the GRM is the main object of this study rather than the theoretical/mathematical aspects of the method. The studies are partly based on the models and field examples presented by Palmer. For comparison, some other refraction interpretation methods and techniques have been employed, namely the ABC method, the ABEM correction method, the mean‐minus‐T method and Hales' method. The comparisons showed that the results, i.e. the depths and velocities determined by Palmer, are partly incorrect due to some errors and misinterpretations when analysing the data from field examples. Due to the limitations of the GRM, some of which are mentioned here, stated by Palmer in his various publications, and other shortcomings of the method (e.g. the erasing of valuable information), the GRM must be regarded as being of limited use for detailed and accurate interpretations of refraction seismics for engineering purposes. 相似文献
17.
18.
在城市活断层勘探中,对于面波干扰大、地震反射波法难以开展工作的区域,可尝试利用地震折射波法进行探测,并对折射波法探测的原始记录采用时间项、差异时距曲线和有限差分成像等方法进行综合计算、分析,以探索折射波法在城市活断层勘探中的应用成效。文中针对郑州市须水断层西段浅层地震折射波法勘探记录,利用时间项、差异时距曲线和有限差分成像等计算方法,获取剖面速度结构与界面构造;综合震相特征、计算结果等资料确定主要地层的界面深度和构造特征,3种方法都取得了相近的结论。后又通过在测线上4个钻孔资料的验证,认为3种方法的计算结果与钻孔资料相吻合,说明折射波法勘探在城市活断层探测中的应用是可行的 相似文献
19.
Fang Shengming 《中国地震研究》2007,21(3):236-244
This paper introduces briefly the basic principles of various seismic prospecting techniques and working methods according to nationwide practices of seismic prospecting of active faults beneath big cities in recent years.Furthermore,it analyzes the application range of different seismic prospecting methods,main achievements and solved problems,and discusses the best combination of seismic exploration methods for detecting crustal structures and locating the faults used in the present stage,that is,to trace faults which are at depths of hundred of meters underground using shallow seismic investigation,to detect the faults which are above basement(at a depth of kilometers) using high resolution refraction sounding,and the deep crustal faults using combined seismic prospecting methods of reflection seismic sounding and wide-angle reflection/refraction sounding,and furthermore,to use the 3-D deep seismic sounding method to obtain 3-D velocity structures beneath urban areas.Thus,we can get information about fault attitude and distribution at different depths and a complete image of faults from their shallow part to deep part using the combined seismic exploration method.Some application examples are presented in the article. 相似文献
20.
A new direction for shallow refraction seismology: integrating amplitudes and traveltimes with the refraction convolution section 总被引：1，自引：0，他引：1
Derecke Palmer 《Geophysical Prospecting》2001,49(6):657-673
The refraction convolution section (RCS) is a new method for imaging shallow seismic refraction data. It is a simple and efficient approach to full‐trace processing which generates a time cross‐section similar to the familiar reflection cross‐section. The RCS advances the interpretation of shallow seismic refraction data through the inclusion of time structure and amplitudes within a single presentation. The RCS is generated by the convolution of forward and reverse shot records. The convolution operation effectively adds the first‐arrival traveltimes of each pair of forward and reverse traces and produces a measure of the depth to the refracting interface in units of time which is equivalent to the time‐depth function of the generalized reciprocal method (GRM). Convolution also multiplies the amplitudes of first‐arrival signals. To a good approximation, this operation compensates for the large effects of geometrical spreading, with the result that the convolved amplitude is essentially proportional to the square of the head coefficient. The signal‐to‐noise (S/N) ratios of the RCS show much less variation than those on the original shot records. The head coefficient is approximately proportional to the ratio of the specific acoustic impedances in the upper layer and in the refractor. The convolved amplitudes or the equivalent shot amplitude products can be useful in resolving ambiguities in the determination of wave speeds. The RCS can also include a separation between each pair of forward and reverse traces in order to accommodate the offset distance in a manner similar to the XY spacing of the GRM. The use of finite XY values improves the resolution of lateral variations in both amplitudes and time‐depths. The use of amplitudes with 3D data effectively improves the spatial resolution of wave speeds by almost an order of magnitude. Amplitudes provide a measure of refractor wave speeds at each detector, whereas the analysis of traveltimes provides a measure over several detectors, commonly a minimum of six. The ratio of amplitudes obtained with different shot azimuths provides a detailed qualitative measure of azimuthal anisotropy and, in turn, of rock fabric. The RCS facilitates the stacking of refraction data in a manner similar to the common‐midpoint methods of reflection seismology. It can significantly improve S/N ratios.Most of the data processing with the RCS, as with the GRM, is carried out in the time domain, rather than in the depth domain. This is a significant advantage because the realities of undetected layers, incomplete sampling of the detected layers and inappropriate sampling in the horizontal rather than the vertical direction result in traveltime data that are neither a complete, an accurate nor a representative portrayal of the wave‐speed stratification. The RCS facilitates the advancement of shallow refraction seismology through the application of current seismic reflection acquisition, processing and interpretation technology. 相似文献