Shallow structure of the Somma–Vesuvius volcano from 3D inversion of gravity data

A gravity investigation was carried out in the Somma–Vesuvius complex area (Campania, Italy) based on a dataset recently enlarged with new measurements. These cover the volcanic top and fill some other important spatial gaps in previous surveys. Besides the new gravity map of the Vesuvius, we also present the results of a 3D inverse modelling, carried out by using constraints from deep well exploration and seismic reflection surveys. The resulting density model provides a complete reconstruction of the top of the carbonate basement. This is relevant mostly on the western side of the survey area, where no significant information was previously available. Other new information regards the Somma–Vesuvius structure. It consists of an annular volume of rocks around the volcanic vent and that extends down to the carbonate basement. It results to be denser with respect to the surrounding sedimentary cover of the Campanian Plain and to the material located just along the central axis of the volcanic structure. The coherence between these features and other geophysical evidences from previous studies, will be discussed together with the other results of this research.     

基于模型空间压缩的大地电磁三维反演研究

本文提出了一种基于模型空间压缩技术的大地电磁三维反演方法.该方法在传统大地电磁三维反演理论的基础上,通过小波变换将待反演的空间域模型参数映射到小波域进行反演,获得小波域更新模型后再通过小波逆变换得到空间域反演模型.由于小波变换具有压缩特性和多尺度分辨能力,本文反演方法可在一定程度上提高反演分辨率.为了提高反演效率,我们针对基于L1范数的模型约束求解不易收敛的反演问题,提出了一种基于模型粗糙度的简单有效的预条件处理技术.为验证本文算法的有效性,本文首先对经典的"棋盘"模型进行三维反演测试.反演结果表明本文算法的反演效率与传统方法相当,但对于深部异常体具有更好的分辨能力.最后,我们通过对实测数据反演进一步验证本文算法的有效性.     

3D resistivity inversion using 2D measurements of the electric field

Field and 'noisy' synthetic measurements of electric-field components have been inverted into 3D resistivities by smoothness-constrained inversion. Values of electrical field can incorporate changes in polarity of the measured potential differences seen when 2D electrode arrays are used with heterogeneous 'geology', without utilizing negative apparent resistivities or singular geometrical factors. Using both the X - and Y -components of the electric field as measurements resulted in faster convergence of the smoothness-constrained inversion compared with using one component alone. Geological structure and resistivity were reconstructed as well as, or better than, comparable published examples based on traditional measurement types. A 2D electrode grid (20 × 10), incorporating 12 current-source electrodes, was used for both the practical and numerical experiments; this resulted in 366 measurements being made for each current-electrode configuration. Consequently, when using this array for practical field surveys, 366 measurements could be acquired simultaneously, making the upper limit on the speed of acquisition an order of magnitude faster than a comparable conventional pole–dipole survey. Other practical advantages accrue from the closely spaced potential dipoles being insensitive to common-mode noise (e.g. telluric) and only 7% of the electrodes (i.e. those used as current sources) being susceptible to recently reported electrode charge-up effects.     

石墨-石英模型电导率的数值模拟研究及地球物理应用

地壳高导异常及各向异性成因是地学研究的热点之一,断裂带区域常呈现出高电导率特征,对于石墨富集的断裂带,石墨形态及分布对导电性的影响已经引起关注.本文利用有限元数值模拟方法,构建了石墨-石英三维模型,探讨了石墨含量、排列方式、孔隙度及温度对模型电导率的影响.结果表明,石墨连通前,温度的增加能使模型电导率显著增大,而在石墨连通后,石墨含量的增加、片状石墨沿导电方向的定向排列以及石墨直径与厚度比值的增大,均能导致模型电导率发生显著增加,孔隙度的增大则导致模型电导率减小.同时,这些参数的变化也会使模型电导率突变时的石墨含量阈值发生改变,当石墨排列方向改变时,这一阈值变化显著,对于石墨富集的一些断裂带区域,石墨定向排列可能是引起高导的一种成因.

     

