Three-dimensional forward modeling for magnetotelluric sounding by finite element method

A finite element algorithm combined with divergence condition was presented for computing three-dimensional(3D) magnetotelluric forward modeling. The finite element equation of three-dimensional magnetotelluric forward modeling was derived from Maxwell's equations using general variation principle. The divergence condition was added forcedly to the electric field boundary value problem, which made the solution correct. The system of equation of the finite element algorithm was a large sparse, banded, symmetric, ill-conditioned, non-Hermitian complex matrix equation, which can be solved using the Bi-CGSTAB method. In order to prove correctness of the three-dimensional magnetotelluric forward algorithm, the computed results and analytic results of one-dimensional geo-electrical model were compared. In addition, the three-dimensional magnetotelluric forward algorithm is given a further evaluation by computing COMMEMI model. The forward modeling results show that the algorithm is very efficient, and it has a lot of advantages, such as the high precision, the canonical process of solving problem, meeting the internal boundary condition automatically and adapting to all kinds of distribution of multi-substances.

Forward modeling for “earth-ionosphere” mode electromagnetic field

A fixed artificial source(greater than 200 kW) was used and the source location was selected at a high resistivity region to ensure high emission efficiency. Some publications used the "earth-ionosphere" mode in modeling the electromagnetic(EM) fields with the offset up to a thousand kilometer, and such EM fields still have a signal/noise ratio of 10-20 dB. This means that a new EM method with fixed source is feasible, but in their calculation, the displacement in air was neglected. In this work, some three-layer modeling results were presented to illustrate the basic EM fields' characteristics in the near, far and waveguide areas under "earth-ionosphere" mode, and a standard is given to distinguish the boundary of near, far and waveguide areas. Due to the influence of the ionosphere and displacement current in the air, the "earth-ionosphere" mode EM fields have an extra waveguide zone, where the fields' behavior is very different from that of the far field zone.

Three-dimensional forward modeling for the SBTEM method using an unstructured fi nite-element method

In this study,we propose a three-dimensional(3D)forward modeling algorithm of surface-to-borehole transient electromagnetic(SBTEM)fields based on an unstructured vector fi nite-element method to analyze the characteristics of SBTEM responses for complex geoelectrical models.To solve the double-curl diff usion equation for the electric fi eld,we use an unstructured tetrahedral mesh to discretize the model domain and select the unconditionally stable backward Euler scheme to discretize the time derivative.In our numerical experiments,we use a grounded wire as a transmitting source.After validating the algorithm’s eff ectiveness,we first analyze the diffusion characteristics and detectability of the electromagnetic field.After that,we focus our attention on the distribution and the cause of zero bands for Ex and dBy/dt components with the hope of guiding future field surveys.Finally,by simulating diff erent models,we analyze the capability of the SBTEM method in detecting typical mineral veins so that we can provide a reference for mineral resource exploration in the deep earth.

Seismic Centrifuge Modelling of Offshore Wind Turbine with Tripod Foundation

The development of offshore wind energy is fast as it is clean, safe and of high efficiency. The harsh marine environment raises high demand on the foundation design of offshore wind turbine. Earthquake loading is one of the most significant factors which should be considered in the design phase. In this paper, a group of earthquake centrifuge tests were conducted on a physical wind turbine model with tripod foundation. The seismic responses of both wind turbine model and foundation soil were analyzed in terms of the recorded accelerations, pore water pressures, lateral displacements and settlements. The results were also compared with those measured in the previous research on mono-pile foundation. It is demonstrated that the tripod foundation can provide better resistance in the lateral displacement and structural settlement under earthquake loading.

FEM Using Adaptive Topography in 2-D MT Forward Modeling

Complex terrain causes great MT noise.This paper puts forward a FEM model method using adaptive topography and quadratic elements based on studies by previous researchers.This method can model all kinds of complicated terrain and geoelectric bodies preferably.The numeric modeling,calculation of the auxiliary field and definition of resistivity are deduced by electromagnetic equations.Lastly,several examples are presented,which show the method is rapid,effective and has highly accurate.

基于跨孔电阻率CT的污染阻隔帷幕性能检测

为克服传统钻孔取样和常规阵列测量的跨孔电阻率法在帷幕检测工作中的不足,采用阵列优化方法对3种常规跨孔电阻率阵列进行优选,构建了高分辨率和高信噪比的阻隔帷幕测量组合阵列。采用数值模拟方法对该组合阵列进行阻隔帷幕检测效果验证,并应用于某垃圾填埋场的阻隔帷幕性能检测。根据垃圾填埋场2处阻隔帷幕的跨孔电阻率CT探测结果得出:膨润土阻隔帷幕主体在固结初期表现为均匀完整的低阻带,电阻率为0.001~0.06Ω·m,2处阻隔帷幕的注浆影响半径分别为0.55,0.95 m,满足帷幕厚度和完整性要求。结合阻隔帷幕芯样的渗透系数测试分析得出:2个测量区域阻隔材料渗透系数均<10^(-6)cm/s,满足污染阻隔的性能要求。利用组合阵列的跨孔电阻率CT技术,实现对污染场地阻隔帷幕完整性和防渗性能的高效准确评价。

Measurement and modeling of the effective thermal conductivity for porous methane hydrate samples

The effective thermal conductivities of gas-saturated porous methane hydrates were measured by a single-sided transient plane source (TPS) technique and simulated by a generalized fractal model of porous media that based on self-similarity.The density of porous hydrate,measured by the volume of the sample in the experimental system,was used to evaluate the porosity of methane hydrate samples.The fractal model was based on Sierpinski carpet,a thermal-electrical analogy technique and one-dimensional heat flow assumption.Both the experimental and computational results show the effective thermal conductivity of methane hydrate decreases with the porosity increase.The porosity of 0.3 can reduce the thermal conductivity of the methane hydrate by 25%.By analysis of the experimental data and the simulative result,the optimized thermal conductivity of the zero-porosity methane hydrate is about 0.7 W m-1K-1.