基于退火BP神经网络的GPS高程转换

阐述模拟退火算法的基本思想和原理,提出并介绍模拟退火算法优化的BP神经网络模型在GPS高程转换中的具体应用,同时编写相应的MATLAB处理程序,结合大量数据进行仿真实验,结果表明文中提出的退火BP神经网络模型具有收敛速度快、精度高、避免陷入局部最小的优良特性。     

混合禁忌搜索算法在湖泊富营养化评价中的应用

本文提出了幂函数加和型湖泊富营养化评价指数公式,同时将禁忌搜索算法与模拟退火、免疫进化算法有机地结合起来,建立了求解优化问题的混合禁忌搜索算法(HTS).将此算法用于评价指数公式的参数优化,从而建立起湖泊富营养化综合评价模型,为湖泊富营养化评价和预测提供了一条有效的途径.将优化好的湖泊富营养化模型用于东湖富营养化评价的实例研究,结果表明此评价模型能较好地反映湖泊富营养化水平,并具有较好的通用性和实用性.     

单纯形—模拟退火混合算法反求水文地质参数及其并行求解

本文利用单纯形法局部搜索速度快和模拟退火算法全局寻优的特点,同时为了克服各自算法的弊病,提出采用单纯形—模拟退火混合算法(SMSA)进行水文地质逆问题的求解。论文详细描述了SMSA算法的具体操作算子的实现,并将该算法应用于一个大型水源地的地下水模拟反演。计算结果表明,SMSA算法在水文地质参数反演计算具有求解速度快,精度高的特点,而且易于实现并行运算。     

单纯形一模拟退火混合算法反求水文地质参数及其并行求解

本文利用单纯形法局部搜索速度快和模拟退火算法全局寻优的特点，同时为了克服各自算法的弊病，提出采用单纯形一模拟退火混合算法（SMSA）进行水文地质逆问题的求解。论文详细描述了SMSA算法的具体操作算子的实现，并将该算法应用于一个大型水源地的地下水模拟反演。计算结果表明，SMSA算法在水文地质参数反演计算具有求解速度快，精度高的特点，而且易于实现并行运算。     

Raman microspectroscopy: A non-destructive tool for routine calibration of apatite crystallographic structure for fission-track analyses

Compositional control on the annealing kinetics of fission-tracks (FT) in apatite requires routine measurement of sample grain composition. However, for practical reasons the bulk composition of analysed grains is not routinely measured and instead grain chlorine content or etch-pit dimensions are used to characterise a samples annealing behaviour. A more desirable approach is to measure crystallographic parameters (i.e. unit cell dimension) of a grain as these represent the summed effect of all substitutions and crystal defects. We show how Raman microspectrometry can be used as a routine non-destructive tool to obtain rapid measurement of the crystallographic structure of apatite grains etched for FT analysis. Variations of unit cell parameter a are found to correspond to a systematic variation of Raman shift in the range of 452–440 cm^− 1 for measurements made on c-parallel sections of apatite where the direction of the polarized incident beam is parallel to the c axis.     

Effect of alpha-damage annealing on apatite (U–Th)/He thermochronology

