电性任意各向异性且分块连续变化CSAMT三维有限元数值模拟

考虑地球介质电导率任意各向异性且随空间位置连续变化的情况,本文实现了直接求解电磁场的可控源音频大地电磁测深（CSAMT）三维有限元数值模拟.首先给出了电导率任意各向异性介质中CSAMT二次电场满足的控制方程及其相应变分问题,然后采用任意六面体单元对研究区域进行剖分,在网格单元中对任意各向异性电导率进行线性插值,解决了实际工作中岩矿石电导率各向异性且连续变化的情况,将变分问题转化为线性代数方程组的求解.电导率各向异性且连续变化一维模型三维有限元数值模拟结果与电导率各向异性且分层均匀渐进模型解析解结果对比验证了方法的有效性;三维地电模型电导率随位置线性变化且各向同性、主轴各向异性、方位各向异性和倾斜各向异性的数值模拟结果表明,电导率各向异性且连续变化对CSAMT视电阻率和相位数据均有明显的影响.     

Calculation of magnetic anomalies caused by 2D bodies of arbitrary shape with consideration of demagnetization

Forward calculations of magnetic anomalies caused by two-dimensional bodies of any shape and magnetic properties may be performed either without considering demagnetization as in the equivalent source technique or taking demagnetization into account as in the volume integral equation (VIE) approach, in which, for this purpose, magnetized bodies are divided into a set of rectangular prismatic cells. Ignoring demagnetization may result in distortion of the shape and the amplitude of an anomaly, whereas rectangular cells may not be an optimal representation of the source. Moreover, an inaccurate form approximation in the VIE technique may lead to inconsistent results in the near-body region. In this paper, a method is proposed, based on the VIE approach but differing by applying triangular elementary cells. The method largely overcomes the above-mentioned limitations of the VIE technique. It allows us to delineate large and complex structures exactly and only requires the source to be divided into a few elementary cells to take demagnetization into account satisfactorily. These improvements have been attained through analytical calculation of the Green's function in the complex plane, using the theory of the Cauchy-type integral. Comparing numerical solutions with analytical solutions for homogeneous elliptic cylinders without remanence, the method is found to be consistent with the theory in the range of relative magnetic permeability of 2–20, not only far from but also at subcell distances from the body. The method is appropriate for modelling highly and inhomogeneously magnetized 2D bodies of any shape. It may be of value in interpreting underground measurements or topographic effects, as well as in modelling regional geomagnetic profiles, and it is also a convenient tool for testing questionable geological hypotheses. In the framework of the method, the gravitational anomaly for the same causative bodies can be easily calculated. However, at higher and geologically uncommon values of relative magnetic permeability, the algorithm may become unstable but may be stabilized with SVD regularization. The fact that discrepancies were found with the method employed is a basis for further research.     

电导率任意各向异性海洋可控源电磁三维矢量有限元数值模拟

现有海洋可控源电磁三维数值模拟方法大多基于电导率各向同性介质理论,不能模拟海底地层电导率各向异性的实际情况.本文给出了电导率各向异性三维介质中电性源海洋可控源电磁二次电场的边值问题以及相应的变分问题,采用长方体单元对研究区域剖分,将场分量定义在剖分单元的边上,利用矢量有限单元法求解变分问题,实现了电导率任意各向异性海洋可控源电磁三维矢量有限元数值模拟.这个新的正演方法可以计算电导率任意各向异性三维地电模型的海洋可控源电磁响应,基于二次场矢量有限元法直接求解电磁场,避免了传统有限元方法可能遇到的伪解问题和难于处理电场法向分量不连续的问题,提高了数值模拟计算精度.一维电导率各向异性模型电磁场数值解与解析解吻合得相当好,无论在源附近还是远离源处相对误差均不超过1%.电导率各向异性二维模型的计算结果与已有文献采用的非结构有限元模拟结果十分吻合.三维地电模型数值模拟结果显示,电导率各向异性张量电导率主轴分量和欧拉角对不同装置海洋可控源电磁响应均有着明显的影响.     

A hybrid optimization scheme for self-potential measurements due to multiple sheet-like bodies in arbitrary 2D resistivity distributions

Self-potential is a passive geophysical method that can be applied in a straightforward manner with minimum requirements in the field. Nonetheless, interpretation of self-potential data is particularly challenging due to the inherited non-uniqueness present in all potential methods. Incorporating information regarding the target of interest can facilitate interpretation and increase the reliability of the final output. In the current paper, a novel method for detecting multiple sheet-like targets is presented. A numerical framework is initially described that simulates sheet-like bodies in an arbitrary 2D resistivity distribution. A scattered field formulation based on finite differences is employed that allows the edges of the sheet to be independent of the grid geometry. A novel analytical solution for two-layered models is derived and subsequently used to validate the accuracy of the proposed numerical scheme. Lastly, a hybrid optimization is proposed that couples linear least-squares with particle-swarm optimization in order to effectively locate the edges of multiple sheet-like bodies. Through numerical and real data, it is proven that the hybrid optimization overcomes local minimal that occurs in complex resistivity distributions and converges substantially faster compared to traditional particle-swarm optimization.     