Numerical simulation of divided-bar thermal conductivity measurements

Summary A three-dimensional numerical model of the divided-bar thermal conductivity measuring device that can include both solid and fluid sample portions has been developed. The model has been employed to investigate the effects of porosity in divided-bar measurements and the effects of the distribution of solid sample and contained fluids in porous material on the thermal conductivity calculated from the bar temperature measurements. The results show that the positions of the temperature sensors and thus the nature of the vertical column between the temperature measuring points is of prime importance, and that the sizes of the solid and fluid portions of the sample affect the sample conductivity estimates. The divided bar measures the thermal conductivity parallel to the axis of the bar. Thermal conductivity variations in the radial direction have little effect on the calculated conductivity. Whether or not the sides of the bar are insulated strongly affects the calculated conductivity.Presented at the International Meeting on Terrestrial Heat Flow and the Structure of Lithosphere, Bechyn Castle, Czech Republic, September 2 – 7, 1991.     

Numerical simulation of organic chemicals in a marine environment using a coupled 3D hydrodynamic and ecotoxicological model

For ecotoxicological risk assessment in a marine ecosystem, we constructed a coupled three-dimensional hydrodynamic and ecotoxicological model (EMT-3D), and applied it to Tokyo Bay. The model was calibrated with field data obtained in 2002. The results of sensitivity analysis for dissolved Bisphenol A showed that biodegradation rate was the most important factor for concentration change. Bioconcentration coefficient was the most important factor for Bisphenol A in phytoplankton. Therefore, the parameters must be carefully considered in the modeling. The mass balance results showed that standing stocks of Bisphenol A in water, in particulate organic carbon and in phytoplankton are 7.85 x 10(4), 1.78 x 10(2) and 3.44 x 10(-1) g, respectively. With respect to flux, biodegradation in the water column had the highest value of 1.06 x 10(3) g/day, and next were effluent to the open sea, partition to particulate organic carbon, and bioconcentration in phytoplankton.     

淀山湖风浪场的数值模拟

对淀山湖的风场和风浪场进行了观测,根据所观测的风场利用SWAN模型计算了湖区的风浪场.通过将计算的波高和观测值进行比较,说明计算结果的变化趋势和观测结果的变化趋势相吻合;观测结果的波动幅度较大,而计算结果比较平滑.总体说来,基于SWAN模型所计算的波高可信、可靠.在此基础上,利用SWAN模型较为系统地计算了不同水位和不...     

三维正交各向异性孔隙介质震电波场数值模拟研究

各向异性等介质物性参数会对震电波场的传播产生影响.为了定量研究介质物性参数对震电波场传播的影响, 本文首先研究了三维正交各向异性介质的震电耦合方程, 然后采用时域有限差分方法模拟了爆炸源在该介质中的震电响应.在数值模拟中, 首先计算地震波场, 然后把地震波场作为输入来计算伴随电场.为了验证数值算法的精度, 我们将数值结果与从Pride方程获得的解析解进行了对比, 结果表明数值解可以很好地近似解析解.我们设计了几种各向异性模型, 分析了介质物性参数对震电波场的影响, 结果表明爆炸源在各向异性孔隙介质中会产生四种地震波, 分别是快慢纵波和两个可分离的横波.产生的地震波都会在均匀各向异性孔隙介质中产生伴随电场.弯曲度对慢纵波的传播有影响, 当弯曲度为各向异性时, 慢纵波的波场快照为椭圆形, 当弯曲度为各向同性时, 慢纵波的波场快照为圆形.

     

Numerical simulation and development of data inversion in borehole ultrasonic imaging

