基于模拟退火算法的饲料配方设计

本文在饲料配方设计中首次引入新颖、高效的组合优化算法--模拟退火算法，并获得理想的数值结果，说明了在传统的饲料配方设计中高技术的应用可以获得更高的效益.     

模拟退火算法在控制网平差中的应用

线性最小二乘估计在对非线性函数进行线性近似的过程中会产生模型误差,而一些非线性参数估计方法可能因为函数复杂而难以求导,法方程系数矩阵秩亏或呈病态矩阵时难以求解,非线性迭代解法有时对初始值的选择存在依赖性,不恰当的初始值会导致迭代无法收敛。针对这些问题,引入了模拟退火算法,介绍了该算法的基本原理、计算步骤和收敛性,并以3个控制网平差应用为例,说明该算法具有无需求导求逆,简洁实用,易于编程等优势,并能实现全局优化,获得高精度的平差结果。     

Characterizing heterogeneous permeable media with spatial statistics and tracer data using sequential simulated annealing

Characterizing heterogeneous permeable media using flow and transport data typically requires solution of an inverse problem. Such inverse problems are intensive computationally and may involve iterative procedures requiring many forward simulations of the flow and transport problem. Previous attempts have been limited mostly to flow data such as pressure transient (interference) tests using multiple observation wells. This paper discusses an approach to generating stochastic permeability fields conditioned to geologic data in the form of a vertical variogram derived from cores and logs as well as fluid flow and transport data, such as tracer concentration history, by sequential application of simulated annealing (SA). Thus, the method incorporates elements of geostatistics within the framework of inverse modeling. For tracer-transport calculations, we have used a semianalytic transit-time algorithm which is fast, accurate, and free of numerical dispersion. For steady velocity fields, we introduce a transit-time function which demonstrates the relative importance of data from different sources. The approach is illustrated by application to a set of spatial permeability measurements and tracer data from an experiment in the Antolini Sandstone, an eolian outcrop from northern Arizona. The results clearly reveal the importance of tracer data in reproducing the correlated features (channels) of the permeability field and the scale effects of heterogeneity.     

改进的模拟退火－单纯形综合反演方法

实际中的大量地球物理反演是一个多参数、非线性优化问题,所采用的目标函数,即度量由参数化的理论模型得出的预测值与观测值的吻合程度,往往具有多个局部极值．针对这类问题,本文综合全局反演方法具有的全域搜索能力强、局部方法收敛速度快和“均匀设计”布点效率高的特点,提出了模拟退火－单纯性综合反演方法,并通过一维声波非线性反演验证了这种综合方法的搜索能力和效率．     

Modeling the seismic response of basin-edge structure by a sequential waveform method

A sequential waveform method is developed to simulate the seismic response of basin-edge structure excited by a plane incident P-wave. The full procedure involves: (a) a previous parameterization of the investigated model using the seismic wave velocities and depths of the sedimentary stratifications; (b) an input motion determined from the records at stations installed on hard rock; (c) forward computation of the P-SV elastic wave field by means of a two-dimensional finite difference (FD) method; (d) the optimization of the model vector using simulated annealing technique and comparing the simulated seismic response of the tested structure with the observed wave field; (e) the correction of the initial model by trial-and-error by testing the differences between synthetics and observed data, and (f) the final solution obtained by iteration using the conjugate gradient algorithm. The search of an optimal basin-edge model has been parallel processed by varying the shapes and velocities of strata on the basis of the fitting of relative timing, amplitude and phase between the output and the observed data. The input motion and sensitivity have been checked and the validity of the method has been demonstrated by numeric analysis. Using the teleseismic records generated by 7 earthquakes recorded at 26 broadband seismic stations, we have studied the seismic velocity structure of the southern edge of the Jiyang depression located in the Bohai Bay basin, northern China. Two cross sections show an agreement between the velocity results and the geological sections available in the region. In addition, we obtain evidence of three hidden faults under the sections and features that suggest major extensions at the Paleogene.     

