基于GM(1,1)优化模型的岩石边坡变形预报

岩石边坡系统是一典型灰色系统，其变形发展过程可以用灰色预测模型完成。尽管传统GM（1，1）模型预测有很多成功的实例，但是也存在一些预测偏差过大的情况，必须对其进行优化。逐步迭代法GM（1，1）模型不仅收敛速度快，而且与原始数据库序列的凹凸性保护一致。利用自编的计算程序对马步坎边坡预测测点G1沉降和开裂进行预测分析，结果表明逐步迭代法GM（1，1）优化模型计算精度较传统GM（1，1）模型和背景构造法GM（1，1）优化模型高，较好地反映了岩石边坡的变形趋势。     

福建福安下白石特大桥海上岩土勘察方法

下白石特大桥是同江至三亚国道主干线福鼎至宁德高速公路建设中的一项控制性工程 ,是交通部指定的三阶段设计工程 ,桥长 110 5 m,主跨采用 2× 2 6 0 m预应力混凝土连续钢构。桥位位于白马河入海口处 ,工程地质、水文地质条件复杂 ,通过采用地表大面积的综合地质调查、海上地震映像物探以及海上钻探、原位测试、岩土测试等一系列行之有效的岩土勘察方法 ,确保了岩土工程勘察质量 ,其结果切实可行 ,为后期施工单位提供科学依据。     

FLAC数值模拟软件及其在地学中的应用

FLAC是以拉格朗日差分法为基础的一种数值模拟应用软件,主要用于岩土力学工程。近年地质学家将其应用到地质构造和块体运动学、动力学研究中,以证实或证伪地球学家们提出的相关地球动力学模型。本文初步介绍了FLAC应用软件的基本原理、使用方法以及它在国内外地质学上的几个应用实例。     

铅锌矿床Rb-Sr定年研究综

本文在总结Rb-Sr等时线定年的基本前提和矿床Rb-Sr等时线定年存在的问题的基础上,综述了铅锌矿床流体包裹体和闪锌矿Rb-Sr等时线定年的研究进展,同时对该定年方法存在的问题进行了简要评述。综合分析表明,Rb-Sr等时线定年(尤其是闪锌矿Rb-Sr等时线定年)为铅锌矿床较理想的直接定年方法,但该定年方法还存在许多理论和技术问题,有待深入研究;铅锌矿床Rb-Sr等时线定年应以深入细致的地质、地球化学研究为基础,定年过程中必须严慎对待样品采集、样品挑选、流体包裹体溶液提取和质谱分析等各个环节。     

遗传算法在边坡数值计算中的应用

改进了进化方向的遗传算法与有限元数值法，结合并研制了相应的软件。应用该软件对多类型岩土边坡进行弹性模量、内聚力、内摩擦角等参数反演分析，显示误差很小，收敛速度也很快，这说明改进进化方向遗传算法这种新型的优化算法在多类型岩土参数优化估计中具有独特的优势。     

软土地基沉降计算的误差分析

对引起路堤荷载下软土地基沉降从竖应力、固结沉降、瞬时沉降和后期沉降的计算四个方面,常见计算方法误差的诸多因素,从受力、沉降原理上进行了详细分析并提出一些建议,以使沉降更为接近实际.     

Analogue study of particle segregation in pyroclastic density currents, with implications for the emplacement mechanisms of large ignimbrites

Abstract Analogue flume experiments were conducted to investigate the transport and sedimentation behaviour of turbulent pyroclastic density currents. The experimental currents were scaled approximately to the natural environment in three ways: (1) they were fully turbulent; (2) they had a very wide range of particle sizes and associated Rouse numbers (the ratio of particle settling velocity to effective turbulent eddy velocity in the current); and (3) they contained particles of two different densities. Two sets of surge‐type experiments were conducted in a 5 m long, water‐filled lock‐exchange flume at five different volumetric particle concentrations from 0·6% to 23%. In one set (one‐component experiments), the currents contained just dense particles; in the other set (two‐component experiments), they contained both light and dense particles in equal volume proportions. In both sets of experiments, the population of each component had a log‐normal size distribution. In the two‐component experiments, the size range of the light particle population was selected in order to be in hydrodynamic equivalence with that of the dense particles. Dense particles were normally graded, both vertically and downstream, in the deposits from both sets of experiments. The mass loading (normalized to the initial mass of the suspension) and grain size of the dense component in the deposits decreased with distance from the reservoir and were insensitive to initial total particle concentration in the currents. On the other hand, in the two‐component experiments, the light particles were extremely sensitive to concentration. They were deposited in hydrodynamic equivalence with the dense particles from dilute currents, but were segregated efficiently at concentrations higher than a few per cent. With increasing particle concentration, the large, light particles were carried progressively further down the flume because of buoyancy effects. Deposits from the high‐concentration currents exhibited reverse vertical grading of the large, light particles. Efficient segregation of the light component was observed even if the bulk density of the current was less than that of the light particles. In both sets of experiments, marked inflexions in the rate of downstream decline in mass loading and maximum grain size of the dense component can be attributed to the presence of two different particle settling regimes in the flow: (1) particles with Rouse numbers >2·5, which did not respond to the turbulence and settled rapidly; and (2) particles with Rouse numbers <2·5, which followed the turbulent eddies and settled slowly. The results are applied to the transport and sedimentation dynamics of pyroclastic density currents that generate large, widespread ignimbrites. Field data fail to reveal significant departures from aerodynamic equivalence between pumice and lithic clasts in three such ignimbrites: the particulate loads of some large ignimbrites are transported principally in turbulent suspensions of low concentration. In some ignimbrites, the well‐developed inflexions in curves of maximum lithic (ML) size vs. distance can be attributed to the existence of distinct high and low Rouse number particle settling regimes that mark the transition from an overcharged state to one in which the residual particulate load is transported more effectively by turbulence.     

