用三维半解析边界元的子结构法分析凸起山包对地震波的散射

本文建立了一种三维半解析边界元子结构法。这种方法可在频域就三维轴对称凸起山包对地震波的散射进行分析。这种方法节省存储空间、计算工作量小。另外还叙述了拐角的处理方法。     

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.     

Shinkansen high-speed train induced ground vibrations in view of viaduct–ground interaction

Train viaduct behavior and nearby ground motion under the high-speed train passage have been studied in this paper. First, the findings from the field measurement alongside the high-speed Shinkansen railway in Japan are interpreted. Then, the computer simulation is made based on the soil-foundation-viaduct interaction analysis under moving axle loads. The solution method is to apply the dynamic substructure method in the frequency domain. The viaduct girders including track structure and pier supports are modeled by the three-dimensional beam-column elements. The supporting pile foundation and nearby field are discretized by the axisymmetric three-dimensional finite elements and analyzed in a semi-analytical way, with a transmitting boundary replacing the far field based on the thin layer element method. Nearby ground motion during train passage on a viaduct have been calculated by superimposing the effects from neighboring pile foundations.The main parameters affecting viaduct vibrations are discussed by taking environmental vibration into consideration. The nearby ground motion along the viaduct is recomputed by applying the above determined forces to the foundation tops. The results from numerical studies are compared with the field test data, thus proving the present simulation to be effective and reliable.     

某人防工程地震响应分析与抗震性能评价

目前为止还没有成熟的地下结构抗震性能评价方法.本文结合地下结构地震响应的特点及已有相关规范,初步建立了地下结构抗震性能评价方法.并以某人防工程为例,建立三维分析模型,计算其在两种不同幅值的人工合成波作用下的三维地震响应,计算了结构的水平横向相对位移,评价了该人防工程的抗震性能.评价结果表明此人防工程具有良好的抗震性能....     

Three-dimensional analysis of spatially varying ground motions around a uniform canyon in a homogeneous half-space

Department of Civil Engineering, University of California, Berkeley, CA 94720, U.S.A. A direct boundary element method to determine the three-dimensional seismic response of an infinitely-long canyon of arbitrary but uniform cross-section cut in a homogeneous viscoelastic half-space is presented. The seismic excitation is represented by P, SV, SH or Rayleigh waves at arbitrary angles with respect to the axis of the canyon. The accuracy of the procedure and implementing computer program is demonstrated by comparison with previous solutions for the limiting case of two-dimensional response, recently obtained three-dimensional response results for infinitely-long canyons, and three-dimensional boundary method solutions presented in this paper for finite canyons.     

三维土-结构动力相互作用的一种时域直接分析方法

本文提出了一种分析三维土-结构动力相互作用的时域直接方法。该方法采用集中质量显式有限元和透射人工边界模拟无限域地基,通过编制的FORTRAN程序实现;采用ANSYS软件对上部结构进行建模分析,并通过FORTRAN程序对ANSYS软件的调用,实现了土与结构系统在地震作用下的整体分析。该方法为显隐式相结合的方法,地基和上部结构可采用不同的时间步距进行分析,可大大提高效率。通过两算例,验证了该方法的可行性。     

三维高密度电法在底板水探测中应用

针对目前采用二维高密度探测矿井底板水, 只能获得巷道底板的富水性特征, 而不能获得工作面内部的底板富水性情况的现状,研制出了通过在工作面周边铺设电缆,采集工作面底板三维数据体的设备及数据处理软件,并结合三维成像技术,实现了工作面底板富水性采前三维探测、三维数据体垂直及水平切割技术.结合新汶煤田协庄煤矿71504W工作面应用实例,介绍了三维高密度电法技术在底板水探测中的应用, 实现了对工作面底板富水性采前三维数据体的采集、采场底板含水性三维电阻率成像, 可方便地进行水平及垂直切片, 将工作面底板富水状态直观地表现出来,说明了该技术的实用性和有效性.     

Impedance functions for three-dimensional foundations supported on an infinitely-long canyon of uniform cross-section in a homogeneous half-space

A direct boundary element procedure is presented to determine the impedance matrix for a three-dimensional foundation supported on an infinitely-long canyon of uniform cross-section cut in a homogeneous half-space. The uniform cross-section of the canyon permits analytical integration along the canyon axis leading to a series of two-dimensional boundary problems involving Fourier transforms of the full-space Green's functions. Solution of these two-dimensional boundary problems leads to a dynamic flexibility influence matrix which is inverted to determine the impedance matrix. The accuracy of the procedure is demonstrated by comparison with previous solutions for a surface-supported, square foundation and results obtained by a three-dimensional boundary element method (BEM) for a foundation of finite-width supported on an infinitely-long canyon. Compared with the three-dimensional BEM, the present method requires less computer storage and is more accurate and efficient. The foundation impedance matrix determined by this procedure can be incorporated in a substructure method for earthquake analysis of arch dams.     

大跨度空间网格结构多维多点随机地震反应分析