任意地球物理模型的三角形和四面体有限单元剖分

系统介绍了适用于复杂地球物理模型的三角形和四面体单元剖分算法,讨论了单元优良性,优良的网格剖分判定规则,建立了C++执行类库与地球物理模型单元剖分库.利用C++剖分类库与模型单元库,可以进行任意地形条件下的地球物理场的正演与反演计算,剖分类的高效性为实际模型的并行有限元法计算提供了基础.     

钢框架结构的非线性静力抗震可靠性分析:Ⅰ有限元反应及其灵敏度

作为一种随机有限元方法,有限元可靠度方法通过有限元反应灵敏度分析将结构可靠度分析的近似解析方法与结构确定性分析的有限元方法结合起来,可以有效地处理结构反应是基本随机变量的隐式函数这一难题,因此成为大型复杂结构可靠度分析的有效工具。采用基于位移的非线性纤维梁柱单元对钢框架结构进行有限元建模,推导了单元与截面的基本方程。针对基于位移的纤维梁柱单元,采用逐级递进方式,分别推导了整体级、单元级、截面级和材料级的有限元反应灵敏度直接微分表达式。研究结果可为采用FORM或SORM等近似解析法进行钢框架结构的非线性静力有限元可靠性分析提供算法支持和编程依据。     

Improved tilt-depth method for fast estimation of top and bottom depths of magnetic bodies

The tilt-depth method can be used to make fast estimation of the top depth of magnetic bodies. However, it is unable to estimate bottom depths and its every inversion point only has a single solution. In order to resolve such weaknesses, this paper presents an improved tilt-depth method based on the magnetic anomaly expression of vertical contact with a finite depth extent, which can simultaneously estimate top and bottom depths of magnetic bodies. In addition, multiple characteristic points are selected on the tilt angle map for joint computation to improve reliability of inversion solutions. Two- and threedimensional model tests show that this improved tilt-depth method is effective in inverting buried depths of top and bottom bodies, and has a higher inversion precision for top depths than the conventional method. The improved method is then used to process aeromagnetic data over the Changling Fault Depression in the Songliao Basin, and inversion results of top depths are found to be more accurate for actual top depths of volcanic rocks in two nearby drilled wells than those using the conventional tilt-depth method.     

Seismic analysis of base-isolated liquid storage tanks using the BE–FE–BE coupling technique

A three-dimensional soil–structure–liquid interaction problem is numerically simulated in order to analyze the dynamic behavior of a base-isolated liquid storage tank subjected to seismic ground motion. A dynamic analysis of a liquid storage tank is carried out using a hybrid formulation, which combines the finite shell elements for structures and the boundary elements for liquid and soil. The system is composed of three parts: the liquid–structure interaction part, the soil–foundation interaction part, and the base-isolation part. In the liquid–structure interaction part, the tank structure is modeled using the finite elements and the liquid is modeled using the internal boundary elements, which satisfy the free surface boundary condition. In the soil–foundation interaction part, the foundation is modeled using the finite elements and the half-space soil media are modeled using the external boundary elements, which satisfy the radiation condition in the infinite domain. Finally, above two parts are connected with the base-isolation system to solve the system's behavior. Numerical examples are presented to demonstrate the accuracy of the developed method, and an earthquake response analysis is carried out to demonstrate the applicability of the developed technique. The properties of a real LNG tank located in the west coast of Korea are used. The effects of the ground and the base-isolation system on the behavior of the tank are analyzed.     

A comparison of continuous mass‐lumped finite elements with finite differences for 3‐D wave propagation

The finite‐difference method on rectangular meshes is widely used for time‐domain modelling of the wave equation. It is relatively easy to implement high‐order spatial discretization schemes and parallelization. Also, the method is computationally efficient. However, the use of finite elements on tetrahedral unstructured meshes is more accurate in complex geometries near sharp interfaces. We compared the standard eighth‐order finite‐difference method to fourth‐order continuous mass‐lumped finite elements in terms of accuracy and computational cost. The results show that, for simple models like a cube with constant density and velocity, the finite‐difference method outperforms the finite‐element method by at least an order of magnitude. Outside the application area of rectangular meshes, i.e., for a model with interior complexity and topography well described by tetrahedra, however, finite‐element methods are about two orders of magnitude faster than finite‐difference methods, for a given accuracy.     