Isothermal annealing of partially metamict zircon: evidence for a three-stage recovery process

 Raman spectroscopy and the powder diffraction technique have been used to monitor the recovery process of two partially metamict zircons (2.6 and 4.8 × 10¹⁸ α-decays g⁻¹) from Sri Lanka during a series of isothermal annealing experiments in the temperature range from 870 to 1622 K. These experiments show for the first time that structural recovery in partially metamict zircon proceeds via three distinct recovery stages, each of which occurs within a distinct time-temperature regime. Whereas the first two stages have previously been recognized (recovery of damaged crystalline remnants and epitaxial recrystallization), the third stage has not yet been identified as a single activated process. It is suggested that anisotropic defect annealing during the first stage at low temperatures, where the structure recovers preferentially along the a(b) plane, produces a geometrical situation where large structural rearrangements are necessary to remove the remaining defects inside the crystalline material. This situation is approximately reached when the amorphous domains start to recrystallize. The reason for anisotropic annealing can be found in a different connectivity between polyhedral linkages in both directions of the zircon lattice. High apparent activation energies, in the range of 6.4 to 7.9 eV, were determined for the third recovery stage from the Raman data, which are interpreted to reflect large structural rearrangements (i.e. polyhedral tilting) associated with the final recovery of the c axis. This explains the occurrence of a distinct recrystallization stage without defect annealing. Finally, it should be mentioned that the first recovery stage is not necessarily expected to occur in less damaged zircon crystals (<∼2 × 10¹⁸ α-decays g⁻¹), since less stable defects along the basal plane might have already been self-annealed during radiation damage accumulation under ambient temperatures. Received: 6 September 2001 / Accepted: 25 February 2002     

Hydration and refractive indices of Holocene tephra glass in Hokkaido, Northern Japan

Refractive indices of incompletely hydrated tephra glasses vary widely. Thus, glass refractive index is not a practical indicator for identifying Holocene tephras. Hydration near the surface of tephra glass shards can be removed by either hydrofluoric acid treatment or annealing for 12 h at 400°C. The annealing procedure is a particularly reliable and simple way to dehydrate tephra glasses. Standard deviations (s) of refractive indices for glasses after 12 h annealing are small (s=0.0014–0.0018), in contrast with untreated glasses (s=0.0034–0.0405). The refractive indices of dehydrated tephra glasses are 0.006–0.014 lower than those of the untreated (hydrated) glasses. Using the 400°C 12-h annealing procedure, values for the refractive indices of eight Holocene tephras in Hokkaido were obtained. These refractive indices are useful indicators for identification and correlation of Holocene tephras.     

地震静校正全局最优化问题的求解

针对地震静校正存在的非线性和多参数性等问题，本文提出了一种全局优化的剩余静校正方法.在模型更新规则上，把模拟退火方法和均匀优化设计方法结合在一起，使模型的更新更为合理，从而加快了搜寻全局最优解的速度.同时提出了模型的不同分量的温度参数和退火过程的选取方式，把温度参数的选取与地震剖面的能量联系起来，使温度参数的选取具有自适应的特点.该方法克服了常规模拟退火方法所具有的寻优空间不均匀以及退火参数需通过多次试验选取的缺陷.通过实际地震资料的计算证明，本文提出的方法对地震静校正问题合理而有效.     

镍铜硫化物矿石中磁黄铁矿固溶体的退火及其选矿意义

磁黄铁矿固治体从硫化物熔体结晶后，在缓慢冷却过程中经历了显著的退火。出治和出治体的租化是固治体退火的两种方式。叶片状的单斜磁黄铁矿和“火焰状”的镍黄铁矿原始出治相在降温过程中均可发生退火和租化。分布于磁黄铁矿等矿物粒间或包于磁黄铁矿粒内的粒状镍黄铁矿，不只是高温出治的直接产物，有一部分可能是由火焰状出治体租化而成的。磁黄铁矿中单斜变体的出治和租化可使矿石的磁性发生改变，镍黄铁矿出治体的租化使含镍矿物的粒度加大。因而，退火作用对矿石的选矿工艺性能有着显著影响。     

模拟退火算法在动态建立基准地价模型中的应用

阐述了传统基准地价评估模型建立流程将地价样本检验与模型建立分离的不足，引入了模拟退火算法思想，对目标函数、初始状态、新状态产生、新状态接受、初温、温度更新函数以及抽样准则、算法终止准则等关键问题进行设计，将地价样本检验、优化组合与模型建立有机结合，并设计了基于模拟退火算法的动态建立基准地价模型的流程，最后应用实例验证了方法的可行性。