Large strain similarity solution for a spherical or circular opening excavated in elastic‐perfectly plastic media

This paper deals with the unloading problem of a spherical or circular opening excavated in elastic‐perfectly plastic media with a nonassociated Mohr–Coulomb yield criterion. A large strain similarity solution, using incremental velocity approach, is presented by replacing partial differential equations from stress equilibrium, constitutive law, consistency condition, and displacement equation with first‐order ordinary differential equations. The classical Runge–Kutta method is used to solve the first‐order ordinary differential equations. Comparisons among small and large strain solutions are made using some data sets of soil and rock. The results show that the displacements by large strain similarity solution are smaller than those by exact small strain solution and somewhat larger than those by large strain solution using total strain approach. Copyright © 2014 John Wiley & Sons, Ltd.     

Free vibrations of simply supported double plate on two models of elastic soils

In this paper the free vibrations of simply supported rectangular plates, resting on two different models of soils, are considered. The first model called Hetényi, by the name of the deviser, assumes a continuous plate, embedded in a Winkler‐type soil to realize the foundation partial continuity. The upper plate rests on a multiple layer characterized by a Winkler‐type soil, a continuous plate and another Winkler‐type soil in sequence. The two Winkler‐type soils have different modules. In the second model, the plate will be embedded in two layers of soils, whose behavior is similar to that of the Pasternak–Kerr‐type soil. The two models have been already used for the study of the double beam system. The free motions, in both cases, are described by a homogeneous set of partial differential equations, based on Kirchhoff–Love theory. Next, the homogeneous equations of motion are solved by using the classical Navier method. The free frequencies and associated vibration mode shapes of double plate system are found and numerical examples are illustrated to compare the two models. Copyright © 2008 John Wiley & Sons, Ltd.     

Proterozoic tectonothermal processes imaged with magnetotellurics and seismic reflection in southern Australia

Over the last two decades,co-located seismic and magnetotelluric(MT) profiles provided fundamental geophysical data sets to image the Australian crust.Despite their complimentary nature,the data are processed and often interpreted separately without common processes in mind.We here qualitatively compare 2 D resistivity inversion models derived from MT and seismic reflection profiles across a region of Archean-Proterozoic Australia to address the causes of variations in seismic response and anomalous conductivity in the crust.We find that there exists a spatial association between regions of low reflectivity in seismic sections and low resistivity in co-located2 D MT modelled sections.These relationships elucidate possible signatures of past magmatic and fluid-related events.Depending on their diffuse or discrete character,we hypothesize these signatures signify fossil melting of the crust due to mafic underplating,magma movement or hydrothermal fluid flow through the crust.The approach discussed herein is a process-oriented approach to interpretation of geophysical images and a significant extension to traditional geophysical methods which are primarily sensitive to a singular bulk rock property or state.     

计算地球流体力学的回顾、进展及展望

简要回顾了计算地球流体力学的发展历史,概括介绍了计算地球流体力学的研究进展及最新发展方向。针对线性发展方程的计算稳定性问题,介绍了判定线性发展方程初值问题的稳定性判别条件：CLF条件。对分析线性偏微分方程差分格式计算稳定性判据的Fourier方法、启发性稳定性分析方法,也做了简要的介绍。同时,重点介绍了判定线性发展方程初边值问题计算稳定性的GKS理论。在非线性发展方程的计算稳定性方面,重点介绍的主要内容包括：非线性发展方程的计算紊乱现象和计算不稳定的原因;克服非线性发展方程计算不稳定的方法;瞬时平方守恒型差分格式的构造;隐式和显式完全平方守恒格式的设计;强迫耗散非线性发展方程的计算稳定性问题。在计算地球流体力学的近期进展方面,重点介绍了非线性发展方程的计算稳定性与初值的关系,强迫耗散非线性发展方程显式准完全平方守恒格式的构造。对计算地球流体力学需要进一步研究的问题,也做了简要介绍。这些研究工作的介绍,无疑对推动计算地球流体力学的研究和大气海洋模式的研制具有一定的指导意义。     

Tephra Fallout Models: The Effect of Different Source Shapes on Isomass Maps