关于用有限元法作磁法勘探正演计算的理论问题

本文提出用有限元法解磁法勘探的正问题,不再需要引入均匀磁化假定。 传统的正演计算方法是在均匀磁化的假定下建立的,它不能适应解释复杂磁异常的需要。近几年来,国外学者提出了用迭代方法作非均匀磁化条件下的正演计算,对形状比较简单的磁性体取得了较好的效果。但随着磁性体的形状变得比较复杂,计算所需时间很快增加,精度和稳定性也明显降低。以求能量函数极小的变分原理为依据,用有限元法作非均匀磁化条件下的正演计算,得到的磁位势函数满足磁性体内外的偏微分方程和全部界面条件（文中对此作了证明）。因此,这种方法可以研究各种复杂因素的影响,宜于计算非均匀磁化磁性体的有效磁化强度和磁异常。     

Digital computation of gravitational and magnetic anomalies and their derivatives for two-dimensional bodies of arbitrary shape

Summary A simple method is presented for computing the gravity and magnetic profiles across two-dimensional bodies of arbitrary shape. The computer programme makes use of expressions for the attraction caused by an infinitely long sheet and it takes only a few seconds to compute a profile both for gravity and magnetic anomalies and their vertical and horizontal derivatives. Accuracy of the method is shown to be as good as that of graticules and owing to the use of simpler formulae, the method appears to be faster than the method of polygons.     

三维高密度电阻率E SCAN法有限元模拟异常特征研究

工程电法勘探中，常需要探测埋深不大、规模较小、分布复杂的目标体，对探测方法的精度要求较高.由于三维高密度电阻率法数据采集密度大，可对目标体进行多方位观测，是目前工程探测中常用的方法之一.本文对三维高密度电阻率E SCAN法的视电阻率异常特征进行研究.设计几种典型地质体模型，利用有限单元法进行正演计算，并分析其异常特征和分布规律.算例结果表明，E SCAN法对低阻体的分辨能力强于高阻体.无需对所有电极进行扫描供电观测，即可分辨地质体电性特征及水平位置，但观测结果难以分辨地质体纵向延伸.     

Numerical modelling of suspended sediments

The turbulent advection-diffusion mathematical model in three-dimensional space is solved by a mixed finite element finite difference method. Linear finite elements in the vertical direction and central finite differences in the horizontal directions are used coupled with the Galerkin error minimization procedure. The integration in time is performed in fractional steps (one explicit one implicit) by splitting the differential operator. The method is illustrated by application to the three-dimensional movement of suspended sediment. Its accuracy is checked by comparison to analytical solutions and its efficiency is gauged relative to finite elements and implicit finite difference solutions for two-dimensional suspended sediment transport over a dredged channel.     

Coupling of boundary and finite elements for soil-structure interaction problems

The investigation of complex soil-structure interaction problems is usually carried out with numerical solution procedures such as the finite element or the boundary element method. It must be noted, however, that the choice of one or the other of these approaches is not just a matter of preferences; depending on the type of the problem under consideration, either boundary or finite elements may be more advantageous. A considerable expansion in the computational power can be obtained, on the other hand, if one resorts to hybrid schemes which retain the main advantages of the two methods and eliminate their respective disadvantages. This paper presents results obtained with a boundary element-finite element coupling procedure, and discusses its applicability to some representative soil-structure interaction problems. The structures considered are elastic systems, such as foundations, tunnels and filled trenches (modelled by finite elements), which are coupled with homogeneous elastic halfspaces (modelled by boundary elements). The examples demonstrate the importance of using a model that includes wave radiation effects. The coupling approach is formulated entirely in the time domain so that an extension of the algorithm to non-linear analyses seems to present no further difficulties.     

Free vibration of frame shear wall structures on flexible foundations

The finite strip method is used to study the effect of an elastic foundation on the natural frequencies of coupled frame shear wall structures. The solid wall in the structure is divided into several strip elements, the column is treated as a line element and the effect of the connecting beams is dealt with through the compatibility matrices which transfer their structural properties to the adjacent strip or line elements. The comparison functions which satisfy the boundary conditions of being free at the top and being spring supported at the bottom are used for the displacement field in the longitudinal direction. A series of numerical examples is presented to show the accuracy and applicability of the proposed method.     

Symmetric modes of vibration for A circular plate using finite dynamic elements

The dynamic element method has been shown previously to provide a computational advantage over the ordinary finite element method for various beam elements. The Taylor expansions are computed here for the dynamic shape functions (two terms) and dynamic stiffness matrix (four terms) for the axisymmetric vibrations of an annular plate element. The complicated matrices which result are made more tractable by expressing them as power series in powers of the aspect ratio. The percentage error in the natural frequencies is then calculated using both the two- and the three-term dynamic stiffness matrix, demonstrating the increased accuracy for a given number of elements.     

