共查询到20条相似文献，搜索用时 31 毫秒
1.
In this work, a hybrid boundary integral equation method (BIEM) is developed, based on both displacement and hypersingular traction formulations, for the analysis of time-harmonic seismic waves propagating through cracked, multi-layered geological regions with surface topography and under plane strain conditions. Specifically, the displacement-based BIEM is used for a multi-layered deposit with interface cracks, while the regularized, traction-based BIEM is used when internal cracks are present within the layers. The standard uni-dimensional boundary element with parabolic shape functions is employed for discretizing the free surface and the layer interfaces, while special discontinuous boundary elements are placed near the crack tips to model the asymptotic behaviour of both displacements and tractions. This formulation yields displacement amplitudes and phase angles on the free surface of a geological deposit, as well as stress intensity factors near the tips of the cracks. Finally, in the companion paper, numerical results are presented which show that both scattered wave and stress concentration fields are sensitive to the incidence seismic wave parameters and to specific site conditions such as surface topography, layering, the presence of cracks and crack interaction. 相似文献
2.
深反射地震剖面法为了获取深部结构特征常常采取大的偏移距采集数据.目前公开发表的相关资料中,鲜有利用深反射地震炮集数据获取近地表的结构特征.为此,本文通过正演测试了相关数据处理流程,即利用有限差分正演了起伏地表模型的大偏移距地震单炮弹性波场特征,通过共检波点域面波信号F-K频谱叠加构建新方法,从深反射地震数据集中提取了高品质的多阶面波频散曲线,再利用多阶面波联合反演获得了近地表的结构特征.在前述正演流程基础上,利用跨越班公湖—怒江缝合带的SinoProbe深反射地震剖面中的实际炮集数据,求取了基阶和一阶瑞利波频散曲线,联合反演后得到近地表横波速度结构.该结果与初至波走时反演获取的纵波速度结构具有较好的一致性,且在近地表的浅层分辨率较纵波速度结构特征更高,而更与已有地质认识相吻合.本文提供的相关数据处理流程表明利用深反射地震炮集数据,也能够获取近地表浅层的横波速度结构. 相似文献
3.
A discrete wave number approach in conjunction with a propagator-based formalism is used to synthesize the Loma Prieta earthquake ground motion at both the near and the far field, taking into account all kinds of seismic waves (body and surface). A bilaterally propagating shear slip over a rectangular fault is used to describe the seismic source mechanism, while the earth model is based on geological profiles of the Santa Cruz mountain area and consists of three layers overlaying a half-space. The synthesized ground motion is first compared with actual records from the Loma Prieta earthquake and the agreement between the two is found to be satisfactory, as far as magnitude, duration and essential wave form characteristics are concerned. Then, ground motions are synthesized and plotted at a dense grid of observer locations over a large area around the epicenter, at different time instants. Using such plots, it is possible to study the generation ana propagation of different kinds of seismic waves, the spatial variability of ground motion, as well as the development of the permanent gound deformation. 相似文献
4.
低速带具有“自由表面、低速度、高吸收”的特征,能引起反射波能量减弱,产生很强的频散,造成近地表波场畸变,严重影响地震勘探的精度。而地震采集过程又不可避免遇到以下低速带问题:①低速带存在与否会对地震波传播造成什么影响;②近地表变化大,怎么才能得到精确的低速带参数;③存在低速带时,激发井深、最佳的药量(与子波频率有关)、最佳的耦合(与低速带速度相关)等采集参数会对采集资料的品质产生什么影响。本文通过数值模拟,集中对低速带与地震波激发的井深、速度、子波主频等核心采集参数的关系进行分析研究,同时在低速带背景下对数值模拟的网格、边界等问题也进行了讨论。模拟得到的波场特征及低速带参数影响的机理分析,对做好野外地震波激发,获取高品质的地震资料,具有理论意义和指导作用。 相似文献
5.
南黄海海相地层的地震波场特征和层位标定一直是困扰地震勘探的重要问题.为了近距离、高精度和高分辨率地观测井周围构造特征和岩石性质引起的波场变化,为地震资料的采集、处理与解释提供地震波衰减规律、速度与层位标定等信息,对大陆架科学钻探CSDP-2井实施了近零偏移距垂直地震剖面（VSP）观测.针对海相地层顶部强反射界面地震波穿透难的问题,采用了大容量气枪震源并设计了气枪阵列组合方式,提高了激发地震波的能量,获得了强反射界面之下清晰的PP、PS下行波和上行波信号.采用了三分量偏振合成、组合滤波和波场分离等处理方法,对VSP观测数据进行处理,获得了海相三叠系—志留系的精细的纵波、横波速度结构和地层吸收因子等物性数据,建立了钻井地层、测井、VSP上行波和多道地震剖面对应关系,实现了不同尺度的地质和地球物理属性资料的有效衔接,标定了钻井地质剖面上各深度地质体的地震反射特性,厘定了过井地震剖面上反射同相轴的地质属性.此次观测取得的纵波、横波速度信息,成为建立南黄海海相地层速度模型主要的资料来源,也是地震资料的岩性反演处理不可缺少的信息. 相似文献
6.
Seismic attenuation in Faroe Islands basalts 总被引：2，自引：1，他引：1
F. Shaw M.H. Worthington R.S. White M.S. Andersen U.K. Petersen the Seifaba Group 《Geophysical Prospecting》2008,56(1):5-20
We analysed vertical seismic profiling (VSP) data from two boreholes at Glyvursnes and Vestmanna on the island of Streymoy, Faroe Islands, to determine the magnitude and causes of seismic attenuation in sequences of basalt flows. The work is part of SeiFaBa, a major project integrating data from vertical and offset VSP, surface seismic surveys, core samples and wireline log data from the two boreholes. Values of effective seismic quality factor (Q) obtained at Glyvursnes and Vestmanna are sufficiently low to significantly degrade the quality of a surface reflection seismic image. This observation is consistent with results from other VSP experiments in the North Atlantic region. We demonstrate that the most likely cause of the low values of effective Q at Glyvursnes and Vestmanna is a combination of 1D scattering and intrinsic attenuation due to seismic wave‐induced fluid flow within pores and micro‐cracks. Tests involving 3D elastic wave numerical modelling with a hypothetical basalt model based on field observations, indicate that little scattering attenuation is caused by lateral variations in basalt structure. 相似文献
7.
正演模拟技术在地震采集设计中的应用 总被引：1，自引：0，他引：1
随着地震勘探开发的不断深入和发展,地震勘探的主战场逐渐向复杂地区转移.复杂地区既指地表条件复杂的地区,也指地下地质构造和地层条件复杂的地区,这些都对地震勘探提出了新的挑战和更高的要求.地震正演模拟正是开展此类问题研究的一个重要手段和方法.目前市场上具有正演模拟功能的软件大多是根据射线理论采用射线追踪的方法来完成正演模拟的,这种方法不能很好地反映地震波的动力学特征,特别是在复杂地区难以得到正确的结果.本文利用高阶有限差分有效克服了常规有限差分算法求解波动方程的频散问题,并以高效的OpenMP并行计算模式进行了并行优化,较大程度上提高了正演计算的速度和精度;同时实现了二阶Higdon边界条件,改善了边界吸收效果;也在一定程度上提高了计算速度.塔中地区主要目的层埋藏深,逆断层发育,地震反射特征复杂,增加了地震勘探的难度.本文依据该地区的地质模型,利用波动方程正演技术论证了该地区的地震采集观测系统,为该区地震采集观测系统的设计提供了科学依据. 相似文献
8.
9.
10.
Monitoring of geophysical conditions of marine sedimentary basins is necessary for predicting seismic events and for adaptation of geothermal technologies for seismically active (as a rule) sea bottom geothermal areas. These conditions are characterized by seismo-hydro-electromagnetic (EM) geophysical field interaction in the presence of gravity. Based on the main physical principles, geophysical and petrophysical data, we formulate a mathematical model of seismo-hydro-EM interaction in a basin of a marginal sea and calculate the transformation of a seismic excitation in the upper mantle under the central part of the sea of Japan into the low-frequency (0.1 to 10 Hz) EM signals at the top of the sea bottom sedimentary layer, at the sea surface and in the atmosphere up to the lower boundary of the ionosphere. Physics of the EM generation and propagation process is shown including: generation of EM waves in the upper mantle layer M by a seismic wave from under M, spatial modulation of diffusive EM waves by a seismic wave, stopping of the EM wave arrived (before the seismic P wave) from the upper mantle at the top of the sediments because of the high electric conductivity of seawater (3.5 S/m), immediate penetration of the EM wave through the seawater thickness after the delayed seismic P wave shock into the sea bottom, and EM emission from the sea surface into the atmosphere. Let us note that the EM signal in the sea bottom sediments is the first measurable signal of a seismic activation of geological structures beneath the seafloor and this signal is protected by seawater from the influence ionosphere disturbances. Amplitude of the computed magnetic signals (300, 200, 50, and 30 pT at the ocean–atmosphere interface and at the height of 10, 30 and 50 km, respectively), their predominant frequency (0.25 Hz), the delay of the seismic P wave in regard to the magnetic signal for the receivers at the shore (20 s), the amplitude of temperature disturbances in sediments (up to 0.02 K), the parameters of the long (150 km) tsunami wave of a small (up to 20 cm) amplitude far from the shore and other values that characterize the seismo-hydro-EM process are of the orders observed. Recommendations for the EM monitoring of dynamic processes beneath seafloor geothermal areas are given. 相似文献
11.
Michel Bouchon 《Pure and Applied Geophysics》1996,148(1-2):3-20
We review the application of the discrete wave number method to problems of scattering of seismic waves formulated in terms of boundary integral equation and boundary element methods. The approach is based on the representation of the diffracting surfaces and interfaces of the medium by surface distributions of sources or by boundary source elements, the radiation from which is equivalent to the scattered wave field produced by the diffracting boundaries. The Green's functions are evaluated by the discrete wave number method, and the boundary conditions yield a linear system of equations. The inversion of this system allows the calculation of the full wave field in the medium. We investigate the accuracy of the method and we present applications to the simulation of surface seismic surveys, to the diffraction of elastic waves by fractures, to regional crustal wave propagation and to topographic scattering. 相似文献
12.
煤层中存在的裂隙会导致介质表现为各向异性,本文以HTI型煤层为例,结合各向异性介质弹性矩阵和各向异性裂隙理论,推导出不同充填物的垂直裂隙中各向异性参数表达式,将其应用于地震波响应分析;通过改进的交错网格差分法和各向异性Christoffel方程波场分解法,得到地震波合成记录和分解后的P波和SV波记录;将Thomsen群速度与相速度公式,经过坐标轴旋转变换,得到HTI型煤层中不同各向异性参数的地震波速度响应表达式;建立不同类型煤层地质模型,分析了裂隙密度、裂隙充填物以及煤层厚度等参数变化时的地震波响应特征.研究结果为分析垂向裂隙各向异性薄煤层地震波传播规律提供工具,为选用相应地震数据进行地震波各向异性参数反演提供依据. 相似文献
13.
14.
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. 相似文献
15.
Cachar area, a part of the Assam–Arakan frontal fold belt in the eastern part of India is characterized by a series of narrow anticlines which are dissected by a number of faults including thrusts, and broad synclines. Out of several exposed structures, Badarpur anticline is most important because of its major share of hydrocarbon strikes in Cachar. However, the seismic image of this area so far has been extremely poor for deciding locations for exploratory drilling. The problems in obtaining a good seismic image in this area include proper spread configuration including far and near offsets to capture reflected signals from the complex subsurface, variation of charge-sizes and depths along a line, large variation in statics due to rapid changes in elevation and near surface properties. The presence of coherent noise due to scattering and proper migration are major processing problems. This paper deals with the problem of spread configuration. For areas of complex subsurface, spread configuration obtained through field experiment at one spot may be quite different from other spots, requiring a large number of costly experiments. In this work, therefore, we used numerical modeling to derive an optimum spread configuration. For this, we chose a line passing through four wells on Badarpur anticline, and prepared a depth model using available seismic, surface geological and well data. Continuous synthetic common shot records were generated through oblique incidence ray tracing and were processed to obtain stacks with various spread configurations. Comparison of these stacks of synthetic data shows that end-on shooting configuration is inadequate for acquiring a good seismic image of the complex subsurface. From the present study it is evident that a split–spread configuration with near and far offsets of 50 m and 2400 m, respectively, can bring out optimum seismic images in Badarpur and other anticlines in this area. 相似文献
16.
A numerical study on the influence that cracks and discontinuities (closed cracks) can have on the seismic response of a hypothetical soil–structure system is presented and discussed. A 2-D finite-difference model of the soil was developed, considering a bilinear failure surface using a Mohr–Coulomb model. The cracks are simulated with interface elements. The soil stiffness is used to characterize the contact force that is generated when the crack closes. For the cases studied herein, it was considered that the crack does not propagate during the dynamic event. Both cases, open and closed cracks, are considered. The nonlinear behavior was accounted for approximately using equivalent linear properties calibrated against several 1-D wave propagation analyses of selected soil columns with variable depth to account for changes in depth to bed rock. Free field boundaries were used at the edges of the 2-D finite-difference model to allow for energy dissipation of the reflected waves. The effect of cracking on the seismic response was evaluated by comparing the results of site response analysis with and without crack, for several lengths and orientations. The changes in the response obtained for a single crack and a family of cracks were also evaluated. Finally, the impact that a crack may have on the structural response of nearby structures was investigated by solving the seismic-soil–structure interaction of two structures, one flexible and one rigid to bracket the response. From the results of this investigation, insight was gained regarding the effect that discontinuities may have both on the seismic response of soil deposits and on nearby soil–structure systems. 相似文献
17.
高密度空间采样地震数据特点分析 总被引：3，自引：0，他引：3
中国陆相沉积盆地的地质结构复杂,储层岩性的多变,需要有高精度的勘探方法。高密度空间采样是提高地震勘探精度的一项新技术。本文简要说明了点激发和点接收技术,分析高密度空间采样的野外工作方法,介绍了Gijs j．o．Vermeer提出的对称采样原理,从波场连续性的角度探讨了高密度空间采样技术。重点分析高密度空间采样数据的特点, 即:高密集的初至波有利于近地表结构的调查,可提高静校正的精度;小偏移距、小点距接收增加了浅层的有效覆盖次数,提高浅层反射的成像精度; 对规则噪声可实现无假频采样,室内的噪声分析和噪声压制的精度随着提高;高密度采集使波场连续性增强,使得各种数学变换的精度提高,有利于不同波场的分离。最后指出高密度空间采样地震勘探技术的难点在于室内数据处理,针对海量数据的分析和处理方法还需要进行深入细致的研究工作。 相似文献
18.
《Journal of Applied Geophysics》2006,58(2):87-98
A geophysical campaign to characterize the subsurface of a contaminated site down to a depth of several tens of meters was carried out under the HYGEIA-CEE project. On this site, seismic techniques were combined to image the geological structures; i.e. seismic reflection, P-wave tomography and spectral analysis of surface waves. Because these techniques consider different wave components in the processing, they can be expected to provide complementary information concerning the site lithology. The special feature of this experiment is the fact that the same seismic acquisition device, consisting of a mobile central unit, a drop-weight seismic source, and a sensor line of gimbal mounted geophones, was used for each of the techniques. Two perpendicular seismic lines were set up in the field for testing two geophone spacings. Three processing procedures, one each for the seismic reflection, P-wave tomography and spectral analysis of surface waves, were developed for producing seismic images from the P-wave reflectivity, the first P-wave arrivals and the dispersion of Rayleigh waves, respectively. The images show good complementarity in terms of investigation depth. The results are also in good agreement with available borehole data: the sandy layers seem to be related to low velocities, since the high velocities are better explained by the presence of clayey and gravelly intervals. The contribution and the limits of this seismic multi-approach method is discussed. 相似文献
19.
本文用位错模式研究了含粗糙面断层的远场地震体波辐射。假设粗糙面大小和整个错动面大小的比值为C,并假设破裂从中心开始,计算了不同C值对远场地震波谱的影响。结果表明,与均匀圆盘状断层的远场振幅谱相比,不同之处是,由于粗糙面的存在,当裂纹前端传播至粗糙面的边缘处时,引起波谱开始部分的分裂,表现在谱的第一个波谷处。粗糙面的大小对振幅谱的影响不显著。另外,比较P波和S波振幅谱,发现S波波谱的拐角频率比P波的低。 相似文献
20.
Iwona Stan-Kłeczek 《Acta Geophysica》2008,56(4):1065-1073
The existence of cracks in a rock mass causes the reduction of seismic wave velocity and the anisotropy of seismic wave velocity, which is characteristic for rocks with preferred orientation of cracks. The present study concerns sedimentary rocks (sandstone, limestone) and igneous rocks (basalt, granite). Studying the relationship between seismic anisotropy and cracks anisotropy in rocks I was only interested in fractures perpendicular to the layering. This allowed me to calculate two-dimensional crack tensors and velocity tensors in planes parallel to the layer surface. An application of tensor calculus enables to take into account both geometry and orientation of cracks. The obtained results confirm that the directions of major axes of second-rank velocity tensor prove the relationship consistent with predictions of theoretical models. Owing to these dependences, seismic methods can be used to study the cracks anisotropy in rocks inaccessible for direct observations. 相似文献