The major function of underground nuclear waste disposal is to avoid the migration of radionucleides towards the biosphere during their active period. This function can be deteriorated by the EDZ (excavated damaged zone) around the excavation. The EDZ analysis is therefore crucial in the performance assessment of the storage. The paper deals with the determination of the EDZ around a nuclear waste storage cavity using borehole ultrasonic imaging (azimuthal tomography). Indeed, before processing experimental data obtained with this tool, it is necessary to establish that data can be satisfactorily inverted. This analysis is based on a method that is able to sound and image the rock mass velocity field. The velocity field is numerically simulated (3D geomechanical modeling) based on an assumption on the relationship between stress and velocity fields. In order to evaluate a radial velocity profile starting from inter-sensor distance and their corresponding traveltimes, different ray tracing algorithms are tested using synthetic data. These tests led to a simple and fast approach (implemented in a Mathematica package) to process a large quantity of data.     

Seismic traveltime inversion of 3D velocity model with triangulated interfaces

Seismic traveltime tomographic inversion has played an important role in detecting the internal structure of the solid earth. We use a set of blocks to approximate geologically complex media that cannot be well described by layered models or cells. The geological body is described as an aggregate of arbitrarily shaped blocks, which are separated by triangulated interfaces. We can describe the media as homogenous or heterogeneous in each block. We define the velocities at the given rectangle grid points for each block, and the heterogeneous velocities in each block can be calculated by a linear interpolation algorithm. The parameters of the velocity grid positions are independent of the model parameterization, which is advantageous in the joint inversion of the velocities and the node depths of an interface. We implement a segmentally iterative ray tracer to calculate traveltimes in the 3D heterogeneous block models. The damped least squares method is employed in seismic traveltime inversion, which includes the partial derivatives of traveltime with respect to the depths of nodes in the triangulated interfaces and velocities defined in rectangular grids. The numerical tests indicate that the node depths of a triangulated interface and homogeneous velocity distributions can be well inverted in a stratified model.     

Study on 3D simulation of wave fields in acoustic reflection image logging

The borehole acoustic reflection imaging logging is a newly developed acoustic logging method that has attracted many interests. These converted and reflected waves for imaging are usually mixed up with borehole guided waves and therefore difficult to be clearly identified. To improve the downhole tool design and develop more sophisticate data processing and interpretation algorithms,studies on precisely numerical modeling of the wave fields in the acoustic reflection imaging logging are neces-sary and critical. This paper developed a parallelized scheme of 3D finite difference (3DFD) with non-uniform staggered grid and PML absorbing boundary to simulate the acoustic wave fields in isotropic and anisotropic formations. Applications of this scheme to the typical cases of isotropic and anisot-ropic formations and comparison with the results from published analytical solutions have demon-strated the validation and efficiency of the scheme. Higher accuracy and lower computation cost (3.5 times faster than the conventional schemes) have been achieved with this scheme for modeling such a complex wave fields of 60 dB dynamic range with higher frequency (10 kHz). This simulating program provides a quantitative analytical means for studying acoustic reflection imaging tool and development of the data processing and interpretation methods.     

Determination of the regularization level of truncated singular-value decomposition inversion: The case of 1D inversion of MT data

The total mean‐square error (MSE) of the estimated model, defined as the sum of the standard model variance and the bias variance, is used to define the truncation level of the singular‐value decomposition to give a reasonable balance between model resolution and model variance. This balance is determined largely by the data and no further assumptions are necessary except that the bias terms are estimated sufficiently well. This principle has been tested on the 1D magnetotelluric inverse problem with special emphasis on high‐frequency radio magnetotelluric (RMT) data. Simulations clearly demonstrate that the method provides a good balance between resolution and variance. Starting from a homogeneous half‐space, the best solution is sought for a fixed set of singular values. The model variance is estimated from the sum of the inverse eigenvalues squared, up to a certain threshold, and the bias variance is estimated from the model projections on the remaining eigenvectors. By varying the threshold, the minimum of the MSE is found for an increasing number of fixed singular values until the number of active singular values becomes greater than or equal to the estimated number. As a side‐effect, the depth of penetration of a given set of measurements can be estimated very efficiently by simply noting at which depth the final model deviates little from the starting homogeneous half‐space model. A suite of synthetic data is inverted and an example of inversion of one site is shown to illustrate how the truncation is carried out as the non‐linear inversion process proceeds. A field example with a profile across a plume of contaminated groundwater in the Netherlands shows good agreement with the electrical resistivity obtained in a nearby borehole.     

