近震资料确定三维地震速度结构的一种方法

我们研究了由近震得到的地震走时资料估计地震波速度三维结构的反演方法。这种方法由同时确定速度界面的二维深度分布、三维速度分布、台站校正和震源来表征。速度界面的深度分布和慢度扰动的分布以纬度、经度和深度的幂级数来拟合。这种三维地震结构的表示不仅有助于减少未知参数的个数,也有助于解析地求取界面深度、速度和走时关于未知参数偏导数的值。连接震源和台站之间的射线轨迹,是在具有确切界面深度分布的简化的速度分布条件下由称为试射法的射线追踪方法确定的。数值试验证明,我们的反演方法可成功地估计三维速度结构。     

用地震层析法确定三维结构和速度

为了能够从反针资杆同时得到体波速度和界面的三维结构,发展了一种用计算机处理的地震县析方法.这种介质由一些任意连续的三维弯曲界面的层组成,而这些界面把具有不同声学特性的均匀物质分隔成层.界面用一个多项式曲面表示. 恨定弹性波在其针线路径满足斯奈尔(Snen)定律的弯曲界面上传播或反针.时每个源检对的射线追综是通过解一个非线性方程组确定的.这种三维封线追综的方法在计算许多地震射钱中是有效的,其中包括转换波和多次波. 用公式来表示一种阵尼最小二乘反演方法,并且通过使观测到的走时与计算出的走时的差减到最小值来重建界面的层速度和三维结构.从综合模型的结果看出,这些解能很快收敛于真实模型. 1984年在温哥华岛曾用这种地震层析方法取得了一条可控震源的地震剖面,这条部面是岩石圈计划的一部分.通过用标称的二维地震反射刚量中的弯曲线,我们发现这种方法讨俯冲板块底面的三维成象是非常有用的.     

动水压力问题的边界单元法

本文介绍了用边界元法求解动水压力的原理与方法。在二维及三维情况下,分别给出了一些数值结果。计算表明,这种方法在分析坝-库水动力相互作用时为考虑动水压力提供了一种方便的途径。     

稀疏测线大地电磁资料三维反演研究:合成算例(英文)

受勘探成本和工区环境等因素的影响,当前大多数大地电磁实际工作采取布置稀疏测线采集数据和使用二维反演方法解释这些稀疏测线数据的方式。然而,二维反演方法在解释三维地电构造数据时存在局限性,有时甚至做出错误的地质解释。本文尝试了使用三维反演方法对大地电磁稀疏测线数据进行反演解释。使用大地电磁全信息资料三维共轭梯度反演程序对理论模型合成稀疏测线数据进行了三维反演。结果表明:这种反演方案是可行与有效的。同时,我们发现在不同数据的三维反演结果中,四个张量阻抗元素和两个倾子数据同时反演的结果相对更为准确,更接近理论模型。     

用地震资料测定岩石层的三维速度结构

本文发展了一种利用远震P波资料反演岩石层三维速度结构的新模型与方法;并将射线追踪方法引进上述模型与方法中。这种模型和方法同样可用于利用近震资料同时反演介质结构与震源参数。 数值实验结果表明,上述模型和方法是可行的。     

用FPS定量解释法实现三维近似成像的探讨

目的 在激发极化数据解释中，过去常采用原始曲线的特征值法，现介绍一种用FPS定量解释法确定异常体空间位置的三维近似成像方法。方法 用FPS定量解释法确定异常体空间位置的三维近似成像方法。结果 该方法可以较准确的确定异常体的上下左右边界以及异常体在垂直测线方向上的延伸情况。结论 算例表明，这种作图解释成像的方法效果较好，值得推广应用。     

基于数据规则化和稀疏反演的三维表面多次波压制方法

表面多次波是海洋地震勘探中的主要问题.目前,二维数据驱动的表面多次波压制技术(SRME)已经比较成熟,并且已经成为工业界压制海洋表面多次波的主流方法.但是由于二维SRME算法没有考虑横测线方向上多次波的贡献,导致在处理实际三维海洋资料时存在比较大的误差.将二维SRME算法扩展到三维空间后可以得到三维SRME算法,但是由于目前实际采集的三维海洋资料的观测系统存在拖缆漂移,而且横测线方向采样过于稀疏,直接应用三维SRME算法无法准确预测表面多次波.本文提出的通过数据规则化配合稀疏反演的三维表面多次波压制方法能够解决这种实际资料和三维SRME算法之间的矛盾.本文通过研究数据规则化与反规则化技术,使得数据分布满足三维SRME的要求;通过研究稀疏反演技术,有效解决了横测线方向采样稀疏对于多次波预测的影响,三维实际海洋资料的应用结果验证了方法的有效性和可行性.     

三维F－K域DMO方法及应用效果

