工程地质力学综合集成理论及其在五强溪水电站船闸边坡上的应用

为了使工程地质学在工程建设和防灾减灾方面发挥更大的作用，作者主张在工程地质力学和系统科学指导下发展工程地质力学综合集成（EGMS）。本文主要论述了EGMS形成的基础、相应的方法和技术以及综合集成中的4个关键问题（包括应变体系、优先监测和加固突破口的设计方法，综合地质信息分析系统、时空综合分析方法等）。另外，还扼要地讨论了如何将EGMS用于五强溪水电站船闸边坡的问题。     

The phytoplankton of some saline lakes in Central Asia

The Aral Sea, Lake Balkhash, and Lake Kamyslybas are closed lakes in Central Asia. They range from oligosaline to metasaline. The salinity of the Aral Sea has increased by more than 30 g L⁻¹ since widespread irrigation began in its catchment area. Few studies of the phytoplankton have been conducted on these lakes since extensive irrigation started. The investigation reported here compares the flora of phytoplankton in these saline lakes. In the Small Aral Sea, phytoplankton density gradually decreased with increasing electrical conductivity (EC) (∼ salinity), but there was no such relation in Lake Balkhash and Lake Kamyslybas. In the Aral Sea, Dinophyceae and Bacillariophyceae were frequently observed in most areas of high EC value, and Cyanophyceae were most conspicuous in the area of medium and lower EC values. In Lake Balkhash, Cyanophyceae were most conspicuous, but Chlorophyceae were also noticeable. Most Cyanophyceae in Aral Sea formed filaments with heterocysts. The distinct characteristic of the phytoplankton of the Lake Balkhash was that all dominant species form colonies covered with a gelatinous film. Siliceousplankton diversity gradually decreased with increasing EC values in the Aral Sea and Lake Balkhash.     

The effects of the spinel—garnet phase transition on the formation of rifted sedimentary basins

We have examined the effects of the spinel-garnet phase transition on subsidence of extensional sedimentary basins. For a constant positive Clapeyron slope ( dP/dT ), the phase boundary moves downwards in the syn-rift and upwards in the post-rift phase. For a non-linear Clapeyron curve ( dP/dT > 0 above 900°C and dP/dT < 0 below 900°C), theory predicts for the reaction of the spinel-garnet phase transition, the direction of phase boundary movement is dependent on the stretching factor, the position of the Clapeyron curve and the lithospheric thickness. A smaller syn-rift and larger post-rift subsidence are predicted for a deeper phase boundary and a thicker lithosphere. The model with a non-linear Clapeyron curve is applied to the subsidence histories of a young extensional basin (Gulf of Lion) and an old continental margin (eastern Canada). The observed syn-rift uplift and the larger post-rift subsidence can be reasonably explained by this model, where the optimum depth of the phase boundary for eastern Canada (˜90 km) is consistent with the estimate from seismic observations and is larger than that for the Gulf of Lion (˜ 50 km). The depth of the spinel-garnet phase boundary is sensitive to the composition of mantle rocks and increases with the extraction of basaltic components from the lithosphere, compatible with our result that the phase boundary is deeper for an older and thicker lithosphere. Thus the surface movement associated with the rifting for these areas may reflect the chemical evolution of the continental lithosphere.     

Shear-wave anisotropy: spatial and temporal variations in time delays at Parkfield, Central California

Shear-wave splitting is analysed on data recorded by the High Resolution Seismic Network (HRSN) at Parkfield on the San Andreas fault, Central California, during the three-year period 1988-1990. Shear-wave polarizations either side of the fault are generally aligned in directions consistent with the regional horizontal maximum compressive stress, at some 70° to the fault strike, whereas at station MM in the immediate fault zone, shear-wave polarizations are aligned approximately parallel to the fault. Normalized time delays at this station are found to be about twice as large as those in the rock mass either side. This suggests that fluid-filled cracks and fractures within the fault zone are elastically or seismically different from those in the surrounding rocks, and that the alignment of fault-parallel shear-wave polarizations are associated with some fault-specific phenomenon.
Temporal variations in time delays between the two split shear-waves before and after a M_L = 4 earthquake can be identified at two stations with sufficient data: MM within the fault zone and VC outside the immediate fault zone. Time delays between faster and slower split shear waves increase before the M_L = 4 earthquake and decrease near the time of the event. The temporal variations are statistically significant at 68 per cent confidence levels. Earthquake doublets and multiplets also show similar temporal variations, consistent with those predicted by anisotropic poroelasticity theory for stress modifications to the microcrack geometry pervading the rock mass. This study is broadly consistent with the behaviour observed before three other earthquakes, suggesting that the build-up of stress before earthquakes may be monitored and interpreted by the analysis of shear-wave splitting.     

Periodicity of magnetic intensities in magnetic anomaly profiles: the Cenozoic of the South Atlantic