本文建立了三维正交地震动多点激励下大跨度空间网格结构的随机地震反应分析方法,依据现行抗震设计规范的有关规定,确定了平稳随机地震动功率谱密度的模型参数。数值仿真分析了一柱距80m的正方形平板网架分别在一维地震动或三维地震动的一致激励、行波激励和考虑部分相干效应的随机激励下的地震反应。结果表明:考虑地震动的空间效应会很大程度地改变结构杆件的内力,其中控制杆件的内力增幅达到30%;地震动的行波效应对结构杆件内力的影响比随机地震动的部分相干效应的影响更大;三维地震作用比一维地震作用下结构杆件的内力大。由此得出结论,对于大跨度空间网格结构,必须进行多维多点地震激励下的随机地震反应分析。     

三维Radon变换的一种快速解析方法

3D Radon变换及其反变换是X-CT三维图像重建理论的核心,它在其他许多学科领域也有广泛应用。3D Radon变换的表达式是一个三重积分,按照定义直接计算相当费时。为此,研究一种新的快速的方法实现3D Radon变换,对X-CT图像重建理论及相关领域的发展有重要意义。本文以算法仿真常用椭球模型为基础,通过求解椭球模型与空间任意平面的面积,实现了用解析的方法快速得到模型的Radon变换,进一步比较了它与传统方法的优缺点,最后根据Radon反变换重建出原物体模型;计算机仿真结果验证了这种方法的正确。     

THE DELINEATION OF THREE-DIMENSIONAL SHALLOW GEOMETRY OF ACTIVE FAULT BASED ON TLS AND GPR: A CASE STUDY OF AN NORMAL FAULT ON THE NORTH MARGIN OF MAOYABA BASIN IN LITANG,WESTERN SICHUAN PROVINCE

It is crucial to reveal the surface traces and activity of active faults by obtaining high-precision microtopography and three-dimensional shallow geometry. However, limited by the traditional geological investigation methods in the field and geological condition factors, the measurement method on microtopography and shallow geometry of active fault is badly insufficient. In this study, the TLS and GPR are firstly used comprehensively to delineate the microtopography and shallow geometry of the normal fault scarp on the north margin of Maoyaba Basin in Litang. Firstly, the vertical displacements of two landforms produced by the latest two periods of normal faulting and the two-dimensional GPR profiles are obtained separately. Secondly, the three-dimensional measurement method of active fault based on TLS and GPR is preliminarily established. On this basis, three-dimensional model of fault scarp and three-dimensional images of subsurface geometry are also obtained. These data all reveal a graben structure at normal fault scarps. Thirdly, the fusion and interpretation of three-dimensional data from the surface and subsurface are realized. The study results show:1)the vertical displacements of the T₁ and T₂ terraces by the normal fault movement is 1.4m and 5.7m, the GPR profile shows a typical fault structure and indicates the existence of small graben structure with a maximum width of about 40m in the shallow layer, which further proves that it is a normal fault. 2)the shallow geometry of the normal fault scarp can be more graphically displayed by the three-dimensional radar images, and it also makes the geometry structure of the fault more comprehensive. The precise location and strike of faults F₁ and F₂ on the horizontal surface are also determined in the three-dimensional radar images, which further proves the existence of small graben structure, indicating the extensional deformation characteristics in the subsurface of the fault scarps. Furthermore, the distribution of small graben structure on the surface and subsurface is defined more precisely. 3)the integrated display of microgeomorphology and shallow geometry of normal fault scarp is realized based on the three-dimensional point cloud and GPR data. The fusion of the point cloud and GPR data has obvious advantages, for the spatial structure, morphological and spectral features from the point cloud can improve the recognition and interpretation accuracy of GPR images. The interpreted results of the GPR profiles could minimize the transformation of the surface topography by the external environment at the most extent, restore the original geomorphology, relocate the position and trend of faults on the surface and constrain the width of deformation zones under the surface, the geological structure, and the fault dislocation, etc. In a word, the TLS and GPR can quickly and efficiently provide the spatial data with multi-level and multi-visual for non-destructive inspection of the microgeomorphology and shallow structure for the active fault in a wide range, and for the detection of active fault in the complex geological environments, and it is helpful to improve the accuracy and understanding of the investigation and research on microtopography and shallow geometry of active faults. What's more, it also offers important data and method for more comprehensive identification and understanding of the distribution, deformation features, the behaviors of active faults and multi-period paleoseismicity. Therefore, to continuously explore and improve this method will significantly enhance and expand the practicability and application prospects of the method in the quantitative and elaborate studies of active faults.     

基于BFGS法融合InSAR和GPS技术监测地表三维形变

虽然InSAR技术具有高精度、大范围和高空间分辨率等优点,但只能监测雷达视线方向上的一维地表形变;而GPS技术虽可以监测地表的三维形变,但其空间分辨率很低.本文针对融合InSAR和GPS技术监测地表高空间分辨率三维形变展开研究.首先证明了简单的局部最优化迭代算法就能求得综合InSAR和GPS监测地表形变速率的能量函数模型的全局最优估值.随后提出了利用BFGS局部最优算法反演最优的地表三维形变速率.该方法既能避免全局最优化算法计算复杂且难以收敛的问题,又能克服传统的解析法中数值计算不稳定的缺点.最后,通过模拟实验和美国南加州真实数据实验表明,该方法能够得到高精度的地表三维形变速率场.而且当观测或插值误差导致解析法误差较大时,BFGS方法仍能得到高精度、稳定的全局最优解.