基于粒子群算法与混合罚函数法的有限元优化反演模型及应用

岩土工程优化反分析是一个典型的复杂非线性函数优化问题,采用全局优化算法是解决这个问题的理想途径。针对常规反演方法应用于岩土工程参数反演时搜索效率低的缺点,结合粒子群算法和遗传算法的特点,充分考虑二者的互补性,提出一种效率较高的全局优化算法,以测点的实测值与计算值建立一种新的评价函数,将多目标优化问题转化为单目标优化问题,用混合罚函数法将约束问题变为无约束问题,构建了一种新的目标函数,将有限元程序ABAQUS作为一个模块嵌入到优化算法程序中,编制了有限元优化反演分析程序。并给出了应用实例验证了该法的有效性和实用性,是一种可行的参数反演方法,可应用于实际工程中复杂岩土介质初始应力场反演、渗流场以及位移反分析     

我国原位地应力测量与地应力场分析研究进展

随着我国交通、能源和水利水电等行业快速发展,在复杂地形地貌和地质条件下兴建大量岩土工程,尤其是深埋长隧道,原位地应力测量及其相关研究逾加重要。我国自20世纪60年代首次成功进行原位地应力实测以来,在地应力测量及相关研究方面已取得了长足进展。文中扼要介绍了我国地应力测量的发展历程,着重总结了在钻孔应力解除法、水压致裂法和其他常用原位地应力测量法及工程区初始地应力场分析研究方面所取得的进展,并结合我国区域构造应力环境和岩土工程的特点,对地应力测量和工程区初始地应力场研究方面的展望阐述了一点认识。     

Influence of shear angle on hangingwall deformation during tectonic inversion

Tectonic inversion is a common phenomenon in island arc settings, especially in back‐arc basins. The reactivation of normal faults as thrusts, triggered by tectonic inversion, produces typical inversion fault‐related folds and thrusts in the hangingwall. These hangingwall inversion geometries are affected by two factors: the geometry of the underlying master fault and the angle of inclined simple shear relative to the regional dip of strata, in the case that the deformation is approximated by simple shear. This study employed numerical simulations to analyse the influence of the antithetic shear angle on the geometry of the hangingwall and displacement along the master fault. The simulation results reveal that a steeply inclined shear vector during extension produces a narrow, steep‐sided half‐graben, whereas a gently inclined shear produces a wide, open basin. After tectonic inversion, a tight anticline is formed under steeply inclined shear, whereas an open anticline is formed under gently inclined shear. Antithetic shear results in reduced total displacement along the master fault, and the greater the angle between the shear direction and the regional dip, the greater the displacement along the master fault. Because the deformation geometry of syn‐extension layers is affected by extension followed by contraction, a change in the shear angle during tectonic inversion produces a wide variety of deformation geometries. Comparison of the simulation results with the results of analogue modelling suggests that the shear angle decreases by 5° during the transition from extension to tectonic inversion and that such a change may be commonly observed in natural geological structures. These results highlight the benefits of numerical simulations, which can be used to readily examine a variety of constraining parameters and thereby lead to a better understanding of the mechanism of hangingwall deformation, avoiding erroneous estimates of the amount of fault displacement.     

地面核磁共振信号相位求取电阻率

椭圆极化现象是电磁场在导电介质中传播时产生的,为得到准确的相位信息,在计算激发磁场时应考虑椭圆极化。本文研究了地面核磁共振(SNMR)信号相位求取电阻率的基本方法,利用SNMR信号通过广义逆反演法得到电阻率未知情况下的含水率及层位分布,作为相位反演电阻率的先验值,通过模拟退火法实现对电阻率的反演。对SNMR的研究解决了通过相位对电阻率的求取,有利于进一步利用相位信息提高含水层和含水率反演的精度。      

湖南保靖X区页岩气储层“甜点”区地震资料非线性反演预测方法

针对湖南保靖地区地表覆盖层复杂、地下高陡构造多变、地震勘探资料信噪比低,以及页岩气储层与周围泥岩层波阻抗对比度小等特点,研究敏感测井曲线重构与优化处理、地震资料相控非线性反演和叠后地震属性流体检测的综合储层评价方法。通过对多条测线的地震资料反演和流体检测并与实钻井结果对比,获得了良好的地质效果,初步验证了所研究方法的有效性和实用性。      

A new series solution method for two-dimensional elastic wave scattering along a canyon in half-space

