与花岗岩有关锡矿成岩成矿作用研究若干问题讨论

文章基于近年来大量典型锡矿床和成矿区带成岩成矿作用及高温高压实验研究成果,对锡矿成矿年代学、成矿花岗岩类型、物质来源、源区物质的熔融条件及锡在岩浆中的存在形式、熔体-流体相间分配行为及热液流体中的存在形式等方面的研究进展及存在问题进行了系统综述,指出:①锡石原位U-Pb测年已成为直接测定锡成矿年龄的有效手段,近年来发展起来的多种矿石矿物U-Pb测年方法与传统蚀变矿物测年方法的结合,将为研究锡多金属矿床多期次成矿作用叠加、矿集区内不同类型矿床成因联系及区域矿床组合模型的建立提供关键的年代学依据;②锡的富集成矿不仅与花岗质岩浆的高程度的结晶分异密切相关,而且受源区组成及部分熔融条件的制约;③全球锡多金属矿床的不均匀分布、与准铝质花岗岩有关锡矿的富集成矿机制以及各类复杂的锡-钨-铜-钼-铅-锌银多金属矿种组合差异及其金属元素富集机制亟待查明;④已有的关于锡在花岗质熔体、热液流体相中的存在形式及两相间的分配行为的高温高压实验结果受锡与贵金属容器间合金效应的影响,人工合成流体包裹体等新的实验方法的开发及XAS同步辐射和激光拉曼等原位分析技术的应用,将为研究花岗岩浆-热液过程锡的迁移、富集及沉淀成矿机制提供更为直接、可靠的实验依据。     

A Discrete Wavenumber Boundary Element Method for study of the 3-D response 2-D scatterers

A mathematical formulation of the 2·5D elastodynamic scattering problem is presented and validated. The formulation is a straightforward extension of the Discrete Wave number Boundary Integral Equation Method (DWBIEM) originally proposed by Kawase¹ for 2D scattering problems and subsequently extended to the 3D problem by Kim and Papageorgiou.² It is demonstrated that the Green's function which is appropriate for a boundary formulation of the 2·5D elastodynamic scattering problem is the one corresponding to a unit force moving on a straight line with constant velocity. Such a Green's function is derived in the present study. The formulation may be used to study the wavefields in models of sedimentary deposits (e.g. valleys) or topography (e.g. canyons or ridges) with a 2D variation in structure but obliquely incident plane waves. The advantage of a 2·5D formulation is that it provides the means for calculations of 3D wavefields in scattering problems by requiring a storage comparable to that of the corresponding 2D calculations. © 1998 John Wiley & Sons, Ltd.     

Groundwater contaminant source identification based on iterative local update ensemble smoother

Identification of the location and intensity of groundwater pollution source contributes to the effect of pollution remediation, and is called groundwater contaminant source identifcation. This is a kind of typical groundwater inverse problem, and the solution is usually ill-posed. Especially considering the spatial variability of hydraulic conductivity field, the identification process is more challenging. In this paper, the solution framework of groundwater contaminant source identification is composed with groundwater pollutant transport model (MT3DMS) and a data assimilation method (Iterative local update ensemble smoother, ILUES). In addition, Karhunen-Loève expansion technique is adopted as a PCA method to realize dimension reduction. In practical problems, the geostatistical method is usually used to characterize the hydraulic conductivity feld, and only the contaminant source information is inversely calculated in the identifcation process. In this study, the identification of contaminant source information under Kriging K-field is compared with simultaneous identification of source information and K-field. The results indicate that it is necessary to carry out simultaneous identification under heterogeneous site, and ILUES has good performance in solving high-dimensional parameter inversion problems.     

一维连续导电剖面的反演问题

基于Schlumberger装置测得的视电阻率,在Coen理论的基础上对一维连续直流导电剖面的反演问题进行研究。分别给出了导电剖面是常函数、线性函数、幂函数、指数函数、分段函数情况下的精确数据和噪音数据反演结果。利用一个野外场地反演的例子,把利用该方法进行反演得到的连续模型和Inman得到的层状模型进行了相应的比较,发现该方法在处理实际场地数据也是有效的。     

Elastoplastic deformation around underground openings under biaxial initial stress field

Estimation of elastoplastic deformation around an underground opening induced by the excavation of it, especially displacement and strain field in plastic region, is presented in this paper, as well as the formulation for calculating the displacement and strain in the plastic region around the underground opening by the coupled Boundary Element Method - Characteristics Method (BEM-CM). In this method, the non-associated flow rule is adopted to calculate the displacement and strain field in the plastic region, which is determined by the integration of the displacement along characteristics lines under the boundary condition of the elastic displacement on an elastoplastic interface analysed. It is shown that this method is one of the accurate and effective methods for estimating not only the shape and extent of the plastic region but also the state of the displacement and strain in the plastic region around the underground opening, comparing the theoretical solution with numerical results by this method for a circular opening under hydrostatic initial stress condition. Furthermore, this method is applied to rectangular and horse-shoe shaped openings and the characteristics of the strain field in the plastic region are discussed. © 1997 by John Wiley & Sons, Ltd.     

Estimating free parameters in 2D inversion: example of gravity inversion in a rifted basement