Numerous tephra dispersion and sedimentation models rely on some abstraction of the volcanic plume to simplify forecasts of tephra accumulation as a function of the distance from the volcano. Here we present solutions to the commonly used advection–dispersion equation using a variety of source shapes: a point, horizontal and vertical lines, and a circular disk. These may be related to some volcanic plume structure, such as a strong plume (vertical line), umbrella cloud (circular disk), or co-ignimbrite plume (horizontal line), or can be used to build a more complex plume structure such as a series of circular disks to represent a buoyant weak plume. Basing parameters upon eruption data, we find that depositions for the horizontal source shapes are very similar but differ from the vertical line source deposition. We also compare the deposition from a series of stacked circular disk sources of increasing radius above the volcanic vent with that from a vertical line source.     

Solving the Vialov Equation of Glaciology in Terms of Elementary Functions

Very few exact solutions are known for the non-linear Vialov ordinary differential equation describing the longitudinal profiles of alpine glaciers and ice caps under the assumption that the ice deforms according to Glen’s constitutive relationship. Using a simple, yet wide, class of models for the accumulation rate of ice and Chebysev’s theorem on the integration of binomial differentials, many new exact solutions of the Vialov equations are obtained in terms of elementary functions.     

大比例尺航空物探在勘查中的应用——以加拿大新斯科舍省悉尼铜钼多金属矿为例

简要介绍了加拿大航空物探发展的现况,并以新斯科舍省悉尼铜钼多金属矿为例,介绍了大比例尺航空物探在勘查中的应用及取得的成果。通过收集矿区区域航磁以及部份地段的地面高磁,以资料的处理、融合、分析为手段,以地质为依据,利用已经掌握的地质规律,建立测区内可能有的几种地质模型及地球物理模型,划分出矿区区域地质构造、岩体,并圈定了有价值成矿的靶区,为地质勘查和找矿提供了帮助。     

EH4物探找矿技术在非洲赞比亚的应用示范

非洲已成为我国企业海外铜矿资源开发的投资主要目的地之一,针对非洲赞比亚铜矿山的深边部地质找矿技术问题开展的EH4物探找矿技术试验性应用研究,具有重要的示范效应.使用EH4连续电导率剖面测量仪完成的物探剖面测量所获取的某铜矿西矿区14线和东南矿区59线物探电性断面图,较好地反映了加丹加沉积盆地基底构造层与盆地沉积构造层之间的角度不整合沉积接触界面的空间形态展布,而且在含铜页岩层底板与基底构造层顶板之间的石英岩层与高阻异常具有较好的对应关系.这些地球物理信息可为赞比亚沉积岩型铜矿床的深部找矿提供十分重要的间接找矿线索.     

The Problem of Viscoelastic Relaxation of the Earth Solved by a Matrix Eigenvalue Approach Based on Discretization in Grid Space

Viscoelasticity is a geophysical process which operates over intermediate timescales between a few years to millions of years, depending on the ambient thermal conditions of the self-gravitating spherical planet. Topography undulations with time represent geological signatures to both the internal and external loading processes, such as post-glacial rebound, volcanic eruptions, sedimentary loading, meteoritic impacts and mountain building. The span of relaxation timescales or relaxation spectrum in a viscoelastic spherical body has traditionally been determined by employing the Laplacian transform method and the correspondence principle relating the elastic solution to its viscoelastic counterpart.We have devised a novel approach based on the method of lines in which the equations at each angular order in the spherical harmonic expansion are discretized in the radial co-ordinate. The finite-dimensional space spanned by the discretized points in the radial direction of the planetary model then forms the basis of a matrix Eigenvalue problem. The Eigenvalues can be computed very fast because of the availability of public domain software. We can, for instance, compute the entire range of viscoelastic relaxation with computational times from 0.05 to 50 sec using only 30 to 300 radial grid points. The models can have both realistic density and elastic parameter profiles, derived from seismology. We show results here for complicated viscosity profiles with an asthenosphere in the upper mantle and a viscosity hill in the middle portion of the lower mantle. Because of the rapidity of the code, we may use this new method for exploring non-linear inversion problems by parameter sweeping.     

Bayesian Gaussian Mixture Linear Inversion for Geophysical Inverse Problems

