A Comparison of Two-Phase Dynamic Upscaling Methods Based on Fluid Potentials

Although a number of methods for calculating dynamic pseudo-functions have been developed over the years, there is still a lack of understanding as to why a certain method will succeed in some cases but fail in others. In this paper, we describe the results of an assessment of several upscaling methods, namely the Kyte and Berry (KB) method, the Stone method, the Hewett and Archer (HA) method and the Transmissibility-Weighted (TW) method. We have analyzed the equations for deriving the methods and investigated the results of numerical simulations of gas displacing oil, in a variety of models to enable us to gain new insights into these, and related, upscaling methods. In particular, some novel observations on methods based on fluid potential are presented and the issue of using predicted fluid mobilities as a criterion of accuracy of an upscaling method is clarified.     

Experimental Study on Mine Spoils During Superficial Geochemistry Processes

Based on the mineral and chemical compositions of spoils taken from the Shilu Cu-Mo deposit and Hetai gold deposit, the leaching and batch experiments were made on spoils taken from these two deposits using the leaching column designed by the authors. The experimental results showed that it is not always true that the mine drainage is acidic. Its acidity depends on acid-buffering capacity of gangue, host rock and its alteration mineral assemblage. The composition of the drainage water is related with the interaction between superficial or underground water and solid materials in the spoils, including minerals, hydroxides and amorphous substances. The leaching extent of the element is related with its occurrence form in the deposit. The preferential flow results in leaching-out of heavy metals in large amounts. So it is important to prevent the generating of preferential flow in the system. The results will provide very important grounds for the comprehensive management of land and ecological rehabilitation of the mine site.     

MODELING THE FATE AND TRANSPORT OF HYDROPHOBIC ORGANIC COMPOUNDS IN AN UNSTEADY RIVER-ESTUARINE SYSTEM

1 INTRODUCTION Extensive literature (Brown et al., 1985; Sawhney et al., 1981; Bierman and Swain, 1982; Connolly, 1980; Lopez-Avila and Hites, 1980; O扖onnor, 1988) described lots of sorbed pollutants or toxic substances in bed sediments of rivers, even after the effluent was halted for a long time. This is particularly true for hydrophobic organic compounds that can be sorbed on the particles and accumulated in the river bed sediments (Karickhoff et al., 1979). Pollution events of…     

Crystal-rich mass flow deposits related to the eruption of a sublacustrine silicic cryptodome (Early Permian Collio Basin, Italian Alps)

Sedimentological, granulometric and petrographic data are presented from a detailed study on a crystal-rich mass flow deposit, which is presumably related to the eruption of a sublacustrine cryptodome. The deposit forms a prominent intercalation in the Lower Permian Collio Formation in the Italian Alps north of Brescia. Outcrops of the 10–20-m-thick volcaniclastic deposit (Dasdana I Beds, DB) can be traced over 12 km from east to west. The DB consists of a thick, crystal-rich, sandy–gravelly lower subunit representing a sequence of amalgamated Bouma-a(b) divisions overlain by a thin, well-bedded, sandy–muddy subunit that is rich in outsize porphyritic silicic fragments. Modal and computer-aided image analyses reveal that the crystal-rich lower subunit contains up to 80% of volcanogenic crystals. Some samples contain up to 60% of porphyritic fragments, which have a phenocryst content of about 20%. The wide textural range from cryptocrystalline, poikilomosaic, to rarer medium-grained granophyric groundmass, the irregular to lensoid shapes of the porphyritic fragments, and the presence of basement and sedimentary clasts suggest that the DB originated from a sublacustrine eruption of a partially extrusive cryptodome (ca. 1.6 km³). Two other porphyritic felsic cryptodomes (Dosso dei Lupi, Dosso del Bue), described briefly here, emplaced into the Collio Formation sometime after the DB event, and expose flat bases and tilted sediments at their sides. Textures observed in these domes are comparable to those found in the DB porphyritic fragments.     

Hydrogeochemistry of the Simav geothermal field, western Anatolia, Turkey

Thermal waters hosted by Menderes metamorphic rocks emerge along fault lineaments in the Simav geothermal area. Thermal springs and drilled wells are located in the Eynal, Çitgöl and Na a locations, which are part of the Simav geothermal field. Studies were carried out to obtain the main chemical and physical characteristics of thermal waters. These waters are used for heating of residences and greenhouses and for balneological purposes. Bottom temperatures of the drilled wells reach 163°C with total dissolved solids around 2225 mg/kg. Surface temperatures of thermal springs vary between 51°C and 90°C. All the thermal waters belong to Na–HCO₃–SO₄ facies. The cold groundwaters are Ca–Mg–HCO₃ type. Dissolution of host rock and ion-exchange reactions in the reservoir of the geothermal system shift the Ca–Mg–HCO₃ type cold groundwaters to the Na–HCO₃–SO₄ type thermal waters. Thermal waters are oversaturated at discharge temperatures for aragonite, calcite, quartz, chalcedony, magnesite and dolomite minerals giving rise to a carbonate-rich scale. Gypsum and anhydrite minerals are undersaturated with all of the thermal waters. Boiling during ascent of the thermal fluids produces steam and liquid waters resulting in an increase of the concentrations of the constituents in discharge waters. Steam fraction, y, of the thermal waters of which temperatures are above 100°C is between 0.075 and 0.119. Reservoir pH is much lower than pH measured in the liquid phase separated at atmospheric conditions, since the latter experienced heavy loss of acid gases, mainly CO₂. Assessment of the various empirical chemical geothermometers and geochemical modelling suggest that reservoir temperatures vary between 175°C and 200°C.     