模拟退火方法在三维速度模型地震波走时反演中的应用

采用块状建模以及三角形拼接的界面描述方式,并通过立方体速度网格线性插值获得块体内部的速度分布。正演过程中采用逐段迭代射线追踪方法计算三维复杂地质模型中的射线走时,并采用模拟退火方法进行了三维模型中的地震波走时反演研究。模型测试结果表明,使用的射线追踪和走时反演算法有效。     

Constraints in 3D gravity inversion

A three-dimensional (3D) inversion program is developed to interpret gravity data using a selection of constraints. This selection includes minimum distance, flatness, smoothness and compactness constraints, which can be combined using a Lagrangian formulation. A multigrid technique is also implemented to resolve separately large and short gravity wavelengths. The subsurface in the survey area is divided into rectangular prismatic blocks and the problem is solved by calculating the model parameters, i.e. the densities of each block. Weights are given to each block depending on depth, a priori information on density and the density range allowed for the region under investigation. The present computer code is tested on modelled data for a dipping dike and multiple bodies. Results combining different constraints and a weight depending on depth are shown for the dipping dike. The advantages and behaviour of each method are compared in the 3D reconstruction. Recovery of geometry (depth, size) and density distribution of the original model is dependent on the set of constraints used. From experimentation, the best combination of constraints for multiple bodies seems to be flatness and a minimum volume for the multiple bodies. The inversion method is tested on real gravity data from the Rouyn-Noranda (Quebec) mining camp. The 3D inversion model for the first 10 km is in agreement with the known major lithological contacts at the surface; it enables the determination of the geometry of plutons and intrusive rocks at depth.     

火山Mogi模型形变反演的虚拟观测法及方差分量估计

针对Mogi模型垂直位移与水平位移联合反演中的病态问题,改进火山形变总体最小二乘（Total Least Squares,TLS）联合反演的虚拟观测法,并使用方差分量估计（Variance Components Estimation,VCE）方法确定病态问题的正则化参数.将附有先验信息的参数作为观测方程,与垂直位移和水平位移的观测方程联合解算,推导了三类观测方程联合反演的求解公式及基于总体最小二乘方差分量估计确定正则化参数的表达式,给出了算法的迭代流程.通过算例实验,研究了总体最小二乘联合反演的虚拟观测法在火山Mogi模型形变反演中的应用;算例结果表明,三类数据的联合平差及方差分量估计方法可以确定权比因子并得到修正后的压力源参数,具有一定的实际参考价值.

     

三维复电阻率模型电磁场正演模拟研究

复电阻率法是上世纪70年代发展起来的一种激电分支方法,自提出伊始就受到了国内外学者的广泛关注.很多学者开展了复电阻率电磁场正演计算,但缺少对正演结果的深入分析.本文采用体积分方程法,进行了复电阻率模型的正演模拟,通过对正演结果分析,得出三维复电阻率体的电磁响应是由累积电荷,电磁感应和激电效应共同产生的,并总结了三种效应的电磁响应规律,为准确分析异常提供了依据.     

含湿孔隙岩石有效热导率的数值分析

本文采用有限元方法研究含湿孔隙岩石的有效热导率,即随机划分网格并指定材料性质,建立三维含湿孔隙岩石的有限元模型,模型的上下表面施加不同的温度,侧面绝热,计算出总热流,然后结合上下表面的温度梯度计算出岩石的有效热导率.考虑到单个随机模型不一定具有代表性,对给定的孔隙率和饱和度均生成了200种矿物、水、空气随机分布的岩石模型,进行Monte Carlo实验和统计分析,统计分析结果与前人实验结果吻合良好.数值分析结果表明,孔隙岩石的有效热导率与岩石的孔隙率、饱和度、固体矿物组分及孔隙的分布情况有关,数值计算的误差随着网格数目的增加而减小.此有限元方法可以用来估算岩石的有效热导率,在已知组分性质的多矿物岩石物性计算方面有广阔应用前景.     

