基于PSR-WSVM模型的边坡位移预测

为建立高精度的边坡位移预测模型，采用相空间重构（PSR）将边坡位移时间序列数据转换为多维数据，同时构造小波核函数改进的支持向量机模型，建立PSR-WSVM模型并应用于边坡位移预测。将PSR-WSVM模型预测结果与传统支持向量机（SVM）模型、小波支持向量机（WSVM）模型和基于相空间重构的支持向量机（PSR-SVM）模型预测结果进行对比，通过平均绝对误差（MAE）、平均绝对误差百分比（MAPE）和均方根误差（RMSE）3个精度评价指标验证PSR-WSVM模型的可行性。工程实例结果表明，PSR-WSVM模型预测结果的3个精度评价指标都优于另外3种模型，边坡位移预测的精度明显提升。     

Analytical study of ground motion caused by seismic wave propagation across faulted rock masses

Studying seismic wave propagation across rock masses and the induced ground motion is an important topic, which receives considerable attention in design and construction of underground cavern/tunnel constructions and mining activities. The current study investigates wave propagation across a rock mass with one fault and the induced ground motion using a recursive approach. The rocks beside the fault are assumed as viscoelastic media with seismic quality factors, Q_p and Q_s. Two kinds of interactions between stress waves and a discontinuity and between stress waves and a free surface are analyzed, respectively. As the result of the wave superposition, the mathematical expressions for induced ground vibration are deduced. The proposed approach is then compared with the existing analysis for special cases. Finally, parametric studies are carried out, which includes the influences of fault stiffness, incident angle, and frequency of incident waves on the peak particle velocities of the ground motions.     

Transient‐state groundwater flow in various geometries with wells at arbitrary positions: analytical methods

We have developed a method for analytically solving the porous medium flow equation in many different geometries for horizontal (two‐dimensional), homogeneous and isotropic aquifers containing impermeable boundaries and any number of pumping or injection wells located at arbitrary positions within the system. Solutions and results are presented for rectangular and circular aquifers but the method presented here is easily extendible to many geometries. Results are also presented for systems where constant head boundary conditions can be emulated internal to the aquifer boundary. Recommendations for extensions of the present work are briefly discussed. Copyright © 2003 John Wiley & Sons, Ltd.     

Metamorphism in Archaean greenstone belts: calculated fluid compositions and implications for gold mineralization

An inescapable consequence of the metamorphism of greenstone belt sequences is the release of a large volume of metamorphic fluid of low salinity with chemical characteristics controlled by the mineral assemblages involved in the devolatilization reactions. For mafic and ultramafic sequences, the composition of fluids released at upper greenschist to lower amphibolite facies conditions for the necessary relatively hot geotherm corresponds to those inferred for greenstone gold deposits (X_CO2= 0.2–0.3). This result follows from the calculation of mineral equilibria in the model system CaO–MgO–FeO–Al₂O₃–SiO₂–H₂O–CO₂, using a new, expanded, internally consistent dataset. Greenstone metamorphism cannot have involved much crustal over-thickening, because very shallow levels of greenstone belts are preserved. Such orogeny can be accounted for if compressive deformation of the crust is accompanied by thinning of the mantle lithosphere. In this case, the observed metamorphism, which was contemporaneous with deformation, is of the low-P high-T type. For this type of metamorphism, the metamorphic peak should have occurred earlier at deeper levels in the crust; i.e. the piezothermal array should be of the ‘deeper-earlier’type. However, at shallow crustal levels, the piezothermal array is likely to have been of ‘deeper-later’type, as a consequence of erosion. Thus, while the lower crust reached maximum temperatures, and partially melted to produce the observed granites, mid-crustal levels were releasing fluids prograde into shallow crustal levels that were already retrograde. We propose that these fluids are responsible for the gold mineralization. Thus, the contemporaneity of igneous activity and gold mineralization is a natural consequence of the thermal evolution, and does not mean that the mineralization has to be a consequence of igneous processes. Upward migration of metamorphic fluid, via appropriate structurally controlled pathways, will bring the fluid into contact with mineral assemblages that have equilibrated with a fluid with significantly lower X_CO2. These assemblages are therefore grossly out of equilibrium with the fluid. In the case of infiltrated metabasic rocks, intense carbonation and sulphidation is predicted. If, as seems reasonable, gold is mobilized by the fluid generated by devolatilization, then the combination of processes proposed, most of which are an inevitable consequence of the metamorphism, leads to the formation of greenstone gold deposits predominantly from metamorphic fluids.     

An Early Investigation Of The Striated Tail Of Comet Hale-Bopp (C/1995 O1)

