确定边坡最危险滑动面的计算机模拟

边坡稳定性分析 ,关键是确定潜在的滑动面 ,目的是有利于工程安全。本文编制了二、三维计算机模拟程序 ,在二维模拟中提出了四种计算模式 ,三维模拟中考虑了空间、荷载、渗流效应 ,结合工程实例对边坡进行模拟并对土层抗剪强度参数进行了反演 ,模拟计算与实测结果吻合 ,模拟出的抗剪强度参数平均值与实测的抗剪强度参数标准值相当     

原油在储层介质中的加水裂解生气模拟实验

采用高温高压热模拟实验方法,开展了原油在砂岩和火成岩储层介质中的加水裂解生气模拟实验研究.结果表明,原油开始大量裂解的温度是400℃,随模拟温度增加,甲烷相对含量增大,乙烷以上重烃气尤其是丙烷相对含量减小.其中砂岩的油水混合物裂解生气主要发生在450～500℃之间,生气窗范围小,对应的烃气产率高,火成岩的油水混合物裂解生气主要发生在450～600℃之间,生气窗范围大,对应的烃气产率小.模拟烃气的组分碳同位素分馏显著,随模拟温度增加呈变重趋势.在裂解生气过程中,水解加氢和催化作用对烃气的组成、产率和碳同位素分布有重要影响。     

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

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

Evaluation of the Reduction in Uncertainty Obtained by Conditioning a 3D Stochastic Channel to Multiwell Pressure Data

A stochastic channel embedded in a background facies is conditioned to data observed at wells. The background facies is a fixed rectangular box. The model parameters consist of geometric parameters that describe the shape, size, and location of the channel, and permeability and porosity in the channel and nonchannel facies. We extend methodology previously developed to condition a stochastic channel to well-test pressure data, and well observations of the channel thickness and the depth of the top of the channel. The main objective of this work is to characterize the reduction in uncertainty in channel model parameters and predicted reservoir performance that can be achieved by conditioning to well-test pressure data at one or more wells. Multiple conditional realizations of the geometric parameters and rock properties are generated to evaluate the uncertainty in model parameters. The ensemble of predictions of reservoir performance generated from the suite of realizations provides a Monte Carlo estimate of the uncertainty in future performance predictions. In addition, we provide some insight on how prior variances, data measurement errors, and sensitivity coefficients interact to determine the reduction in model parameters obtained by conditioning to pressure data and examine the value of active and observation well data in resolving model parameters.     

Combination of Dependent Realizations Within the Gradual Deformation Method

Gradual deformation is a parameterization method that reduces considerably the unknown parameter space of stochastic models. This method can be used in an iterative optimization procedure for constraining stochastic simulations to data that are complex, nonanalytical functions of the simulated variables. This method is based on the fact that linear combinations of multi-Gaussian random functions remain multi-Gaussian random functions. During the past few years, we developed the gradual deformation method by combining independent realizations. This paper investigates another alternative: the combination of dependent realizations. One of our motivations for combining dependent realizations was to improve the numerical stability of the gradual deformation method. Because of limitations both in the size of simulation grids and in the precision of simulation algorithms, numerical realizations of a stochastic model are never perfectly independent. It was shown that the accumulation of very small dependence between realizations might result in significant structural drift from the initial stochastic model. From the combination of random functions whose covariance and cross-covariance are proportional to each other, we derived a new formulation of the gradual deformation method that can explicitly take into account the numerical dependence between realizations. This new formulation allows us to reduce the structural deterioration during the iterative optimization. The problem of combining dependent realizations also arises when deforming conditional realizations of a stochastic model. As opposed to the combination of independent realizations, combining conditional realizations avoids the additional conditioning step during the optimization process. However, this procedure is limited to global deformations with fixed structural parameters.     

Numerical Calculation of Equivalent Cell Permeability Tensors for General Quadrilateral Control Volumes

A new method for upscaling fine scale permeability fields to general quadrilateral-shaped coarse cells is presented. The procedure, referred to as the conforming scale up method, applies a triangle-based finite element technique, capable of accurately resolving both the coarse cell geometry and the subgrid heterogeneity, to the solution of the local fine scale problem. An appropriate averaging of this solution provides the equivalent permeability tensor for the coarse scale quadrilateral cell. The general level of accuracy of the technique is demonstrated through application to a number of flow problems. The real strength of the conforming scale up method is demonstrated when the method is applied in conjunction with a flow-based gridding technique. In this case, the approach is shown to provide results that are significantly more accurate than those obtained using standard techniques.     

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…     

Study on dynamics of tectonic evolution in the Fushun Basin, Northeast China