A model of equilibrium bed topography for meander bends with erodible banks

Channel curvature produces secondary currents and a transverse sloping channel bed, along which the depth increases towards the outer bank. As a result deep pools tend to form adjacent to the outer bank, promoting bank collapse. The interaction of sediment grains with the primary and secondary flow and the transverse sloping bed also causes meanders to move different grain sizes in different proportions and directions, resulting in a consistent sorting pattern. Several models have been developed to describe this process, but they all have the potential to over‐predict pool depth because they cannot account for the influence of erodible banks. In reality, bank collapse might lead to the development of a wider, shallower cross‐section and any resulting flow depth discrepancy can bias associated predictions of flow, sediment transport, and grain‐size sorting. While bed topography, sediment transport and grain sorting in bends will partly be controlled by the sedimentary characteristics of the bank materials, the magnitude of this effect has not previously been explored. This paper reports the development of a model of flow, sediment transport, grain‐size sorting, and bed topography for river bends with erodible banks. The model is tested via intercomparison of predicted and observed bed topography in one low‐energy (5·3 W m^?2 specific stream power) and one high‐energy (43·4 W m^?2) study reach, namely the River South Esk in Scotland and Goodwin Creek in Mississippi, respectively. Model predictions of bed topography are found to be satisfactory, at least close to the apices of bends. Finally, the model is used in sensitivity analyses that provide insight into the influence of bank erodibility on equilibrium meander morphology and associated patterns of grain‐size sorting. The sensitivity of meander response to bank cohesion is found to increase as a function of the available stream power within the two study bends. Copyright © 2002 John Wiley & Sons, Ltd.     

流动成像测井研究进展

流动成像测井属于对非均匀介质动态的非线性测量，通过对油气井内多相流体信息的投影扫描和反演处理，实时成像显示流体分布及流动状况，国外研制的电导法，电容法仪器采用电导探针或电容元件构成阵列测量探头，分别利用油气与水的导电特性和介电特性差异辨识井内流体，对流动截面的测量局限于个别点上，已在水平井中见到应用效果，但对物场信息投影测量的数据量和分辨率未能满足成像要求。我国提出的电磁法采用环状阵列电极构成测量探头，综合利用油气与水的导电和介电特性差异辨识井内流体，可以测量获得多相流体流动截面的清晰图像，为研究开发全新的流动成像测井系统奠定了良好基础。     

Efficient methods for solving water flow in variably saturated soils under prescribed flux infiltration

The non-linear solvers in numerical solutions of water flow in variably saturated soils are prone to convergence difficulties. Many aspects can give rise to such difficulties and in this paper we address the gravity term and the prescribed-flux boundary in the Picard iteration. The problem of the gravity term in the Picard iteration is iteration-to-iteration oscillation as the gravity term is treated, by analogy with the time-step advance technique, ‘explicitly’ in the iteration. The proposed method for the gravity term is an improvement of the ‘implicit’ approach of Zhang and Ewen [Water Resour. Res. 36 (2000) 2777] by extending it to heterogeneous soil and approximating the inter-nodal hydraulic conductivity in the diffusive term and the gravity term with the same scheme. The prescribed-flux boundary in traditional methods also gives rise to iteration-to-iteration oscillation because there is no feedback to the flux in the solution at the new iteration. To reduce such oscillation, a new method is proposed to provide such a feedback to the flux. Comparison with traditional Picard and Newton iteration methods for a wide range of problems show that a combination of these two proposed methods greatly improves the stability and consequently the computational efficiency, making the use of small time step and/or under-relaxation solely for convergence unnecessary.     

克孜尔水库大坝的面波检测

1999年 3月 15日克孜尔地震之后对克孜尔水库坝体出现的裂缝进行了面波勘探。首先介绍了面波勘探法的基本原理 ,然后通过典型剖面检测与资料处理找出异常区、带并进行解释。试验表明面波勘探法适用于克孜尔水库大坝监测     

不同损失条件下的泥石流预报模型

泥石流预报是泥石流减灾的重要手段之一，然而泥石流形成的复杂性使泥石流预报准确度低，误报和漏报率较高。泥石流误报和漏报都会造成损失，但二者造成的损失有很大的差别。为了减少泥石流误报或漏报造成的损失，应当考虑两种错报造成损失的不同。根据使总平均损失达到最小的原则，建立了不同损失条件下的泥石流预报模型，并将该模型应用到云南东川蒋家沟。