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

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

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

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

Mesozoic-Cenozoic tectonics and geodynamics of Altai, Tien Shan, and Northern Kazakhstan, from apatite fission-track data

Apatite fission-track (AFT) thermochronological modeling as a diagnostic tool for periods of stability (peneplanation) and tectonic activity (orogeny) has been broadly used in tectonic studies of Central Asia in recent years. We discuss more than 100 AFT ages of samples from the Kyrgyz Tien Shan and Altai and compare them with AFT data from northern Kazakhstan. Geological, geomorphological, and AFT data indicate intense activity in the Late Cenozoic Eurasian continental interior. The impact from the India-Eurasia collision on the northern Tien Shan, Altai, and northern Kazakhstan regions showed up at 11, 5, and 3 Ma, respectively, as a result of stress propagation into the continent, with the ensuing reactivation and mountain growth. We hypothesize that a distant effect of the Late Cenozoic India-Eurasia collision was to rejuvenate Paleozoic fault zones and to deform the Mesozoic sedimentary cover north of the collision front as far as the West Siberian Plate. The reactivation facilitated formation of tectonic oil and gas traps. The activity in northern Central Asia under the effect of the Indian indentation into Eurasia appears to continue and may evolve to include uplift of southern West Siberian plate with uplift.     

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.     

青藏高原东缘龙门山晚新生代剥蚀厚度与弹性挠曲模拟

龙门山是青藏高原东缘边界山脉,具有青藏高原地貌、龙门山高山地貌和山前冲积平原三个一级地貌单元。利用数字高程模式图像和裂变径迹年代测定方法研究和计算龙门山晚新生代剥蚀厚度与剥蚀速率,结果表明:3.6 Ma以来龙门山的剥蚀厚度介于1.91-2.16 km之间,剥蚀速率介于0.53-0.60 mm/a之间。在此基础上,开展了该地区岩石圈的弹性挠曲模拟,结果表明龙门山的隆升机制具有以构造缩短隆升和剥蚀卸载隆升相叠合的特点。3.6 Ma之前,龙门山的隆升与逆冲推覆构造负载有关,以构造缩短驱动的构造隆升为特色;3.6 Ma之后,龙门山的隆升与剥蚀卸载驱动的抬升有关,并以剥蚀卸载隆升为特色,进而提出了龙门山晚新生代以来的隆升机制以剥蚀成山作用为主的认识。     

Denudation and cooling of the Lake Teletskoye Region in the Altai Mountains (South Siberia) as revealed by apatite fission-track thermochronology

Lake Teletskoye occupies a narrow graben located in the northwestern sector of the Altai fold belt in South Siberia. The lake basin is thought to have formed during the Pleistocene as a distant result of the Cenozoic collision of India and Eurasia that caused a tectonic reactivation of the Palaeozoic Gorny–Altai (GA) and West Sayan (WS) blocks.The present work reports of a pilot fission-track study performed on 13 apatite separates collected from rocks that were sampled along two profiles in close proximity of the lake. The age–length data and AFT thermochronological modelling reveal two important phases of cooling in the Altai Mountains, a first one during the Late Jurassic–Early Cretaceous and a second one that started in the Miocene–Pliocene and that persists until today. The first event is interpreted to result from uplift-induced denudation probably related to the closure of the Mongol–Okhotsk Ocean; the second event can be linked to the young Cenozoic movements that lie at the origin of the formation of the Lake Teletskoye basin.     

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.