岩体水力学基础(三)──岩体渗流场与应力场耦合的集中参数模型及连续介质模型

本文在分析了岩体系统的结构性和水力学特征后，提出了岩体渗流场与应力场耦合数学模型的机理分析、混合分析及系统辨识建模方法。运用系统辨识方法建立了岩体渗流场与应力场耦合的集中参数模型，并应用于解决实际问题；运用机理分析和混合分析方法建立了岩体渗流场与应力场双场耦合及与温度场三场耦合的连续介质分布参数模型。     

Inferring the relative measure of principal stress components

The phase velocity and the attenuation coefficient of compressional seismic waves, propagating in poroelastic, fluid-saturated, laminated sediments, are computed analytically from first principles. The wavefield is found to be strongly affected by the medium heterogeneity. Impedance fluctuations lead to poroelastic scattering; variations of the layer compressibilities cause inter-layer flow (a 1-D macroscopic local flow). These effects result in significant attenuation and dispersion of the seismic wavefield, even in the surface seismic frequency range, 10–100 Hz. The various attenuation mechanisms are found to be approximately additive, dominated by inter-layer flow at very low frequencies. Elastic scattering is important over a broad frequency range from seismic to sonic frequencies. Biot's global flow (the relative displacement of solid frame and fluid) contributes mainly in the range of ultrasonic frequencies. From the seismic frequency range up to ultrasonic frequencies, attenuation due to heterogeneity is strongly enhanced compared to homogeneous Biot models. Simple analytical expressions for the P -wave phase velocity and attenuation coefficient are presented as functions of frequency and of statistical medium parameters (correlation lengths, variances). These results automatically include different asymptotic approximations, such as poroelastic Backus averaging in the quasi-static and the no-flow limits, geometrical optics, and intermediate frequency ranges.     

Wide-angle seismic velocities in heterogeneous crust

Seismic velocities measured by wide-angle surveys are commonly used to constrain material composition in the deep crust. Therefore, it is important to understand how these velocities are affected by the presence of multiscale heterogeneities. The effects may be characterised by the scale of the heterogeneity relative to the dominant seismic wavelength (λ); what is clear is that heterogeneities of all scales and strengths bias wide-angle velocities to some degree. Waveform modelling was used to investigate the apparent wide-angle P -wave velocities of different heterogeneous lower crusts. A constant composition (50 per cent felsic and 50 per cent ultramafic) was formed into a variety of 1- and 2-D heterogeneous arrangements and the resulting wide-angle seismic velocity was estimated. Elastic, 1-D models produced the largest velocity shift relative to the true average velocity of the medium (which is the velocity of an isotropic mixture of the two components). Thick (width > λ) horizontal layers, as a result of Fermat's Principle, provided the largest increase in velocity; thin (width ≪λ) vertical layers produced the largest decrease in velocity. Acoustic 2-D algorithms were shown to be inadequate for modelling the kinematics of waves in bodies with multiscale heterogeneities. Elastic, 2-D modelling found velocity shifts (both positive and negative) that were of a smaller magnitude than those produced by 1-D models. The key to the magnitude of the velocity shift appears to be the connectivity of the fast (and/or slow) components. Thus, the models with the highest apparent levels of connectivity between the fast phases, the 1-D layers, produced the highest-magnitude velocity shifts. To understand the relationship between measured seismic velocities and petrology in the deep crust it is clear that high-resolution structural information (which describes such connectivity) must be included in any modelling.     

P-SH conversions in a flat-layered medium with anisotropy of arbitrary orientation

P-SH conversion is commonly observed in teleseismic P waves, and is often attributed to dipping interfaces beneath the receiver. Our modelling suggests an alternative explanation in terms of flat-layered anisotropy. We use reflectivity techniques to compute three-component synthetic seismograms in a 1-D anisotropic layered medium. For each layer of the medium, we prescribe values of seismic velocities and hexagonally symmetric anisotropy about a common symmetry axis of arbitrary orientation. A compressional wave in an anisotropic velocity structure suffers conversion to both SV -and SH -polarized shear waves, unless the axis of symmetry is everywhere vertical or the wave travels parallel to all symmetry axes. The P-SV conversion forms the basis of the widely used 'receiver function' technique. The P-SH conversion occurs at interfaces where one or both layers are anisotropic. A tilted axis of symmetry and a dipping interface in isotropic media produce similar amplitudes of both direct ( P ) and converted ( Ps ) phases, leaving the backazimuth variation of the P-Ps delay as the main discriminant. Seismic anisotropy with a tilted symmetry axis leads to complex synthetic seismograms in velocity models composed of just a few flat homogeneous layers. It is possible therefore to model observations of P coda with prominent transverse components with relatively simple 1-D velocity structures. Successful retrieval of salient model characteristics appears possible using multiple realizations of a genetic-algorithm (GA) inversion of P coda from several backazimuths. Using GA inversion, we determine that six P coda recorded at station ARU in central Russia are consistent with models that possess strong (> 10 per cent) anisotropy in the top 5 km and between 30 and 43 km depth. The symmetry axes are tilted, and appear aligned with the seismic anisotropy orientation in the mantle under ARU suggested by SKS splitting.     