Spectral analyses of several published magnetic anomaly profiles from Candé & Kent (1992a) were undertaken prior to analysing, in the same way, raw magnetic anomaly data from similar parts of the South Atlantic. It was found that similar and distinct medium and short wavelengths were present in both the published and raw data. When these are converted into the time domain using the average rate of spreading for each profile, these periodicities appear similar, possibly identical, to those expected from the long-term eccentricity orbital parameters (Fischer, DeBoer & Premoli Silva 1990). While such correlations are not necessarily causative, they suggest that magnetohydro-dynamical processes near the core-mantle boundary may be affected by gravitational changes due to planetary orbital perturbations.     

Three-dimensional Rayleigh hysteresis of oriented core samples from the German Continental Deep Drilling Program: susceptibility tensor, Rayleigh tensor, three-dimensional Rayleigh law

Rayleigh hysteresis, as defined by the well-known Rayleigh relations, has been observed not only when magnetization of pyrrhotite-bearing KTB-samples is measured in parallel to a weak dc magnetic field, but also in experiments where field and measuring directions have been adjusted strictly perpendicularly to each other. Nine-tupels of independent Rayleigh hysteresis loops could thus be compiled. Their characteristic coefficients X ^ijk of initial susceptibility together with the Rayleigh loss coefficients α^jk have been proved to determine completely the samples' weak-field magnetic anisotropy. Interpreting the coefficient matrices ( X ^ijk) and (α^jk) as the tensor of initial susceptibility and the Rayleigh tensor, respectively, generalization of the isotropic Rayleigh relations in terms of corresponding tensor relationships has been suggested for the anisotropic case. Application to the KTB samples showed 3-D Rayleigh hysteresis measurements to be an excellent tool for rock magnetic analysis in terms of ore content and crystalline texture. In particular, a magnetocrystalline double texture of the basal planes of pyrrhotite precipitates and their [1120] directions of easy magnetization have been clearly detected. Surprisingly, the welt-known theorem α= const. X ²_I, formulated by Néel (1942) for the isotropic case, has been found to hold true even in tensor generalization (α_jk) = const ( X ²_jk). To reach sufficient resolution for the measurements performed, a sensitive vibrating coil magnetometer (VCM) has been developed.     

The S reflector west of Galicia: the seismic signature of a detachment fault

The so-called S reflector is a group of bright, continuous reflections underlying the landward-tilted fault blocks of the west Galicia rifted margin, S has been interpreted as the brittle-ductile transition, the lop of an intrusion, a detachment fault, and the crust-mantle boundary. To constrain the internal structure of the reflector, we have carried out detailed analyses of these reflections. We compare the waveforms of the seafloor reflection and its first multiple, both to determine the amplitude of the seafloor reflection and to show (hat the seafloor is in effect a spike in the reflectivity series so that the seafloor reflection can be used as the far-field wavelet, including both source and receiver ghosts. We compare (he waveform of the seafloor and 5 and show that, within the resolution of our data, S is a reflection from a step increase in acoustic impedance. This result is confirmed through complex trace analysis, and in particular the determination of the apparent polarity of S, and the examination of the instantaneous frequency function: S is consistently positive polarity, and shows no significant frequency anomaly. Simple modelling shows that S is very unlikely to be a reflection from a thin layer. We thus conclude that S is probably a single steplike interface. From the varying frequency content of the data, we determine a value for the effective Q between S and the seafloor, and use this to assess the amplitude loss due to attenuation and scattering. We use a comparison between the seafloor and the S reflection to constrain the amplitude of S, and estimate a reflection coefficient for S of at least 0.2 in places, decreasing landwards. By analogy with structures developed in the highly extended regions of the western United States, we consider that the most likely interpretation of S is as a sharp west-dipping detachment fault separating a 'granitic' upper plate from a higher-velocity lower plate, locally probably serpentinized mantle.     

Comment on 'Estimation of resolution and covariance for large matrix inversions' by J. Zhang and G. A. McMechan

When inverting large matrices, iterative techniques are necessary because of their speed and low memory requirements, as opposed to singular value decomposition (SVD). Recently, there have been attempts to obtain information on the quality of the solutions calculated using conjugate gradient (CG) methods such as LSQR. The purpose of this note is to comment on the paper titled "Estimation of resolution and covariance for large matrix inversions' by Zhang & McMechan (1995), who extend Paige and Saunders' LSQR algorithm to obtain an orthonormal basis used to approximate resolution and covariance. We show that for larger problems, where the number of orthogonal vectors is several orders of magnitude smaller than the number of model parameters, the vectors obtained do not adequately span the range of the model space. We use a synthetic borehole experiment to illustrate the differences between the singular value spectrum obtained through the more complete method of SVD and the Ritz value spectrum that results from a simple extension of LSQR, We also present a trivial numerical example to illustrate the differences between Zhang & McMechan's approximate resolution matrix and the true resolution.