Joint parameter estimation from magnetic resonance and vertical electric soundings using a multi‐objective genetic algorithm

Magnetic resonance sounding (MRS) has increasingly become an important method in hydrogeophysics because it allows for estimations of essential hydraulic properties such as porosity and hydraulic conductivity. A resistivity model is required for magnetic resonance sounding modelling and inversion. Therefore, joint interpretation or inversion is favourable to reduce the ambiguities that arise in separate magnetic resonance sounding and vertical electrical sounding (VES) inversions. A new method is suggested for the joint inversion of magnetic resonance sounding and vertical electrical sounding data. A one‐dimensional blocky model with varying layer thicknesses is used for the subsurface discretization. Instead of conventional derivative‐based inversion schemes that are strongly dependent on initial models, a global multi‐objective optimization scheme (a genetic algorithm [GA] in this case) is preferred to examine a set of possible solutions in a predefined search space. Multi‐objective joint optimization avoids the domination of one objective over the other without applying a weighting scheme. The outcome is a group of non‐dominated optimal solutions referred to as the Pareto‐optimal set. Tests conducted using synthetic data show that the multi‐objective joint optimization approximates the joint model parameters within the experimental error level and illustrates the range of trade‐off solutions, which is useful for understanding the consistency and conflicts between two models and objectives. Overall, the Levenberg‐Marquardt inversion of field data measured during a survey on a North Sea island presents similar solutions. However, the multi‐objective genetic algorithm method presents an efficient method for exploring the search space by producing a set of non‐dominated solutions. Borehole data were used to provide a verification of the inversion outcomes and indicate that the suggested genetic algorithm method is complementary for derivative‐based inversions.     

A Dynamical Statistical Model for Prediction of a Tropical Cyclone

A dynamical statistical method is applied for operational forecasting of the Bay of Bengal tropical cyclone “Nargis” of April–May 2008. The method consists of three forecast components, namely (a) analysis of Genesis Potential Parameter (GPP) and maximum potential intensity, (b) track prediction, and (c) 12 hourly intensity prediction for forecasts up to 72 hours. The results of the study showed that GPP could provide necessary predictive signal at early stages of development on the further intensification of the low pressure system into a tropical cyclone. The landfall forecast position errors by different operational numerical models (NWP) showed landfall position errors ranging from 10 km to 150 km and landfall time error ranges from 6 hours early to 6 hours delay. The dynamical statistical model is capable to provide 12 hourly nearly realistic intensity forecasts up to 60 hours of forecast.     

Moveout approximation for horizontal transversely isotropic and vertical transversely isotropic layered medium. Part I: 1D ray propagation‡

Anisotropy in subsurface geological models is primarily caused by two factors: sedimentation in shale/sand layers and fractures. The sedimentation factor is mainly modelled by vertical transverse isotropy (VTI), whereas the fractures are modelled by a horizontal transversely isotropic medium (HTI). In this paper we study hyperbolic and non‐hyperbolic normal reflection moveout for a package of HTI/VTI layers, considering arbitrary azimuthal orientation of the symmetry axis at each HTI layer. We consider a local 1D medium, whose properties change vertically, with flat interfaces between the layers. In this case, the horizontal slowness is preserved; thus, the azimuth of the phase velocity is the same for all layers of the package. In general, however, the azimuth of the ray velocity differs from the azimuth of the phase velocity. The ray azimuth depends on the layer properties and may be different for each layer. In this case, the use of the Dix equation requires projection of the moveout velocity of each layer on the phase plane. We derive an accurate equation for hyperbolic and high‐order terms of the normal moveout, relating the traveltime to the surface offset, or alternatively, to the subsurface reflection angle. We relate the azimuth of the surface offset to its magnitude (or to the reflection angle), considering short and long offsets. We compare the derived approximations with analytical ray tracing.     

Definition and interests of reciprocity and reciprocity gap principles for groundwater flow problems

We introduce the reciprocity and reciprocity gap principles for flow problems in hydrogeology and illustrate their interest in addressing identification problems. The reciprocity principle is derived from mechanics and establishes for flow problems a relationship between different sets of forcing terms, including sources, sinks and boundary conditions, and the resulting head fields. The reciprocity gap principle compares different head fields resulting from the same forcing terms applied to different structures. We give general 2D expressions of the reciprocity and reciprocity gap principles for transient flow problems and give two examples of applications for the identification of transmissivity values and interfaces between different transmissivities. Identification capacities of the reciprocity and reciprocity gap principles yielding direct inversion methods could be used as initial guesses for more advanced inverse problem methodologies.     

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

