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

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

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.     

Bistable plant–soil dynamics and biogenic controls on the soil production function

Soil formation results from opposite processes of bedrock weathering and erosion, whose balance may be altered by natural events and human activities, resulting in reduced soil depth and function. The impacts of vegetation on soil production and erosion and the feedbacks between soil formation and vegetation growth are only beginning to be explored quantitatively. Since plants require suitable soil environments, disturbed soil states may support less vegetation, leading to a downward spiral of increased erosion and decline in ecosystem function. We explore these feedbacks with a minimal model of the soil–plant system described by two coupled nonlinear differential equations, which include key feedbacks, such as plant‐driven soil production and erosion inhibition. We show that sufficiently strong positive plant–soil feedback can lead to a ‘humped’ soil production function, a necessary condition for soil depth bistability when erosion is assumed to vary monotonically with vegetation biomass. In bistable plant–soil systems, the sustainable soil condition engineered by plants is only accessible above a threshold vegetation biomass and occurs in environments where the high potential rate of erosion exerts a strong control on soil production and erosion. Vegetation removal for agriculture reduces the stabilizing effect of vegetation and lowers the system resilience, thereby increasing the likelihood of transition to a degraded soil state. Copyright © 2016 John Wiley & Sons, Ltd.     

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,与实测值基本一致,取得了较好的工程实践效果。     

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

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

LF-Y油田动态石油地质储量评价

LF-Y油田S区在ODP开发调整方案实施中,钻井表明构造边部变陡,含油面积变小,静态方法重新计算地质储量为1 328.30×10⁴ m³,但该储量规模与油田生产动态不相符。结合生产动态法、储量反算法及油藏数值模拟法分析油田储量规模及合理性,指出了S区存在储量计算不足的问题并分析了原因,通过动态信息反算地质储量的方法,反算储量规模在1 450.00×10⁴~1 700.00×10⁴之间,因此该油田西边无井控区域储量潜力可能性较大。结合重新开展的孔隙度反演、精细地质模型、储量计算及油藏数值模拟等研究,再通过静态方法最终计算S区的地质储量为1 566.03×10⁴ m³,储量规模增加了237.73×10⁴ m³。以动态反算储量,并反馈到静态储层评价与储量计算之中,使两者达到较为一致的结果。研究成果不但解决了S区储量不足的矛盾,使模型更加合理,剩余油分布预测也更加可靠,为ODP调整方案成功实施提供了储量基础;而且验证了油藏动态方法对储量问题研究的反证作用,从而实现了油藏动、静态结合解决地质难点问题。