Recent data on He diffusion challenge the temperature sensitivity of apatite (U–Th)/He thermochronology: the damage induced by recoil of U and Th decay series during emission of α particles (α-recoil damage) has been proposed to modify He-diffusion properties through time. However, we propose that annealing of these irradiation defects may be an important phenomenon and may be significant in case of slowly-cooled or reheated basement rocks. To test this hypothesis, we developed a quantitative model including an explicit treatment of α-recoil damage, annealing, and their effect on He-diffusion kinetics, and calibrate it against literature data. Our model is based on two hypotheses: (1) helium is in equilibrium between an apatite crystal and its defects and (2) alpha-recoil damage annealing can be described analogously to fission-track annealing. This model has been embedded into a Monte Carlo simulation of helium production/ejection/diffusion and applied to data from the French Massif Central; a complex slowly-cooled terrain with burial reheating, where the thermal history has been constrained by previous fission-track (FT) data including FT length distributions. (U–Th)/He ages are close to the FT ages from the same samples and are generally reproducible among replicates, but some samples present He-age dispersion that is not correlated with crystal size. Our model reproduces the Massif Central data very well except for three samples where He ages are older than corresponding FT ages. We show that annealing of irradiation damage has an important impact on retentivity of helium and that the He content, [He] is only a rough approximation of the damage level. In particular our results show that independence of He ages on crystal sizes, in case of reheated samples, is a clear indication of the higher He retentivity induced by α-recoil defects and that an explicit treatment of defect annealing is required for a correct interpretation of (U–Th)/He ages in such a case. More generally a correlation with the crystal size can bring information on the thermal path only if the age of defects, well represented by the fission-track age, is available, due to the dependence of the partial retention zone on damages. Conversely, in case of rapid cooling or for samples having low U and Th contents, damage effects can be ignored without significant effects on He ages.     

模拟退火算法在江苏及邻区地震定位中的应用

使用模拟退火算法,以2019年4月6日江苏溧水M_L 3.3地震为典型震例,进行地震定位反演试算,确定适用于江苏及邻区地震定位的最佳初始温度、降温策略和扰动函数。设置初始温度为25,降温策略为T = 0.96T,选用效率较高的局部收敛加强型方法产生随机模型,对2019年江苏及周边地区105次M_L≥1.8地震进行重新定位,结果表明,76%的事件震中经、纬度定位误差绝对值在0.1°以内,但深度定位误差较大。鉴于该算法对P波到时信息及地震台站空间分布的依赖性,认为模拟退火算法在江苏及邻区陆地地震定位中效果较好,对于海域地震,则定位误差较大。     

模拟退火方法在三维速度模型地震波走时反演中的应用

采用块状建模以及三角形拼接的界面描述方式,并通过立方体速度网格线性插值获得块体内部的速度分布。正演过程中采用逐段迭代射线追踪方法计算三维复杂地质模型中的射线走时,并采用模拟退火方法进行了三维模型中的地震波走时反演研究。模型测试结果表明,使用的射线追踪和走时反演算法有效。     

Thomsen's parameters from p‐wave measurements in a spherical sample‡

The aim of this paper is to understand the seismic anisotropy of the overburden shale in an oilfield in the North West Shelf of Western Australia. To this end, we first find the orientation of the symmetry axis of a spherical shale sample from measurements of ultrasonic P‐wave velocities in 132 directions at the reservoir pressure. After transforming the data to the symmetry axis coordinates, we find Thomsen's anisotropy parameters δ and ? using these measurements and measurements of the shear‐wave velocity along the symmetry axis from a well log. To find these anisotropy parameters, we use a very fast simulated re‐annealing algorithm with an objective function that contains only the measured ray velocities, their numerical derivatives and the unknown elasticity parameters. The results show strong elliptical anisotropy in the overburden shale. This approach produces smaller uncertainty of Thomsen parameter δ than more direct approaches.     

基于混合优化算法的MT阻抗张量畸变分解方法

在前人研究基础上,对Groom-Bailey(GB)张量分解畸变因子和区域阻抗的求解方法进行了改进.首先,通过Swift旋转与GB分解的扭变和剪切矩阵的求逆变换,利用变换后区域阻抗主对角元素为0的条件获得关于扭变因子和剪切因子的超定方程组,采用模拟退火全局优化算法进行求解.其次,由得到的扭变因子和剪切因子,结合Swift旋转确定的走向角和区域阻抗元素的估计,作为非线性最小二乘局部优化算法的初始值,对GB分解定义式的超定方程组进行求解,得到各畸变参数和区域阻抗的解.通过模型试验验证了方法的正确,对方法的稳定性进行了比较与评价,并通过与已有结果的对比和实际资料的应用,表明了方法实际应用的效果.