The updated study shows that the taphrogenesis of basement of the Fushun Basin is not a kind of instantaneous process. It intensified gradually and went to extreme in the sedimentary stage of the Guchengzi formation, and then, it weakened rapidly and stopped soon afterwards; the depression did not take place after the taphrogenesis. On the contrary, it almost happened simultaneously with the taphrogenesis. The depression went at a high speed from the beginning of the sedimentary period of the Xilutian formation, and then weakened gradually in the sedimentary period of the Gengjiajie formation. The evolution course of the synsedimentary structure of the Fushun Basin can be summarized as the following six stages: slow taphrogenesis and high speed depression to accelerated taphrogenesis and high speed depression to high speed taphrogenesis and high speed depression to retarded taphrogenesis and high speed depression to gradual halt of taphrogenesis and reduced depression to slow depression and gradual halt of depression. The tectonic evolution resulted in the formation of the "lower taphrogenesis and upper depression" structure. The formation of the binary structure might be due to the suspension of taphrogenesis and the change of the regional structure stress field, but the depression kept going. The result of calculation combining the analysis of the synsedimentary structural frame, the back-stripping method of the subsidence history of the basin basement and the simulation of thermo-settlement history indicates that the great sedimentary space required by the "upper depression part" consists of two parts, namely, 40% from compaction of sediments and 60% from slow depression of the basin basement during a long period of time. Gradual halt of the depression in the Fushun Basin may be attributed to the reversal of the lithosphere hot-recession and gravity isostasy adjustment which may be the result of new hot-events in the depths and accompanied invasion of extremely thick diabase sill, thus revealing a new forming mechanism of "fault subsidence at the base and depression on the top" structure.     

A quantitative investigation into the temporal and spatial variations in the emergence of adult blackflies (Diptera: Simuliidae) from the Breitenbach, a small upland stream in Germany

As a contribution to the long-term emergence studies carried out as part of the “Breitenbach ecosystem project”, this paper presents the results obtained on emergence patterns and population dynamics of blackflies from 1984 to 1988. The Breitenbach is a small first order stream near Schlitz/Hesse, running into the Fulda river. Adult blackflies were caught in 4 greenhouse emergence traps, each spanning the whole width of the stream for a length of 6 m. A manual and then a partly automated method for collecting trapped insects was applied for three years and one year, respectively. Sixteen blackfly-species of different origin were found in the traps: a) autochthonous species: Prosimulium tomosvaryi, Simulium vernum s.l., S. cryophilum s.l., S. ornatum s.l. (common); S. costatum, S. trifasciatum, S. monticola, S. argyreatum (rare); b) species of doubtful origin: S. lundstromi, S. angustitarse, S. angustipes, S. aureum (very rare); c) allochthonous species: S. lineatum, S. equinum, S. erythrocephalum, S. noelleri. It was shown that the last 4 species had not emerged from the Breitenbach but had flown into the traps as adults. In addition, females of autochthonous species with blood or with mature eggs were trapped, which were also considered to have flown in. Besides revealing a limitation of the trap construction, they supported the detailed interpretation of some intricated patterns of appearance. Variations in emergence patterns and specimen numbers between years and traps were pronounced, but only in a few cases could they be attributed to changes in abiotic factors such as water temperature or discharge. P. tomosvaryi had the most simple and regular life cycle, with one well synchronised emergence peak annually from April to May (or even to June). No gradient of specimen numbers along the stream was evident. The two closely related species S. vernum and S. cryophilum had quite similar emergence patterns: There were two broad peaks per year, extending mainly from March to June and from July to October or November. In some cases the number and separation of consecutive generations was not clear. The abundance of both species clearly decreased downstream, more so for S. cryophilum than for S. vernum. In a trap closest to a tributary spring, S. cryophilum was the dominant species during three of the four years examined. The adults of S. ornatum displayed an intricate pattern of appearance, with very low specimen numbers in spring and medium to very high numbers in July/August and September/October. S. ornatum is the only blackfly species that inhabits both the Breitenbach and the adjacent section of the Fulda river. It was shown that females emerging from the Fulda river regularly invade the Breitenbach valley in greatly varying numbers and oviposit there. This leads to overlapping larval cohorts with corresponding emergence peaks. Although S. ornatum was the most abundant species in one year in the lower traps (60 to 80% of all individuals), it remains uncertain whether it is a long-term, permanent member of the autochthonous blackfly fauna of the Breitenbach. Estimates of total numbers of flown-in adults, actual emergence, dry weight biomass, the ecological separation of the species and their life cycle strategies are discussed.     

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

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