A Bayesian linear inversion methodology based on Gaussian mixture models and its application to geophysical inverse problems are presented in this paper. The proposed inverse method is based on a Bayesian approach under the assumptions of a Gaussian mixture random field for the prior model and a Gaussian linear likelihood function. The model for the latent discrete variable is defined to be a stationary first-order Markov chain. In this approach, a recursive exact solution to an approximation of the posterior distribution of the inverse problem is proposed. A Markov chain Monte Carlo algorithm can be used to efficiently simulate realizations from the correct posterior model. Two inversion studies based on real well log data are presented, and the main results are the posterior distributions of the reservoir properties of interest, the corresponding predictions and prediction intervals, and a set of conditional realizations. The first application is a seismic inversion study for the prediction of lithological facies, P- and S-impedance, where an improvement of 30% in the root-mean-square error of the predictions compared to the traditional Gaussian inversion is obtained. The second application is a rock physics inversion study for the prediction of lithological facies, porosity, and clay volume, where predictions slightly improve compared to the Gaussian inversion approach.     

粘弹性波动方程的系数反演

波动方程的系数反演是一种用于识别地下介质物理力学参数的重要方法。过去这方面的研究一般都建立在弹性模型的基础上,本文则提出了粘弹性反演模型,这种模型更能真实地反映波在地下传播的实际情形。在粘弹性反演模型的基础上,我们还提出了一整套行之有效的数值反演方法,并在频域中完成了反演计算。最后的数值计算结果是令人满意的,证明了本文的模型和方法的合理性。     

抗拔桩轴力分布函数和位移分布函数间的关系模型

为了从荷载传递微分方程导出具有实际意义的非线性解析解,建立一种简单的沿桩长轴力分布函数和位移分布函数之间的关系模型U(z)-?(z)是非常重要的。采用建立经验公式类数学模型的方法,从几何作图法获得的U(z)-?(z)曲线的形状入手,通过对若干个数学模型的试算,确定了指数形式的模型,并对模型中的参数意义和确定进行了讨论,认为参数?,b只与桩顶和桩底位移有关。     

凝聚态地质问题

在地质学背景后面隐藏着众多凝聚态问题。对它们需要运用凝聚态科学工具进行探索 ,才能发现蕴含其中的对称、相变、周期和突变等现象 ,以抓住其共有的科学本质。文中对结晶过程中的固熔体相变 ,上地幔部分熔融和地核差异旋转等典型地学问题 ,从凝聚态角度重新作了分析 ,指出尽管它们具有不同的时空范畴和地学内涵 ,曾分别被划归于不同的地学分支学科 ,但却共同体现出了大自然中模式形成的一个基本范式 ,即 :无序定态的空间对称破缺可以导致诸多宏观有序结构的形成。这也是其它一些类似的地学现象得以产生的共同原因。由此可认为 ,从各种地质分支学科的单一性探索 ,过渡到凝聚态地质学的统一研究 ,是地质科学发展过程中的一个必由阶段。而凝聚态物理学的成熟与深化 ,恰好为之提供了发展的契机。因此 ,就有可能从凝聚态地质学这一新的观点出发 ,将传统的地质研究领域进一步拓深 ,并且将其拓宽、推广到水利 ,环境 ,材料和生命科学中去。     

Mortar coupling and upscaling of pore-scale models

Pore-scale models are becoming increasingly useful as predictive tools for modeling flow and transport in porous media. These models can accurately represent the 3D pore-structure of real media. Currently first-principles modeling methods are being employed for obtaining qualitative and quantitative behavior. Generally, artificial, simple boundary conditions are imposed on a model that is used as a stand-alone tool for extracting macroscopic parameters. However, realistic boundary conditions, reflecting flow and transport in surrounding media, may be necessary for behavior that occurs over larger length scales or including pore-scale models in a multiscale setting. Here, pore-scale network models are coupled to adjacent media (additional pore-scale or continuum-scale models) using mortars. Mortars are 2D finite-element spaces employed to couple independent subdomains by enforcing continuity of pressure and flux at shared boundary interfaces. While mortars have been used in the past to couple subdomains of different models, physics, and meshes, they are extended here for the first time to pore-scale models. The approach is demonstrated by modeling single-phase flow in coupled pore-scale models, but the methodology can be utilized to model dynamic processes and perform multiscale modeling in 3D continuum simulators for flow and transport.     

