地下水数值计算中断层处理的"切割-导通法"

从对断层导水性能的分析出发 ,提出了地下水数值计算中断层处理的“切割—导通法” ,既考虑断层侧向导水性 ,又考虑了其垂向导水性 ,更符合断层导水的实际情况 ,也易于描述较大断裂在走向上的导水性差异。     

Propriétés fractales de la fragmentation et processus stochastiques de fracturation : approche géométrique 3D à l'aide du modèle ObsifracFractal properties of fragmentation and stochastic fracturing processes: geometrical 3D approach by OBSIFRAC model

A numerical rock fragmentation model was elaborated, producing a 3D puzzle of convex polyhedra, geometrically described in a database. In the first scenario, a constant proportion of blocks are fragmented at each step of the process and leads to fractal distribution. In the second scenario, division affects one random block at each stage of the process, and produces a Weibull volume distribution law. Imposing a minimal distance between the fractures, the third scenario reveals a power law. The inhibition of new fractures in the neighbourhood of existing discontinuities could be responsible for fractal properties in rock mass fragmentation. To cite this article: L. Empereur-Mot, T. Villemin, C. R. Geoscience 334 (2002) 127–133.     

2-D viscoelastic FEM simulation on stress state in the deep part of a subducted slab

Introduction It is found that there are some relationships between the thermal structures of subduction zones and the deep seismicity, while the mechanism relates the thermal structure and the deep seismicity is still unsure (Helffrich, Brodholt, 1991; Furukawa, 1994; Kirby, et al, 1996). From 1980s, geoscientists have constituted a series of numerical simulations on the stress states of subduction slabs. Based on the kinetic computation of Sung and Burns (1976a, b), Goto, et al (1983, 1987…     

含预制软弱带的岩石破裂过程的数值模拟

考虑到岩石脆性破坏过程中介质的不均匀性之特点，对岩石样品中含预制弱介质条带的岩石样品破坏过程进行二维有限元数值模拟，并对弱介质带的破坏贯通过程、新的断层的产生和有关的地震活动进行了研究．数值模型展示了岩石从变形、微观破坏到整体破坏的全过程以及微震活动的时空分布特征．应力、应变和微震活动的时空分布形象地描绘了岩石变形的局部化和时空迁移等现象，这与实际地壳中所观测到的现象是一致的．此外，模拟结果与实验观测结果也是一致或相似的．      

Improved treatments for evaluating horizontal magnetic components through the 3-D FDM in E -polarization induction problems

To improve the accuracy of the numerical evaluation through the 3-D finite difference method, the surface boundary conditions are added to modify the old program. The author has tested the new program by making calculations for the model constructed by Wanamaker, et al (1984). The comparison between the numerical results obtained from this paper and those by Wannamaker, et al (1984) indicates that a pronounced improvement is realized in the evaluation of the horizontal magnetic components. Moreover, better calculations for the vertical magnetic components are also obtainable by using the new program.     

地震次生毒气泄漏与扩散的数值模拟与动态仿真

文章对破坏性地震中可能发生的次生毒气泄漏与扩散灾害进行了数值模拟与动态仿真研究，设定了几种有实际意义的数值分析模型，分别进行了结构地震破坏，有毒有害气体泄漏，在一定泄漏和气象条件下的向周围空间扩散的动态仿真研究，给出了典型毒气扩散过程中可致人员死亡，严重危险和有感区域的范围的动态变化结果，作者认为地震次生毒气泄漏与扩散研究对大地震应急、城市和大型企业防灾都有重要的意义和广阔的应用前景。     

Comparison of Seismic Dispersion and Attenuation Models

The frequency-dependent attenuation of seismic waves causes decreased resolution of seismic images with depth, and the difference in transmission losses induces amplitude variations with offset. Transmission losses may occur due to friction or fluid movement, or may result from scattering in thin-layer. Whatever the physical mechanism, they can often be conveniently described using an empirical formulation wherein the elastic moduli and propagation velocity are complex functions of frequency.We have compiled and compared algebraically and numerically eight different models involving complex velocity: the Kolsky-Futterman model, the power-law model, Kjartansson's model, Müller's model, Azimi's second and third model, the Cole-Cole model, and the standard linear-solid model.For two different parameter sets, the attenuation and phase velocity are computed in the seismic frequency band, and the plane-wave propagation of a Ricker wavelet for the other models is compared with that for the Kolsky-Futterman model. The first parameter set consists of parameters for each of the models calculated from expressions given in the appendix. These expressions make the different models behave similarly to the KF model. The second parameter set consists of model parameters that are numerically adapted to the KF model.By selecting proper parameters, all models, except the standard linear-solid model, show behavior similar to that of the Kolsky-Futterman model. The SLS model behaves differently from the other models as the frequency goes to zero or infinity. Broadband measurement data is needed to select a specific model for a given seismic experiment.     

Study on the characteristics of rock failure strain and acoustic emission field for two parallelling faults with the same slip direction including asperities

Introduction Concerning about research on grouped-earthquake seismogenic and occurrence, at present it is only limited on the stage of the theoretical analysis and numerical model. For instance, ZHANG (1987) simulated the grouped-earthquake generation and physical mechanism making use of spring-block principle. Applying for non-linear dynamic model, referring to Chinese continental plate as object, LI, et al (1997) performed the numerical simulation on the seismogenic and occurring of group…     

Natural hydraulic cracking: numerical model and sensitivity study

Natural hydrofracturing caused by overpressure plays an important role in geopressure evolution and hydrocarbon migration in petroliferous basins. Its mechanism is quite well understood in the case of artificial hydraulic fracturing triggered by high-pressure fluid injection in a well. This is not so for natural hydraulic fracturing which is assumed to initiate as micro-cracks with large influence on the permeability of the medium. The mechanism of natural hydraulic cracking, triggered by increasing pore pressure during geological periods, is studied using a fracturing model coupled to the physical processes occurring during basin evolution. In this model, the hydraulic cracking threshold is assumed to lie between the classical failure limit and the beginning of dilatancy. Fluid pressure evolution is calculated iteratively in order to allow dynamic adjustment of permeability so that the fracturing limit is always preserved. The increase of permeability is interpreted on the basis of equivalent fractures. It is found that fracturing is very efficient to keep a stress level at the rock’s hydraulic cracking limit: a fracture permeability one order of magnitude larger than the intrinsic permeability of the rock would be enough. Observations reported from actual basins and model results strongly suggest that natural hydraulic cracking occurs continuously to keep the pressure at the fracturing limit under relaxed stress conditions.     

Models of passive margin inversion: implications for the Rhenohercynian fold-and-thrust belt, Belgium and Germany

Positive tectonic inversion is related to the transmission of compressional stresses along a décollement into the foreland of an orogenic zone. This stress and strain concentration in regions remote from the main orogenic front is commonly related to the presence of pre-existing rheological heterogeneities such as normal syn-depositional faults. During inversion, these pre-existing normal faults are reactivated as reverse faults. Tectonic inversion in the Rhenohercynian fold-and-thrust belt during the Variscan Orogeny shows that inversion is likely synchronous with the onset of collision in the hinterland. Here, we present the results of a simplified thermo-mechanical model (STM) which allows one to study strain partitioning between two orogenic zones. We show that, if the two orogenic zones have the same mechanical properties, the viscosity of the décollement, which links them, controls the initial strain partitioning. During subsequent finite shortening, erosional processes determine the partitioning of strain rate. The presence of a weak structure in the inverted zone and of a low-viscosity décollement leads to initial strain concentration in the inverted track rather than in the collision zone and a progressive decrease in strain partitioning between the two orogenic zones. The STM results are in good agreement with results of a 2D finite-element model. We conclude that, in the western part of the Rhenohercynian Massif, simultaneous uplift and deformation within the Mid-German Crystalline Rise (the main collision zone) and the Ardenne Anticlinorium (the inverted zone) lead to interpreting this orogenic event as a case of vice tectonic rather than the propagation of a ‘wave of folding’ towards the Variscan front, as suggested by previous authors.