在常速情况下，从三维时距方程出发，通过倾角代换，将三维倾角时差分裂为两个正交方向上的二维倾角时差，然后在这两个正交方向上分别做一次二维倾角时差校正，便可达到三维倾角时差校正的目的．应用实例表明这种三维ＤＭＯ分裂方法是切实有效的．     

三维F－K域DMO方法及应用效果

在常速情况下，从三维时距方程出发，通过倾角代换，将三维倾角时差分裂为两个正交方向上的二维倾角时差，然后在这两个正交方向上分别做一次二维倾角时差校正，便可达到三维倾角时差校正的目的．应用实例表明这种三维ＤＭＯ分裂方法是切实有效的．     

栅状地震剖面的投影制图

<正> 无论是正常测网的地震勘探,还是三维地震勘探,都可以组成栅状剖面。这种剖面能帮助解释人员想像地下地层的形态,确定断点的组合,进行测网内的追踪及闭合。轴测投影是机械制图中的常用方法,下面介绍如何用这种方法在计算机内制成栅状剖面。     

Substructure Method of Three-Dimensional Semi-Analytic Boundary Element and Its Application to Seismic Scattering of Hill Topography

In this paper,a substructure method of three-dimensional semi-analytic boundary element is established.The seismic scattering by three-dimensional topography of a hill can be analyzed by the method in frequency domain.Using this method,the computational effort and storage space are reduced considerably.Finally,analytic results are given.     

Magnetotelluric responses of three-dimensional conductive and magnetic anisotropic anomalies

A magnetotelluric finite-element modelling algorithm is developed, which is capable of handling three-dimensional conductive and magnetic anisotropic anomalies. Different from earlier three-dimensional magnetotelluric anisotropic modelling methods, the algorithm we presented has taken the magnetic anisotropy into consideration. The variational equations are produced by the Galerkin method and the governing equations are solved using a hexahedral vector edge finite-element method. The accuracy of this algorithm is firstly validated by comparing its solutions with the results of finite-difference method for a three-dimensional conductive arbitrary anisotropic model, and then validated by comparing with analytical solutions for a one-dimensional magnetic model. The responses of four kinds of models under different conditions are studied, and some conclusions are obtained. It shows that for materials with a high magnetic permeability, its influence on magnetotelluric responses cannot be ignored in some circumstances. Especially, if the magnetic susceptibility is exceptionally high, it may really distort the apparent resistivities of lower resistive anomalies. These conclusions are also beneficial for magnetotelluric survey.     

塔里木地区勘探地震正演模拟研究

以塔里木前陆盆地地质背景为模型,采用三维任意差分精细积分方法并行算法实现了库车地区三维正演模拟.三维任意差分精细积分方法通过时间域采用局部积分半解析方法求得波动方程的递推算子,与常规的差分法相比计算精度有较大提高;文中计算稳定性根据实际算例进行分析,采用稳定因子约束,得到较好的稳定性;边界条件采用改进的自适应吸收边界,并通过串行程序并行化,大大减少了三维正演模拟的耗时,完成了大数据量三维正演模拟.     

Three-dimensional dynamic soil–structure interaction analysis in the time domain

A new numerical procedure is proposed for the analysis of three-dimensional dynamic soil–structure interaction in the time domain. In this study, the soil is modelled as a linear elastic solid, however, the methods developed can be adapted to include the effects of soil non-linearities and hysteretic damping in the soil. A substructure method, in which the unbounded soil is modelled by the scaled boundary finite-element method, is used and the structure is modelled by 8–21 variable-number-node three-dimensional isoparametric or subparametric hexahedral curvilinear elements. Approximations in both time and space, which lead to efficient schemes for calculation of the acceleration unit-impulse response matrix, are proposed for the scaled boundary finite-element method resulting in significant reduction in computational effort with little loss of accuracy. The approximations also lead to a very efficient scheme for evaluation of convolution integrals in the calculation of soil–structure interaction forces. The approximations proposed in this paper are also applicable to the boundary element method. These approximations result in an improvement over current methods. A three-dimensional Dynamic Soil–Structure Interaction Analysis program (DSSIA-3D) is developed, and seismic excitations (S-waves, P-waves, and surface waves) and externally applied transient loadings can be considered in analysis. The computer program developed can be used in the analysis of three-dimensional dynamic soil–structure interaction as well as in the analysis of wave scattering and diffraction by three-dimensional surface irregularities. The scattering and diffraction of seismic waves (P-, S-, and Rayleigh waves) by various three-dimensional surface irregularities are studied in detail, and the numerical results obtained are in good agreement with those given by other authors. Numerical studies show that the new procedure is suitable and very efficient for problems which involve low frequencies of interest for earthquake engineering. Copyright © 1999 John Wiley & Sons Ltd     

基于重力卫星数据监测地表质量变化的三维点质量模型法