This paper presents a semi-analytical method for studying the two-dimensional problem of elastic wave scattering by surface irregularities in a half-space. The new method makes use of the member of a c-completeness family of wave functions to construct the scattering fields, and then applies equal but opposite tractions to those of the foregoing constructed scattering fields on the horizontal surface of the half-space to produce additional scattering fields. These additional scattering fields are a series of Lamb's solutions. Thus the whole scattering field constructed in the series automatically satisfies the Navier equations, the condition of zero traction on the half-space surface, and the radiation boundary conditions at infinity. Using the traction-free conditions along the canyon surface, the coefficients of the series solutions are determined via a least-squares method. For incident P, SV, and Rayleigh waves, the numerical results are presented for the scattering displacements in the vicinity of a semi-circular canyon in the half-space.     

On the spatial distribution of seismicity and the 3D tectonic stress field in western Greece

We analyzed a large number of focal mechanisms and relocated earthquake hypocenters to investigate the geodynamics of western Greece, the most seismically active part of the Aegean plate-boundary zone. This region was seismically activated multiple times during the last decade, providing a large amount of enhanced quality new information that was obtained by the Hellenic Unified Seismological Network (HUSN). Relocated seismicity using a double-difference method appears to be concentrated above ∼35 km depth, exhibiting spatial continuity along the convergence boundary and being clustered elsewhere. Earthquakes are confined within the accreted sediments escarpment of the down-going African plate against the un-deformed Eurasian hinterland. The data arrangement shows that Pindos constitutes a seismic boundary along which large stress heterogeneities occur. In Cephalonia no seismicity is found to be related with the offshore Cephalonia Transform Fault (CTF). Onshore, NS crustal extension dominates, while in central and south Peloponnesus the stress field appears rotated by 90°. Shearing-stress obliquity by 30° is indicated along the major strike-slip faults, consistent with clockwise crustal rotation. Within the lower crust, the stress field appears affected by plate kinematics and distributed deformation of the lower crust and upper mantle, which guide the regional geodynamics.     

Crustal attenuation characteristics of S-waves beneath the Eastern Tohoku region,Japan

An inversion method was applied to crustal earthquakes dataset to find S-wave attenuation characteristics beneath the Eastern Tohoku region of Japan.Accelerograms from 85 shallow crustal earthquakes up to25 km depth and magnitude range between 3.5 and 5.5were analyzed to estimate the seismic quality factor Q_s. A homogeneous attenuation model Q_s for the wave propagation path was evaluated from spectral amplitudes, at 24 different frequencies between 0.5 and 20 Hz by using generalized inversion technique. To do this, non-parametric attenuation functions were calculated to observe spectral amplitude decay with hypocentral distance. Then, these functions were parameterized to estimate Q_s. It was found that in Eastern Tohoku region, the Q_s frequency dependence can be approximated with the function 33 f~(1.22) within a frequency range between 0.5 and 20 Hz. However,the frequency dependence of Q_s in the frequency range between 0.5 and 6 Hz is best approximated by Q_s(f) = 36f~(0.94) showing relatively weaker frequency dependence as compared to the relation Q_s(f) = 6 f~(2.09) for the frequency range between 6 and 15 Hz. These results could be used to estimate source and site parameters for seismic hazard assessment in the region.     

山区谷底沉陷预测模型及其参数反演

我国西部地区地势复杂,沟壑纵横,地下开采极易导致边坡失稳,引发采动滑坡。在地下采动沉降与滑坡体挤压上升的叠加影响下,谷底区域地表沉降值明显小于类似地质采矿条件下的平原地区。为准确预测山区谷底区域地表沉降值,基于简支梁的弹性变形理论,并借助概率密度函数建立了山区谷底区域地表沉陷预计修正模型,明确模型参数物理意义及其取值方法。依据修正模型,以实测值和预测值之差平方和最小为原则构建适应值函数,基于模拟退火粒子群算法提出新的模型参数反演方法,借助MATLAB语言编制了相应的参数反演程序。最后将研究成果应用于山西某矿,得到谷底区域预测结果中误差为73 mm,与实测值基本一致,取得了较好的工程实践效果。     

基于响应面法的高堆石坝瞬变-流变参数反演方法

如何准确确定原级配堆石体的力学参数是高堆石坝建设中亟待解决的一个关键问题,参数反演是解决这一问题的可行方法之一,传统的参数反演方法因需要进行大量的有限元正分析,其计算工作量大,反演效率较低。响应面法可以有效克服以上问题,但已有方法仅针对堆石体的瞬变参数,未考虑流变参数,不能满足堆石坝长期变形预测的需要。综合考虑瞬变和流变参数,通过构造更加合理的响应面函数,提出了基于响应面法的高堆石坝瞬变-流变材料参数反演方法,大大提高了反演的效率和精度。以水布垭面板堆石坝为例,采用该方法对坝体瞬变和流变参数进行了反演分析。反演结果表明,计算值与实测值在数值和变化规律上总体符合较好,反演结果合理可靠且更加高效。