The Sekanina-Farrell particle fragmentation model for the striated tails of dust comets is successfully applied to two images of comet Hale-Bopp to study the motions of 12 striae in a time span of March 12–15, 1997. There is evidence for recurring outbursts with a periodicity of 11^h21^m, consistent with results based on analysis of dust jets. The ejecta in all the striae appear to have been released from one source on the nucleus between the end of January and the second half of February 1997, some 60 to 40 days before perihelion. The parent particles were subjected to a radiation pressure acceleration of β_p ≃ 0.55 and their fragmentation lifetimes in 11 of the 12 striae were practically constant and equal to 13–15 days, when normalized to 1 AU from the Sun. Brief analysis of Watanabe et al.'s measurements of striae on their images from March 5–9, 1997 shows even shorter fragmentation lifetimes for the parent particles, mostly about 7–11 days at1 AU. This revised version was published online in July 2006 with corrections to the Cover Date.     

青藏块体东北缘地壳水平运动状态

应用青藏块体东北缘1999～2003年多期GPS观测资料，计算了不同时段GPS点水平运动速率。通过分析发现：甘青块体可分为东部块体和西部块体，东、西部块体的运动状态存在明显的差异；受2001年11月14日昆仑Ms8．1地震的影响，震后地壳运动状态发生了明显的改变。     

Structures in the phase space of a four dimensional symplectic map

We numerically investigate the projections of non periodic orbits in a 4-dimensional (4-D) symplectic map composed of two coupled 2-dimensional (2-D) maps. We describe in detail the structures that are produced in different planes of projection and we find how the morphology of the 4-D orbits is influenced by the features of the 2-D maps as the coupling parameter increases. We give an empirical law that describes this influence.     

Numerical prediction of blast‐induced stress wave from large‐scale underground explosion

This paper presents a numerical model for predicting the dynamic response of rock mass subjected to large‐scale underground explosion. The model is calibrated against data obtained from large‐scale field tests. The Hugoniot equation of state for rock mass is adopted to calculate the pressure as a function of mass density. A piecewise linear Drucker–Prager strength criterion including the strain rate effect is employed to model the rock mass behaviour subjected to blast loading. A double scalar damage model accounting for both the compression and tension damage is introduced to simulate the damage zone around the charge chamber caused by blast loading. The model is incorporated into Autodyn3D through its user subroutines. The numerical model is then used to predict the dynamic response of rock mass, in terms of the peak particle velocity (PPV) and peak particle acceleration (PPA) attenuation laws, the damage zone, the particle velocity time histories and their frequency contents for large‐scale underground explosion tests. The computed results are found in good agreement with the field measured data; hence, the proposed model is proven to be adequate for simulating the dynamic response of rock mass subjected to large‐scale underground explosion. Extended numerical analyses indicate that, apart from the charge loading density, the stress wave intensity is also affected, but to a lesser extent, by the charge weight and the charge chamber geometry for large‐scale underground explosions. Copyright © 2004 John Wiley & Sons, Ltd.     

激光天文动力学弱光锁相的研究

近年来激光物理与应用的进展促成了激光天文动力学空间任务概念的提出。在激光天文动力学任务概念研究方面 ,必须对由远程航天器上传回的激光进行讯号的测量与处理。激光经过长距离的传输后 ,功率大幅下降 ,因此在空间任务概念的考虑上 ,弱光锁相是计划中关键性的技术。由航天器射来的激光 ,经过望远镜聚光后与本地激光进行锁相 ,由本地激光承载及传达太空激光的相位信息。实验中 ,我们使用 2支半导体激光泵浦非平面环形共振腔钇镏石激光 (Diodelaserpumpednon -planarringcavityNd :YAGlaser) ,分别代表远程的弱光及代表本地的强光 ,建立弱光锁相环路系统 (weak -lightphase -lockedloopsystem)。以中性光强度滤光器 (ND -filter;neutraldensityfilter)减弱光讯号来仿真远程激光传来的弱光。在相位探测部分使用均衡探测法 ,消除激光强度噪声 ,以提高讯噪比。同时配合适当的环路滤波器 ,控制激光频率 ,提高锁相的能力。对 2nW的弱光与 2mW的强光可长时间锁相 ,其均方根相位误差为 57mrad。 2 0 0pW的弱光与 2 0 0 μW的强光锁相时间可达 2h以上 ,其相位误差为 2 0 0mrad。 2 0pW的弱光与 2 0 0 μW的强光锁相时间亦可达 2h以上 ,其相位误差为 1 60mrad。最后 ,我们对 2pW的弱光与 2 0 0 μW的强光锁相 ,