We advance a principle directed to assigning numerical values to free parameters usually present in inversion methods. It may be formulated as: ‘Optimum estimates of free parameters in an inversion procedure must lead, in tests using synthetic data, to solutions whose geometrical expression reflects the main qualitative or semiquantitative geological characteristic of the study area.’ To this end, the interpreter should (i) specify a typical anomalous source geometry which incorporates the most relevant geological information for the study area, (ii) compute the corresponding gravity anomaly and (iii) invert the anomaly for the source geometry finding the numerical values of the free parameters that lead to a solution closest to the typical source. Application of the above methodology to synthetic and real data from the basement relief of a rift basin has asserted its efficacy.     

Estimation of groundwater storage from seismic data using deep learning

Convolutional neural networks can provide a potential framework to characterize groundwater storage from seismic data. Estimation of key components, such as the amount of groundwater stored in an aquifer and delineate water table level, from active-source seismic data are performed in this study. The data to train, validate and test the neural networks are obtained by solving wave propagation in a coupled poroviscoelastic–elastic media. A discontinuous Galerkin method is applied to model wave propagation, whereas a deep convolutional neural network is used for the parameter estimation problem. In the numerical experiment, the primary unknowns estimated are the amount of stored groundwater and water table level, while the remaining parameters, assumed to be of less of interest, are marginalized in the convolutional neural network-based solution. Results, obtained through synthetic data, illustrate the potential of deep learning methods to extract additional aquifer information from seismic data, which otherwise would be impossible based on a set of reflection seismic sections or velocity tomograms.     

Some remarks on Q-compensated sparse deconvolution without knowing the quality factor Q

The subsurface media are not perfectly elastic, thus anelastic absorption, attenuation and dispersion (aka Q filtering) effects occur during wave propagation, diminishing seismic resolution. Compensating for anelastic effects is imperative for resolution enhancement. Q values are required for most of conventional Q-compensation methods, and the source wavelet is additionally required for some of them. Based on the previous work of non-stationary sparse reflectivity inversion, we evaluate a series of methods for Q-compensation with/without knowing Q and with/without knowing wavelet. We demonstrate that if Q-compensation takes the wavelet into account, it generates better results for the severely attenuated components, benefiting from the sparsity promotion. We then evaluate a two-phase Q-compensation method in the frequency domain to eliminate Q requirement. In phase 1, the observed seismogram is disintegrated into the least number of Q-filtered wavelets chosen from a dictionary by optimizing a basis pursuit denoising problem, where the dictionary is composed of the known wavelet with different propagation times, each filtered with a range of possible values. The elements of the dictionary are weighted by the infinity norm of the corresponding column and further preconditioned to provide wavelets of different values and different propagation times equal probability to entry into the solution space. In phase 2, we derive analytic solutions for estimates of reflectivity and Q and solve an over-determined equation to obtain the final reflectivity series and Q values, where both the amplitude and phase information are utilized to estimate the Q values. The evaluated inversion-based Q estimation method handles the wave-interference effects better than conventional spectral-ratio-based methods. For Q-compensation, we investigate why sparsity promoting does matter. Numerical and field data experiments indicate the feasibility of the evaluated method of Q-compensation without knowing Q but with wavelet given.     

SEQUENTIAL STIFFNESS DESIGN FOR SEISMIC DRIFT RANGES OF A SHEAR BUILDING–PILE–SOIL SYSTEM

A shear building supported by a prescribed pile–soil system is subjected to bedrock earthquake input. A new design procedure is presented for generating a sequence of stiffness designs satisfying the constraints on interstorey drifts. The mean peak interstorey drifts of the shear building subjected to a set of spectrum-compatible ground motions at the bedrock are evaluated by a modal combination rule. Tuning of the fundamental natural period of a shear building with a fixed base with that of a shear beam ground results in a non-monotonic sequence of stiffness designs with respect to a ground stiffness parameter and previous approaches cannot be applied to such a problem. This difficulty in finding such a non-monotonic sequence is overcome by utilizing the ground stiffness parameter and the superstructure stiffness parameter alternately in multiple design phases and by developing a new multi-phase perturbation technique. Fundamental characteristics of this sequence of stiffness designs and the effect of ground stiffnesses on the design of the shear building are disclosed. It is further shown that the stiffness contour method is also useful for the design procedure such that a scattering effect in the estimates of ground stiffnesses is taken into account. The usefulness of the proposed procedure of sequential stiffness design and contour line method is demonstrated through several sequential design examples.     

THE SETTLEMENT OF A RIGID CIRCULAR FOUNDATION RESTING ON A HALF-SPACE EXHIBITING A NEAR SURFACE ELASTIC NON-HOMOGENEITY

The present paper examines the elastostatic problem pertaining to the axisymmetric loading of a rigid circular foundation resting on the surface of a non-homogeneous elastic half-space. The non-homogeneity corresponds to a depth variation in the linear elastic shear modulus according to the exponential form G(z)=G₁+G₂e^-ζz. The equations of elasticity governing this type of non-homogeneity are solved by employing a Hankel transform technique. The mixed boundary value problem associated with the indentation of the half-space by the rigid circular foundation is reduced to a Fredholm integral equation which is solved via a numerical technique. The numerical results presented in the paper illustrate the influence of the near-surface elastic non-homogeneity on the settlement of the foundation.