Geophysical evidence for a large impact structure on the Falkland (Malvinas) Plateau

A large, roughly circular structural basin is recognised on the Falkland (Malvinas) Plateau to the NW of West Falkland (Gran Malvina) Island (S 51°00′, W 62°00′). The basin, seen in seismic‐reflection profiles and evident as a large negative gravity anomaly, has a diameter of ~250 km. The age of the basin is estimated to be Late Palaeozoic. It is completely buried by younger sediments and has no topographic expression on the sea floor. We propose that the basin and geophysical anomalies, especially the combination of a large circular negative gravity anomaly with a rim of positive anomalies, and a marked circular series of positive magnetic anomalies in the same area, may be best explained by the presence of a large buried impact structure.     

基于水合物指示因子的地震识别方法

速度是水合物勘查的重要参数,但是速度的影响因素众多,给高精度识别水合物带来一定困难.水合物储层岩石物理是研究多种因素对速度影响的有效手段之一.考虑到水合物储层未固结等特征,本文采用SCA-DEM模型,重点对比了孔隙度、饱和度等对纵横波速度以及AVO特征的影响,并构建了一种新的水合物指示因子,消除了由于孔隙变化带来的假异...     

Comparing electromagnetic and seismic geophysical methods: Estimating the depth to water in geologically simple and complex arid environments

An integrated electromagnetic (EM) and seismic geophysical study was performed to evaluate non-invasive approaches to estimate depth to shallow groundwater (i.e., < 5 m) in arid environments with elevated soil salinity, where the installation of piezometers would be limited or prohibited. Both methods were tested in two study areas, one serving as a control site with relatively simple hydrogeology and the other serving as the experimental site with complex hydrogeology. The control site is located near the shore of Utah Lake (Palmyra, Utah, USA) where groundwater is shallow and unconfined in relatively homogeneous lacustrine sediments. The experimental site is in Carson Slough, Nevada, USA near the Ash Meadows National Wildlife Refuge in Amargosa Valley. Carson Slough is underlain by valley fill, with variable shallow depths to water beneath an ephemeral braided stream system. The geophysical methods used include frequency domain electromagnetic induction with multiple antenna–receiver spacings. High-resolution P-wave seismic profiles using a short (0.305 m) geophone spacing for common depth-point reflection stacking and first arrival modeling were also acquired. Both methods were deployed over several profiles where shallow piezometer control was present. EM results at both sites show that water surfaces correspond with a drop in conductivity. This is due to elevated concentrations of evaporative salts in the vadose zone immediately above the water table. EM and seismic profiles at the Palmyra site accurately detected the depth to groundwater in monitoring wells, as well as interpolated depths between them. This demonstrates that an integrated approach is ideal for relatively homogeneous aquifers. On the other hand, interpreting the EM and seismic profiles at Carson Slough was challenging due to the laterally and vertically variable soil types, segmented perched water surfaces, and strong salinity variations. The high-resolution images and models provided by the geophysical profiles confirm the simple soil and hydrological structure at the Palmyra site as well as the laterally complex structure at Carson Slough. The integrated approach worked well for determining depth to water in the geologically simple site, but was less effective in the geologically complex site where multiple water tables appear to be present.     

上地幔顶部电导率异常的起因

地球物理深部探测研究表明,上地幔顶部个别地区具有显著的电导率异常,表现为高电导率和各向异性。造成上地幔顶部电导率异常的机理,一直以来是地球科学中极具争议的问题之一。目前比较流行的解释是由部分熔融产生的熔体和橄榄石中耦合的结构水所致,其它的解释包括由颗粒边界的石墨或硫化物以及其它一些导电性强的矿物所致。这些不同的模型,对于认识上地幔的结构、组成和性质有直接的影响。综合评述了近年来基于岩石学、地球化学、地球物理学、数值模拟和高温高压实验研究等方面的进展,对已有模型伴随的各种问题进行了探讨。局部地区的电导率异常可能主要由固态导电机制造成,但其主要载体可能不是构成地幔主体的橄榄岩,而很可能与上地幔岩石学组成上的宏观不均一性有关。不同地区上地幔顶部的电导率异常,可能由不同的因素造成,可能并不存在所谓的"单一"导电机制。