利用卫星重力测量手段监测全球质量变化取得了巨大成功,本文基于牛顿万有引力定律在三维空间直角坐标系中导出利用重力卫星观测数据监测全球质量变化的三维点质量模型法,该方法可直接利用重力卫星的轨道和星间观测数据或时变重力场模型计算全球质量变化,由于利用卫星观测数据计算地表质量变化的向下延拓过程以及观测数据噪声的影响,需要采用合适的空间约束方程或正则化技术对解算结果进行约束或平滑处理.利用合成全球质量变化模型模拟一个月的GRACE双星轨道和星间距离变率数据计算全球质量变化,对三维点质量模型法进行分析验证,采用零阶Tikhonov正则化技术处理病态问题.结果表明,三维点质量模型法可有效用于重力卫星观测数据监测全球质量变化,为利用重力卫星观测数据监测全球质量变化提供一种可选的途径.     

One-dimensional magnetotelluric inversion using an adaptation of Zohdy's resistivity method

An iterative refinement method for determining a layered resistivity model from a Schlumberger or Wenner sounding curve is adapted to determine a layered resistivity model by using apparent resistivity and phase derived from the magnetotelluric impedance. Magnetotelluric observations presented as a function of period are first converted to an approximate resistivity–depth profile using Schmucker's transformation and this is used to construct an initial guess (starting) model. A two-stage procedure is then invoked. Keeping resistivities constant, layer boundaries are first adjusted to give a minimum misfit between measured data and responses and this is followed by resistivity adjustments with fixed layer boundaries to reduce the misfit further. The method is illustrated by application to some synthetic data both exact and with added noise, to a real field data set and to some magnetotelluric profile data obtained in a survey over the Carnmenellis granites in south Cornwall. The method is validated by recovering conductivity models from the exact and noisy 1D synthetic data. For complicated three-dimensional data at a single site and along a profile of stations, the method is shown to produce acceptable solutions which may be used as starting models in further two- or three-dimensional studies.     

3-D seismic response analysis of long lined tunnels in half-space

The dynamic response of infinitely long lined tunnels with a uniform cross-section buried into an elastic or viscoelastic half-space to body and surface harmonic seismic waves is numerically determined by a special direct boundary element method in the frequency domain. The waves have an arbitrary direction of propagation with respect to the axis of the tunnel and this renders the problem three-dimensional. However, this problem is effectively reduced to a two-dimensional one by a coordinate transformation and appropriate integration of the full-space dynamic fundamental solution along the direction of the tunnel axis. Quadratic isoparametric boundary line elements and advanced numerical integration techniques for the treatment of singular integrals produce results of high accuracy. Numerical results are presented for the case of an infinitely long lined tunnel of circular cross-section and compared against those of a full three-dimensional boundary element analysis, as well as those of other methods. Thus the proposed method is illustrated and its merits demonstrated.     

西安咸阳国际机场污水排放口潜在滑坡稳定性分析

对西安咸阳国际机场污水排放口潜在滑坡的工程地质勘察结果表明,该滑坡处于极限平衡状态。采用B ish-op法参数反演和直快剪试验结果,分别用二维、三维极限平衡法和二维、三维有限元法,对该滑坡进行稳定性计算。结果表明：有限元法（二维、三维）的计算结果偏大,二维极限平衡法的计算结果最小;极限平衡法（二维、三维）与二维有限元法宜选取试验结果的大值平均值,三维有限元法则取试验结果的小值平均值较为合理。最后,根据三维有限元法的计算结果,提出了削坡压脚的防护措施。     

饱和多孔介质三维时域黏弹性人工边界与动力反应分析的显式有限元法

本文基于Biot的饱和多孔介质本构方程,考察具有辐射阻尼的外行球面波,推导了饱和多孔介质三维黏弹性人工边界的法向和切向边界方程;在已有的饱和多孔介质二维显式有限元数值计算方法基础上,提出该理论的三维方法,并开发了实现该三维方法的有限元程序.算例表明饱和多孔介质三维时域黏弹性人工边界与动力反应分析的显式有限元法具有较好的精度和稳定性.     

Feasibility of reducing three-dimensional wave energy by introducing periodic foundations

A new seismic isolation foundation called the periodic foundation is proposed, which can be described as a three-dimensional typical cell consisting of a high density core, a soft coating and a concrete matrix. Utilizing the attenuation zones (AZs) resulted by the periodicity of the foundation, the mechanism of this new seismic isolation system is different from the traditional ones. The seismic waves with the frequencies in the AZs cannot propagate across the foundation. Thus, the seismic responses on the top surface of the periodic foundation can be reduced significantly. In this paper, the dispersion curves of the three-dimensional three-component (3D–3C) periodic foundations are analyzed by the finite element method (FEM). The influencing factors such as physical and geometrical parameters of the typical cell are investigated. In order to verify its feasibility, a three-dimensional (3D) soil-foundation finite element model is analyzed, and the frequency zones of seismic attenuations are consistent with the AZs. The periodic foundation can greatly reduce the seismic response, which implies that the periodic foundation has a great potential application on seismic isolation.