将Cole—Cole模型表示的地下异常体各剖分单元复电阻率,引入到体积分方程求解均匀大地三维电磁场计算中,实现电偶源地面激发、接收复电阻率三维电磁场正演模拟。     

Identification of diffusion parameters in a nonlinear convection–diffusion equation using the augmented Lagrangian method

Numerical identification of diffusion parameters in a nonlinear convection–diffusion equation is studied. This partial differential equation arises as the saturation equation in the fractional flow formulation of the two-phase porous media flow equations. The forward problem is discretized with the finite difference method, and the identification problem is formulated as a constrained minimization problem. We utilize the augmented Lagrangian method and transform the minimization problem into a coupled system of nonlinear algebraic equations, which is solved efficiently with the nonlinear conjugate gradient method. Numerical experiments are presented and discussed. This work was partially supported by the Research Council of Norway (NFR), under grant 128224/431.     

Robust Regulation of Alkalinity in Highly Uncertain Continuous Anaerobic Digestion Processes

In this paper, a multiple‐input multiple‐output (MIMO) model‐based robust control scheme is proposed for the indirect control of both total alkalinity and the ratio (intermediate alkalinity)/(total alkalinity) by regulating volatile fatty acid concentrations and strong ions concentration, while guaranteeing the so‐called operational stability in anaerobic digestion (AD) processes. The proposed MIMO regulator is an adaptive controller derived from an AD model which incorporates the physicochemical equilibrium of the system as well as the use of a robust interval observer to estimate key process bounds that are used in the computation of the control efforts. Numerical simulations were carried out for a number of operating conditions under the most uncertain scenarios. Results showed that the proposed multivariable control law is able to recover the system stability around a pre‐determined set point in the face of parameter uncertainty and load disturbances.     

集合预报误差在GRAPES全球四维变分同化中的应用研究Ⅱ：参数选取与数值试验

研究的第一部分讨论了如何有效应用集合预报误差的科学方案,确定了集合预报误差在GRAPES（Global Regional Assimilation and PrEdiction System）全球4DVar（four dimensional variational data assimilation）中应用的分析框架。在此基础上研究了针对集合预报误差实际应用于GRAPES全球4DVar,解决接近或超过100个集合样本数时高效生成的计算效率问题,以及与GRAPES全球4DVar匹配的同化关键参数确定问题。选择基于4DVar的集合资料同化方法生成集合样本,通过将第1个样本极小化迭代过程中产生的预调节信息用于其他样本极小化做预调节,将计算效率提高了2倍。通过时间错位扰动方法增加集合样本数,实现集合样本增加到3倍。对集合方差进行膨胀,并选择水平局地化相关尺度为流函数背景误差水平相关的1.4倍。通过批量数值试验方法确定背景误差与集合预报误差的权重系数,对60个集合样本当集合预报误差权重为0.7时预报效果最好。对北半球夏、冬两季各52 d的批量试验表明,对于南、北半球En4DVar （ensemble 4DVar）较4DVar的改进在冬季主要集中在700—30 hPa,而在夏季主要集中在400—150 hPa。赤道地区受季节影响较小,En4DVar对位势高度、风场与温度的改进都较为明显,且经向风场的改进最为显著。文中研发的集合预报误差在GRAPES全球4DVar中应用的方法合理可行。      

Application of Proper Orthogonal Decomposition (POD) to inverse problems in saturated groundwater flow

We develop a new Proper Orthogonal Decomposition (POD) reduced order model for saturated groundwater flow, and apply that model to an inverse problem for the hydraulic conductivity field. We use sensitivities as the POD basis. We compare the output when the optimizer uses the reduced order model against results obtained with a full PDE based model. The solutions generated using the POD reduced model are comparable in residual norm to the solutions formed using only the full-scale model. The material parameters are similarly comparable. The time to solution when using the reduced model is reduced by at least an order of magnitude, as are the number of calls to the full model.     

钢管水泥土围护墙等效弹性模量的有限元计算

分析了水泥土围护墙弹性模量对于控制基坑变形的重要性、根据上海某钢管水泥土墙基坑支护的工程实例建立ABAQUS平面二维有限元模型,计算钢管水泥土围护墙等效弹性模量,量化了内插钢管对提高水泥土重力式围护结构的支护强度所起到的作用,并对钢管长度、外径、间距和插入位置进行了参数分析。