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.     

Hydrogeochemical stratification of the unconfined Samobor aquifer (Zagreb,Croatia)

In the western part of the city of Zagreb, Croatia, hydrogeochemical and isotopic investigations of the Samobor aquifer were carried out with the aim of determining the differences in hydrogeochemical characteristics at increasing aquifer depths. The aquifer comprises 40-m thick gravelly–sandy deposits, with lenses and interlayers of silt and clay. The analyses have proven that with increasing aquifer depth, there are decreases in groundwater temperature and the values of electrical conductivity and increases in the sodium, iron and manganese concentrations. The δ¹³C distribution shows an evident increase in biogenic carbon concentrations with increasing aquifer depth. The measured specific ¹⁴C activities showed that the deeper part of the aquifer is characterized by slow water exchange, while the shallower part is influenced by current recharge, although the pumping wells located on the well-field downstream penetrate the aquifer fully. A direct exchange of water from the Sava River and groundwater occurs in the near vicinity of the river. This exchange weakens further away, while the difference in hydrogeochemical characteristics between the Sava River water and groundwater increases.     

弹性层状半空间中无限长洞室对斜入射平面SH波的三维散射（Ⅰ）——方法及验证

利用间接边界元法, 求解了弹性层状半空间中无限长洞室对斜入射平面SH波的三维散射问题. 通过与文献结果进行对比, 验证了本方法的正确性. 与工程中常用的二维模型比较表明, 工程中将三维散射问题简单地分解为平面内问题和平面外问题的做法存在较大误差. 文中并研究了斜入射角度对散射的影响, 表明斜入射角度对地表位移幅值有着重要影响.      

Amplification of in-plane seismic ground motion by group cavities in layered half-space (Ⅰ)

Amplification of in-plane seismic ground motion by underground group cavities in layered half-space is studied both in frequency domain and time domain by using indirect boundary element method (IBEM), and the effect of cavity interval and spectrum of incident waves on the amplification are studied by numerical examples. It is shown that there may be large interaction between cavities, and group cavities with certain intervals may have significant amplification to seismic ground motion. The amplification of PGA (peak ground acceleration) and its PRS (peak response spectrum) can be increased up to 45.2% and 84.4%, for an example site in Tianjin, under the excitation of Taft wave and El Centro wave; and group cavities may also affect the spectra of the seismic ground motion. It is suggested that the effect of underground group cavities on design seismic ground motion should be considered.     

Evolution of borate minerals from contact metamorphic to hydrothermal stages: Ludwigite-group minerals and szaibélyite from the Vysoká – Zlatno skarn,Slovakia

Mineralogy and Petrology - Borate minerals of the ludwigite group (LGM) and szaibélyite in association with hydroxylclinohumite, clinochlore, a serpentine mineral, magnesian magnetite, spinel,...     

深埋电极的地电阻率观测研究

首先研究了四极观测系统装置系数与电极埋深的关系;然后给出了点电流源在3层地壳模型的地表和第二层时,电源所在层的电位的解析表达式;最后将天津宝坻地区的电性结构简化成一个3层模型,计算给出了当地表层和基岩中的电阻率出现变化时,在地表和基岩上层开展四极地电阻率观测结果与供电极距和深度的关系.     

1975-2008 年西藏色林错湖面变化对气候变化的响应

根据1975 年地形图、80 年代至2008 年的TM、CBERS 卫星遥感资料和近34 年(1975-2008 年) 的气温、降水量、蒸发量、最大冻土深度等气候资料分析得出,西藏那曲地区西部的色林错及其周围的错鄂、雅根错的面积在近30 年来呈较显著的扩大趋势,到2008 年面积分别为2196.23 km2、279.24 km²、103.07 km²,与1975 年分别增长了574.46 km²、11.59 km2和68.13 km²,增长速度分别为35.4%、4.3%和195%。色林错从1999-2008 年湖面扩大速度为20%,平均上涨了420 km²/10a,已超过纳木错面积,成为西藏第一大咸水湖。冰雪融水量的增加是湖泊上涨的根本原因,其次与降水量的增加和蒸发量的减少、冻土退化等暖湿化的气候变化存在很大关系。     

Height transformation models from ellipsoidal into the normal orthometric height system for the territory of the City of Zagreb

The paper presents the testing of the possibility of determining the heights of GPS points in the homogeneous field in the new Croatian Height Reference System (HVRS71) by using the method of height transformation. The testing was made in the area of Zagreb. As part of the field works, normal orthometric heights of 27 GPS points were determined according to the new height system, by transferring the benchmark heights using the geometric levelling method, thus obtaining GPS/levelling points of known ellipsoidal and normal orthometric heights. The GPS/levelling points served as the basis for determining the transformation models that enabled the computation of normal orthometric heights from ellipsoidal heights of any GPS point in the observed area. The empirical data used for modelling were reduced undulation dN values of GPS/levelling points. As part of the dN modelling with parametric functions, the approximation surfaces were obtained on the basis of three polynomials: FN310, FN312 and FN318. The transformation models were also tested using non-parametric Watson and Loess algorithms. The FN318 and Loess models yielded the best results.     

基于动力学震源模型的三维沉积盆地直下型断层地震动模拟

沉积盆地与近断层地震共同作用会增加地震破坏的风险水平,尤其是盆地下方直下型断层发震情况。采用动力学震源模型刻画断层破裂发震过程,开展沉积盆地直下型断层谱元法地震动模拟研究,探讨不同断层面初始剪应力和成核区位置下三维沉积盆地地表响应规律。研究结果表明,断层面应力降对盆地地表地震动的影响显著,在断层面强度一定的情况下,随着初始剪应力的增大,即应力降增大,盆地地表峰值响应增大,原因在于应力降的改变影响了断层破裂释放能量,进而引起断层破裂速度改变,最终导致盆地地表响应发生变化;改变断层面成核区位置会对盆地内部地震动分布规律产生影响,当成核区位置从断层中间向断层左侧移动时,盆地左侧地震动逐渐减小,而右侧地震动逐渐增大,最终表现为盆地右侧地震动显著高于盆地左侧,原因在于改变成核区位置后,导致近断层地震动的方向性效应发生变化。