共查询到20条相似文献,搜索用时 312 毫秒
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.
定向设计的震源、检波器组合可增强特定方向传播的地震波能量,同时压制其他方向噪声,从而改善低信噪比数据的资料品质.本文提出了一种基于地质模型的定向组合激发-接收技术方法,该技术方法包含三个步骤:定向方向计算、定向组合参数设计、定向组合数据处理.首先对地质模型做波动方程正演并统计地震波出射到地表的能量和角度(方向),然后优化设计组合参数将波场能量定向聚焦到该方向,在进行激发-接收定向组合处理后,最终实现提升地下弱照明地区成像质量的目的.基于此思路,本文首先分析了地震波场产生方向性的原因,以模型照明结果为例展示了地质结构与波场传播方向的相关性;然后定义了表征波场传播方向和能量大小的能流密度矢量,在此基础上推导出声学介质条件下地表出射波场的角度和能量计算方法;最后给出了组合叠加波场的表达式和求解各组合参数的目标函数.经理论模型和实际资料验证,本文提出的定向组合方法能显著改善资料信噪比,在低信噪比地区应用前景广阔. 相似文献
3.
滑坡的地震波反演是确定滑坡动力学过程的有效途径.然而,目前的滑坡地震波反演大多基于点源模型,将滑坡面近似为点源,仅考虑滑体作用于滑床的摩擦力及其变化.但实际滑坡过程中,由于重力和滑动距离的存在,滑床还受到一个逐渐增加的正压力力偶作用.为定量探讨这一正压力力偶对滑坡地震波场的贡献,我们基于有限滑坡面模型,分别正演模拟了滑坡过程中摩擦力和正压力力偶的地震波场,得到了正压力力偶地震波场在整个滑坡地震波场中的能量比例及其随距离的变化.研究发现,正压力力偶地震波场具备近场项衰减特征,相较于摩擦力波场随距离衰减更快.正演模拟显示,正压力力偶地震波场在滑坡尺度2~3倍距离范围内较为显著;但在更远距离上,由于其幅度大致随距离的平方衰减,其贡献几可忽略.此外,正压力力偶波场在波形特征上与摩擦力波场也存在明显区别,二者的相对大小与摩擦系数和坡度等参数密切相关.本文模拟分析表明,完整的滑坡近场地震波场应该同时考虑摩擦力和正压力力偶的贡献;滑坡动力学参数反演中,利用近场地震波记录,考虑有限尺度的滑坡面,原则上我们不仅可以确定传统滑坡动力学参数,还有可能确定摩擦系数和流体压力等参数在空间上的变化分布. 相似文献
4.
深反射地震剖面法为了获取深部结构特征常常采取大的偏移距采集数据.目前公开发表的相关资料中,鲜有利用深反射地震炮集数据获取近地表的结构特征.为此,本文通过正演测试了相关数据处理流程,即利用有限差分正演了起伏地表模型的大偏移距地震单炮弹性波场特征,通过共检波点域面波信号F-K频谱叠加构建新方法,从深反射地震数据集中提取了高品质的多阶面波频散曲线,再利用多阶面波联合反演获得了近地表的结构特征.在前述正演流程基础上,利用跨越班公湖—怒江缝合带的SinoProbe深反射地震剖面中的实际炮集数据,求取了基阶和一阶瑞利波频散曲线,联合反演后得到近地表横波速度结构.该结果与初至波走时反演获取的纵波速度结构具有较好的一致性,且在近地表的浅层分辨率较纵波速度结构特征更高,而更与已有地质认识相吻合.本文提供的相关数据处理流程表明利用深反射地震炮集数据,也能够获取近地表浅层的横波速度结构. 相似文献
5.
深反射地震剖面法为了获取深部结构特征常常采取大的偏移距采集数据.目前公开发表的相关资料中,鲜有利用深反射地震炮集数据获取近地表的结构特征.为此,本文通过正演测试了相关数据处理流程,即利用有限差分正演了起伏地表模型的大偏移距地震单炮弹性波场特征,通过共检波点域面波信号F-K频谱叠加构建新方法,从深反射地震数据集中提取了高品质的多阶面波频散曲线,再利用多阶面波联合反演获得了近地表的结构特征.在前述正演流程基础上,利用跨越班公湖—怒江缝合带的SinoProbe深反射地震剖面中的实际炮集数据,求取了基阶和一阶瑞利波频散曲线,联合反演后得到近地表横波速度结构.该结果与初至波走时反演获取的纵波速度结构具有较好的一致性,且在近地表的浅层分辨率较纵波速度结构特征更高,而更与已有地质认识相吻合.本文提供的相关数据处理流程表明利用深反射地震炮集数据,也能够获取近地表浅层的横波速度结构. 相似文献
6.
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. 相似文献
7.
低速带具有“自由表面、低速度、高吸收”的特征,能引起反射波能量减弱,产生很强的频散,造成近地表波场畸变,严重影响地震勘探的精度。而地震采集过程又不可避免遇到以下低速带问题:①低速带存在与否会对地震波传播造成什么影响;②近地表变化大,怎么才能得到精确的低速带参数;③存在低速带时,激发井深、最佳的药量(与子波频率有关)、最佳的耦合(与低速带速度相关)等采集参数会对采集资料的品质产生什么影响。本文通过数值模拟,集中对低速带与地震波激发的井深、速度、子波主频等核心采集参数的关系进行分析研究,同时在低速带背景下对数值模拟的网格、边界等问题也进行了讨论。模拟得到的波场特征及低速带参数影响的机理分析,对做好野外地震波激发,获取高品质的地震资料,具有理论意义和指导作用。 相似文献
8.
南黄海海相地层的地震波场特征和层位标定一直是困扰地震勘探的重要问题.为了近距离、高精度和高分辨率地观测井周围构造特征和岩石性质引起的波场变化,为地震资料的采集、处理与解释提供地震波衰减规律、速度与层位标定等信息,对大陆架科学钻探CSDP-2井实施了近零偏移距垂直地震剖面(VSP)观测.针对海相地层顶部强反射界面地震波穿透难的问题,采用了大容量气枪震源并设计了气枪阵列组合方式,提高了激发地震波的能量,获得了强反射界面之下清晰的PP、PS下行波和上行波信号.采用了三分量偏振合成、组合滤波和波场分离等处理方法,对VSP观测数据进行处理,获得了海相三叠系—志留系的精细的纵波、横波速度结构和地层吸收因子等物性数据,建立了钻井地层、测井、VSP上行波和多道地震剖面对应关系,实现了不同尺度的地质和地球物理属性资料的有效衔接,标定了钻井地质剖面上各深度地质体的地震反射特性,厘定了过井地震剖面上反射同相轴的地质属性.此次观测取得的纵波、横波速度信息,成为建立南黄海海相地层速度模型主要的资料来源,也是地震资料的岩性反演处理不可缺少的信息. 相似文献
9.
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. 相似文献
10.
本文论述准各向异性粘弹波动方程有限元数值解地震波响应的原理与计算方法,并给出地震衰减因子与方程中粘滞系数的计算关系式。列出了所研制的地面记录地震剖面与井中记录垂直地震剖面,以及波场时间切片图、散度与旋度图、应变能与动能的时间切片图。本方法适于任意分层形态且层间物性参数为强间断的实际地震剖面的仿真,旨在研究正问题以实现反问题。 相似文献
11.
随着地震勘探开发的不断深入和发展,地震勘探的主战场逐渐向复杂地区转移.复杂地区既指地表条件复杂的地区,也指地下地质构造和地层条件复杂的地区,这些都对地震勘探提出了新的挑战和更高的要求.地震正演模拟正是开展此类问题研究的一个重要手段和方法.目前市场上具有正演模拟功能的软件大多是根据射线理论采用射线追踪的方法来完成正演模拟的,这种方法不能很好地反映地震波的动力学特征,特别是在复杂地区难以得到正确的结果.本文利用高阶有限差分有效克服了常规有限差分算法求解波动方程的频散问题,并以高效的OpenMP并行计算模式进行了并行优化,较大程度上提高了正演计算的速度和精度;同时实现了二阶Higdon边界条件,改善了边界吸收效果;也在一定程度上提高了计算速度.塔中地区主要目的层埋藏深,逆断层发育,地震反射特征复杂,增加了地震勘探的难度.本文依据该地区的地质模型,利用波动方程正演技术论证了该地区的地震采集观测系统,为该区地震采集观测系统的设计提供了科学依据. 相似文献
12.
本文论述准各向异性粘弹波动方程有限元数值解地震波响应的原理与计算方法,并给出地震衰减因子与方程中粘滞系数的计算关系式。列出了所研制的地面记录地震剖面与井中记录垂直地震剖面,以及波场时间切片图、散度与旋度图、应变能与动能的时间切片图。本方法适于任意分层形态且层间物性参数为强间断的实际地震剖面的仿真,旨在研究正问题以实现反问题。 相似文献
13.
对近年来有关地震波衰减的物理机制进行评述和初步研究认为,对于地壳介质,造成地震波衰减的主要原因是地壳内存在大量裂隙,裂纹中饱含水或部分含水,地震波传播时引起裂隙中的流体运动,从而造成地震波的衰减,对于上地幔,扩散控制的位错阻尼机制可能是造成地幔地震波衰减的一个重要原因,而上地幔软流层部分熔融的存在也是造成地震波衰减不可忽视的原因. 相似文献
14.
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. 相似文献
15.
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. 相似文献
16.
煤层中存在的裂隙会导致介质表现为各向异性,本文以HTI型煤层为例,结合各向异性介质弹性矩阵和各向异性裂隙理论,推导出不同充填物的垂直裂隙中各向异性参数表达式,将其应用于地震波响应分析;通过改进的交错网格差分法和各向异性Christoffel方程波场分解法,得到地震波合成记录和分解后的P波和SV波记录;将Thomsen群速度与相速度公式,经过坐标轴旋转变换,得到HTI型煤层中不同各向异性参数的地震波速度响应表达式;建立不同类型煤层地质模型,分析了裂隙密度、裂隙充填物以及煤层厚度等参数变化时的地震波响应特征.研究结果为分析垂向裂隙各向异性薄煤层地震波传播规律提供工具,为选用相应地震数据进行地震波各向异性参数反演提供依据. 相似文献
17.
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. 相似文献
18.
利用时域有限差分方法进行地震波模拟能够得到丰富的波场信息, 对研究复杂地质构造条件下的地震波波场十分有效。 该文在对时域有限差分方法进行概述的基础上, 以2002年4月施测的安新—宽城地震测深剖面资料为基础, 采用时域有限差分方法对该剖面进行二维全波数值模拟, 同时以理想匹配层作为模拟的吸收边界条件。 将模拟得到的合成地震图, 通过与野外实验得到的记录截面和射线追踪合成记录相比较和讨论, 验证了时域有限差分方法在模拟地震波传播过程的有效性。 相似文献
19.
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. 相似文献
20.
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. 相似文献
|