Crustal structure of the Ruby Mountains and southern Carlin Trend region, Nevada, from magnetotelluric data

We have collected about 150 magnetotelluric (MT) soundings in northeastern Nevada in the region of the Ruby Mountains metamorphic core complex uplift and southern Carlin mineral trend, in an effort to illuminate controls on core complex evolution and deposition of world-class gold deposits. The region has experienced a broad range of tectonic events including several periods of compressional and extensional deformation, which have contributed to the total expression of electrical resistivity. Most of the soundings reside in three east–west profiles across increasing degrees of core uplift to the north (Bald Mountain, Harrison Pass, and Secret Pass latitudes). One short cross-line was also taken to assess an east–west structure to the north of the northern profile. Model resistivity cross-sections were derived from the MT data using a 2-D inversion algorithm, which damps departures of model parameters from an a priori structure. Geological interpretation of the resistivity combines previous seismic, potential field and isotope models, structural and petrological models for regional compression and extension, and detailed structural/stratigraphic interpretations incorporating drilling for petroleum and mineral exploration. To first order, the resistivity structure is one of a moderately conductive, Phanerozoic sedimentary section fundamentally disrupted by intrusion and uplift of resistive crystalline rocks. Late Devonian and early Mississippian shales of the Pilot and Chainman Formations together form an important conductive marker sequence in the stratigraphy and show pronounced increases in conductance (conductivity–thickness product) from east to west. These increases are attributed to graphitization caused by Elko–Sevier era compressional shear deformation and possibly by intrusive heating. The resistive crystalline central massifs adjoin the host stratigraphy across crustal-scale, steeply dipping fault zones. The zones provide pathways to the lower crust for heterogeneous, upper crustal induced, electric current flow. Resistive core complex crust appears steeply bounded under the middle of the neighboring grabens and not to deepen at a shallow angle to arbitrary distances to the west. The numerous crustal breaks imaged with MT may contribute to the low effective elastic thickness (T_e) estimated regionally for the Great Basin and exemplify the mid-crustal, steeply dipping slip zones in which major earthquakes nucleate. An east–west oriented conductor in the crystalline upper crust spans the East Humboldt Range and northern Ruby Mountains. The conductor may be related to nearby graphitic metasediments, with possible alteration by middle Tertiary magmatism. Lower crustal resistivity everywhere under the profiles is low and appears quasi one-dimensional. It is consistent with a low rock porosity (<1 vol.%) containing hypersaline brines and possible water-undersaturated crustal melts, residual to the mostly Miocene regional extension. The resistivity expression of the southern Carlin Trend (CT) in the Pinon Range is not a simple lineament but rather a family of structures attributed to Eocene intrusion, stratal deformation, and alteration/graphitization. Substantial reactivation or overprinting by core complex uplift or Basin–Range extensional events seems likely. We concur with others that the Carlin Trend may result in part from overlap of the large Eocene Northeast Nevada Volcanic Field with Precambrian–Paleozoic deep-water clastic source rocks thickening abruptly to the west of the Pinon Range, and projecting to the north–northwest.     

Probability Density Functions as Botanical-Climatological Transfer Functions for Climate Reconstruction

We present a new procedure, the pdf method (pdf=probability density function), for reconstructing Quaternary climate utilizing botanical data. The procedure includes the advantages of the indicator species method by considering the fossil and modern presence and absence of taxa rather than their frequencies, thus avoiding the need for modern analog plant communities. Overcoming the problematic use of absolute limits to describe climate response ranges is the main progress of the pdf method in comparison to the indicator species method. This advantage results from estimating probability density functions (pdfs) for monthly mean January and July temperature conditional on the present day occurrence of single taxa. Gaussian distributions sufficiently approximate pdfs of many, although not all, studied taxa. On the assumption of statistical independence, the procedure calculates a joint pdf as the product of the pdfs of the individual taxa. This algorithm weights each taxon according to the extent of its climate response range expressed by its covariance structure. We interpret the maximum of the resulting pdf as the most likely climate and its confidence interval as the uncertainty range. To avoid an artificial reduction of uncertainty arising from the use of numerous similar pdfs, a preselection method is proposed based on the Mahalanobis distance between pdfs. The pdf method was applied to the Carpinus phase of a profile from Gröbern, Germany, that spans the last interglaciation (Eemian). The reconstructed most probable January and July temperatures of about 0.0°C and 18.4°C barely differ from the modern values of −0.5°C and 18.3°C.     

有限弹性变形地壳中的地震波

地壳由半无限大的基岩上一层厚度为H^-的表土层组成，入射地震波为垂直的SH波，产生水平地面运动。当浅源大地震发生时，在极震区以外行波传播产生地面运动将使地壳介质有非线性的有限弹性变形。用小参数摄动法使非线性控制方程为线性化的小参数各阶控制方程，得出头两阶线性控制方程的解析解。