渗流场反分析及其在多孔介质中的应用

本文讨论了复杂地质环境中渗流场的反分析问题，利用地下水流数值模拟结合最优化技术，通过对地下渗流场的模拟和寻优，求得描述地下含水介质非稳定渗流场的重要参位──渗透系数和给水度的模型解，并应用于各向同性的多孔介质渗流场反分析实例中。     

多媒体创作工具Compel的应用

本文介绍了一种多媒体应用的创作工具。利用这种工具可以编辑和组织文字、图象、图表、声音以及视频等多种媒体元素。并通过媒体链接以及超级链接，完成多媒体应用的制作。同时，也介绍了多媒体应用的播放。     

A coherency function model of ground motion at base rock corresponding to strike-slip fault

IntroductionEarthquakedamagesurveyandresearchresultshavedemonstratedthatspatialdistributiondifferenceofgroundmotionisoneoftheimportantreasonswhichcausedlongstructure(eglongspanbridge,undergroundpipe)destroy.Thathowtoprovideareasonableinputofgroundmotionfieldforaseismicdesignoflongstructureisaurgentprobleminearthquakeengineeringfield.Atpresent,themethodtostudyspatialvariationofgroundmotionsisadoptingstatisticanalysisbasedondensearrayrecordssuchasSMART-1array,etc,togetcoherencyfunctionofground…     

Dynamic response of partially sealed circular tunnel in viscoelastic saturated soil

Based on Biot's wave equation, dynamic response of a circular tunnel with partially sealed liner in viscoelastic saturated soil is investigated. By introducing two scalar potential functions, the analytical solutions of stresses, displacements and pore pressure induced by axisymmetric gradually applied step load are derived in Laplace transform domain. Numerical results are obtained by inverting Laplace transform presented by Durbin and used to analyze the influences of partial permeable property of boundary and viscoelastic damping coefficient of soil on dynamic response of the tunnel. It is shown that the attenuation of radial displacement appeared with the increase of viscoelastic damping coefficient of soil, and relative rigidity of liner and soil, and the influence of partial sealing property of boundary on stresses, displacements and pore pressure is remarkable. The available solutions of permeable and impermeable boundary conditions are only two extreme cases of this paper.     

多孔介质非饱和导水率预测的分形模型

多孔介质非饱和导水率是地下水污染预测与评价的重要参数。根据分形几何的基本原理和方法,推导出了与Campbell经验公式在形式上完全一致的多孔介质非饱和导水率的预测公式。公式中的幂指数为介质孔隙分维和随机行走分维的函数,分别体现了多孔介质的静态性质与动态性质对其中水分运动的影响,但静态性质的影响是主要的,即导水率主要受多孔介质的结构控制。根据文献中报道的大量数据,利用笔者推导的预测公式计算得到的幂指数的统计值与试验测定的幂指数的统计值基本一致,说明推导的理论公式预测多孔介质非饱和导水率是较为可靠的。     

Nuclear Waste Disposal Simulations: Couplex Test Cases

We give some results obtained for the Couplex test cases proposed by the ANDRA. In this paper our aim is twofold. Firstly, to compute the release of nuclides out of the repository by concentrating on the 3D near field (Couplex 2). The simulation of the transport phenomena takes into account the dissolution of the glass containers and congruent emissions of the radio-nuclides including filiation chains and some simplified chemistry. Secondly, it is to use the near field computations in order to simulate the nuclide migrations in a 2D far field (Couplex 3). Coupling in between the two simulations takes into consideration the periodicity of the disposal modules and the geometry of the repository described in Couplex 1. The mixed finite element and discontinuous Galerkin methods are used to solve the convection–diffusion equations. In order to handle the nonlinear precipitation/dissolution term, we developed a new iterative technique that combines Picard and Newton–Raphson methods.