Behavior of electric potential fields over random conductivity,layered earth models

Solutions are derived for the potential distributions over one-layer and two-layer random conductivity earth models for the direct current resistivity method. Monopole and Wenner arrays are considered. It is shown that the random potential due to a spatially stationary Gaussian random conductivity is also Gaussian but not spatially stationary. The ensemble and sample statistics of the random kernel and apparent resistivity are examined. Representative curves of these functions are presented to show the effects of a random conductivity profile on them. It is shown that the effect of the random nature of the conductivity profile on an apparent resistivity sounding curve is greatest for small electrode spacings.     

斜长石、辉石混合模型的电导率有限元数值计算研究

以岩石实验中矿物的几何形态及空间分布为建模依据,以实验条件及单矿物电导率的测量结果为约束条件,用有限元数值方法模拟了不同微观结构的斜长石、辉石混合物在施加电压后电势及电流的分布情况,并计算了混合模型在不同温度条件下的电导率.研究结果显示,数值模型网格数及矿物颗粒数的选取对电导率计算结果的精度有较大影响,在体导电情况下,模型电导率因矿物比例含量和排列结构而异.当斜长石及辉石随机分布时,随着辉石含量的增加,混合模型电导率在不同温度下均有所增加,且温度越高,增加幅度越大,电导率的有限元模拟计算结果接近于有效介质渗透理论模型,且位于并、串联模型之间以及HS模型的上、下边界范围内;在斜长石及辉石含量一定的情况下,各矿物的排列分布对电导率计算结果也有一定的影响,当矿物颗粒大小接近且分布均匀时,模型中电势沿电流传导方向变化较为均匀,模拟计算得出的电导率相对较高,当矿物颗粒大小差别较大及分布不均匀时,电势分布受到一定的扰动,电导率计算结果也较低.将混合模型电导率有限元计算结果与辉长岩、辉绿岩及玄武岩实验测量结果进行比较,显示这3种岩石样品电导率与温度变化关系的斜率均与混合模型计算结果的斜率相接近,表明这些岩石在所选温度段导电机制与斜长石、辉石混合模型相似,用斜长石、辉石混合模型的电导率研究玄武岩、辉长岩及辉绿岩的导电性具有适用性.将混合模型有限元计算结果与玄武岩、辉长岩、辉绿岩覆盖区地壳大地电磁实测结果对比,发现大地电磁电导率结果位于混合模型计算结果范围内,用斜长石、辉石混合模型模拟玄武岩、辉长岩等岩石地壳具有一定的可行性.

     

采用有限元等效体力方法数值分析火山区岩浆压力变形源—以长白山火山为例

火山区岩浆压力变形源的反演计算采用解析方法存在难以考虑地形的限制,采用传统有限元方法则存在网格依赖和计算量大的问题,反演过程中每一次正演由于岩浆房位置和大小变化都需要重新生成一次网格,耗费巨大的计算量和网格生成时间.为了克服上述问题,首次在长白山火山区使用"有限元等效体力"方法考虑地形影响反演地下岩浆压力变形源,计算岩浆应力扰动对周边断层稳定性的影响.在火山区地下压力变形源引起的地表形变计算中,地表地形影响不可忽略.埋深越浅,地表最大径向位移ur所在的位置越靠近岩浆囊中心.当坡度达到30°时,最大垂向位移uz所在位置不再位于岩浆囊正上方.椭球状岩浆囊压力源可以较好地模拟长白山火山地区2002—2003年间的GPS和水准测量.岩浆房扰动应力场和区域构造应力场的叠加有可能造成天池西部近EW向,天池北部以NW-NNW向为主的现今应力方向.岩浆房压力源引起的库仑应力变化有利于天池火山口NW向震群在空间上主要分布于火山口的西南和东北部.