Finite element-fractional steps solution of the 3-D coastal circulation model

The mathematical model for the nearly horizontal circulation due to wind, tides and density gradients in 3-D coastal areas is solved by a combined use of the method of finite elements and the integration in fractional steps. The discretisation of the flow domain is achieved through a system of 1-D finite elements over the depth, z, and 2-D finite elements in x?y space. The differential operators of the momentum equations in x and y, are split and integrated separately in z and x?y dimensions. The method is an extension of a previously presented approach combining finite differences and expansion in series. The application refers to the wind induced circulation in the 3-D coastal basin of Thessaloniki Bay.     

Rapid computation of magnetic anomalies and demagnetization effects caused by bodies of arbitrary shape

Summary With the aid of formulae derived for the field caused by a finite rectangular prism, a machine method is developed for rapid computation of magnetic anomalies due to a body of any shape. Conversely, the method could be utilised for determination of the magnetization vector from the observed anomalies. An example is given to demonstrate the applicability, accuracy and speed of the method. The method is shown to be suitable also for the evaluation of magnetic terrain effects.In addition, a new method is presented which enables high speed calculation of demagnetization effect caused by a body of arbitrary shape where the magnetization in general may be inhomogeneous. Applicability and accuracy of the method are discussed. The method is applied to study the effect of susceptibility on inhomogeneity of the induced magnetization in a cube specimen.For a body of any shape with low susceptibility, a new concept of average demagnetization factor is introduced. Average demagnetization factors so calculated for a cylindrical specimen are tabulated.

---

Zusammenfassung Unter Anwendung der Formeln für das Feld eines Quaders wird eine Methode entwickelt um die durch einen Körper beliebiger Form erzeugten magnetischen Anomalien mit Hilfe eines Rechenautomaten zu berechnen. Umgekehrt kann die Methode verwendet werden, um aus den beobachteten Anomalien den Magnetisierungsvektor zu bestimmen. An einem Beispiel wird die Anwendbarkeit, die Genauigkeit und die Schnelligkeit der Methode erläutert. Die Methode kann auch für die Berechnung von magnetischen Terraineffekten verwendet werden.Zudem wird eine neue Methode dargelegt, welche eine rasche Berechnung des Entmagnetisierungseffektes eines Körpers beliebiger Gestalt und inhomogener Magnetisierung ermöglicht. Anwendbarkeit und Genauigkeit der Methode werden besprochen. Die Methode wurde verwendet, um den Einfluss der Suszeptibilität auf die Inhomogenität der induzierten Magnetisierung in einem Würfel zu untersuchen.Für Körper beliebiger Form mit kleiner Suszeptibilität wls neuer Begriff der mittlere Entmagnetisierungsfaktor eingeführt. Mittlere Entmagnetisierungsfaktoren für eine zylindrische Probe werden tabelliert.

---

---

Part I of the dissertation Theretical study of the magnetic attraction due to rock bodies and experimental investigation of the stability of rock magnetism submitted to the Swiss Federal Institute of Technology (ETH), Zurich, for the degree of Doctor of Natural Sciences.     

Computation of shallow water waves with hybrid finite elements

The numerical computation of two-dimensional incompressible long-period shallow water waves using the method of finite elements is presented. It is shown that the set of equations having been solved in each time step is reduced to a third of the usual size by satisfying the equations of motion on the element level only (hybrid finite elements). The water levels in all nodes are the sole parameters computed from the set of equations. The velocities are subsequently calculated on the element level in each time step. Special emphasis is placed on reproducing eddies behind corners using an anisotropic turbulent exchange tensor.     

Love numbers of a generalized Maxwell sphere

By elementary methods, we study the Love numbers of a homogeneous, incompressible, self-gravitating sphere characterized by a generalized Maxwell rheology, whose mechanical analogue is represented by a finite or infinite system of classical Maxwell elements disposed in parallel. Analytical, previously unknown forms of the complex shear modulus for the generalized Maxwell body are found by algebraic manipulation, and studied in the particular case of systems of springs and dashpots whose strength follows a power-law distribution. It will be shown that the sphere is asymptotically stable for any choice of the mechanical parameters that define the generalized Maxwell body and analytical forms of the Love numbers are always available for generalized bodies composed by less than five classical Maxwell bodies. For the homogeneous sphere, real Laplace inversion methods based on the Post-Widder formula can be applied without performing a numerical discretization of the n-th derivative, which can be computed in a closed-form with the aid of the Faà di Bruno formula.